[bzoj1018] [SHOI2008]堵塞的交通

时间:2024-11-23 20:05:49

题目描述

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

输出格式:

对于每个查询,输出一个YN

输入输出样例

输入样例#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_1\qquad a_2\\
-\\
a_3\qquad 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_1\qquad a_2\qquad b_1\qquad b_2\\
-\\
a_3\qquad a_4 \qquad b_3 \qquad b_4
\]

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

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

\[a_1\qquad a_2\qquad b_1\qquad b_2\qquad c_1\qquad c_2\\
-\\
a_3\qquad a_4 \qquad b_3 \qquad b_4 \qquad c_3 \qquad c_4
\]

注意一些情况,比如询问\(b_1\)到\(b_4\),可能存在一条路线

\[b_1\rightarrow a_2\rightarrow a_4\rightarrow b_3 \rightarrow b_2
\rightarrow c_1 \rightarrow c_3 \rightarrow b_4\\
\text {或者} b_1\rightarrow a_2\rightarrow a_4\rightarrow b_3 \rightarrow b_4\text{等}
\]

其余几种也类似。

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

#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\n",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;
}