思路：从一维扩展到三维。可以看看poj2155的解法。
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#define Maxn 102
#define lowbit(x) (x&(-x))
using namespace std;
int C[Maxn][Maxn][Maxn],n;
int Sum(int i,int j,int k)
{
int sum=;
int y,z;
y=j,z=k;
while(i)
{
j=y;
while(j)
{
k=z;
while(k)
{
sum+=C[i][j][k];
k-=lowbit(k);
}
j-=lowbit(j);
}
i-=lowbit(i);
}
return sum;
}
void update(int i,int j,int k)
{
int y,z;
y=j,z=k;
while(i<=n)
{
j=y;
while(j<=n)
{
k=z;
while(k<=n)
{
C[i][j][k]++;
k+=lowbit(k);
}
j+=lowbit(j);
}
i+=lowbit(i);
}
}
int main()
{
int x1,x2,y1,y2,z1,z2,m,x;
while(scanf("%d%d",&n,&m)!=EOF)
{
memset(C,,sizeof(C));
while(m--)
{
scanf("%d",&x);
if(x==)
{
scanf("%d%d%d",&x1,&y1,&z1);
printf("%d\n",Sum(x1,y1,z1)&);
}
if(x==)
{
scanf("%d%d%d%d%d%d",&x1,&y1,&z1,&x2,&y2,&z2);
update(x1,y1,z1);
update(x1,y2+,z1);
update(x2+,y2+,z1);
update(x2+,y1,z1);
update(x1,y1,z2+);
update(x1,y2+,z2+);
update(x2+,y2+,z2+);
update(x2+,y1,z2+);
}
}
}
return ;
}