TYVJ 1476 三维树状数组

裸题~

注意三维的前缀体积怎么算就好

 

#include <cstdio>
#include <cstring>
#include <cstdlib>

#define N 120

using namespace std;

int n,c[N][N][N];

inline int lowbit(int x)
{
    return x&-x;
}

void updata(int x,int y,int z,int dt)
{
    while(x<=n)
    {
        int y1=y;
        while(y1<=n)
        {
            int z1=z;
            while(z1<=n)
            {
                c[x][y1][z1]+=dt;
                z1+=lowbit(z1);
            }
            y1+=lowbit(y1);
        }
        x+=lowbit(x);
    }
}

int getsum(int x,int y,int z)
{
    int rt=0;
    while(x)
    {
        int y1=y;
        while(y1)
        {
            int z1=z;
            while(z1)
            {
                rt+=c[x][y1][z1];
                z1-=lowbit(z1);
            }
            y1-=lowbit(y1);
        }
        x-=lowbit(x);
    }
    return rt;
}

void read()
{
    scanf("%d",&n);
    int pd,x,y,z,x1,y1,z1,p,s;
    while(scanf("%d",&pd))
    {
        if(pd==3) break;
        if(pd==1)
        {
            scanf("%d%d%d%d",&x,&y,&z,&p);
            x++; y++; z++;
            updata(x,y,z,p);
        }
        else
        {
            scanf("%d%d%d%d%d%d",&x1,&y1,&z1,&x,&y,&z);
            x++; y++; z++;
            s=getsum(x,y,z)-getsum(x,y,z1)-getsum(x,y1,z)-getsum(x1,y,z);
            s=s+getsum(x1,y1,z)+getsum(x1,y,z1)+getsum(x,y1,z1);
            s=s-getsum(x1,y1,z1);
            printf("%d\n",s);
        }
    }
}

int main()
{
    read();
    return 0;
}