题目:http://www.lydsy.com/JudgeOnline/problem.php?id=1018
分析:首先我们能想到要从c1列到c2列,必须要先从c1列到k列(c1<=k<=c2), 再从k列到c2, 于是问题变为c1->k是否可达,k->c2是否可达。 所以我们想到用线段树来处理区间的可达的情况,对于一个区间,我们要考虑左上角到右上角的可达情况,左上角到右下角的可达情况,左下角到右上角,左下角到右下角, 可能你想到这里就觉得够了,很可惜这样是不够的。 因为假设我们要问的区间是L的左上角,到R的右下角, 但是我们可能会不得不先往左走,然后绕到L的左下角,然后再到R的右下角。 同样也可能L到R的右上角后,继续往右,再绕回来R的右下角......所以我们对于一个区间我们还要维护左上角和左下角的连通性(只需要在所在区间内走,不需要考虑走出去的)同样维护右上角和右下角的。 然后是区间的合并。。。。分类讨论吧。。。应该能分清。 更新的处理没什么好说的,线段书从叶子(1个或者是2个)开始往上更新就行了。
重头戏是询问连通性。。。假设询问的是[L,R] ,我们先把[1,L]的区间从线段树中抽出来,这是用来判断左上角和左下角连通性的,然后用一个值表示哪个点能作为起点往右走,同样把[R,n]的区间按线段树分离出来,然后看右上角和右下角连通性,用一个值表示哪个点能作为终点。再把[L,R]的区间分出来,然后从左边走到右边,看能不能走到该走的点。
总的来说,写起来很麻烦,因为很多细节可能会因为脑子不清晰而写错。。。
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <stack>
#include <string.h>
#include <queue>
#include <vector>
#include <algorithm>
#include <cassert>
#include <set>
#include <map>
#include <cmath>
#include <ctime>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<(b);++i)
#define rrep(i,a,b) for(int i=(a);i>=(b);--i)
#define clr(a,x) memset(a,(x),sizeof(a))
#define eps 1e-9
#define LL long long
#define mp make_pair
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
const int maxn=100000+5;
bool lk[maxn<<2][2][2];
bool lr[maxn<<2][2];
bool open[maxn][3];
void PushUp(int l,int m,int r,int rt)
{
int ls=rt<<1,rs=rt<<1|1;
rep(i,0,2) rep(j,0,2) {
lk[rt][i][j]=false;
rep(k,0,2) {
lk[rt][i][j] |= lk[ls][i][k] && open[m][k+1] && lk[rs][k][j];
if(lk[rt][i][j]) break;
}
}
lr[rt][0] = lr[rt][1] = false;
if(open[l][0]) {
lr[rt][0] = true;
rep(i,0,2) rep(j,0,2)
lk[rt][i][j] |= lk[rt][i^1][j];
}
if(open[r][0]) {
lr[rt][1] = true;
rep(i,0,2) rep(j,0,2)
lk[rt][i][j] |= lk[rt][i][j^1];
}
lr[rt][0] |= lr[ls][0] || (lr[rs][0] && open[m][1] && open[m][2] && ((lk[ls][0][0]&&lk[ls][1][1]) || (lk[ls][0][1] && lk[ls][1][0])));
lr[rt][1] |= lr[rs][1] || (lr[ls][1] && open[m][1] && open[m][2] && ((lk[rs][0][0]&&lk[rs][1][1]) || (lk[rs][0][1] && lk[rs][1][0])));
}
void build(int l,int r,int rt)
{
if(l==r) {
rep(i,0,2) rep(j,0,2) lk[rt][i][j]=false;
rep(i,0,2) lk[rt][i][i]=true;
lr[rt][0]=lr[rt][1]=false;
return;
}
int m=(l+r)>>1;
build(lson); build(rson);
PushUp(l,m,r,rt);
}
void update(int L,int R,int l,int r,int rt)
{
if(l==r) {
lr[rt][0]=lr[rt][1]=open[l][0];
rep(i,0,2) lk[rt][i][i^1]=open[l][0];
return;
}
int m=(l+r)>>1;
if(L<=m) update(L,R,lson);
if(m<R) update(L,R,rson);
PushUp(l,m,r,rt);
}
int n;
char op[10];
void DoUpdate(bool on)
{
int r1,r2,c1,c2;
scanf("%d%d%d%d",&r1,&c1,&r2,&c2);
if(r1>r2) swap(r1,r2),swap(c1,c2);
else if(r1==r2&&c1>c2) swap(c1,c2);
if(r1!=r2) open[c1][0]=on;
else open[c1][r1]=on;
update(c1,c2,1,n,1);
}
struct Segment
{
int l,r;
int lk[2][2];
int lr[2];
}sgt[500];
int tot;
void GetSegment(int L,int R,int l,int r,int rt)
{
if(L<=l&&r<=R) {
Segment & st = sgt[tot];
rep(i,0,2) rep(j,0,2) st.lk[i][j]=lk[rt][i][j];
st.l=l,st.r=r;
rep(i,0,2) st.lr[i]=lr[rt][i];
++tot;
return;
}
int m=(l+r)>>1;
if(L<=m) GetSegment(L,R,lson);
if(m<R) GetSegment(L,R,rson);
}
bool DoQuery()
{
int r1,c1,r2,c2; scanf("%d%d%d%d",&r1,&c1,&r2,&c2);
if(r1==r2&&c1==c2) return true;
if(c1>c2) swap(r1,r2),swap(c1,c2);
else if(c1==c2&&r1>r2) swap(r1,r2);
tot = 0;
GetSegment(1,c1,1,n,1);
bool t1=false;
rrep(i,tot-1,0) {
if(sgt[i].lr[1]) { t1=true; break; }
if(i==0) break;
int l=sgt[i].l-1;
if(!open[l][1]||!open[l][2]) break;
if(sgt[i].lk[0][0]&&sgt[i].lk[1][1]) continue;
if(sgt[i].lk[0][1]&&sgt[i].lk[1][0]) continue;
break;
}
tot = 0;
GetSegment(c2,n,1,n,1);
bool t2=false;
rep(i,0,tot) {
if(sgt[i].lr[0]) { t2=true; break; }
if(i==tot-1) break;
int r=sgt[i].r;
if(!open[r][1]||!open[r][2]) break;
if(sgt[i].lk[0][0]&&sgt[i].lk[1][1]) continue;
if(sgt[i].lk[0][1]&&sgt[i].lk[1][0]) continue;
break;
}
int s1=1<<(r1-1),s2=1<<(r2-1);
if(t1) s1=3; if(t2) s2=3;
tot = 0;
GetSegment(c1,c2,1,n,1);
rep(i,0,tot) {
int ret=0;
rep(j,0,2) {
if((s1&(1<<j))==0) continue;
rep(k,0,2) ret |= (sgt[i].lk[j][k]<<k);
}
s1=ret;
if(i<tot-1) {
int r=sgt[i].r;
if(open[r][1]&&open[r][2]) s1 &= 3;
else if(open[r][1]) s1 &= 1;
else if(open[r][2]) s1 &= 2;
else s1 &= 0;
}
if(s1==0) return false;
}
return (s1&s2)>0;
}
void solve()
{
while(scanf("%s",op)==1) {
if(op[0]=='E') break;
if(op[0]=='O') DoUpdate(true);
else if(op[0]=='C') DoUpdate(false);
else puts(DoQuery()?"Y":"N");
}
}
int main()
{
// freopen("in.txt","r",stdin);
while(scanf("%d",&n)==1) {
clr(open,0); clr(lk,0); clr(lr,0);
build(1,n,1);
solve();
}
}