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