UVA 10755 Garbage Heap 最大子立方体和

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define ll long long
#define inf 1ll<<60//加ll
using namespace std;
ll sum[27][27][27];
ll getsum(int x1,int x2,int y1,int y2,int z)//求从x1到x2，y1到y2在z轴上的数字之和
{
return sum[x2][y2][z]-sum[x1-1][y2][z]-sum[x2][y1-1][z]+sum[x1-1][y1-1][z];
}
int main()
{
int t,cas=1,flag;
scanf("%d",&t);
flag=t;
while(t--)
{
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
memset(sum,0,sizeof(sum));
for(int i=1;i<=a;i++)
for(int j=1;j<=b;j++)
for(int k=1;k<=c;k++)
{
scanf("%lld",&sum[i][j][k]);
sum[i][j][k]+=sum[i-1][j][k]+sum[i][j-1][k]-sum[i-1][j-1][k];//预处理得到在k这个平面上前矩阵（i，j）的和
}
ll maxx=-inf;
for(int sx=1;sx<=a;sx++)//x方向开始位置
for(int ex=sx;ex<=a;ex++)//x方向结束位置
for(int sy=1;sy<=b;sy++)//y方向开始位置
for(int ey=sy;ey<=b;ey++)//y方向结束位置
{
ll ans=0;
for(int sz=1,ez=1;ez<=c;ez++)//遍历z轴
{
ans+=getsum(sx,ex,sy,ey,ez);
maxx=max(maxx,ans);
while(ans<0&&sz<=ez)
{
ans-=getsum(sx,ex,sy,ey,sz);
sz++;
if(sz<=ez&&ans>maxx)maxx=ans;
}
}
}
printf("%lld\n",maxx);
if(cas!=flag)printf("\n");
cas++;
}
}