[SHOI2008]堵塞的交通

题目描述

有一天，由于某种穿越现象作用，你来到了传说中的小人国。小人国的布局非常奇特，整个国家的交通系统可
以被看成是一个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；

题解

这么**的题我居然能一A。。。

我写的比较麻烦，用二元组[0/1][0/1]来表示当前矩形的四个角，再去维护每个小矩形，线段树维护合并后的大矩形。

讨论比较恶心。。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 100009
using namespace std;
char s[10];
int n;
inline int rd(){
    int x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return f?-x:x;
}
struct node{
    bool a[2][2][2][2];
    node operator +(const node &b)const{
        node c;
        c.a[0][0][0][1]=(a[0][0][1][1]&&b.a[1][0][0][1])||(a[0][0][0][1]&&b.a[0][0][0][1]);
        c.a[0][0][1][0]=(a[0][0][1][0])||(a[0][0][0][1]&&b.a[0][0][1][0]&&a[1][0][1][1]);
        c.a[0][0][1][1]=(a[0][0][1][1]&&b.a[1][0][1][1])||(a[0][0][0][1]&&b.a[0][0][1][1]);
        c.a[1][0][0][1]=(a[1][0][1][1]&&b.a[1][0][0][1])||(a[1][0][0][1]&&b.a[0][0][0][1]);
        c.a[1][0][1][1]=(a[1][0][1][1]&&b.a[1][0][1][1])||(a[1][0][0][1]&&b.a[0][0][1][1]);
        c.a[0][1][1][1]=(b.a[0][1][1][1])||(a[0][1][1][1]&&b.a[0][0][0][1]&&b.a[1][0][1][1]);
        return c;
    }
    inline void clear(){memset(a,0,sizeof(a));}
    inline void print(){
      cout<<a[0][0][0][1]<<" "<<a[1][0][1][1]<<" "
      <<a[0][0][1][1]<<" "<<a[0][0][1][0]<<" "<<a[0][1][1][1]<<" "<<a[1][0][0][1]<<endl;
    }
}now[N];
struct seg{
    node val;int ls,rs;
}tr[N<<1];
int tot=1;
void build(int cnt,int l,int r){
    if(l==r)return;
    int mid=(l+r)>>1;tr[cnt].ls=++tot;tr[cnt].rs=++tot;
    build(tr[cnt].ls,l,mid);build(tr[cnt].rs,mid+1,r);
}
void upd(int cnt,int l,int r,int tag,int x,int y){
    if(l==r){
        if(tag<2)now[l].a[tag][0][tag][1]=y;
        else now[l].a[0][tag-2][1][tag-2]=y;
        tr[cnt].val.clear();
        tr[cnt].val.a[0][0][1][0]=now[l].a[0][0][1][0]||(now[l].a[0][0][0][1]&&now[l].a[0][1][1][1]&&now[l].a[1][0][1][1]);
        tr[cnt].val.a[0][0][0][1]=now[l].a[0][0][0][1]||(now[l].a[0][0][1][0]&&now[l].a[1][0][1][1]&&now[l].a[0][1][1][1]);
        tr[cnt].val.a[0][0][1][1]=(now[l].a[0][0][0][1]&&now[l].a[0][1][1][1])||(now[l].a[0][0][1][0]&&now[l].a[1][0][1][1]);
        tr[cnt].val.a[1][0][0][1]=(now[l].a[0][0][1][0]&&now[l].a[0][0][0][1])||(now[l].a[1][0][1][1]&&now[l].a[0][1][1][1]);
        tr[cnt].val.a[1][0][1][1]=(now[l].a[1][0][1][1])||(now[l].a[0][0][1][0]&&now[l].a[0][0][0][1]&&now[l].a[0][1][1][1]);
        tr[cnt].val.a[0][1][1][1]=(now[l].a[0][1][1][1])||(now[l].a[0][0][0][1]&&now[l].a[0][0][1][0]&&now[l].a[1][0][1][1]);
    //    cout<<l<<" ";tr[cnt].val.print();
        return;
    }
    int mid=(l+r)>>1;
    if(mid>=x)upd(tr[cnt].ls,l,mid,tag,x,y);
    else upd(tr[cnt].rs,mid+1,r,tag,x,y);
    tr[cnt].val=tr[tr[cnt].ls].val+tr[tr[cnt].rs].val;
}
node _query(int cnt,int l,int r,int x,int tag){
    if(l==r){
      node x;
      int xx=tr[cnt].val.a[0][tag][1][tag];
      x.a[0][0][0][1]=x.a[1][0][1][1]=1;
      x.a[0][0][1][1]=x.a[0][0][1][0]=x.a[0][1][1][1]=x.a[1][0][0][1]=xx;
      return x;
    }
    int mid=(l+r)>>1;
    if(mid>=x)return _query(tr[cnt].ls,l,mid,x,tag);
    else return _query(tr[cnt].rs,mid+1,r,x,tag);
}
node query(int cnt,int l,int r,int L,int R){
    if(L>R)return node();
    if(l>=L&&r<=R)return tr[cnt].val;
    int mid=(l+r)>>1;
    if(mid>=L&&mid<R)return query(tr[cnt].ls,l,mid,L,R)+query(tr[cnt].rs,mid+1,r,L,R);
    else if(mid>=L)return query(tr[cnt].ls,l,mid,L,R);
    else return query(tr[cnt].rs,mid+1,r,L,R); 
}
int main(){
    n=rd();
    build(1,1,n-1);
    int r1,r2,x,y;
    while(1){
        scanf("%s",s);if(s[0]=='E')break;
        r1=rd();x=rd();r2=rd();y=rd();r1--;r2--;
        if(x>y)swap(x,y),swap(r1,r2);
        if(s[0]=='C'){
            if(x==y){if(x!=n)upd(1,1,n-1,2,x,0);if(x>1)upd(1,1,n-1,3,x-1,0);}
            else upd(1,1,n-1,r1,min(x,y),0);
        }
        else if(s[0]=='O'){
            if(x==y){if(x!=n)upd(1,1,n-1,2,x,1);if(x>1)upd(1,1,n-1,3,x-1,1);}
            else upd(1,1,n-1,r1,min(x,y),1);
        }
        else{
            bool ans=0;
            node xx=query(1,1,n-1,x,y-1),x1=query(1,1,n-1,1,x-1),x2=query(1,1,n-1,y,n-1);
            if(x==y)xx=_query(1,1,n-1,x==1?x:x-1,x==1?0:1);
    //        xx.print();x1.print();x2.print();
            if(r1){
                if(r2){
                    if(xx.a[1][0][1][1])ans=1;
                    if(x1.a[0][1][1][1]&&x2.a[0][0][1][0]&&xx.a[0][0][0][1])ans=1;
                }
                else{
                    if(xx.a[1][0][0][1])ans=1;
                    if(x1.a[0][1][1][1]&&xx.a[0][0][0][1])ans=1;
                    if(x2.a[0][0][1][0]&&xx.a[1][0][1][1])ans=1;
                }
            }
            else{
                if(r2){
                    if(xx.a[0][0][1][1])ans=1;
                    if(x1.a[0][1][1][1]&&xx.a[1][0][1][1])ans=1;
                    if(x2.a[0][0][1][0]&&xx.a[0][0][0][1])ans=1;
                }
                else{
                    if(xx.a[0][0][0][1])ans=1;
                    if(x1.a[0][1][1][1]&&x2.a[0][0][1][0]&&xx.a[1][0][1][1])ans=1;
                }
            }
            if(ans)puts("Y");else puts("N"); 
        }
    }
    return 0;
}