题目
给定一个数字序列A1,A2,…,An,求i,j(1<=i<=j<=n),使得Ai+…+Aj最大,输出这个最大和。
样例
-2 11 -4 13 -5 -2
显然11+(-4)+13=20为和最大的选取情况
分析
- 步骤1:令状态dp[i]表示以A[i]作为末尾的连续序列的最大和,以样例为例:序列-2 11 -4 13 -5 -2,下标分别记为0 1 2 3 4 5,那么
dp[0]=-2
dp[1]=11
dp[2]=7(11+(-4)+13=20)
dp[4]=15(因为由dp数组的含义,A[4]=-5必须作为连续序列的结尾,于是最大和就是11+(-4)+13+(-5)=15,而不是20)
dp[5]=13(11+(-4)+13+(-5)+(-2)=13)
通过设置一个dp数组,要求的最大和其实就是dp[0],dp[1],…dp[n-1]中的最大值
2. 步骤2:因为dp[i]要求是必须以A[i]结尾的连续序列,那么只有两种情况
1)这个最大和的连续序列只有一个元素,即以A[i]开始,以A[i]结尾
2)这个最大和的连续序列有多个元素,即从前面某处A[p]开始(p<i),一直到A[i]结尾。
对第一种情况,最大和就是A[i]本身
对第二种情况,最大和是dp[i-1]+A[i],即A[p]+…+A[i-1]+A[i]=dp[i-1]+A[i]。
只有这两种情况,于是得到状态转移方程:
dp[i] = max{A[i],dp[i-1]+A[i]}
这个式子只和i与i之前的元素有关,且边界为dp[0]=A[0],由此从小到大枚举i,即可得到整个dp数组。接着输出dp[0],dp[1],…,dp[n-1]中最大值即为最大了连续子序列的和。
#include <cstdio>
#include <algorithm>
using namespace std;
const int maxn = 10010;
int A[maxn],dp[maxn]; //A[i]存放序列,dp[i]存放以A[i]结尾的连续序列的最大和
int main()
{
int n;
scanf("%d", &n);
for(int i = 0; i < n; i++) //读入序列
{
scanf("%d", &A[i]);
}
//边界
dp[0] = A[0];
for(int i = 1; i < n; i++)
{
//状态转移方程
dp[i] = max(A[i], dp[i-1]+A[i]);
}
//dp[i]存放以A[i]结尾的连续序列的最大和,需要遍历i得到最大的才是结果
int k = 0;
for(int i = 1; i < n; i++)
{
if(dp[i] > dp[k])
{
k = i;
}
}
printf("%d\n", dp[k]);
return 0;
}
-2 11 -4 13 -5 -2
20
状态的无后效性
当前状态记录了历史信息,一旦当前状态确定,就不会再改变,且未来的决策只能在已有的一个或若干个状态的基础上进行,历史信息只能通过已有的状态去影响未来的决策。
如何设计状态和状态转移方程,才是动态规划的核心,而它们也是动态规划最难的地方。