[HDU3584] Cube - 三维树状数组

---

题目描述

给定一个体积为N*N*N立方体，每个单位小立方体A[x,y,z]里有一个值，初始值全部为0，我们可以对立方体进行一下两种操作：
操作“Not”：改变A[i,j,k]=!A[i,j,k]。意思是改变A[i,j,k]的值，从0->1或者1->0。 (x1<=i<=x2,y1<=j<=y2,z1<=k<=z2)。
操作“Query”：询问A[i,j,k]的值。

---

输入格式

多组测试数据。对于每组测试数据：
第一行包含两个整数N和M，接下来M行，每行首先一个整数X，若X=1表示“Not”操作；X=0表示“Query”操作；如果X=1，接下来x1, y1, z1, x2, y2, z2；如果X=0，接下来x,y,z。

---

输出格式

对于每个“Query”操作输出一行。(1<=n<=100,m<=10000)

---

样例数据

样例输入

2 5
1 1 1 1 1 1 1
0 1 1 1
1 1 1 1 2 2 2
0 1 1 1
0 2 2 2

样例输出

1
0
1

---

题目分析

三维树状数组
区间修改点查询

---

源代码

#include<algorithm>#include<iostream>#include<iomanip>#include<cstring>#include<cstdlib>#include<vector>#include<cstdio>#include<cmath>#include<queue>using namespace std;inline const int Get_Int() {    int n=0,bj=1;    char x=getchar();    while(x<'0'||x>'9') {        if(x=='-')bj=-1;        x=getchar();    }    while(x>='0'&&x<='9') {        n=n*10+x-'0';        x=getchar();    }    return n*bj;}const int maxn=205;struct BIT { //树状数组    int n,c[maxn][maxn][maxn];    inline int Lowbit(int x) { //低位操作        return x&(-x);    }    void init(int n) {        this->n=n;        memset(c,0,sizeof(c));    }    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))                    c[i][j][k]+=v;    }    int sum(int x,int y,int z) {        int s=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))                    s+=c[i][j][k];        return s;    }} bit;int n,m;int main() {    while(scanf("%d%d",&n,&m)!=EOF) {        bit.init(n);        for(int i=1; i<=m; i++) {            int order=Get_Int();            if(order==1) {                int x1=Get_Int(),y1=Get_Int(),z1=Get_Int(),x2=Get_Int(),y2=Get_Int(),z2=Get_Int();                bit.add(x1,y1,z1,1);                bit.add(x2+1,y1,z1,-1);                bit.add(x1,y2+1,z1,-1);                bit.add(x1,y1,z2+1,-1);                bit.add(x2+1,y2+1,z1,1);                bit.add(x1,y2+1,z2+1,1);                bit.add(x2+1,y1,z2+1,1);                bit.add(x2+1,y2+1,z2+1,-1);            } else {                int x=Get_Int(),y=Get_Int(),z=Get_Int();                printf("%d\n",bit.sum(x,y,z)%2);            }        }    }    return 0;}

---