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,导致考虑不全的代码也能通过
下次复习这道题的时候解决这个问题,更新考虑更全面的代码。