【BZOJ】【1018】【SHOI2008】堵塞的交通traffic

线段树

---

　　这题的线段树+分类讨论蛮神奇的……我以前学的线段树简直就是渣渣QAQ

　　看了下ydc题解里的思想>_>用线段树维护连通性！那么就自己写吧……每个节点表示一段区间的连通性（我的叶子节点表示的是一个方块型的四个点之间的连通性，所以我直接n--了）对线段树上每个节点维护6个信息，即四个端点中任意一对点之间的连通性。

　　维护连通性的时候要进行信息的整合，也就是说间接连通的要全部找到并标记上连通，举个例子：如果左边连通且下边两端点连通，则左上到右下连通。这是一个简单的讨论我就不细说了，自己想一下真的很容易<_<。

　　但是这时候我遇到了一个蛋疼的问题：对于叶子节点来说，如果我上下左右四条边都连着，那么肯定是完全连通的……但是如果这时候我删一条边（Close）该怎么办？因为连通性并没有任何变化，但是图跟原来不一样了……所以我选择单独开一个a数组存叶子的直接连通情况（也就是存原图）然后用线段树t维护直接间接连通。

　　

　　好的这个时候这个线段树我们就维护好了，那么查询的时候就方便啦~我们可以查询一下L=[1,c1-1]这个区间的右端点是否连通

（保证了这种情况：aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAAAJYAAABUCAIAAADvQ1kKAAAABmJLR0QA/wD/AP+gvaeTAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADJUlEQVR4nO2dP3KjMByFpZ29imc8aVxYvoQ926TbbeEEzgEyW24BJ7DbLdOZSyyNm4wHDkMKxxPMykLoT+A576syRAh++pBB44Qny7JcLpeCwPJt7BMgvlAhPFQIDxXCY6WwSKWUq7wOdMw6X0kp0yJQd3EJXLsWzwEpy7Ixc0iESA7drVWmNFt1aFvqO50cvrVrCTwgvQqrTAmVVZ2Dnek7pKlllanJS/SpXUuUAelTePPisL9qbrWc/ET0q/1DV7dx4AHpuRcWL3uRPK4dP6TNrB8TsX+Z7h3Rp/YilZvjZf5WmTjZ3EsdB8SssD4dhXqYDe3UktmDEker4sbAo/Y6f96r7O/2fefZdre16sdtQPqfSBfzaArni8uPRSo/uDyajbvRo/b6tXTatTUgA/jusE8E1rum2U1rIwz9szDeR119OkbqORSOtbveIdwGxKxwNl+I8tX+XIatg10/bz4Hj9pn299J+fTrMg51ntoNieOADF9UtJ+WhRCtldMhEVfrqNstPZ6hPw+v2rWLiigD4rC0Nza1PgXMpb2xqV810Zb2jf3FUWXKst4BnY5MlNq9jvQ/Fgrf57/nKbaoMiUgBDZN8Nq1+A2I5B9eoMPvC+GhQnioEB4qhIcK4aFCeKgQHiqEhwrhoUJ4qBAeKoSHCuGhQnioEB4qhCewwiKN/Y9cpAtnITxUCA8VwkOF8FAhPFQIDxXCQ4XwUCE8VAgPFcJDhfBQITyB33ix3jVx3jNEbsJZCA8VwkOF8FAhPFQIDxXCQ4XwUCE8VAgPFcJDhfBQITxUCA8VwkOF8FAhPExO6+Eek9OuXxFu8SrbO0pOG1y7lrGT0w5J60g9rwS/t+S0IbVrmVZymuXvDU0mPxH9akdJTvPJ5UJPTjPVjpOcVvx5KtXPH46BL9jJaYbaYZLT6ny1aZ/pYICT08y1YySnFanc7FVW/XMWeAVUclrg2n1xSk4rUrnZJ4fGswbI5DSb2ieenFbnq/Odeqe50d95cppt7dNOTjtn03S4rJ3uPDltUO1MTpsETE4L0enIMDlND5PTOjA57YvD7wvhoUJ43gCCj3eO9OkLOgAAAABJRU5ErkJggg==" alt="" />）

　　查询一下R=[c2,n]这个区间左端点是否连通，跟上图类似，是在右边连通的

　　最后再查询一下mid=[c1,c2-1]这个区间的连通情况，如果L的右端点连通，则mid的左端点是连通的，mid的右端点同理（R的左端点）

　　然后重新维护下mid的连通情况……OK啦~

　　这题写的好辛苦……分类讨论总是想不清……

　　最后给个样例图吧~（红色的为查询点）

aaarticlea/png;base64,iVBORw0KGgoAAAANSUhEUgAAAjcAAACKCAIAAADZgfx7AAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAAArlSURBVHhe7dVRjty2FkXRDOR9Zv4ze2NwAGef6JpJdeuKYllE7QV9baurjYNm8Y8fkiQ9lbeU9NH+/78///3wb9IDeEtJn2i4ll49vC39Pt5S0mcZ7qEzDz8p/Q7eUtrJ8O35z8M/6zvDbucffl56O28pbWD4xvzi4Qf0L8NQ1x4+S3ojbyk93fBF+e3Dj+lXw0rXHj5LeiNvKT3X8BXZevgI/TSMM/PwidK7eEvpoYYvxwsPH/TxhlnmHz5XegtvKT3R8LV4+eHjPtuwyfzD50pv4S2lxxm+EycfPvRTDWvc9fDp0nreUnqc4Qtx8uFDP9Wwxl0Pny6t5y2lZxm+DW95+OjPM+xw78PvkBbzltKzDF+Ftzx89OcZdrj34XdIi3lL6UGG78EbH37BhxlGuPfhd0iLeUvpQYbvwRsffsGHGUa49+F3SIsdt9Qf0u82fA/e+PALvsRJ2Ar/9ReGEe59+B3SrfjLLryl9CDD9+CND7/gS5yErfBff2EY4d6H3yHdir/swltKDzJ8D9748Au+xEnYCv/1F4YR7n34HdKt+MsuvKX0IMP34I0Pv+BLnISt8F9/YRjh3offId2Kv+zCW0oPMnwP3vjwC77ESdgK//XXhh3uevh06W78ZRcvbymq5rBmUPXa8G14y8NHf6RhirsePv2zcaqDqibmC2rhLbUWawZVrw3fhrc8fPSnGtaYf/jcj8epDqqamC+ohbfUWqwZVH1p+E6cfPjQDzYMMv/wuR+PUx1UNTFfUAtvqbVYM6j60vCdOPnwoZ9t2GTm4RPl6b4J8wW18JZaizWDqu8M34yXHz5OTroApzqoamK+oBbeUmuxZlB1wvDleOHhg/TTMM61h8/ST5zqoKqJ+YJaeEutxZpB1TnDV2Tr4SP0q2Gl8w8/r4JTHVQ1MV9QC2+ptVgzqOoYvi6/ffgxvTDMdebhJ/UrTnVQ1cR8QS28pdZizaCqafjS/OLhB/SdYbdXD2/rv3Cqg6om5gtq4S21FmsGVVcN36H/PPyz+oYl/374N32JUx1UNTFfUAtvqbVYM6iawJRB1VXsGFSdwGRBVRPzBbU4Eq8EVXNYM6iawJRB1VXsGFSdwGRBVRPzBbU4Eq8EVXNYM6iawJRB1VXsGFSdwGRBVRPzBbU4Eq8EVXNYM6iawJRB1VXsGFSdwGRBVRPzBbU4Eq8EVXNYM6iawJRB1VXsGFSdwGRBVRPzBbU4Eq8EVXNYM6iawJRB1VXsGFSdwGRBVRPzBbU4Eq8EVXNYM6iawJRBld6OP8Ggqon5glociVeCqjmsGVRNYMqgSm/Hn2BQ1cR8QS2OxCtB1RzWDKomMGVQpbfjTzCoamK+oBZH4pWgag5rBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQJe2PUx1UNTFfUAtvqbVYM6iS9sepDqqamC+ohbfUWqwZVEn741QHVU3MF9TCW2ot1gyqpP1xqoOqJuYLauEttRZrBlXS/jjVQVUT8wW18JZaizWDKml/nOqgqon5glp4S63FmkGVtD9OdVDVxHxBLbyl1mLNoEraH6c6qGpivqAW3lJrsWZQNYEpgyq9HX+CQVUT8wW1OBKvBFVzWDOomsCUQZXejj/BoKqJ+YJaHIlXgqo5rBlUTWDKoEpvx59gUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrQdUc1gyqJjBlUHUVOwZVJzBZUNXEfEEtjsQrkiT9JlxIhbeUJOkpuJAKbylJ0lNwIRXeUpKkp+BCKrylJElPwYVU/EeSJOkhvKUkSc/lLSVJei5vKUnSU/348Rc66nPre/ELGwAAAABJRU5ErkJggg==" alt="" width="350" height="85" />

 /**************************************************************
     Problem: 1018
     User: Tunix
     Language: C++
     Result: Accepted
     Time:1068 ms
     Memory:4208 kb
 ****************************************************************/
 
 //BZOJ 1018
 #include<vector>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 #define pb push_back
 using namespace std;
 inline int getint(){
     int v=,sign=; char ch=getchar();
     while(ch<''||ch>''){ if (ch=='-') sign=-; ch=getchar();}
     while(ch>=''&&ch<=''){ v=v*+ch-''; ch=getchar();}
     return v*sign;
 }
 const int N=1e5+,INF=~0u>>;
 typedef long long LL;
 /******************tamplate*********************/
 int n;
 struct node{
     bool z,y,s,x,zs,zx;
     node(bool z=,bool y=,bool s=,bool x=,bool zs=,bool zx=):
         z(z),y(y),s(s),x(x),zs(zs),zx(zx){}
     //z 左端两结点是否连通
     //y 右端
     //s 上方
     //x 下方
     //zs 左上->右下
     //zx 左下->右上
 }a[N],t[N<<],ans;
 #define L (o<<1)
 #define R (o<<1|1)
 #define mid (l+r>>1)
 void maintain(node &a){
     a.zs|=(a.z&a.x) | (a.s&a.y);
     a.zx|=(a.z&a.s) | (a.x&a.y);
     a.z|=(a.zs&a.x) | (a.zx&a.s);
     a.y|=(a.zs&a.s) | (a.zx&a.x);
     a.s|=(a.zs&a.y) | (a.zx&a.z);
     a.x|=(a.zs&a.z) | (a.zx&a.y);
 }
 node Union(node x,node y){
     node tmp;
     tmp.z =x.z; tmp.y=y.y;
     tmp.s =(x.s&y.s)  | (x.zs&y.zx) ;
     tmp.x =(x.x&y.x)  | (x.zx&y.zs) ;
     tmp.zs=(x.s&y.zs) | (x.zs&y.x) ;
     tmp.zx=(x.x&y.zx) | (x.zx&y.s) ;
     return tmp;
 }
 void update(int o,int l,int r,int pos,int num,bool v){
 //将pos位置的num号连通情况改为v
     if (l==r){
         if (num==) a[l].z=v; if (num==) a[l].y=v;
         if (num==) a[l].s=v; if (num==) a[l].x=v;
         t[o]=a[l];
     }else{
         if (pos<=mid) update(L,l,mid,pos,num,v);
         else update(R,mid+,r,pos,num,v);
         t[o]=Union(t[L],t[R]);
     }
     maintain(t[o]);
 }
 int ql,qr;bool sign=;
 void query(int o,int l,int r){
     if (ql<=l && qr>=r){
         if (!sign) {ans=t[o];sign=;}
         else ans=Union(ans,t[o]);
     }
     else{
         if (ql<=mid) query(L,l,mid);
         if (qr>mid) query(R,mid+,r);
     }
 }
 int main(){
 #ifndef ONLINE_JUDGE
     freopen("1018.in","r",stdin);
     freopen("1018.out","w",stdout);
 #endif
     n=getint()-;
     int r1,r2,c1,c2;
     char cmd[];
     while(scanf("%s",cmd)!=EOF && cmd[]!='E'){
         r1=getint(); c1=getint(); r2=getint(); c2=getint();
         if (c1>c2){swap(c1,c2); swap(r1,r2);}
         if (cmd[]=='A'){
             ans=node();
             if (c1==c2){
                 if (c1<n+ && c1>){//不是右端点
                     node t1,t2;
                     ql=; qr=c1-; sign=;
                     query(,,n); t1=ans;
                     ql=c1; qr=n; sign=;
                     query(,,n);
                     ans.z|=t1.y; ans.z|=t2.z;
                     puts(ans.z?"Y":"N");
                 }else if(c1==n+){
                     ans=t[];
                     puts(ans.y?"Y":"N");
                 }else{
                     ans=t[];
                     puts(ans.z?"Y":"N");
                 }
             }else{
                 node t1,t2;
                 if (c1!=){ql=; qr=c1-; sign=; query(,,n); t1=ans;}
                 if (c2!=n+){ql=c2; qr=n; sign=; query(,,n); t2=ans;}
                 ql=c1; qr=c2-; sign=;
                 query(,,n);
                 ans.z|=t1.y | (t1.s & t1.z & t1.x);
                 ans.y|=t2.z | (t2.s & t2.y & t2.x);
                 maintain(ans);
                 if (r1== && r2==) puts(ans.s ?"Y":"N");
                 if (r1== && r2==) puts(ans.x ?"Y":"N");
                 if (r1== && r2==) puts(ans.zs?"Y":"N");
                 if (r1== && r2==) puts(ans.zx?"Y":"N");
             }
         }else{
             bool type=;
             if (cmd[]=='O') type=;
             else type=;
             if (r1==r2){
                 if (r1==) update(,,n,c1,,type);
                 else update(,,n,c1,,type);
             }else{
                 if (c1<n+) update(,,n,c1,,type);
                 if (c1>  ) update(,,n,c1-,,type);
             }
         }
     }
     return ;
 }
 

1018: [SHOI2008]堵塞的交通traffic

Time Limit: 3 Sec  Memory Limit: 162 MB
Submit: 2030  Solved: 638
[Submit][Status][Discuss]

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”作为结束。我们假设在一开始所有的道路都是堵塞的。
对30%测试数据，我们保证C小于等于1000，信息条数小于等于1000； 对100%测试数据，我们保证
C小于等于100000，信息条数小于等于100000。

Output

对于每个查询，输出一个“Y”或“N”。

Sample Input

2
Open 1 1 1 2
Open 1 2 2 2
Ask 1 1 2 2
Ask 2 1 2 2
Exit

Sample Output

Y
N

HINT

Source

[Submit][Status][Discuss]