8601 最大长方体问题（优先做）

8601 最大长方体问题（优先做）
时间限制:800MS 代码长度限制:10KB
提交次数:950 通过次数:383

题型: 编程题 语言: G++;GCC;VC;JAVA
Description
一个长,宽,高分别是m,n,p的长方体被分割成mnp个小立方体。每个小立方体内含一个整数。

试着设计一个算法，计算所给长方体的最大子长方体。子长方体的大小由它内部所含所有整数之和确定。

约定：当该长方体所有元素均为负数时，输出的最大子长方体为0。

输入格式
第一行3个正整数m,n,p，其中 1<=m,n,p<=50

接下来的m*n行中每行p个整数，表示小立方体中的数。

输出格式
第一行中的数是计算出的最大子长方体的大小。

输入样例
3 3 3

0 -1 2

1 2 2

1 1 -2

-2 -1 -1

-3 3 -2

-2 -3 1

-2 3 3

0 1 3

2 1 -3

输出样例
14

提示

---

这里，很多同学不会动态申请数组，还和静态定义数组一样写成如下：
int n;
n=10000; //或cin>>n;
int c[n];

最好不要这样写，可能会出错的啊，且不报错，可以运行，但结果可能会不对。
其他的题目如果你想“动态申请数组空间”都同理如下这样来表达。

//一维数组动态申请，c数组大小为： n
int *c=new int[n];

//二维数组动态申请，b数组大小为： n*p
int *b=new int[n];
for(int i=0;i<n;i++)
b[i]=new int[p];

//三维数组动态申请, a数组大小为： mnp:
int ***a=new int **[m];
for(i=0;i<m;i++)
{
a[i]=new int *[n];
for(j=0;j<n;j++)
a[i][j]=new int [p];
}

另外，当不再需要一个动态分配的多维数组时，可按以下步骤释放它所占用的空间。首先释放在for循环中为每一行所分
配的空间。然后释放为行指针分配的空间。

//一维空间释放：
delete []c;
c=0; //可在释放空间后将c置为0，以防继续访问已被释放的空间。这句可以不要。

//二维空间释放：
for (int i=0;i<n;i++)
delete []b[i];
delete []b;
b=0; //可在释放空间后将b置为0，以防继续访问已被释放的空间。这句可以不要。

//三维空间释放：
for (int i=0;i<m;i++)
{
for(j=0;j<n;j++)
delete []a[i][j];
delete []a[i];
}
delete []a;
a=0; //可在释放空间后将a置为0，以防继续访问已被释放的空间。这句可以不要。

数据输入格式为：

for (i=0;i<m;i++)
for (j=0;j<n;j++)
for (k=0;k<p;k++)
cin >> a[i][j][k];

此题思路：
1，先编写一维的“最大字段和”的解法。
2，基于“最大字段和”，编写二维的“最大子矩阵和”的解法。
3，基于“最大子矩阵和”，编写三维的“最大子长方体和”的解法。
三维的最大子长方体和求解过程：对原三维长方体，枚举在某个方向上切各种厚薄的片（如切吐司面包一样，切一段），
这个各种厚薄的片，得用两重循环来枚举（0<= i

// 8601 最大长方体问题.cpp : 定义控制台应用程序的入口点。
//

#include <iostream>
using namespace std;

//长方体：m行,n列,p层
//[m+1,n+1,p+1]


int MaxSum(int n,int *a){
	int sum=0;
	int b=0;
	for(int i=1;i<=n;i++){
		if(b>0)b+=a[i];
		else b=a[i];
		if(sum<b)sum=b;
	}
	return sum;
}

int MaxSum2(int m,int n, int **a){
	int sum=0;
	int *b = new int[n+1];
	for(int i=1;i<=m;i++){
		for(int k=1;k<=n;k++){
			b[k]=0;
		}

		for(int j=i;j<=m;j++){
			for(int k=1;k<=n;k++){
				b[k] += a[j][k];
			}
			int max=MaxSum(n,b);
			if(max>sum)sum=max;
		}
	
	}
	return sum;
}

int MaxSum3(int m, int n, int p, int ***a){
	int sum=0;
	int **b = new int*[n+1];	//动态分配数组[m+1,n+1]
	for(int i=1;i<=n;i++){
		b[i]=new int[p+1];
	}

	for(int i=1;i<=m;i++){//从第1层开始的所有平面  第2层、第3层...

		for(int j=1;j<=n;j++){
			for(int k=1;k<=p;k++){
				b[j][k] =0;  
			}
		}
		for(int t=i;t<=m;t++){//第1层，第1层到第2层，第1层到第3层...  即第i层开始到第t层

			for(int j=1;j<=n;j++){
				for(int k=1;k<=p;k++){
					b[j][k] += a[t][j][k];  //3维压缩为2维,即相加
				}
			}
			int max = MaxSum2(n,p,b);  //粗心 把这2行代码放到for循环外面了
			if(max>sum)sum=max;     //粗心 把这2行代码放到for循环外面了
									//放里面才是正确的
		}



	}
	return sum;

}



int main()
{

	int m,n,p;
	cin>>m>>n>>p;
	int ***a = new int **[m+1];
	for(int i=0;i<m+1;i++){
		a[i] = new int *[n+1];
		for(int j=0;j<n+1;j++){
			a[i][j] = new int [p+1];
		}
	}

	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			for(int k=1;k<=p;k++){
				cin>>a[i][j][k];
			}
		}
	}
	cout<<MaxSum3(m,n,p,a)<<endl;

	return 0;
}