bzoj3132上帝造题的七分钟 二维树状数组

具体的话，这是一个实用的技巧，通过差分来做到区间加和区间查询的作用，好像二维线段树或者树套树会挂。
http://blog.csdn.net/qq_35866453/article/details/77199794

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=2105;
struct bit
{
    int a[N][N];
    inline int lowbit(int x)
    {
        return x&(-x);
    }
    inline void add(int x,int y,int t)
    {
        int i,j;
        for(i=x;i<N;i+=lowbit(i))
        {
            for(j=y;j<N;j+=lowbit(j))a[i][j]+=t;
        }
    }
    inline int get(int x,int y)
    {
        int ans=0;
        int i,j;
        for(i=x;i>0;i-=lowbit(i))
        {
            for(j=y;j>0;j-=lowbit(j))ans+=a[i][j];
        }
        return ans;
    }
}a,b,c,d;

int n,m;
inline void add(int x1,int y1,int x2,int y2,int t)
{
    a.add(x1,y1,t),a.add(x1,y2+1,-t);
    a.add(x2+1,y1,-t),a.add(x2+1,y2+1,t);

    b.add(x1,y1,t*x1); b.add(x2+1,y1,-t*(x2+1));
    b.add(x1,y2+1,-t*x1); b.add(x2+1,y2+1,t*(x2+1));

    c.add(x1,y1,t*y1); c.add(x2+1,y1,-t*y1);
    c.add(x1,y2+1,-t*(y2+1)); c.add(x2+1,y2+1,t*(y2+1));

    d.add(x1,y1,t*x1*y1); d.add(x2+1,y1,-t*(x2+1)*y1);
    d.add(x1,y2+1,-t*x1*(y2+1)); d.add(x2+1,y2+1,t*(x2+1)*(y2+1));
}
inline int get(int x,int y)
{
    return a.get(x,y)*(x+1)*(y+1)-b.get(x,y)*(y+1)-(x+1)*c.get(x,y)+d.get(x,y);
}
inline int get(int x1,int y1,int x2,int y2)
{
    return get(x2,y2)-get(x2,y1-1)-get(x1-1,y2)+get(x1-1,y1-1);
}
int main()
{
    char op[2];
    scanf("%s",op);
    scanf("%d%d",&n,&m);
    while (scanf("%s",op)!=-1)
    {
        int x1,y1,x2,y2,k;
        if (op[0]=='L')
        {
            scanf("%d%d",&x1,&y1);
            scanf("%d%d%d",&x2,&y2,&k);
            add(x1,y1,x2,y2,k);
        }
        else
        {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            printf("%d\n",get(x1,y1,x2,y2));
        }
    }
}