Matrix POJ

//#include<bits/stdc++.h>
#include<stdio.h>
#include<algorithm>
#include<queue>
#include<string.h>
#include<iostream>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<iomanip>
using namespace std;
#define ll long long
#define pb push_back
const int maxn=1e3+3;

int ans;
int n;
int X1,X2,Y1,Y2;
bool ST[maxn<<2][maxn<<2];

void updatey(int yl,int yr,int xrt,int yrt){
	if(Y1<=yl&&Y2>=yr){ST[xrt][yrt]=!ST[xrt][yrt];return;}
	int m=yl+yr>>1;
	if(Y1<=m)updatey(yl,m,xrt,yrt<<1);
	if(Y2>m)updatey(m+1,yr,xrt,yrt<<1|1);
}
void updatex(int xl,int xr,int xrt){
	if(X1<=xl&&X2>=xr){updatey(1,n,xrt,1);return;}
	int m=xl+xr>>1;
	if(X1<=m)updatex(xl,m,xrt<<1);
	if(X2>m)updatex(m+1,xr,xrt<<1|1);
}
void queryy(int yl,int yr,int xrt,int yrt){
	ans+=ST[xrt][yrt];
	if(yl==yr)return;
	int m=yl+yr>>1;
	if(Y1<=m)queryy(yl,m,xrt,yrt<<1);
	else queryy(m+1,yr,xrt,yrt<<1|1);
}
void queryx(int xl,int xr,int xrt){
	queryy(1,n,xrt,1);
	if(xl==xr)return;
	int m=xl+xr>>1;
	if(X1<=m)queryx(xl,m,xrt<<1);
	else queryx(m+1,xr,xrt<<1|1);
}
int main(){
	int T;
	scanf("%d",&T);while(T--){
		memset(ST,0,sizeof ST);int q;
		scanf("%d%d",&n,&q);
		char op[5];
		while(q--){
			scanf("%s",op);
			if(op[0]=='C'){scanf("%d%d%d%d",&X1,&Y1,&X2,&Y2);updatex(1,n,1);}
			else{ans=0;scanf("%d%d",&X1,&Y1);queryx(1,n,1);
				printf("%d
",ans&1);
			}
		}
		if(T)puts("");
	}
}

这种做法没有区间合并，子一维区间不向上传递信息

所以在查找的时候，自顶向下，对于每一个x区间都要做一次y区间的query，说实话感觉挺浪费时间的

应该会有更好的树套树做法吧，或许是这题特殊，暂时先这样