HDU3584 Cube

分析

三维树状数组，跟二维的差不多，只不过容斥的不一样。
更新区间\((x_1,y_1,z_1)\rightarrow(x_2,y_2,z_2)\)时，需要利用容斥原理。须要更新的节点更新如下：
\[ (x_1,y_1,z_1),\\ (x_2+1,y_1,z_1), (x_1,y_2+1,z_1), (x_1,y_1,z_2+1),\\ (x_1,y_2+1,z_2+1), (x_2+1,y_1,z_2+1), (x_2+1,y_2+1,z_1),\\ (x_2+1,y_2+1,z_2+1) \]
更高维的分类讨论就不容易了，可以找规律，利用容斥原理。其实有一点组合的味道。

代码

#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<ctime>
#include<iostream>
#include<string>
#include<vector>
#include<list>
#include<deque>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<complex>
#pragma GCC optimize ("O0")
using namespace std;
template<class T> inline T read(T&x)
{
    T data=0;
    int w=1;
    char ch=getchar();
    while(!isdigit(ch))
    {
        if(ch=='-')
            w=-1;
        ch=getchar();
    }
    while(isdigit(ch))
        data=10*data+ch-'0',ch=getchar();
    return x=data*w;
}
typedef long long ll;
const int INF=0x7fffffff;

const int MAXN=110;
int a[MAXN][MAXN][MAXN],n;

int lowbit(int x)
{
    return x&-x;
}

void add(int x,int y,int z,int v)
{
    for(int i=x;i<=n;i+=lowbit(i))
        for(int j=y;j<=n;j+=lowbit(j))
            for(int k=z;k<=n;k+=lowbit(k))
                a[i][j][k]+=v;
}

int sum(int x,int y,int z)
{
    int res=0;
    for(int i=x;i;i-=lowbit(i))
        for(int j=y;j;j-=lowbit(j))
            for(int k=z;k;k-=lowbit(k))
                res+=a[i][j][k];
    return res;
}

int main()
{
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);
    int m;
    while(~scanf("%d %d",&n,&m))
    {
        memset(a,0,sizeof(a));
        while(m--)
        {
            int opt;
            read(opt);
            if(opt==1)
            {
                int x1,y1,z1,x2,y2,z2;
                read(x1);read(y1);read(z1);read(x2);read(y2);read(z2);
                add(x1,y1,z1,1);
                add(x2+1,y1,z1,1);
                add(x1,y2+1,z1,1);
                add(x1,y1,z2+1,1);
                add(x1,y2+1,z2+1,1);
                add(x2+1,y1,z2+1,1);
                add(x2+1,y2+1,z1,1);
                add(x2+1,y2+1,z2+1,1);
            }
            else
            {
                int x,y,z;
                read(x);read(y);read(z);
                printf("%d\n",sum(x,y,z)%2);
            }
        }
    }
//  fclose(stdin);
//  fclose(stdout);
    return 0;
}