poj1195（二维树状数组）

题目链接：https://vjudge.net/problem/POJ-1195

题意：有s*s的矩阵，初始化为全0，有两种操作，单点修改(x,y)的值，区间查询(x,y)的值（l<=x<=r,b<=y<=t）。

思路：二维树状数组裸应用查询区间(l,b)~(r,t)的值可转换为tr[r][t]-tr[l-1][t]-tr[r][b-1]+tr[l-1][b-1]。要注意的是输入的x,y是从0开始的，所以要加1。

AC代码：　　

#include<cstdio>
#include<cctype>
using namespace std;

inline int read(){
    int x=,f=;char c=;
    while(!isdigit(c)){f|=c=='-';c=getchar();}
    while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
    return f?-x:x;
}

int op,s,tr[][];

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

void update(int x,int y,int k){
    while(x<=s){
        int tmp=y;
        while(tmp<=s){
            tr[x][tmp]+=k;
            tmp+=lowbit(tmp);
        }
        x+=lowbit(x);
    }
}

int query(int x,int y){
    int ans=;
    while(x>){
        int tmp=y;
        while(tmp>){
            ans+=tr[x][tmp];
            tmp-=lowbit(tmp);
        }
        x-=lowbit(x);
    }
    return ans;
}

int main(){
    op=read(),s=read();
    while(op=read(),op!=){
        if(op==){
            int x=read()+,y=read()+,k=read();
            update(x,y,k);
        }
        else{
            int l=read()+,b=read()+,r=read()+,t=read()+;
            printf("%d\n",query(r,t)-query(l-,t)-query(r,b-)+query(l-,b-));
        }
    }
    return ;
}