HDU 1892 See you~ 【 二维树状数组 】

题意：二维的树状数组
注意的有三个地方，
输入进去的坐标都加1，防止lowbit(0) + 0造成死循环
还有就是询问矩形面积的时候，输入进去的x1,x2,y1,y2，可能不是正对角线，要转化成正对角线

初始化的时候，是每个点的值为1

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include <cmath>
 #include<stack>
 #include<vector>
 #include<map>
 #include<set>
 #include<queue>
 #include<algorithm>
 using namespace std;

 typedef long long LL;
 const int INF = (<<)-;
 const int mod=;
 const int maxn=;

 int a[][],c[][];

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

 int sum(int x,int y){
     int ret=;
     for(int i=x;i>;i-=lowbit(i))
      for(int j=y;j>;j-=lowbit(j))
      ret+=c[i][j];

      return ret;
 }

 void add(int x,int y,int d){
     for(int i=x;i<;i+=lowbit(i))
      for(int j=y;j<;j+=lowbit(j))
      c[i][j]+=d;
 }

 void init(){
     memset(c,,sizeof(c));
     for(int i=;i<;i++){
         for(int j=;j<;j++){
             add(i,j,);
             a[i][j] = ;
         }
     }
 }

 int main(){
     int T;
     int kase = ;
     scanf("%d",&T);
     while(T--){
         printf("Case %d:\n",++kase);
         init();
         int m;
         scanf("%d%*c",&m);
         while(m--){
             char cmd;
             scanf("%c",&cmd);
             if(cmd == 'S') {
                 int x,xx,y,yy,x1,x2,y1,y2;
                 scanf("%d%d%d%d%*c",&x,&y,&xx,&yy);
                 x++;y++;xx++;yy++;
                 x1 = min(x,xx);x2 = max(x,xx);
                 y1 = min(y,yy); y2 = max(y,yy);

                 int ans=;
                 ans += sum(x1-,y1-);
                 ans -= sum(x1-,y2);
                 ans -= sum(x2,y1-);
                 ans += sum(x2,y2);
                 printf("%d\n",ans);
             }

             if(cmd == 'A'){
                 int x1,y1,n1;
                 scanf("%d%d%d%*c",&x1,&y1,&n1);
                 x1++;y1++;
                 a[x1][y1]+=n1;
                 add(x1,y1,n1);
             }
             if(cmd == 'D'){
                 int x1,y1,n1;
                 scanf("%d%d%d%*c",&x1,&y1,&n1);
                 x1++;y1++;
                 n1 = min(n1,a[x1][y1]);
                 a[x1][y1]-=n1;
                 add(x1,y1,-n1);
             }
             if(cmd == 'M'){
                 int x1,y1,x2,y2,n1;
                 scanf("%d%d%d%d%d%*c",&x1,&y1,&x2,&y2,&n1);
                 x1++;y1++;x2++;y2++;
                 n1=min(n1,a[x1][y1]);
                 a[x1][y1]-=n1;a[x2][y2]+=n1;
                 add(x1,y1,-n1);
                 add(x2,y2,n1);
             }
         }
     }
     return ;
 }