poj 2155 二维树状数组

/*
题意： 给定一个n*n的0-1矩阵，执行一些如下的两种操作：
C a,b,c,d ：将矩阵中（a，b）到（c，d）的格子中的0-1进行“非”操作
Q a b ：询问当前（a，b）位置的元素的值
思路： 二维树状数组
本题给出的要求正好与树状数组的操作"相反"，但是还是可以转换成树状数组的操作。
本题的要求是对二维数组的某一个矩形范围进行异或操作，而是询问单个元素的值。
一种思路是用一个二维数组c[][]来标记数组中值的变换次数，通过变换的次数就可以得出当前数组元素的值了。
我们可以先分析一下一维时候的情况，假设一个一维的数组array[] ，
现在要求给区间[a,b]内的每个元素增加一个值add，我们可以这样利用树状数组，
将a元素增加add值，将（b+1）元素增加一个-add值，接下去每次询问k位置的值的时候，
只要将1-k位置的所有元素的值相加就是最后k位置的元素的值（PS：一开始元素的所有值都是0）。
将此思想推广到二维上就可以得出本题的算法。此算法的复杂度为：O( Q*logn*logn) 。
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1002;
int t,n,m,c[maxn][maxn];
char op[2];
int lowbit(int i)
{
    return i&(-i);
}
void add(int x,int y,int d)
{
    int i=x,j=y;
    while(i<=n)
    {
        j=y;
        while(j<=n)
        {
            c[i][j]+=d;
            j+=lowbit(j);
        }
        i+=lowbit(i);
    }
}
int sum(int x,int y)
{
    int i=x,j=y,ret=0;
    while(i)
    {
        j=y;
        while(j)
        {
            ret+=c[i][j];
            j-=lowbit(j);
        }
        i-=lowbit(i);
    }
    return ret;
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        int i,j,x,y,xx,yy;
        memset(c,0,sizeof(c));
        while(m--)
        {
            scanf("%s",&op);
            if(op[0]=='Q')
            {
                scanf("%d%d",&x,&y);
                printf("%d\n",(sum(x,y)&1));
            }else
            {
                scanf("%d%d%d%d",&x,&y,&xx,&yy);
                xx++,yy++;
                add(x,y,1);
                add(x,yy,-1);
                add(xx,y,-1);
                add(xx,yy,1);
            }
        }
        printf("\n");
    }
    return 0;
}
/*
1
2 10
C 2 1 2 2
Q 2 2
C 2 1 2 1
Q 1 1
C 1 1 2 1
C 1 2 1 2
C 1 1 2 2
Q 1 1
C 1 1 2 1
Q 2 1
*/