杭电1024 Max Sum Plus Plus

<span style="font-family: Arial, Helvetica, sans-serif; white-space: normal; ">题目内容:</span>

Problem Description
Now I think you have got an AC in Ignatius.L's "Max Sum" problem. To be a brave ACMer, we always challenge ourselves to more difficult problems. Now you are faced with a more difficult problem.

Given a consecutive number sequence S 1, S 2, S 3, S 4 ... S x, ... S n (1 ≤ x ≤ n ≤ 1,000,000, -32768 ≤ S x ≤ 32767). We define a function sum(i, j) = S i + ... + S j (1 ≤ i ≤ j ≤ n).

Now given an integer m (m > 0), your task is to find m pairs of i and j which make sum(i 1, j 1) + sum(i 2, j 2) + sum(i 3, j 3) + ... + sum(i m, j m) maximal (i x ≤ i y ≤ j x or i x ≤ j y ≤ j x is not allowed).

But I`m lazy, I don't want to write a special-judge module, so you don't have to output m pairs of i and j, just output the maximal summation of sum(i x, j x)(1 ≤ x ≤ m) instead. ^_^
 

Input
Each test case will begin with two integers m and n, followed by n integers S 1, S 2, S 3 ... S n.
Process to the end of file.
 

Output
Output the maximal summation described above in one line.
 

Sample Input
  
  
1 3 1 2 3 2 6 -1 4 -2 3 -2 3
 

Sample Output
  
  
6 8
Hint
Huge input, scanf and dynamic programming is recommended.


题目大意:

首先给予变量m、n,即对于给定的长度为n的数字序列,将其分为m段子序列,使得这m段子序列加起来取得最大值


题目分析:

这是一道典型的dp问题,可以认为是从最长子序列这一问题变形而来,因此其核心算法中必然包含关于单段最长子序列的内容,所不同的是可以进行m段的分割,根据dp的解题思路寻找状态转移方程是解决问题的方案;

根据题目进行分析可以得到如下状态转移方程 f(i, j) = max { f(i - 1, j) + v[i], max {f(k, j - 1) + v[i]}(k = j - 1 ... i - 1) }

f(i,j)意为序列的前i项分为j段,而其中max所包含的是f在状态转移过程中的两种方案,即前i-1项分为j段,同时将第i项归入第j段中,而另一种方案则是前k项分为j-1段,同时将第i项单独成段

但可以看出该算法的时间开支非常之大,因而为了优化时间结构,引入了辅助数组g,g(i,j)代表前i项分为j项时的最大值 (第i项不一定被包含在内)

而g(i,j)所使用的即是常规的最长子序列算法,其dp的转移项是是否加入第i项

所以此时总的状态转移方程可以重新写为

f(i, j) = max{f(i - 1, j), g(i - 1, j - 1)} + v[i],转移变成了O(1) 

g[i][j]与g[i-1][j]或f[i][j]相等

f[i][j]-v[i]与g[i-1][j-1]或f[i-1][j]相等

将其空间优化为一位数组有

在第i行时

f[j]=max(f[j],g[j-1])+v[i];

g[j]=max(g[j],f[j]));


题目代码:

#include<cstdio>

#define min -99999999;

int num[1000005];
int f[1000005];
int g[1000005];

int max(int x,int y)
{
	return x>y?x:y;
}

int minn(int x,int y)
{
	return x<y?x:y;
}

int main()
{
	int m,n;
	while(scanf("%d%d",&m,&n)==2&&m&&n)
	{
		int sum=0;
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&num[i]);
			f[i]=g[i]=min;
		}
		f[1]=g[1]=num[1];
		int i,j;
		for(i=2;i<=n;i++)
		{
			int t=minn(i,m);
			for(j=1;j<=t;j++)
			{
				f[j]=max(f[j]+num[i],g[j-1]+num[i]);
				g[j-1]=max(g[j-1],f[j-1]);
				//printf("%d %d %d %d\n",f[j],g[j],i,j);
			}
			g[j-1]=max(g[j-1],f[j-1]);
			//printf("\n");
		}
			sum=g[m];
			printf("%d\n",sum);
	}
	return 0;
}
//http://www.2cto.com/kf/201308/234556.html
//http://blog.csdn.net/mylovepanning/article/details/8121397
其中最后两行为所参考的博文


解题评价:

这道题属于dp方面的经典题目,由于之前并未接触过类似题目,因而所花时间较长,最后还是依靠两篇博文略有理解,掌握多段最长子序列算法对于今后解题会有较大帮助

当前水平评分:9


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值