[SHOI2008]堵塞的交通

Posted 2021-01-21 hbyer

题目描述

有一天，由于某种穿越现象作用，你来到了传说中的小人国。小人国的布局非常奇特，整个国家的交通系统可以被看成是一个22行CC列的矩形网格，网格上的每个点代表一个城市，相邻的城市之间有一条道路，所以总共有2C个城市和3C?2条道路。

小人国的交通状况非常槽糕。有的时候由于交通堵塞，两座城市之间的道路会变得不连通，直到拥堵解决，道路才会恢复畅通。初来咋到的你决心毛遂自荐到交通部某份差事，部长听说你来自一个科技高度发达的世界，喜出望外地要求你编写一个查询应答系统，以挽救已经病入膏肓的小人国交通系统。 小人国的交通部将提供一些交通信息给你，你的任务是根据当前的交通情况回答查询的问题。交通信息可以分为以下几种格式：

• Close r1 c1 r2 c2：相邻的两座城市((r_1, c_1)(r1,c1))和((r_2, c_2)(r2,c2))之间的道路被堵塞了；
• Open r1 c1 r2 c2：相邻的两座城市((r_1, c_1)(r1,c1))和((r_2, c_2)(r2,c2))之间的道路被疏通了；
• Ask r1 c1 r2 c2：询问城市((r_1, c_1)(r1,c1))和((r_2, c_2)(r2,c2))是否连通。如果存在一条路径使得这两条城市连通，则返回Y，否则返回N。

注：(r_i)表示行数，(c_i)表示列数，(1 leq r_i leq 2), (1 leq c_i leq C)。

输入输出格式

输入格式：

第一行只有一个整数CC，表示网格的列数。接下来若干行，每行为一条交通信息，以单独的一行Exit作为结束。我们假设在一开始所有的道路都是堵塞的。我们保证CC小于等于100000100000，信息条数小于等于100000100000.

输出格式：

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

输入输出样例

输入样例#1：

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

输出样例#1：

Y
N

说明

数据范围：

对于(100\%)的数据，(1 leq C leq 100000)，$1 leq (信息条数)leq 100000$。

题解

使用线段树维护。

考虑节点信息的合并，

对于每个点，图大概长这样
[ a_1qquad a_2-a_3qquad a_4 ]

其中1-4为四个端点，然后记录6个信息 ： (a _{1,2},a_{1,3},a_{1,4},a_{3,2},a_{3,4},a_{2,4})

其中(a_{i,j}?)代表i，j两点联不联通，(a_{i,j} in {0,1}?)。

然后合并：
[ a_1qquad a_2qquad b_1qquad b_2-a_3qquad a_4 qquad b_3 qquad b_4 ]

然后类似的讨论下，具体可以看代码的update部分。

然后询问((l,r))，先处理出((1,l-1),(l,r),(r+1,n))三段的联通情况，然后讨论。

[ a_1qquad a_2qquad b_1qquad b_2qquad c_1qquad c_2-a_3qquad a_4 qquad b_3 qquad b_4 qquad c_3 qquad c_4 ]

注意一些情况，比如询问(b_1)到(b_4),可能存在一条路线
[ b_1 ightarrow a_2 ightarrow a_4 ightarrow b_3 ightarrow b_2 ightarrow c_1 ightarrow c_3 ightarrow b_4\text {或者} b_1 ightarrow a_2 ightarrow a_4 ightarrow b_3 ightarrow b_4 ext{等} ]
其余几种也类似。

具体细节还是挺多的，不懂看看代码吧。

#include<bits/stdc++.h>
using namespace std;

#define maxn 500050

void read(int &x) {
    x=0;int f=1;char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;
    for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';x*=f;
}

void print(int x) {
    if(!x) return ;if(x<0) x=-x,putchar('-');
    print(x/10),putchar(x%10+48);
}
void write(int x) {if(!x) putchar('0');else print(x);puts("");}

#define ls p<<1
#define rs p<<1|1
#define mid ((l+r)>>1)

int e[maxn][4],n;

struct segment_tree {
    struct node {
        int a12,a13,a14,a32,a24,a34;
        node() {a12=a13=a14=a32=a24=a34=0;}
    }tr[maxn<<2];
    node merge(int pos,node x,node y) {
        node a;
        if(e[pos][1]) {
            a.a14|=x.a12&y.a14;
            a.a12|=x.a12&y.a12;
            a.a32|=x.a32&y.a12;
            a.a34|=x.a32&y.a14;
        }
        if(e[pos][2]) {
            a.a34|=x.a34&y.a34;
            a.a12|=x.a14&y.a32;
            a.a14|=x.a14&y.a34;
            a.a32|=x.a34&y.a32;
        }
        a.a13=x.a13|(x.a12&e[pos][1]&e[pos][2]&y.a13&x.a34);
        a.a24=y.a24|(e[pos][1]&e[pos][2]&y.a12&y.a34&x.a24);
        return a;
    }
    void build(int p,int l,int r) {
        if(l==r) return tr[p].a12=tr[p].a34=1,void();
        build(ls,l,mid),build(rs,mid+1,r);
        tr[p]=merge(mid,tr[ls],tr[rs]);
    }
    void modify(int p,int l,int r,int x,int y) {
        if(l==r) {
            if(y==1) tr[p].a13=tr[p].a24=tr[p].a14=tr[p].a32=1;
            else if(y==2) tr[p].a13=tr[p].a24=tr[p].a14=tr[p].a32=0;
            return ;
        }
        if(x<=mid) modify(ls,l,mid,x,y);
        else modify(rs,mid+1,r,x,y);
        tr[p]=merge(mid,tr[ls],tr[rs]);
    }
    node query(int p,int l,int r,int x,int y) {
        if(x<=l&&r<=y) return tr[p];
        if(x<=mid&&y<=mid) return query(ls,l,mid,x,y);
        if(x>mid&&y>mid) return query(rs,mid+1,r,x,y);
        return merge(mid,query(ls,l,mid,x,y),query(rs,mid+1,r,x,y));
    }
    /*
      1   2  1   2  1   2
      3   4  3   4  3   4
     */
    int query(int x,int y,int xx,int yy) {
        if(y>yy) swap(x,xx),swap(y,yy);int ans=0;
        int b1=e[y-1][1],b2=e[y-1][2],b3=e[yy][1],b4=e[yy][2];
        node a1,a2=query(1,1,n,y,yy),a3;//print(a2);
        if(y>1) a1=query(1,1,n,1,y-1);
        if(yy<n) a3=query(1,1,n,yy+1,n);
        if(x==1&&xx==1)
            ans=a2.a12|(b1&b2&a1.a24&a2.a32)|(b3&b4&a2.a14&a3.a13)|(b1&b2&b3&b4&a1.a24&a2.a34&a3.a13);
        if(x==1&&xx==2)
            ans=a2.a14|(b1&b2&a1.a24&a2.a34)|(b3&b4&a2.a12&a3.a13)|(b1&b2&b3&b4&a1.a24&a2.a32&a3.a13);
        if(x==2&&xx==1)
            ans=a2.a32|(b1&b2&a1.a24&a2.a12)|(b3&b4&a2.a34&a3.a13)|(b1&b2&b3&b4&a1.a24&a2.a14&a3.a13);
        if(x==2&&xx==2)
            ans=a2.a34|(b1&b2&a1.a24&a2.a14)|(b3&b4&a2.a32&a3.a13)|(b1&b2&b3&b4&a1.a24&a2.a12&a3.a13);
        return ans;
    }
    void print(node i) {printf("%d %d %d %d %d %d
",i.a12,i.a13,i.a14,i.a32,i.a24,i.a34);}
    void debug() {
        for(int i=1;i<=3;i++) print(tr[i]);
    }
}T;

int main() {
    read(n);int a,b,c,d;char s[10];
    T.build(1,1,n);
    while(1) {
        scanf("%s",s+1);if(s[1]=='E') return 0;
        read(a),read(b),read(c),read(d);
        if(s[1]=='A') {puts(T.query(a,b,c,d)?"Y":"N");continue;}
        if(b>d) swap(a,c),swap(b,d);
        if(b==d)
            if(s[1]=='O') T.modify(1,1,n,b,1);
            else T.modify(1,1,n,b,2);
        else {
            if(s[1]=='O') e[b][a]=1;
            else e[b][a]=0;
            T.modify(1,1,n,b,3);
        }//T.debug();
    }
    return 0;
}