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;

}