uva10755(容斥）

蓝书P56

这个题目求前缀和的时候，利用了容斥的思想；

后来求子长方体的时候也用了容斥的思想，并且用了一个数来保存当前出现过的同x同y的子长方体的体积，这样减少了一个数量级的运算；

;

#include <iostream>
#include<algorithm>
#include<cstdio>
#include<string>
#include<cmath>
#include<vector>
#include<cstdio>
#include<cstring>
#include<map>

using namespace std;
const int maxn=25;
const int maxm=505;
const long long  inf= 1LL<<60;
int A,B,C;

long long s[maxn][maxn][maxn];
void get(int i,int &a,int &b,int &c) //二进制枚举
{
    a=(i&1);  i>>=1;
    b=(i&1);  i>>=1;
    c=(i&1);
}
int chara(int a,int b,int c) //奇加偶减
{
    return (a+b+c)%2==1?1:(-1);
}
long long sum(int x1,int x2,int y1,int y2,int z1,int z2)
{
    int dx=x2-x1+1,dy=y2-y1+1,dz=z2-z1+1;
    int a,b,c;
    long long res=0;
    for(int i=0;i<8;i++)
    {
        get(i,a,b,c);
        res-=s[x2-a*dx][y2-b*dy][z2-c*dz]*chara(a,b,c);
    }
    return res;
}
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        memset(s,0,sizeof(s));
        cin>>A>>B>>C;
        for(int i=1;i<=A;i++)
            for(int j=1;j<=B;j++)
                for(int k=1;k<=C;k++)
                   cin>>s[i][j][k];

        int a,b,c;
        for(int x=1;x<=A;x++)
        {
            for(int y=1;y<=B;y++)
            {
                for(int z=1;z<=C;z++)
                {
                    for(int i=1;i<(1<<3);i++)
                    {
                        get(i,a,b,c);
                        s[x][y][z]+=(s[x-a][y-b][z-c])*chara(a,b,c);
                    }
                }
            }
        }

        long long ans=-inf;
        for(int x1=1;x1<=A;x1++)
          for(int x2=x1;x2<=A;x2++)
            for(int y1=1;y1<=B;y1++)
              for(int y2=y1;y2<=B;y2++)
            {
               long long M=0;
               for(int z=1;z<=C;z++)
               {
                  long long m=sum(x1,x2,y1,y2,1,z);
                  ans=max(ans,m-M);
                  M=min(M,m);
               }
           }
        cout<<ans<<endl;
        if(T)
            cout<<endl;
    }
    return 0;
}