Hdu3584_Cube(三维树状数组)

题意：

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

思路：

区间更新，点查询的三维树状数组。记录改变的次数，次数为奇数则最终为1，偶则为0。
更新(x1,y1,z1)—(x2,y2,z2)的立方体，如果只更新(x1,y1,z1)则会把更大区间的也更新了，所以利用容斥原理，再去掉多余的，因为更新是+1，而这里只需判断奇偶，所以再加1就去掉更新效果了。
列出需要更新的点。

#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
typedef long long LL;
const int INF = 0x3f3f3f3f;
const int maxn = 100+5;
using namespace std;
int n,m;
int C[maxn][maxn][maxn];

inline int lowbit(int x){return x & -x; }
void add(int x1, int y1, int z1, int v){
	for(int i = x1; i <= n; i+= lowbit(i)){
		for(int j = y1; j <= n; j+= lowbit(j)){
			for(int z = z1; z <= n; z+= lowbit(z)){
				C[i][j][z]+= v;
			}
		}
	}
}
int sum(int x1, int y1, int z1){
	int s = 0;
	for(int i = x1; i > 0; i-= lowbit(i)){
		for(int j = y1; j > 0; j-=lowbit(j)){
			for(int z = z1; z > 0; z-=lowbit(z)){
				s+= C[i][j][z];
			}
		}
	}
	return s;
}

int main()
{
   // freopen("in.txt","r",stdin);
    int x1,y1,z1,x2,y2,z2, op;
	while(scanf("%d%d",&n,&m) == 2){
		memset(C, 0, sizeof(C));
		for(int i = 0; i < m; ++i){
			scanf("%d", &op);
			//printf("op = %d\n", op);
			if(1 == op){
				scanf("%d%d%d%d%d%d",&x1,&y1,&z1,&x2,&y2,&z2);
				add(x1, y1, z1, 1); 
				add(x1, y1, z2+1, 1); add(x2+1, y1, z1, 1); add(x1, y2+1, z1, 1);
				add(x2+1, y2+1, z1, 1); add(x2+1, y1, z2+1, 1);add(x1, y2+1, z2+1, 1); 
				add(x2+1, y2+1, z2+1, 1);
			}
			else{
				scanf("%d%d%d",&x1,&y1,&z1);
				if(sum(x1,y1,z1)&1) printf("1\n");
				else printf("0\n");
			}
		} 
	}
    fclose(stdin);
	return 0;
}