Codeforces 521E

考虑图中的一个连通块，对于它的一棵dfs树，显然必须要有返祖边。
我们容易用tarjan算法证明，如果不存在某条树边被两条返祖边覆盖的话一定无解。反之我们选出两条在树上覆盖部分有交的返祖边，取出它们对应的树边，形成了两个相交的环，我们取相交部分的首尾为端点，就得到了三条路径。
时间复杂度$\mathcal{O}(n)$。

#include <bits/stdc++.h>
#define FR first
#define SE second
#define end end2
 
using namespace std;
 
typedef pair<int,int> pr;
 
vector <int> e[200005];
 
int fa[200005],dep[200005];
int d[200005],dfn[200005],end[200005],dfs_cnt;
 
bool vis[200005];
pr num[200005];
int e_cnt;
 
void dfs(int x) {
  dfn[x]=++dfs_cnt;
  vis[x]=1;
  for(int i=0;i<e[x].size();i++)
    if (e[x][i]!=fa[x]) {
    	int u=e[x][i];
    	if (!vis[u]) {
    		fa[u]=x;dep[u]=dep[x]+1;
    		dfs(u);
    		d[x]+=d[u];
		}
		else if (dep[u]<dep[x]) {
			d[x]++;d[u]--;
			num[++e_cnt]=pr(x,u);
		}
	}
  end[x]=dfs_cnt;
}
 
vector <int> vt1,vt2;
int col[200005];
 
vector <int> ans1,ans2,ans3;
 
void solve(int x,int n) {
  int id1=0,id2=0;
  for(int i=1;i<=e_cnt;i++)
    if (dfn[num[i].FR]>=dfn[x]&&dfn[num[i].FR]<=end[x]&&dfn[num[i].SE]<dfn[x]) {
    	if (!id1) id1=i;
    	else {
    		id2=i;
    		break;
		}
	}
  for(int i=num[id1].FR;i!=num[id1].SE;i=fa[i]) vt1.push_back(i);
  vt1.push_back(num[id1].SE);
  for(int i=num[id2].FR;i!=num[id2].SE;i=fa[i]) vt2.push_back(i);
  vt2.push_back(num[id2].SE);
  for(int i=0;i<vt1.size();i++) col[vt1[i]]++;
  for(int i=0;i<vt2.size();i++) col[vt2[i]]++;
  int u=0,v=0;
  for(int i=1;i<=n;i++)
    if (col[i]==2) {
    	if (!u||dep[u]<dep[i]) u=i;
    	if (!v||dep[v]>dep[i]) v=i;
	}
  for(int i=u;i!=v;i=fa[i]) ans1.push_back(i);
  ans1.push_back(v);
  int l=vt1.size()-1;
  while (vt1[l]!=u) l--;
  for(int i=l;i>=0;i--) ans2.push_back(vt1[i]);
  int r=0;
  while (vt1[r]!=v) r++;
  for(int i=vt1.size()-1;i>=r;i--) ans2.push_back(vt1[i]);
  l=vt2.size()-1;
  while (vt2[l]!=u) l--;
  for(int i=l;i>=0;i--) ans3.push_back(vt2[i]);
  r=0;
  while (vt2[r]!=v) r++;
  for(int i=vt2.size()-1;i>=r;i--) ans3.push_back(vt2[i]);
}
 
int main() {
  int n,m;
  scanf("%d%d",&n,&m);
  for(int i=1;i<=m;i++) {
  	int x,y;
  	scanf("%d%d",&x,&y);
  	e[x].push_back(y);
  	e[y].push_back(x);
  }
  for(int i=1;i<=n;i++)
    if (!vis[i]) dfs(i);
  for(int i=1;i<=n;i++)
    if (d[i]>=2) {
    	solve(i,n);
    	puts("YES");
    	printf("%d ",(int)ans1.size());
    	for(int j=0;j<ans1.size();j++) printf("%d ",ans1[j]);
    	printf("\n");
    	printf("%d ",(int)ans2.size());
    	for(int j=0;j<ans2.size();j++) printf("%d ",ans2[j]);
    	printf("\n");
    	printf("%d ",(int)ans3.size());
    	for(int j=0;j<ans3.size();j++) printf("%d ",ans3[j]);
    	printf("\n");
    	return 0;
	}
  puts("NO");
  return 0;
}