Matrix POJ

来源：互联网发布：saa7130监控软件编辑：程序博客网时间：2024/05/19 14:55

//#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_backconst 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\n",ans&1);}}if(T)puts("");}}

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

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

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

阅读全文

0 0