HDU 3584 Cube(三维树状数组)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3584

题意：三维的空间中有两个操作，初始时每个空间元素均为0，更新操作是0变1，1变0，是一个空间内的所有元素都更新，查询是问这个点的元素是0还是1。

题解：三维树状数组的IUPQ模型(Interval Update Point Query段更新，点求值)。

关键在更新的时候的边界处理，好像是容斥原理，并不懂，先记住吧。

二维的情况：

add(x1, y1, v);add(x1,y2+1, -v);add(x2+1,y1,-v);add(x2+1,y2+1,v);

三维的情况：

Add(x1,y1,z1,1);Add(x1,y1,z2+1,-1);Add(x1,y2+1,z1,-1);Add(x2+1,y1,z1,-1);Add(x2+1,y2+1,z1,1);Add(x2+1,y1,z2+1,1);Add(x1,y2+1,z2+1,1);<p>Add(x2+1,y2+1,z2+1,-1);</p>

这道题因为只需转变0和1，所以一直加1，求mod2就行.

代码：

#include <iostream>#include <cstdio>#include <algorithm>#include <vector>#include <cstring>using namespace std;const int MAX=100+5;typedef long long LL;int bit[MAX][MAX][MAX];int N;int sum(int x,int y,int z){    int s=0;    for(int i=x; i>0; i-=(i&(-i)))        for(int j=y; j>0; j-=(j&(-j)))            for(int k=z; k>0; k-=(k&(-k)))                s+=bit[i][j][k];    return s;}void add(int x,int y,int z){    for(int i=x; i<=N; i+=(i&(-i)))        for(int j=y; j<=N; j+=(j&(-j)))            for(int k=z; k<=N; k+=(k&(-k)))                bit[i][j][k]++;}int main(){    int n,m;    while(scanf("%d%d",&n,&m)==2)    {        N=n;        memset(bit,0,sizeof(bit));        for(int i=0; i<m; i++)        {            int op;            scanf("%d",&op);            if(op==1)            {                int x1,y1,z1,x2,y2,z2;                scanf("%d%d%d%d%d%d",&x1,&y1,&z1,&x2,&y2,&z2);                add(x1,y1,z1);                add(x1,y1,z2+1);                add(x1,y2+1,z1);                add(x1,y2+1,z2+1);                add(x2+1,y1,z1);                add(x2+1,y1,z2+1);                add(x2+1,y2+1,z1);                add(x2+1,y2+1,z2+1);            }            else if(op==0)            {                int x,y,z;                scanf("%d%d%d",&x,&y,&z);                printf("%d\n",sum(x,y,z)%2);            }        }    }    return 0;}

关键在更新的时候的边界处理，好像是容斥原理，并不懂，先记住吧。