POJ（1195）（单点修改，区间查询）（二维）

题目大意

给定一个N*N的网格，刚开始每个网格的值都是0，接下来会对这些网格进行操作，有一下两种操作：

1、”X Y A“对网格C[x][y]增加A

2、”L B R T“ 查询所有(L<=X<=R,B<=Y<=T)的网格]，并返回它们的总和

如果是针对于单点修改，那是比区间修改好打很多，因为这个只需要将区间拆分成许多个树状数组就可以了，

然后修改的时候将树状数组修改，十分好理解的。

查询的话就是二维前缀和的查询，然后就可以了。

 #include<cstdio>
 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cmath>
 #define N 1057
 using namespace std;

 int n,k;
 int a[N][N];

 int lowbit(int x)
 {
     return x&(-x);
 }
 void add(int x,int y,int z)
 {
     for (int i=x;i<=n;i+=lowbit(i))
         for (int j=y;j<=n;j+=lowbit(j))
             a[i][j]+=z;
 }
 int query(int x,int y)
 {
     int res=;
     for (int i=x;i>=;i-=lowbit(i))
         for (int j=y;j>=;j-=lowbit(j))
             res+=a[i][j];
     return res;
 }
 int main()
 {
     while(scanf("%d",&k))
     {
         if (k==)
         {
             scanf("%d",&n);
             for (int i=;i<=n;i++)
                 for (int j=;j<=n;j++)
                     a[i][j]=;
         }
         if (k==)
         {
             int x,y,z;
             scanf("%d%d%d",&x,&y,&z);
             x++,y++;
             add(x,y,z);
         }
         if (k==)
         {
             int x1,x2,y1,y2;
             scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
             x1++,x2++,y1++,y2++;
             cout<<query(x2,y2)+query(x1-,y1-)-query(x1-,y2)-query(x2,y1-)<<endl;
         }
         if (k==) break;
     }
 }