HDU-1296 迷宫城堡（Tarjan）

题意

一个 $n$ 个节点， $m$ 条边的单向图，求是不是强连通。
$1 \leq n \leq 10000$
$1 \leq m \leq 100000$

思路

一张强连通图的定义是任意两个节点都可以互相到达。这里使用代码量短，流程清晰的 $Tarjan$ 算法，具体流程如下。

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$ 数组。它的含义是子树中通过一条反祖边能到达 $dfs$ 序最小的点的编号。通过维护这个数组，就能最大范围的找到一块强连通块并记录。其中 $dfn,bel$ 数组要保证初值为 $0$ ，$B,ord,tp$ 的值显然也要赋为零。

代码

#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;
}