一道经典的DP问题------求最大子序列和
算法思想:首先归纳出递归表达式 f[i+1] = f[i]>0?(f[i]+a[i+1]):a[i+1];
1、对这个式子解释一下,f[i]是以a[i]为子序列区间的右端点的最大子序列和,比如(6,-1,5,4,-7),那么f[4] = (6-1+5) + 4 = 14。
2、注意对每一个i 属于 [1...n],f[i]都以a[i]作为右端点,即计算每一个f[i],都要加上a[i]。看上去似乎a[i]对f[i]的贡献只是普通的增加一项而已,但如果我们反过来看这个问题的逻辑,即把a[i]看作是f[i]的主项,考察f[i-1]对主项a[i]的增减影响,那么
if f[i-1]>=0,then f[i-1] + a[i] >= a[i];
else f[i-1] + a[i] < a[i];
这里为什么把a[i]看作计算f[i]的主项呢,原因很简单:从f[i]的定义看,从1...n,依次计算f[i],都是要把a[i]看作是子序列右端点,也即子序列的一部份,这既保证了主项a[i]是f[i]中必然要出现的项,也保证了f[i-1]+a[i]计算出来的f[i]是连续子序列的和。
#include<stdio.h>
#include<string.h>
int a[1000001]; //numbers of the sequence
int begin,end,max;
void max_sum(int n)
{
int start,tmax;
int i;
start = begin = end = 0;
tmax = max = a[0];//初始化当前子序列和最大子序列为a[0]
for(i=1;i<n;++i){
if(tmax >= 0){
tmax += a[i];
}
else{
tmax = a[i]; //抛弃当前子序列
start = i; //开始新的子序列搜索
}
if(max < tmax){ //更新最大子序列
max = tmax;
begin = start;
end = i;
}
}
return;
}
int main()
{
freopen("in","r",stdin);
int T,i;
scanf("%d",&T);
for(i=0;i<T;i++){
int n,j;
scanf("%d",&n);
for(j=0;j<n;j++){
scanf("%d",&a[j]);
//f[j] = min;
}
max_sum(n);
printf("Case %d:\n",i+1);
printf("%d %d %d",max,begin+1,end+1);
if(i!=T-1) printf("\n\n");
else printf("\n");
}
return 0;
}