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