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;
}