题意
一个
n
n
个节点, 条边的单向图,求是不是强连通。
1≤n≤10000
1
≤
n
≤
10000
1≤m≤100000
1
≤
m
≤
100000
思路
一张强连通图的定义是任意两个节点都可以互相到达。这里使用代码量短,流程清晰的 Tarjan T a r j a n 算法,具体流程如下。
void tarjan(int u)
{
dfn[u]=low[u]=++ord;stk[++tp]=u; //初始化
EOR(i,G,u)
{
int v=G.to[i];
if(!dfn[v]) //没遍历过的点
{
tarjan(v);
low[u]=min(low[u],low[v]);
}else if(!bel[v])low[u]=min(low[u],dfn[v]); //在栈里的点
}
if(dfn[u]==low[u]) //找到一个强连通分量
{
int v;B++;
do bel[v=stk[tp--]];
while(u!=v);
}
}
最难理解的是 low l o w 数组。它的含义是子树中通过一条反祖边能到达 dfs d f s 序最小的点的编号。通过维护这个数组,就能最大范围的找到一块强连通块并记录。其中 dfn,bel d f n , b e l 数组要保证初值为 0 0 , 的值显然也要赋为零。
代码
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<stack>
#define FOR(i,x,y) for(int i=(x);i<=(y);i++)
#define DOR(i,x,y) for(int i=(x);i>=(y);i--)
#define lowbit(x) ((x)&-(x))
#define N 100003
typedef long long LL;
using namespace std;
template<const int maxn,const int maxm>struct Linked_list
{
int head[maxn],to[maxm],nxt[maxm],tot;
void clear(){memset(head,-1,sizeof(head));tot=0;}
void add(int u,int v){to[++tot]=v,nxt[tot]=head[u];head[u]=tot;}
#define EOR(i,G,u) for(int i=G.head[u];~i;i=G.nxt[i])
};
Linked_list<N,N>G;
stack<int>stk;
int dfn[N],belong[N],low[N],ord,ans,bcnt,n,m;
void tarjan(int u)
{
dfn[u]=low[u]=++ord;
stk.push(u);
EOR(i,G,u)
{
int v=G.to[i];
if(!dfn[v])
{
tarjan(v);
low[u]=min(low[u],low[v]);
}
else if(!belong[v])
low[u]=min(low[u],dfn[v]);
}
if(dfn[u]==low[u])
{
bcnt++;
int cnt=0,v;
do
{
v=stk.top();stk.pop();
belong[v]=bcnt;
cnt++;
}while(u!=v);
ans=max(ans,cnt);
}
}
void clear()
{
memset(dfn,0,sizeof(dfn));
memset(belong,0,sizeof(belong));
G.clear();
ord=ans=0;
}
int main()
{
while(scanf("%d%d",&n,&m),n||m)
{
clear();
FOR(i,1,m)
{
int u,v;
scanf("%d%d",&u,&v);
G.add(u,v);
}
FOR(i,1,n)if(!dfn[i])tarjan(i);
printf(ans==n?"Yes\n":"No\n");
}
return 0;
}