POJ2155 Matrix

　　原题传送：http://poj.org/problem?id=2155

　　楼教神题，二维树状数组。

　　弄懂一维树状数组，二维的也不难懂。可以参考NOCOW。

　　需要注意的是，修改的四个点分别是(x1, y1), (x2+1, y1), (x1, y2+1), (x2+1, y2+1)，而不是(x1, y1), (x2, y1), (x1, y2), (x2, y2)，因为如果是后者的话，(x2, y2)会被4次覆盖（也就是点(x2, y2)不会被改变）。

#include <stdio.h>
#include <string.h>
#define N 1001
typedef long long LL;

LL c[N][N];

LL lowbit(LL x)
{
    return x & (-x);
}

void modify(int x, int y)
{
    for(int i = x; i < N; i += lowbit(i))
    for(int j = y; j < N; j += lowbit(j))
        c[i][j] ++;
    return ;
}

LL query(int x, int y)
{
    LL 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;
}

int main()
{
    int cas, n, tt, t, x1, y1, x2, y2;
    char op[3];
    scanf("%d", &tt);
    for(cas = 1; cas <= tt; cas ++)
    {
        if(cas != 1)    
            putchar('\n');
        memset(c, 0, sizeof c);
        scanf("%d%d", &n, &t);
        while(t --)
        {
            scanf("%s", op);
            if(op[0] == 'C')
            {
                scanf("%d%d%d%d", &x1, &y1, &x2, &y2);
                modify(x1, y1);
                modify(x2 + 1, y1);
                modify(x1, y2 + 1);
                modify(x2 + 1, y2 + 1);
            }
            else if(op[0] == 'Q')
            {
                scanf("%d%d", &x1, &y1);
                printf("%I64d\n", query(x1, y1) & 1);
            }
        }
    }
    return 0;
}