POJ 2155 Matrix（树状数组+容斥原理）

#include <cstring>
#include <cstdio>
const int N=1010;
using namespace std;
int c[N][N],n,m,T;
char op[2];
void add(int x,int y){
    int i,j,k;
    for(i=x;i<N;i+=i&-i)
    for(j=y;j<N;j+=j&-j)
    c[i][j]^=1;
}
int sum(int x,int y){
    int i,j,ret=0;
    for(i=x;i;i-=i&-i)
    for(j=y;j;j-=j&-j)
    ret^=c[i][j];
    return ret;
}
int main(){
    scanf("%d",&T);
    for(int cas=0;cas<T;cas++){
        if(cas)puts("");
        memset(c,0,sizeof(c));
        scanf("%d%d",&n,&m);
        while(m--){
            int x,y,x1,y1;
            scanf("%s%d%d",op,&x,&y);
            if(op[0]==‘Q‘)printf("%d\n",sum(x,y));
            else{
                scanf("%d%d",&x1,&y1);
                add(x,y); add(x1+1,y);
                add(x,y1+1); add(x1+1,y1+1);
            }
        }
    }return 0;
}