[SHOI2008]堵塞的交通（线段树维护联通性）

题目

2行c列个点，开始时互不联通，支持给同一列或着同一行相邻的两个点连边，和询问两个点能否在一个联通块里。

1≤C≤100000 1<=操作数<=100000;

题解

线段树的又一个骚操作。

我们把2*2的4个点看作线段树上的叶子结点。其他节点就是其儿子的合并（叶子结点的父亲表示2*4八个点，然后是2*8，2*16）。

然后在节点上维护节点表示的点的联通信息。

                                             

信息维护六种，就是用同种颜色连接的点是否联通。

所以信息合并时要写一堆。

然后发现这样叶子节点的信息不好修改（你写写就知道了），所以我们在叶子结点上额外记录实际连边的情况。

查询大体思路是把所有点分成三分，然后枚举所有情况。（假如两个点在同一列要特判）

（这些线代表的是连通情况，并不是实际路线）

然后就可以通过了

这个题告诉我：对拍是个好东西

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 const int N=;
 int n;
 struct jz{
     int a[];
 }hhh;
 struct tree{
     int l,r;
     jz z,tru;
 }tr[N*];
 void build(int l,int r,int now){
     tr[now].l=l;tr[now].r=r;
     if(l==r){
         return;
     }
     int mid=(tr[now].l+tr[now].r)>>;
     build(l,mid,now*);
     build(mid+,r,now*+);
 }
 void update(int now){
     tr[now].z.a[]=(tr[now*].z.a[]&tr[now*+].z.a[])|(tr[now*].z.a[]&tr[now*+].z.a[]);
     tr[now].z.a[]=(tr[now*].z.a[]&tr[now*+].z.a[])|(tr[now*].z.a[]&tr[now*+].z.a[]);
     tr[now].z.a[]=(tr[now*].z.a[]&tr[now*+].z.a[])|(tr[now*].z.a[]&tr[now*+].z.a[]);
     tr[now].z.a[]=(tr[now*].z.a[]&tr[now*+].z.a[])|(tr[now*].z.a[]&tr[now*+].z.a[]);
     tr[now].z.a[]=tr[now*].z.a[]|(tr[now*].z.a[]&tr[now*].z.a[]&tr[now*+].z.a[]);
     tr[now].z.a[]=tr[now*+].z.a[]|(tr[now*+].z.a[]&tr[now*+].z.a[]&tr[now*].z.a[]);
 }
 void update(int x,int k,int c,int now){
     if(tr[now].l==tr[now].r){
         tr[now].tru.a[k]=c;
         for(int i=;i<=;i++){
             tr[now].z.a[i]=tr[now].tru.a[i];
         }
         tr[now].z.a[]=(tr[now].tru.a[]&tr[now].tru.a[])|(tr[now].tru.a[]&tr[now].tru.a[]);
         tr[now].z.a[]=(tr[now].tru.a[]&tr[now].tru.a[])|(tr[now].tru.a[]&tr[now].tru.a[]);
         tr[now].z.a[]=tr[now].z.a[]|(tr[now].z.a[]&tr[now].z.a[])|(tr[now].z.a[]&tr[now].z.a[]);
         tr[now].z.a[]=tr[now].z.a[]|(tr[now].z.a[]&tr[now].z.a[])|(tr[now].z.a[]&tr[now].z.a[]);
         tr[now].z.a[]=tr[now].z.a[]|(tr[now].z.a[]&tr[now].z.a[])|(tr[now].z.a[]&tr[now].z.a[]);
         tr[now].z.a[]=tr[now].z.a[]|(tr[now].z.a[]&tr[now].z.a[])|(tr[now].z.a[]&tr[now].z.a[]);
 //        for(int i=1;i<=6;i++){
 //            cout<<tr[now].z.a[i]<<" ";
 //        }
 //        cout<<endl;
         return;
     }
     int mid=(tr[now].l+tr[now].r)>>;
     if(x<=mid)update(x,k,c,now*);
     else update(x,k,c,now*+);
     update(now);
 }
 jz hb(jz a,jz b,jz c){
     c.a[]=(a.a[]&b.a[])|(a.a[]&b.a[]);
     c.a[]=(a.a[]&b.a[])|(a.a[]&b.a[]);
     c.a[]=(a.a[]&b.a[])|(a.a[]&b.a[]);
     c.a[]=(a.a[]&b.a[])|(a.a[]&b.a[]);
     c.a[]=a.a[]|(a.a[]&a.a[]&b.a[]);
     c.a[]=b.a[]|(b.a[]&b.a[]&a.a[]);
     return c;
 }
 jz check(int l,int r,int now){
     if(tr[now].l==l&&tr[now].r==r){
         return tr[now].z;
     }
     int mid=(tr[now].l+tr[now].r)>>;
     if(l>mid)return check(l,r,now*+);
     else if(r<=mid)return check(l,r,now*);
     else {
         return hb(check(l,mid,now*),check(mid+,r,now*+),hhh);
     }
 }
 int main(){
     scanf("%d",&n);
     build(,n+,);
     string s;
     while(cin>>s){
         if(s[]=='E')break;
         int x1,y1,x2,y2;
         scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
         if(y1>y2)swap(x2,x1),swap(y1,y2);
         if(s[]=='O'){
             if(y1==y2){
                 update(y1,,,);
                 update(y1+,,,);
             }
             else {
                 if(x1==)update(y1+,,,);
                 else update(y1+,,,);
             }
         }
         else if(s[]=='C'){
             if(y1==y2){
                 update(y1,,,);
                 update(y1+,,,);
             }
             else {
                 if(x1==)update(y1+,,,);
                 else update(y1+,,,);
             }
         }
         else {
             if(y1==y2){
                 jz x=check(,y1,);
                 jz y=check(y1+,n+,);
                 if(x.a[]|y.a[])printf("Y\n");
                 else printf("N\n");
             }
             else{
                 jz z=check(y1+,y2,);
                 jz x=check(,y1,);
                 jz y=check(y2+,n+,);
                 if(x1==x2){
                     if(x1==){
                         if(z.a[]|(x.a[]&z.a[]&y.a[])|(x.a[]&z.a[])|(z.a[]&y.a[]))printf("Y\n");
                         else printf("N\n");
                     }
                     else if(z.a[]|(x.a[]&z.a[]&y.a[])|(x.a[]&z.a[])|(z.a[]&y.a[]))printf("Y\n");
                          else printf("N\n");
                 }
                 else{
                     if(x1==){
                         if(z.a[]|(x.a[]&z.a[])|(y.a[]&z.a[])|(x.a[]&z.a[]&y.a[]))printf("Y\n");
                         else printf("N\n");
                     }
                     else if(z.a[]|(x.a[]&z.a[])|(y.a[]&z.a[])|(x.a[]&z.a[]&y.a[]))printf("Y\n");
                          else printf("N\n");
                 }
             }
         }
     }
     return ;
 }