bzoj1018 堵塞的交通traffic(线段树)

题目传送门

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

Output

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

题解

这道题是几天前的考试题了，今天花了一下午才把它打下来。真讨厌数据结构这种代码长又繁琐的题了。

此题用每个节点表示区间内的连通性（不用区间外的边），edge数组储存每个节点的三条边，之后进行维护就好了。

#include <cstdio>#include <iostream>#include <algorithm>#include <cstring>using namespace std;const int Maxn = 100005;#define ls (u<<1)#define rs (u<<1|1)#define MID(a,b) int mid = (a+b) >> 1;struct node {bool a[2], b[2], c[2];<span style="white-space:pre"></span>// a  x  b  <span style="white-space:pre"></span>// c  y  d  <span style="white-space:pre"></span>//a[] means a-c b-d<span style="white-space:pre"></span>//b[] means a-b c-d<span style="white-space:pre"></span>//c[] means a-d b-cnode() {}}N[Maxn << 2];int C, r1, r2, c1, c2;bool edge[Maxn][3];void build(int u,int l,int r) {if (l == r) {N[u].b[0] = N[u].b[1] = true; return;}MID(l,r);build(ls,l,mid); build(rs,mid+1,r);}node pluse(node a, node b,bool x, bool y) {node c;c.a[0] = (a.a[0])||(a.b[0]&&a.b[1]&&x&&y&&b.a[0]);c.a[1] = (b.a[1])||(b.b[0]&&b.b[1]&&x&&y&&a.a[1]);c.b[0] = (a.b[0]&&x&&b.b[0])||(a.c[0]&&y&&b.c[1]);c.b[1] = (a.b[1]&&y&&b.b[1])||(a.c[1]&&x&&b.c[0]);c.c[0] = (x&&a.b[0]&&b.c[0])||(y&&a.c[0]&&b.b[1]);c.c[1] = (y&&a.b[1]&&b.c[1])||(x&&a.c[1]&&b.b[0]);return c;}void update(int u,int l,int r,int val) {if (l == r) {N[u].b[0] = N[u].b[1] = true;N[u].a[0] = N[u].a[1] = N[u].c[0] = N[u].c[1] = edge[val][2];return;}MID(l,r);if (val <= mid) update(ls,l,mid,val);else update(rs,mid+1,r,val);N[u] = pluse(N[ls],N[rs],edge[mid][0],edge[mid][1]);} node query(int u,int l,int r,int x, int y) {if (l >= x && r <= y) return N[u];MID(l,r);if (y <= mid) return query(ls,l,mid,x,y);else if (x > mid) return query(rs,mid+1,r,x,y);return pluse(query(ls,l,mid,x,mid),query(rs,mid+1,r,mid+1,y),edge[mid][0],edge[mid][1]);}int main() {//freopen("B.in","r",stdin);scanf("%d", &C);build(1,1,C);char S[10] = {0};node tmp1, tmp2, tmp3; bool ans;while (scanf("%s",S)) { if (*S == 'E') break;scanf("%d%d%d%d",&r1,&c1,&r2,&c2);--r1, --r2;if (c1 > c2) swap(r1,r2), swap(c1,c2);if (*S == 'O') {if (r1 > r2) swap(r1,r2), swap(c1,c2);if (r1 < r2) {edge[c1][2] = 1; update(1,1,C,c1);} else {edge[c1][r1] = 1; update(1,1,C,c1);}} else if (*S == 'C') {if (r1 > r2) swap(r1,r2), swap(c1,c2);if (r1 < r2) {edge[c1][2] = 0; update(1,1,C,c1);} else {edge[c1][r1] = 0; update(1,1,C,c1);}} else  {tmp1 = query(1,1,C,1,c1), tmp2 = query(1,1,C,c1,c2), tmp3 = query(1,1,C,c2,C);if (r1 && r2) {ans=tmp2.b[1]||(tmp1.a[1]&&tmp2.c[0])||(tmp2.c[1]&&tmp3.a[0]);ans = ans || (tmp1.a[1]&&tmp2.b[0]&&tmp3.a[0]);} else if (! r1 ^ r2) {ans=tmp2.b[0]||(tmp1.a[1]&&tmp2.c[1])||(tmp2.c[0]&&tmp3.a[0]);ans = ans || (tmp1.a[1]&&tmp2.b[1]&&tmp3.a[0]);} else if (r1) {ans=tmp2.c[1]||(tmp1.a[1]&&tmp2.b[0])||(tmp2.b[1]&&tmp3.a[0]);ans = ans || (tmp1.a[1]&&tmp2.c[0]&&tmp3.a[0]);} else {ans=tmp2.c[0]||(tmp1.a[1]&&tmp2.b[1])||(tmp2.b[0]&&tmp3.a[0]);ans = ans || (tmp1.a[1]&&tmp2.c[1]&&tmp3.a[0]);}if (ans) puts("Y"); else puts("N");}}return 0;}

ps：这道题教练给的标程有问题。反对网上随便找题解！！！