【poj2155】Matrix 树状数组

题目大意

给定一个N*N的矩阵a，a矩阵的初值设为0。T次操作，涉及区间翻转0和1以及单点查询。
2 <= N <= 1000, 1 <= T <= 50000

题解

我的做法好像跟其他人不太一样qwq
用二维树状数组维护差分，每次(x1,y1)+1,(x2+1,y1)-1,(x1,y2+1)，因为右下部分多减了了一次，所以(x2+1,y2+1)+1，然后询问的时候只要加起来然后看是偶数还是奇数即可。

代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cctype>
using namespace std;
int c[1010][1010];
int T,n,m;
int lowbit(int x){return x&(-x);}
void update(int x,int y,int z)
{
	for (int i=x;i<=n;i+=lowbit(i))
		for (int j=y;j<=n;j+=lowbit(j)) c[i][j]+=z;
}
int query(int x,int y)
{
	int 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&1?1:0;
}
int main()
{
	for (scanf("%d",&T);T;T--)
	{
		memset(c,0,sizeof(c));
		scanf("%d%d",&n,&m);;
		for (int i=1;i<=m;i++)
		{
			int x1,x2,y1,y2;
			char ch;scanf(" %c",&ch);
			if (ch=='C')
			{
				scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
				update(x1,y1,1);update(x2+1,y2+1,1);
				update(x2+1,y1,-1);update(x1,y2+1,-1);
			}
			else
				scanf("%d%d",&x1,&y1),printf("%d\n",query(x1,y1));
		}
		printf("\n");
	}
	return 0;
}