poj 2155 Matrix（二维树状数组）

链接：http://poj.org/problem?id=2155

题意：给出一个n*n的矩阵，初始化为0，给出q个操作，C x1 y1 x2 y2表示将（x1,y1)至(x2,y2)的矩阵的数翻转（0变成1,1变成0），Q x y表示求（x，y）位置的点的值。

分析：这是hdu1556的二维版本。把（1，1）到四个角的矩阵内的数全翻转1次即update()，然后sum（x,y)即可。

代码如下：

#include<cstdio>
#include<cstring>
const int N=1010;
int c[N][N],n;
int lowbit(int x)
{
    return x&(-x);
}
void update(int x,int y,int num)
{
    int i,j;
    for(i=x;i<=n;i+=lowbit(i))
        for(j=y;j<=n;j+=lowbit(j))
            c[i][j]^=num;
}
int sum(int x,int y)
{
    int i,j,s=0;
    for(i=x;i>0;i-=lowbit(i))
        for(j=y;j>0;j-=lowbit(j))
            s+=c[i][j];
    return s%2;
}
int main()
{
    int t,q,x1,x2,y1,y2,x,y;
    char s[10];
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&q);
        memset(c,0,sizeof(c));
        while(q--)
        {
            scanf("%s",s);
            if(s[0]=='C')
            {
                scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
                update(x1,y1,1);
                update(x2+1,y1,1);
                update(x1,y2+1,1);
                update(x2+1,y2+1,1);
            //    update(x1-1,y1-1,1);      //超时陷进，树状数组的下标不能为0，所以+1变成上面的代码
            //    update(x2,y1-1,1);
            //    update(x1-1,y2,1);
            //    update(x2,y2,1);
            }
            else
            {
                scanf("%d%d",&x,&y);
                printf("%d\n",sum(x,y));
            }
        }
        printf("\n");
    }
    return 0;
}