poj1050 To the Max dp

题意：求一个矩阵的最大子矩阵和。

思路：首先预处理出sum[ i ][ j ][ k ]，即以（i，j）为最上端向下连续k个数的和。接着设dp[ i ][ j ][ k ]为以（i，j）为左上端向下连续k行

的最大矩阵和，那么我们可以很容易写出状态转移方程：dp[ i ][ j ][ k ]=max( 0,dp[ i ][ j + 1 ][ k ] ) + sum[ i ][ j ][ k ]。详见代码：

// file name: poj1050.cpp //
// author: kereo //
// create time:  2014年09月04日 星期四 16时45分31秒 //
//***********************************//
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<stack>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXN=100+5;
const int inf=0x3fffffff;
#define L(x) (x<<1)
#define R(x) (x<<1|1)
int n;
int a[MAXN][MAXN];
int dp[MAXN][MAXN][MAXN],sum[MAXN][MAXN][MAXN];
int main()
{
	while(~scanf("%d",&n)){
		memset(dp,0,sizeof(dp));
		memset(sum,0,sizeof(sum));
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++){
				scanf("%d",&a[i][j]); 
				sum[1][j][i]=sum[1][j][i-1]+a[i][j];
		}
		for(int j=1;j<=n;j++)
			for(int i=1;i<=n;i++)
				for(int k=1;k<=n-i+1;k++)
					sum[i][j][k]=sum[1][j][i+k-1]-sum[1][j][i-1];
		for(int i=n;i>=1;i--)
			for(int j=n;j>=1;j--)
				for(int k=1;k<=n-i+1;k++)
					dp[i][j][k]=max(0,dp[i][j+1][k])+sum[i][j][k];
		int ans=0;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				for(int k=1;k<=n-i+1;k++)
					ans=max(ans,dp[i][j][k]);
		cout<<ans<<endl;
	}
	return 0;
}