1018: [SHOI2008]堵塞的交通traffic
1018: [SHOI2008]堵塞的交通traffic
Description
有一天，由于某种穿越现象作用，你来到了传说中的小人国。小人国的布局非常奇特，整个国家的交通系统可
以被看成是一个2行C列的矩形网格，网格上的每个点代表一个城市，相邻的城市之间有一条道路，所以总共有2C个
城市和3C-2条道路。 小人国的交通状况非常槽糕。有的时候由于交通堵塞，两座城市之间的道路会变得不连通，
直到拥堵解决，道路才会恢复畅通。初来咋到的你决心毛遂自荐到交通部某份差事，部长听说你来自一个科技高度
发达的世界，喜出望外地要求你编写一个查询应答系统，以挽救已经病入膏肓的小人国交通系统。 小人国的交通
部将提供一些交通信息给你，你的任务是根据当前的交通情况回答查询的问题。交通信息可以分为以下几种格式：
Close r1 c1 r2 c2：相邻的两座城市(r1,c1)和(r2,c2)之间的道路被堵塞了；Open r1 c1 r2 c2：相邻的两座城
市(r1,c1)和(r2,c2)之间的道路被疏通了；Ask r1 c1 r2 c2：询问城市(r1,c1)和(r2,c2)是否连通。如果存在一
条路径使得这两条城市连通，则返回Y，否则返回N；
Input
第一行只有一个整数C，表示网格的列数。接下来若干行，每行为一条交通信息，以单独的一行“Exit”作为
结束。我们假设在一开始所有的道路都是堵塞的。我们保证 C小于等于100000，信息条数小于等于100000。
Output
对于每个查询，输出一个“Y”或“N”。
Sample Input
Open 1 1 1 2
Open 1 2 2 2
Ask 1 1 2 2
Ask 2 1 2 2
Exit
Sample Output
N
HINT
题解：
由于只有两行，我们可以在这两行里用线段树维护连通性。把一列的两个点看为整体，维护(i,j)即为维护i~j之间的连通性（注意，只在左右端点之间）。
维护的过程我就不说了，左上端点到右上端点，左上端点到右下端点等等，也够麻烦的。。。
不过事情还没完，一个点还可以通过其他路径到达目标点，比如左侧的点可以向左再向右走，但是这时它肯定经过它这一列的另外一个点，这个点能不能到它这一列的另外一个点可以在(1,L)中知道答案。
到此问题也就基本解决了（提一下，c++用函数返回数组真麻烦。。）
#include<stdio.h> #include<iostream> #include<algorithm> #include<string.h> using namespace std; #define p1 (p<<1) #define p2 (p<<1|1) char s[15]; int n,i,x,y,x1,y1,t[400005][7],f[100005][3],ans[7],ansl[7],ansr[7]; int fans; void renew(int a[7],int b[7],int x,int p,int mid) { int c[7],i; c[1]=(a[1]|(a[3]&f[mid][1]&a[6]&f[mid][2]&b[1])|(a[4]&f[mid][1]&f[mid][2]&b[1]&a[5])); c[2]=((b[2]|(b[3]&f[mid][1]&f[mid][2]&a[2]&b[6])|(b[5]&f[mid][1]&f[mid][2]&a[2]&b[4]))); c[3]=((a[3]&f[mid][1]&b[3])|(a[4]&f[mid][2]&b[5])); c[6]=((a[6]&f[mid][2]&b[6])|(a[5]&f[mid][1]&b[4])); c[4]=((a[3]&f[mid][1]&b[4])|(a[4]&f[mid][2]&b[6])); c[5]=((a[6]&f[mid][2]&b[5])|(a[5]&f[mid][1]&b[3])); for(i=1;i<=6;i++) { if(x==1) t[p][i]=c[i];else if(x==2) ansl[i]=c[i];else if(x==3) ans[i]=c[i];else ansr[i]=c[i]; } } void build(int l,int r,int p) { if(l==r) { t[p][3]=1; t[p][6]=1; return; } int mid=(l+r)>>1; build(l,mid,p1); build(mid+1,r,p2); renew(t[p1],t[p2],1,p,mid); } void update(int l,int r,int x,int y,int z,int k,int p) { if(l==r) { if(z==0) { t[p][1]=k; t[p][2]=k; t[p][4]=k; t[p][5]=k; } else { f[l][y]=k; } return; } int mid=(l+r)>>1; if(x<=mid) update(l,mid,x,y,z,k,p1);else update(mid+1,r,x,y,z,k,p2); renew(t[p1],t[p2],1,p,mid); } void solve(int l,int r,int x,int y,int z,int p) { int i; if(x<=l&&r<=y) { if(l==x) { for(i=1;i<=6;i++) { if(z==2) ansl[i]=t[p][i];else if(z==3) ans[i]=t[p][i];else ansr[i]=t[p][i]; } } else { if(z==2) renew(ansl,t[p],z,0,l-1);else if(z==3) renew(ans,t[p],z,0,l-1);else renew(ansr,t[p],z,0,l-1); } return; } int mid=(l+r)>>1; if(x<=mid) solve(l,mid,x,y,z,p1); if(y>mid) solve(mid+1,r,x,y,z,p2); } int main() { scanf("%d",&n); build(1,n,1); while(true) { scanf("%s",&s); if(s[0]=='E') break; if(s[0]=='O') { scanf("%d%d%d%d",&x,&y,&x1,&y1); if(y>y1) { swap(x,x1); swap(y,y1); } if(y==y1) update(1,n,y,0,0,1,1);else update(1,n,y,x,1,1,1); } else if(s[0]=='C') { scanf("%d%d%d%d",&x,&y,&x1,&y1); if(y>y1) { swap(x,x1); swap(y,y1); } if(y==y1) update(1,n,y,0,0,0,1);else update(1,n,y,x,1,0,1); } else { scanf("%d%d%d%d",&x,&y,&x1,&y1); if(y>y1) { swap(x,x1); swap(y,y1); } for(i=1;i<=6;i++) { ansl[i]=0; ans[i]=0; ansr[i]=0; } solve(1,n,1,y,2,1); solve(1,n,y,y1,3,1); solve(1,n,y1,n,4,1); fans=ans[x*2+x1]; fans=fans|(ansl[2]&ans[(3-x)*2+x1])|(ansr[1]&ans[x*2+3-x1])|(ansl[2]&ansr[1]&ans[(3-x)*2+3-x1]); if(fans) printf("Y\n");else printf("N\n"); } } return 0; }
