CF555E Case of Computer Network

题面:https://www.luogu.com.cn/problem/CF555E
题意:给定一张(n)个点(m)条边的无向图。
给定(q)组有向点对((s,t))
询问是否存在使得所有(s)都能到达(t)的无向图中每条边的定向方案。
n,m,q (leq) 2e5
题解:
看到这种关于无向图连通性的问题,我们可以想到边双。
如果两点在边双中,我们可以将边双里的边定向为一些环,
这样这些点一定是可以到达的。
考虑边双缩点,将整张图缩成一个森林。
接下来,对于每对((s,t)),我们树上差分记录一下
每条边的方向,最后DFS判断一下就行了。
时间复杂度(O(n+m))
代码:

#include<bits/stdc++.h>
using namespace std;
#define re register int
#define F(x,y,z) for(re x=y;x<=z;x++)
#define FOR(x,y,z) for(re x=y;x>=z;x--)
typedef long long ll;
#define I inline void
#define IN inline int
#define STS system("pause")
template<class D>I read(D &res){
	res=0;register D g=1;register char ch=getchar();
	while(!isdigit(ch)){
		if(ch=='-')g=-1;
		ch=getchar();
	}
	while(isdigit(ch)){
		res=(res<<3)+(res<<1)+(ch^48);
		ch=getchar();
	}
	res*=g;
}
struct E{
	int to,nt;
}e[404000];
#define T e[k].to
int n,m,Q,tot,X,Y,sn,head[202000];
int dfn[202000],low[202000],clr[202000],color,mill,vis[202000];
int dep[202000],fa[202000],f[202000][20],lg[202000],top[202000];
int up[202000],down[202000];
stack<int>s;
vector<int>v[202000];
vector<int>::iterator it;
I tarjan(int x,int pre){
	dfn[x]=low[x]=++mill;
	s.push(x);vis[x]=1;
	for(re k=head[x];k!=-1;k=e[k].nt){
		if(k==(pre^1))continue;
		if(!dfn[T]){
			tarjan(T,k);
			low[x]=min(low[x],low[T]);
		}
		else if(vis[T])low[x]=min(low[x],dfn[T]);
	}
	if(dfn[x]==low[x]){
		++color;
		while(s.top()!=x){clr[s.top()]=color;vis[s.top()]=0;s.pop();}
		clr[x]=color;vis[x]=0;s.pop();
	}
}
I D_1(int x,int fat,int depth,int topi){
	//cout<<x<<":";
	f[x][0]=fat;top[x]=topi;dep[x]=depth;
	F(i,1,lg[dep[x]])f[x][i]=f[f[x][i-1]][i-1];
	for(auto k:v[x]){
		//cout<<k<<" ";
		if(k==fat)continue;
		D_1(k,x,depth+1,topi);
	}
	//cout<<endl;
}
I add(int x,int y){v[x].emplace_back(y);}
IN ques_lca(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	re len=dep[x]-dep[y];
	FOR(i,lg[len],0){
		if((len>>i)&1)x=f[x][i];
	}
	if(x==y)return x;
	FOR(i,lg[dep[x]],0){
		if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
	}
	return f[x][0];
}
I D_3(int x,int fat){
	if(!sn)return;
	for(auto k:v[x]){
		if(k==fat)continue;
		D_3(k,x);
		if(!sn)return;
		up[x]+=up[k];down[x]+=down[k];
	}
	//cout<<x<<" "<<up[x]<<" "<<down[x]<<endl;
	if(up[x]&&down[x])sn=0;
}
int main(){
	read(n);read(m);read(Q);tot=-1;
	memset(head,-1,sizeof(head));
	F(i,1,m){
		read(X);read(Y);
		e[++tot].to=Y;
		e[tot].nt=head[X];
		head[X]=tot;
		e[++tot].to=X;
		e[tot].nt=head[Y];
		head[Y]=tot;
	}
	F(i,1,n){
		if(!dfn[i])tarjan(i,-1);
	}
	F(i,1,n){
		for(re k=head[i];k!=-1;k=e[k].nt){
			if(clr[i]!=clr[T])add(clr[i],clr[T]);
		}
	}	
	lg[0]=-1;
	F(i,1,n)lg[i]=lg[i>>1]+1,top[i]=i;
	F(i,1,color){
		sort(v[i].begin(),v[i].end());
		it=unique(v[i].begin(),v[i].end());
		v[i].erase(it,v[i].end());
		//cout<<i<<":";
		//for(auto k:v[i])cout<<k<<" ";
		//cout<<endl;
	}
	F(i,1,color){
		if(top[i]==i)dep[i]=1,D_1(i,0,1,i);
	}//STS;
	re P;sn=1;
	while(Q--){
		read(X);read(Y);X=clr[X];Y=clr[Y];
		if(top[X]!=top[Y]){cout<<"No";return 0;}
		P=ques_lca(X,Y);//cout<<X<<" "<<Y<<" "<<P<<endl;
		up[X]++;down[Y]++;up[P]--;down[P]--;
	}
	//F(i,1,color)cout<<up[i]<<" "<<down[i]<<endl;
	F(i,1,color){
		if(top[i]==i)D_3(i,0);
		if(!sn)break;
	}
	if(sn)cout<<"Yes";
	else cout<<"No";
	return 0;
}
/*
8 13 4
5 3
2 7
5 8
2 4
6 4
2 6
5 8
1 4
6 7
3 1
6 7
3 1
8 3
8 7
8 7
5 1
6 8
*/
原文地址:https://www.cnblogs.com/Purple-wzy/p/12091897.html