题目描述
给定一个体积为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;
}