[poj 2155] Matrix（二维树状数组）

假设原数组设为num， 树状数组设为c。
一维的树状数组操作：
    区间[x, y]取反，那么就在原数组上num[x]++，num[y+1]++，即update(x, 1), update(y+1, 1)。
    求num[x]的状态，就re = getsum(x)，判断re%2的值即可。
推广到二维也是一样的：
    矩阵[x1, y1] - [x2, y2]取反：就在num[x1, y1]、num[x2+1, y2+1]、num[x1, y2+1]、num[x2+1, y1]上各自+1。
    求[x, y]的状态：re = getsum(x, y)，判断re%2的值。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

#define maxn 1010
int c[maxn][maxn];
int n, m;

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

void update(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] ++;
		}
	}
}

int getsum(int x, int y)
{
	int re = 0;
	for(int i = x; i > 0; i -= lowbit(i))
	{
		for(int j = y; j > 0; j -= lowbit(j))
		{
			re += c[i][j];
		}
	}
	return re;
}

void init(int n)
{
	for(int i = 0; i <= n+1; i++)
	{
		for(int j = 0; j <= n+1; j++)
		{
			c[i][j] = 0;
		}
	}
}

int main()
{
	int tot;
	scanf("%d", &tot);
	char op[5];
	int a, b, c, d;
	while(tot--)
	{
		scanf("%d%d", &n, &m);
		init(n);
		while(m--)
		{
			scanf("%s%d%d", op, &a, &b);
			if(op[0] == 'C')
			{
				scanf("%d%d", &c, &d);
				update(a, b);
				update(c+1, d+1);
				update(a, d+1);
				update(c+1, b);
			}
			else if(op[0] == 'Q')
			{
				int re = getsum(a, b);
				printf("%d\n", re&1);
			}
		}
		putchar(10);
	}
	return 0;
}