本文转载自kuangbin博客:http://www.cnblogs.com/kuangbin/archive/2011/08/04/2127085.html
题目链接:http://acm.split.hdu.edu.cn/showproblem.php?pid=1024 (近期HDU域名更新)
本题的大致意思为给定一个数组,求其分成m个不相交子段和最大值的问题。
设Num为给定数组,n为数组中的元素总数,Status[i][j]表示前i个数在选取第i个数的前提下分成j段的最大值,其中1<=j<=i<=n && j<=m,状态转移方程为:
Status[i][j]=Max(Status[i-1][j]+Num[i],Max(Status[0][j-1]~Status[i-1][j-1])+Num[i])
乍看一下这个方程挺吓人的,因为题中n的限定范围为1~1,000,000而m得限定范围没有给出,m只要稍微大一点就会爆内存。但仔细分析后就会发现Status[i][j]的求解只和Status[][j]与Status[][j-1]有关所以本题只需要两个一维数组即可搞定状态转移。
在进行更进一步的分析还会发现其实Max(Status[0][j-1]~Status[i-1][j-1])根本不需要单独求取。在求取now_Status(保存本次状态的数组)的过程中即可对pre_Status(保存前一次状态的数组)进行同步更新。
#include<iostream>
using namespace std;
#include<cstdio>
#include<algorithm>
#define maxn 1000010
#define INF 0x7fffff
int dp[maxn];
int mmax[maxn];
int a[maxn];
int main()
{
int N,M;
int i,j,mmmax;
while (scanf("%d%d",&M,&N)!=EOF)
{
for (int i=1;i<=N;i++)
{
cin>>a[i];
mmax[i]=0;
dp[i]=0;
}
dp[0]=0,mmax[0]=0;
for (int i=1;i<=M;i++)
{
mmmax=-INF;
for (int j=i;j<=N;j++)
{
dp[j]=max(dp[j-1]+a[j],mmax[j-1]+a[j]); //dp[i][j]=Max(dp[i][j-1]+a[j] , max( dp[i-1][k] ) + a[j] )
mmax[j-1]=mmmax; ///我们可以在每次计算dp[i][j]的时候记录下前j个
///的最大值 用数组保存下来
mmmax=max(mmmax,dp[j]);
}
/*
原来的状态转移方程是这样的 Status[i][j]=Max(Status[i-1][j]+Num[i],Max(Status[0][j-1]~Status[i-1][j-1])+Num[i])
实际上 Status[i][j]的求解只和Status[*][j]与Status[*][j-1]有关所以本题只需要两个一维数组即可搞定状态转移
对于程序的r29 r30:
更进一步的分析还会发现其实Max(Status[0][j-1]~Status[i-1][j-1])根本不需要单独求取。在求取now_Status(保存本次状态的数组)的过程中即可对pre_Status(保存前一次状态的数组)进行同步更新
*/
}
cout<<mmmax<<endl;
}
return 0;
}