二维树状数组

题目看这里：http://blog.csdn.net/stone41123/article/details/78848061
然后发一张图，对比一下二维树状数组和四分树：（第一个是二维树状数组）

差的真不是一点半点，差距太大了
首先二维树状数组的写法：

inline void update(int x,int y){
    int tmp=y;
    while(x<=n){
        y=tmp;
        while(y<=m){
            //修改
            y+=y&-y;
        }
        x+=x&-x;
    }
}
inline ll query(int x,int y){
    ll ans=0;int tmp=y;
    while(x){
        y=tmp;
        while(y){
            //累加
            y-=y&-y;
        }
        x-=x&-x;
    }
    return ans;
}

是不是很简单？
然后就是如何用二维树状数组来做区域修改，区域查询
首先建议先用树状数组做一下改段求段，关于公式推导可以查一下
我们来讨论如何用二维树状数组来维护区域和
首先肯定要差分，比如修改$T_{x,y}$就代表修改$(x,y)$到$(n,m)$
然后修改解决了，来讨论如何查询
类似树状数组改段求段，我们要统计每一个值对答案的贡献，也就是它被算了几次
看这个图就差不多懂了：（画的有点儿丑，不要在意，能看就行）

然后把询问拆成四个区域来回答就好了
具体的可以看我的代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#define ll long long
using namespace std;
inline int read(){
    int x=0;char ch=' ';int f=1;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
const int N=805;
int n,m,q;
ll a[N][N],sum[4][N][N];
inline void update(int x,int y,ll val){
    int tmp=y;
    ll val1=val*(ll)(x-1)*(ll)(y-1);
    ll val2=val*(ll)(x-1);
    ll val3=val*(ll)(y-1);
    while(x<=n){
        y=tmp;
        while(y<=m){
            sum[0][x][y]+=val1;
            sum[1][x][y]+=val2;
            sum[2][x][y]+=val3;
            sum[3][x][y]+=val;
            y+=y&-y;
        }
        x+=x&-x;
    }
}
inline ll query(int x,int y,int opt){
    ll ans=0;int tmp=y;
    while(x){
        y=tmp;
        while(y){
            ans+=sum[opt][x][y];
            y-=y&-y;
        }
        x-=x&-x;
    }
    return ans;
}
inline ll Query(int x,int y){
    ll ans=query(x,y,3)*(ll)x*(ll)y+query(x,y,0);
    ans-=query(x,y,1)*(ll)y+query(x,y,2)*(ll)x;
    return ans;
}
int main(){
    n=read();m=read();q=read();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            a[i][j]=read();
            update(i,j,a[i][j]);
            update(i+1,j+1,a[i][j]);
            update(i+1,j,-a[i][j]);
            update(i,j+1,-a[i][j]);
        }
    while(q--){
        int opt=read(),X1=read(),Y1=read(),X2=read(),Y2=read();
        if(opt==1){
            int val=read();
            update(X1,Y1,val);
            update(X2+1,Y2+1,val);
            update(X2+1,Y1,-val);
            update(X1,Y2+1,-val);
        }
        else{
            ll ans=Query(X2,Y2)+Query(X1-1,Y1-1)-Query(X1-1,Y2)-Query(X2,Y1-1);
            printf("%lld\n",ans);
        }
    }
    return 0;
}