HDUOJ-1003 Max Sum(最大子序列和)(动态规划)

O(n)算法

首先,引用一个前人造好的轮子,此轮子来自《数据结构与算法分析-C语言描述》图2.8:

int MaxSubsequenceSum(const int A[],int N)
{
    int ThisSum,MaxSum,j;

    ThisSum = MaxSum = 0;
    for(j=0;j<N;j++)
    {
        ThisSum += A[j];
        if(ThisSum > MaxSum)
            MaxSum = ThisSum;
        else if(ThisSum < 0)
            ThisSum = 0;    
    }
    return MaxSum;
}

对于此算法可以作此理解:
对于(0-N)的每个(A[j]-A[N]),都做一次累加,将每次累加的结果做一次判断(小于零则舍弃),大于MaxSum则统计到MaxSum里,则N次循环内,必然会有一个最优解被更新为MaxSum的值。

状态转移方程为:

dp[i] = max{ dp[i-1]+A[i] , A[i] }

考虑到结尾并不需要每一个dp[i] , 其实用两个变量代替就可以:

ThisSum = max{LastSum+A[i] , A[i]}

本题的要求同时输出最大子序列的上界和下界,需要几个变量来统计。
改造后的完整代码为:

/* * Created by zsdostar in 2016/4/30 */
#include<iostream>
using namespace std;

int start,endd;
int startm,endm;
int ThisSum,LastSum;
int maxx = -2147483648;

int  MaxSubSequence(int arr[],int n)
{
      for(int i=0;i<n;i++)
      {
          if(LastSum>=0)
          {
              ThisSum = LastSum + arr[i];endd = i;
          }
          else
          {
              ThisSum = arr[i];start = endd = i;
          }
          if(maxx<=ThisSum)
          {
              maxx = ThisSum;
              startm = start;endm = endd;
          }
          LastSum = ThisSum;
      }
      return maxx;
}
int main()
{
  int count,len;
  int array[100001];
  cin>>count;
  for(int flag=1;flag<=count;flag++)
  {
      cin>>len;
      for(int i=0;i<len;i++)
          cin>>array[i];
      start = endd = 0,maxx = -1001;
      startm = endm = 0;
      ThisSum = LastSum = 0;
      cout<<"Case "<<flag<<":"<<endl<<MaxSubSequence(array,len);
      cout<<" "<<startm+1<<" "<<endm+1<<endl;
      if(flag<count)
          cout<<endl;
  }
  return 0;
}

关于本题有两个问题:

目前的轮子通过大于小于0来判定是否更新状态,那么如果整个数组全都是负数,或者正数数量等于1,那么结果是错的
HDU1003的测试数据也没有考虑问题1,导致考虑不全的代码也能通过

下次复习这道题的时候解决这个问题,更新考虑更全面的代码。

    原文作者:动态规划
    原文地址: https://blog.csdn.net/zsdostar/article/details/51286482
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞