POJ1050 dp入门

最新推荐文章于 2018-09-25 22:41:45 发布

jiaowopan

最新推荐文章于 2018-09-25 22:41:45 发布

阅读量478

点赞数

分类专栏： OJ

本文链接：https://blog.csdn.net/jiaowopan/article/details/8757400

版权

OJ 专栏收录该内容

17 篇文章 0 订阅

订阅专栏

先来看最大子序列和。有一串数，有正有负，如2，-1，5，4，-9，7，0，3，-5。求：这一串数中，和最大的一段。比如说，从第一个数2开始，发现下一个为-1，加下-1后和显然会变小。再往后看，第三个数是5，所以上一个-1还是要选的，这样才能加上5。哎，不看了，这样求最大和还不得累死。嘿嘿，这时DP就派上用场了。

设这串数为X1 X2 X3 … Xn, 用dp(i,j)表示从Xi…Xj的最大子序列和。

按照DP的思路，想办法减小问题的规模。有n个数，怎样能减少到n-1数？想办法把最后一个数Xn去掉，问题规模就能减少到n-1。

通过观察可以发现：X1…Xn的最大子序列可以分为两类：以Xn结尾、不以Xn结尾。不以Xn结尾的最大子序列，其实就是X1…Xn-1的最大子序列，发现这点很重要。

这样就有：dp( i, j ) = Max( dp( i, j-1 ), Last( j ) ).其中Last( j )表示以Xj结尾的最大子序列的和。

功夫不负有心人，终于把问题规模减少了。但是，一波未平一波又起，新的问题又出现了。Last( j )如何求？即，求以Xj结尾的最大子序列的和。再用DP求解。

Last( j )和Last( j-1 )之间的关系比较简单。Last(j )的值里面必然会包括Xj的值，到底有没有Last( j-1 )也很简单，主要取决于Last( j-1 )是正还是负。

这样就有：Last( j ) = Max( Xj, Last( j-1) + Xj );

状态转换方程：

dp( i, j ) = Max( dp( i, j-1 ), Last( j ) )其中：dp(i,j)表示从Xi…Xj的最大子序列和

Last( j ) = Max( Xj, Last( j-1 ) + Xj ); 其中：Last( j )表示以Xj结尾的最大子序列的和

现在，回到POJ1050。想想能不能利用上面的结果？求最大子矩阵，那么只要确定了子矩阵有几行、几列即可。这样，可以枚举子矩阵的行数和列数。

比如，当子矩阵只要一行时，那么只关心它的列从哪开始到那结束就行。哦，这其实就是一个最大子序列和的问题。这一行就是这一串数，求和最大的一段。那么当子矩阵有两行时，怎么办？如何把两行变为一行？一个聪明的想法就是：把这两行按照对应的列加起来。

好了问题已经漂亮的解决了：在原矩阵中任意画出一部分，然后按照对应的列加起来，问题就转变为一个最大子序列和的问题

以上分析转载自 http://blog.csdn.net/sj13051180/article/details/6667577

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;
#define NUM 105
#define Best( x , y )  ( ( (x) < (y) ) ? (y) : (x) )
int Arr[NUM][NUM];
int N;

int dp[NUM];
int Last[NUM];

int data[NUM];

int getMaxSum(int* d) //对一维数组求最大子串和 杭电1003
{
	memset(dp,0,sizeof(dp));
	memset(Last,0,sizeof(Last));
	dp[0]=d[0];
	Last[0]=d[0];
	int i;
	for(i=1;i<N;i++)
	{
		Last[i]=Best(Last[i-1]+d[i],d[i]);
		dp[i]=Best(dp[i-1],Last[i]);
	}
	return dp[N-1];
}
void Solve()
{
	int maxsum=-100000;
	memset(data,0,sizeof(data));
	int i,j,h,k;
	for(i=0;i<N;i++)
	{
		for(j=i;j<N;j++)
		{
			//compute from line i to line j
			memset(data,0,sizeof(data));
			for(h=0;h<N;h++)
			{
				for(k=i;k<=j;k++)
					data[h]+=Arr[k][h];
			}
			int tmpsum=getMaxSum(data);//求最大子串和
			maxsum=tmpsum>maxsum ? tmpsum : maxsum;
		}
	}
	printf("%d\n",maxsum);
}
int main()
{
	//freopen("poj1050.txt","r",stdin);
	memset(Arr,0,sizeof(Arr));
	while(scanf("%d",&N)!=EOF)
	{
		int i,j;
		for(i=0;i<N;i++)
			for(j=0;j<N;j++)
				scanf("%d",&Arr[i][j]);
		Solve();
	}
	return 0;
}