[IOI2018] werewolf 狼人 - 解题报告

题意：
有一张无向图，一个人他需要从s走到t。
他有两种形态，第一种形态可以走点编号为 \([ l_i，n ]\) ，第二种形态可以走点编号为 \([ 1，r_i ]\)，可以点编号为 \([ l_i，n ]\)切换形态（恰好一次），在起点时为第一种形态。求他是否能从s走到t。
多组询问。

题解：
首先我们可以维护与起点联通的点集合（形态一），同时也可以维护维护与终点联通的点集合（形态二），显然可以用Kruskal重构树维护。
用主席树维护两棵重构树其中子树是否有并集（二位数点）。

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=400010;
int n,m,q,fa[N],w[N];
vector <int> e[N];
int tot=0,rt[N],lf[N*25],rf[N*25],sum[N*25];

inline int read() {
	register int tmp=0;register char c=getchar(); while(c<'0'||c>'9')	c=getchar();
	while(c>='0'&&c<='9')	tmp=(tmp<<1)+(tmp<<3)+(c^48),c=getchar();return tmp;
}
//inline int Max(int x,int y) { return x>y?	x:y; }
//inline int Min(int x,int y) { return x<y?	x:y; }
int Find(int x) { return x==fa[x]?	x:fa[x]=Find(fa[x]); }
struct Kruskal {
	int type;
	int cnt,hed[N],to[N],nxt[N];
	int tot,w[N],idx,dfn[N],bg[N],ed[N];
	int st[N][20];
	inline void add(int x,int y) { to[++cnt]=y,nxt[cnt]=hed[x],hed[x]=cnt; }
	void dfs(int u) {
		dfn[bg[u]=++idx]=u; for(int i=1;i<20;i++)	st[u][i]=st[st[u][i-1]][i-1];
		for(int i=hed[u];i;i=nxt[i])	dfs(to[i]);
		ed[u]=idx;
	}
	void build() {
		tot=n; for(int i=1;i<=n;++i)	fa[i]=i;
		if(type) {
			for(int i=n;i;--i)
				for(int j=0;j<e[i].size();++j)
					if(i<e[i][j]) {
						int v=Find(e[i][j]); if(i==v)	continue;
						add(i,v);fa[v]=st[v][0]=i;
					}
			dfs(1);
		}
		else {
			for(int i=1;i<=n;++i)
				for(int j=0;j<e[i].size();++j)
					if(i>e[i][j]) {
						int v=Find(e[i][j]); if(i==v)	continue;
						add(i,v);fa[v]=st[v][0]=i;
					}
			dfs(n);
		}
	}
	inline int find(int u,int x) {
		for(int i=19;i>=0;i--)
			if(st[u][i]&&((type&&x<=st[u][i])||(!type&&x>=st[u][i])))
				u=st[u][i];
		return u;
	}
};Kruskal mn,mx;// mn:people   mx:wolf

void update(int &u,int v,int l,int r,int x) {
	sum[u=++tot]=sum[v]+1,lf[u]=lf[v],rf[u]=rf[v]; if(l>=r)	return ;
	int mid=(l+r)>>1; x<=mid?	update(lf[u],lf[v],l,mid,x):update(rf[u],rf[v],mid+1,r,x);
}
int ask(int u,int v,int l,int r,int L,int R) {
	if(l==L&&r==R)	return sum[v]-sum[u];
	int mid=(l+r)>>1;
	if(R<=mid)	return ask(lf[u],lf[v],l,mid,L,R);
	if(L>mid)	return ask(rf[u],rf[v],mid+1,r,L,R);
	return ask(lf[u],lf[v],l,mid,L,mid)+ask(rf[u],rf[v],mid+1,r,mid+1,R);
}
int main() {
	n=read(),m=read(),q=read();
	for(int i=1,x,y;i<=m;++i)	x=read()+1,y=read()+1,e[x].push_back(y),e[y].push_back(x);
	mn.type=1,mn.build(),mx.build();
	for(int i=1;i<=n;++i)	w[mn.bg[i]]=mx.bg[i];
	for(int i=1;i<=n;++i)	update(rt[i],rt[i-1],1,n,w[i]);
	for(;q;--q) {
		int s=read()+1,t=read()+1,l=read()+1,r=read()+1;
		s=mn.find(s,l),t=mx.find(t,r);
		if(ask(rt[mn.bg[s]-1],rt[mn.ed[s]],1,n,mx.bg[t],mx.ed[t]))	printf("1\n");
		else	printf("0\n");
	}
	return 0;
}