POJ - 2155 Matrix 二维树状数组 （区间更新，单点查询）

POJ - 2155 Matrix 二维树状数组

没学过的可以看这篇博客：http://www.cnblogs.com/RabbitHu/p/BIT.html

矩阵中每个元素初始化为 0 ，每次更新，相应矩阵中的 0 变为 1 ，1 变为 0 。 更新时只需加 1 ，最后对 2 取模即为答案。 

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
int c[1010][1010];  //差分数组
int 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))
            c[i][j] += z;
}

void range_add(int x1,int y1,int x2,int y2,int z)
{
    add(x1,y1,z);
    add(x1,y2+1,-z);
    add(x2+1,y1,-z);
    add(x2+1,y2+1,z);
}

int ask(int x,int y)
{
    int ans=0;
    for(int i=x;i>0;i -= lowbit(i))
        for(int j=y;j>0;j -= lowbit(j))
            ans += c[i][j];
    return ans%2;
}

int main()
{
    int t,m;
    int x,y,x1,y1;
    char s;
    scanf("%d",&t);
    while(t--)
    {
        memset(c,0,sizeof(c));
        scanf("%d%d",&n,&m);
        for(int i=1;i<=m;i++)
        {
            scanf(" %c",&s);
            if(s=='C')
            {
                scanf("%d%d%d%d",&x,&y,&x1,&y1);
                range_add(x,y,x1,y1,1);
            }
            else
            {
                scanf("%d%d",&x,&y);
                printf("%d\n",ask(x,y));
            }
        }
        printf("\n");
    }
    return 0;
}