HDU3584 Cube

题目大意

给定一个体积为N*N*N立方体，每个单位小立方体A[x][y][z]里有一个值,初始值全部为0，我们可以对立方体进行一下两种操作：

1、查询A[x][y][z]的奇偶性

2、对子立方体的每个元素的值进行增减

题解

三维树状数组，就是POJ2155 Matrix的三维推广形式

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#define MAXN 105
using namespace std;
int c[MAXN][MAXN][MAXN];
int n;
int lowbit(int x)
{
    return x&-x;
}
void add(int x,int y,int z)
{
    int i,j,k;
    for(i=x;i<MAXN;i+=lowbit(i))
    for(j=y;j<MAXN;j+=lowbit(j))
    for(k=z;k<MAXN;k+=lowbit(k))
    c[i][j][k]++;
}
int sum(int x,int y,int z)
{
    int i,j,k,ret=0;
    for(i=x;i>0;i-=lowbit(i))
    for(j=y;j>0;j-=lowbit(j))
    for(k=z;k>0;k-=lowbit(k))
    ret+=c[i][j][k];
    return ret;
}
int main()
{
    int x1,x2,y1,y2,z1,z2,i,m,t;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        memset(c,0,sizeof(c));
        for(i=1;i<=m;i++)
        {
            scanf("%d",&t);
            if(t)
            {
                scanf("%d%d%d%d%d%d",&x1,&y1,&z1,&x2,&y2,&z2);
                add(x1,y1,z1);
                add(x2+1,y1,z1);
                add(x1,y2+1,z1);
                add(x1,y1,z2+1);
                add(x2+1,y2+1,z1);
                add(x2+1,y1,z2+1);
                add(x1,y2+1,z2+1);
                add(x2+1,y2+1,z2+1);
            }
            else
            {
                scanf("%d%d%d",&x1,&y1,&z1);
                printf("%d\n",sum(x1,y1,z1)&1);
            }
        }
    }
    return 0;
}