POJ-2155 Matrix 二维树状数组

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#define maxn 1005
using namespace std;
int c[maxn][maxn],n;
int bit(int n)
{
	return n&(-n);//二进制 
}
int sum(int n,int m)
{
	int ans=0;
	for(int i=n;i>0;i-=bit(i))
	{
		for(int j=m;j>0;j-=bit(j))
		{
			ans+=c[i][j];
		}
	}
	return ans;//计算区间和 
}
int up(int x,int y,int k)
{
	for(int i=x;i<=n;i+=bit(i))
	{
		for(int j=y;j<=n;j+=bit(j))
		{
			c[i][j]+=k;//更新区间值 
		}
	}
}
int main()
{
	int t;
	int x1,y1,x2,y2;
	char s;
	int a;
	cin>>t;
	while(t--)
	{
		memset(c,0,sizeof(c));
		scanf("%d%d",&n,&a);
		getchar();
		for(int i=0;i<a;i++)
		{
			scanf("%c",&s);
			if(s=='C')
			{
				scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
				getchar();
				up(x1,y1,1); 
				up(x1,y2+1,1);
				up(x2+1,y1,1);
				up(x2+1,y2+1,1);//更新矩阵  +1  结果mod 2 
			} 
			else
			{
				scanf("%d%d",&x1,&y1);
				getchar();
				cout<<sum(x1,y1)%2<<endl;//mod 2 用来翻转 
			}
		}
		cout<<endl;
	}
	return 0;
}