有向图强连通分量Tarjan模板

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn=10010;
const int maxm=1000010;
int cnt,p[maxn];
struct node
{
    int v,next;
}E[maxm];
int n,m;
int low[maxn],dfn[maxn];
int Stack[maxn];///数组模拟栈
int num[maxn];///记录第i个连通块点的个数
int Scc[maxn];///第i个点属于哪个连通块
int scc;///连通块计数
int Time;///dfs的时间戳
bool instack[maxn];///点是否在栈里
int top;///栈里点的个数
void init()
{
    cnt=0;
    memset(p,-1,sizeof(p));
}
void add(int u,int v)
{
    E[cnt].v=v;
    E[cnt].next=p[u];
    p[u]=cnt++;
}
void Tarjan(int u)
{
    dfn[u]=low[u]=++Time;
    Stack[++top]=u;
    instack[u]=1;
    int v;
    for(int i=p[u];i!=-1;i=E[i].next)
    {
        v=E[i].v;
        if(!dfn[v])///由于dfn初始化为0，所以可以判断是否访问
        {
            Tarjan(v);
            if(low[u]>low[v]) low[u]=low[v];
        }
        else if(instack[v]&&low[u]>dfn[v]) low[u]=dfn[v];
    }
    if(low[u]==dfn[u])
    {
        scc++;
        while(1)
        {
            v=Stack[top--];
            instack[v]=0;
            num[scc]++;
            Scc[v]=scc;
            if(u==v) break;
        }
    }
}
void solve()
{
    for(int i=1;i<=n+1;i++)
    {
        dfn[i]=num[i]=Scc[i]=0;
        instack[i]=0;
    }
    Time=scc=top=0;
    for(int i=1;i<=n;i++)
        if(!dfn[i]) Tarjan(i);
}
int main()
{
    while(scanf("%d%d",&n,&m),n||m)
    {
        init();
        for(int i=1;i<=m;i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            add(x,y);
        }
        solve();
        if(scc==1) puts("Yes");
        else puts("No");
    }
    return 0;
}