poj 1195 Mobile phones（树状数组）

刚接触树状数组，感觉是个很简练的数据结构，虽然我还没怎么明白它的原理。本来想用浙大模版里的树状数组，不过看半天不知道代码什么意思，而且求和的时候坐标向后错了一位，只能在网上找一个了，速度比浙大模版要快一点，不过浙大模版的简练程度让我真心佩服……

#include <cstdio>
#include <iostream>
#include <string.h>
using namespace std;
#define lowbit(x) ((x)&(-x))
#define MAXN 1050
typedef int elem_t;
int n,m,t;
struct sum
{
    elem_t a[MAXN][MAXN],c[MAXN][MAXN]ret;
    int m,n,t;
    void init(int i,int j)
    {
        memset(a,0,sizeof(a));
        memset(c,0,sizeof(c));
        m=i,n=j;
    }
    //浙大模版，简练的半天没看懂……
    /*void update(int i,int j,elem_t v){
		for (v-=a[i][j],a[i++][j++]+=v,t=j;i<=m;i+=lowbit(i))
			for (j=t;j<=n;c[i-1][j-1]+=v,j+=lowbit(j));
	}
	elem_t query(int i,int j){
		for (ret=0,t=j;i;i^=lowbit(i))
			for (j=t;j;ret+=c[i-1][j-1],j^=lowbit(j));
		return ret;
	}*/
    void modify(int i, int j, int w)
    {
        while(i<=m)
        {
            t=j;
            while(t<=n)
            {
                a[i][t]+=w;
                t+=lowbit(t);
            }
            i+=lowbit(i);
        }
    }
    elem_t query(int i,int j)
    {
        int ans=0;
        while(i>0)
        {
            t=j;
            while(t>0)
            {
                ans+=a[i][t];
                t-=lowbit(t);
            }
            i-=lowbit(i);
        }
        return ans;
    }
}tree;
int main()
{
    //freopen("input.txt","r",stdin);
    int op;
    while(scanf("%d",&op)!=EOF)
    {
        if(op==0)
        {
            scanf("%d",&n);
            tree.init(n,n);
        }
        else if(op==1)
        {
            int x,y,v;
            scanf("%d%d%d",&x,&y,&v);
            modify(x+1,y+1,v);
        }
        else if(op==2)
        {
            int l,b,r,t;
            scanf("%d%d%d%d",&l,&b,&r,&t);
            printf("%d\n",query(r+1,t+1)-query(r+1,b)-query(l,t+1)+query(l,b));
        }
        else continue;
    }
}