bzoj 1018: [SHOI2008]堵塞的交通traffic
由于只有两行，对一个区间可以维护四个角上的连通性 一共6个量，满足区间可加性，用线段树维护，查询时找出3段的量，看是否满足题解上的三种情况之一。
1 #include<cstdio> 2 #include<iostream> 3 #include<cstring> 4 #include<cstdlib> 5 #include<cmath> 6 #include<queue> 7 #include<algorithm> 8 #include<vector> 9 #define M 1000009 10 #define EPS 1e-10 11 #define MO 19650827 12 #define ll long long 13 using namespace std; 14 ll read() 15 { 16 char ch=getchar(); 17 ll x=0,f=1; 18 for(;ch<'0'||ch>'9';ch=getchar()) 19 if(ch=='-') 20 f=-1; 21 for(;ch>='0'&&ch<='9';ch=getchar()) 22 x=x*10+ch-'0'; 23 return x*f; 24 } 25 int n,mp[M]; 26 struct data 27 { 28 bool f[6]; 29 }s[10],a[M]; 30 void build(int x,int l,int r) 31 { 32 if(l==r) 33 { 34 a[x]=s[0]; 35 return; 36 } 37 int mid=(l+r)>>1; 38 build(x*2,l,mid); 39 build(x*2+1,mid+1,r); 40 } 41 int calc(int x,int y) 42 { 43 return x*(n-1)+y; 44 } 45 data he(data a,data b,int x,int y) 46 { 47 data c; 48 c.f[0]=((a.f[0]&&x&&b.f[0])||(a.f[4]&&y&&b.f[5])); 49 c.f[1]=((a.f[1]&&y&&b.f[1])||(a.f[5]&&x&&b.f[4])); 50 c.f[2]=(a.f[2]||(a.f[0]&&x&&b.f[2]&&y&&a.f[1])); 51 c.f[3]=(b.f[3]||(b.f[0]&&x&&a.f[3]&&y&&b.f[1])); 52 c.f[4]=((a.f[4]&&y&&b.f[1])||(a.f[0]&&x&&b.f[4])); 53 c.f[5]=((a.f[5]&&x&&b.f[0])||(a.f[1]&&y&&b.f[5])); 54 return c; 55 } 56 void gai(int x,int l,int r,int x1,int y1,int x2,int c) 57 { 58 int mid=(l+r)>>1; 59 if(x1==x2&&mid==y1) 60 { 61 mp[calc(x1,y1)]=c; 62 a[x]=he(a[x*2],a[x*2+1],mp[calc(0,y1)],mp[calc(1,y1)]); 63 return; 64 } 65 if(l==r) 66 { 67 mp[calc(2,y1)]=c; 68 a[x]=s[c]; 69 return; 70 } 71 if(y1<=mid) 72 gai(x*2,l,mid,x1,y1,x2,c); 73 else 74 gai(x*2+1,mid+1,r,x1,y1,x2,c); 75 a[x]=he(a[x*2],a[x*2+1],mp[calc(0,mid)],mp[calc(1,mid)]); 76 } 77 data query(int x,int l,int r,int L,int R) 78 { 79 int mid=(l+r)>>1; 80 if(L<=l&&R>=r) 81 return a[x]; 82 if(R<=mid) 83 return query(x*2,l,mid,L,R); 84 if(L>mid) 85 return query(x*2+1,mid+1,r,L,R); 86 return he(query(x*2,l,mid,L,R),query(x*2+1,mid+1,r,L,R),mp[calc(0,mid)],mp[calc(1,mid)]); 87 } 88 void xun(int x1,int y1,int x2,int y2) 89 { 90 s[2]=query(1,1,n,1,y1); 91 s[3]=query(1,1,n,y1,y2); 92 s[4]=query(1,1,n,y2,n); 93 int ans; 94 if(x1==x2) 95 ans=(s[3].f[x1]||(s[2].f[3]&&s[3].f[4+x1^1])||(s[3].f[4+x1]&&s[4].f[2])||(s[2].f[3]&&s[3].f[x1^1]&&s[4].f[2])); 96 else 97 ans=(s[3].f[4+x1]||(s[2].f[3]&&s[3].f[x1^1])||(s[3].f[x1]&&s[4].f[2])||(s[2].f[3]&&s[3].f[4+x1^1]&&s[4].f[2])); 98 if(ans) 99 printf("Y\n"); 100 else 101 printf("N\n"); 102 } 103 int main() 104 { 105 s[0]=(data){1,1,0,0,0,0}; 106 s[1]=(data){1,1,1,1,1,1}; 107 n=read(); 108 build(1,1,n); 109 for(;;) 110 { 111 char ch[10]; 112 scanf("%s",ch); 113 if(ch[0]=='E') 114 return 0; 115 int x1=read(),y1=read(),x2=read(),y2=read(); 116 x1--; 117 x2--; 118 if(y1>y2) 119 { 120 swap(y1,y2); 121 swap(x1,x2); 122 } 123 if(ch[0]=='O') 124 gai(1,1,n,x1,y1,x2,1); 125 if(ch[0]=='C') 126 gai(1,1,n,x1,y1,x2,0); 127 if(ch[0]=='A') 128 xun(x1,y1,x2,y2); 129 } 130 return 0; 131 } 132
