smp 1117/1238 最大矩形区域/三维最大子块和（动态规划）

Description

在一个矩形的区域中每个格点(整数坐标点)都有一个价值，请你在这区域中找出价值和最大的矩形区域。
Input

输入有若干种情况，每种情况的第一行是两个正整数m、n，表示矩形的长和宽（1<=m、n<=100）。接着是m行、n列的整数（每个数的绝对值不超过300）表示对应格点的价值。
Ouput

对应输出找出矩形的最大价值和。
Sample Input

4 5
-1 -3 1 1 1
0 -2 3 2 -5
-2 5 3 4 1
-6 3 7 5 0
Sample Output

30

• 思路：
  dp[i][j]+=(dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]);
  记录以左上角（1,1）到右下角（i，j）该子矩阵的和
  例如：
  dp[3][4]+=(dp[2][4]+dp[3][3]-dp[2][3]);
  dp[2][4]为：
  0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzE4ODQxNzYx,size_16,color_FFFFFF,t_70)
  dp[3][3]为：
  
  dp[2][3]为：
  
  显而易见：dp[3][4]可以用(dp[3][4]+dp[2][4]+dp[3][3]-dp[2][3])表示;
  

然后
ans=max(ans,dp[i][j]-dp[i-a][j]-dp[i][j-b]+dp[i-a][j-b]);
是计算 以（i，j）为右下角的所有小矩阵的和，ans记录其最大值

代码：

#include<iostream>
#include<cstring>
#define max(a,b) a>b?a:b
using namespace std;
int main()
{
	ios::sync_with_stdio(false);cin.tie(0);
	int m,n,i,j,dp[105][105];
	
	memset(dp,0,sizeof(dp));
	while(cin>>m>>n)
	{
		int ans=-999;
		for(i=1;i<=m;i++)
		{
			for(j=1;j<=n;j++)
			{
				cin>>dp[i][j];
				dp[i][j]+=(dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1]);
				
				
				for(int a=1;a<=i;a++)
				{
					for(int b=1;b<=j;b++)
					{
						ans=max(ans,dp[i][j]-dp[i-a][j]-dp[i][j-b]+dp[i-a][j-b]);
					}
				}
			}
		}
		
		cout<<ans<<endl;
	}
}

同理，解决三维最大子块和：

Description

给出一个三维的矩阵,它的子矩阵的值定义为子矩阵中所有值之和,求它的最大子矩阵和。
Input

先输入一个正整数T，表示有T种情况。每种情况的第一行输入三个正整数t、n、m，表示方体的高、长、宽，然后有t个n*m的数，以矩阵的阵列方式输入。
Ouput

每种情况一行输出，最大块的和。

Sample Input

1
2 2 2
7 -4
3 7
-9 -6
4 7
Sample Output

21

方法：三维压缩成二维，二维直接计算

其中三维压二维代码：

	for(i=1;i<=t;i++)
	{
		memset(dp,0,sizeof(dp));
		for(j=i;j<=t;j++)
		{				
			for(k=1;k<=m;k++)
			{
				for(w=1;w<=n;w++)
				{
					dp[k][w]+=d[j][k][w];
				}
			}
		
			ans=max(ans,work(dp,m,n));
	}

AC代码：

#include<iostream>
#include<cstring>
#define max(a,b) a>b?a:b
using namespace std;

int work(int dp[50][50],int m,int n)
{
	int i,j,temp[50][50];
	
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
			temp[i][j]=dp[i][j];	
	
	
	int res=-999;
	for(i=1;i<=m;i++)
	{
		for(j=1;j<=n;j++)
		{
			temp[i][j]+=(temp[i-1][j]+temp[i][j-1]-temp[i-1][j-1]);
				
					
			for(int a=1;a<=i;a++)
			{
				for(int b=1;b<=j;b++)
				{
					res=max(res,temp[i][j]-temp[i-a][j]-temp[i][j-b]+temp[i-a][j-b]);
				}
			}
		}
	}
	return res;
}

int main()
{
	int d[50][50][50],dp[50][50];
	int T;
	cin>>T;
	while(T--)
	{
		memset(d,0,sizeof(d));
		int t,m,n,i,j,k,w;
		while(cin>>t>>m>>n)
		{
			for(k=1;k<=t;k++)
			{
				for(i=1;i<=m;i++)
				{
					for(j=1;j<=n;j++)
					{
						cin>>d[k][i][j];
					}
				}		
			}
			
			int ans=-9999;
			for(i=1;i<=t;i++)
			{
				memset(dp,0,sizeof(dp));
				for(j=i;j<=t;j++)
				{
					
					for(k=1;k<=m;k++)
					{
						for(w=1;w<=n;w++)
						{
							dp[k][w]+=d[j][k][w];
						}
					}
					
							
					ans=max(ans,work(dp,m,n));
				}
			}
			cout<<ans<<endl;
		}
	}
}