HDU 3594 Cactus 【强连通】

题目：点击打开链接

题意：求一个图是否为仙人掌图

分析：

仙人掌图的详细说明见链接：点击打开链接

总结下来仙人掌图有三点性质：

①没有横向边

②对于父节点u和其所有子节点v，low[v]<=dfn[u]

③对于父节点u和其所有子节点v，low[v]<dfn[u]的个数和u的反向边的个数之和不超过1

对于性质①，可以在dfs时，对于每次引出的一个v，判断dfn[v]>dfn[u]是否成立，若成立，则表示点v是之前已经遍历过的一个点且v是u的一个后继节点，因为u可以通过其他路径到达v，则u->v这条边是一条横向边。

对于性质②，只要判断强连通分支个数是否为一

对于性质③，开两个数组分别保存对于父节点u和其所有子节点v，low[v]<dfn[u]的个数和u的反向边的个数。容易在计算low[u]的值是累加low[v]<dfn[u]的个数，当遍历到的v点在栈内，说明u->v是一条反向边，累加u的反向边的个数。

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <stdio.h>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
#include <set>
#include <vector>
#include <map>
#define PR pair<int,int>
#define MP make_pair
#define fi first
#define se second
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define sqr(x) ((x)*(x))
#define ll __int64
const ll INF=1e18;
const int inf=0x3f3f3f3f;
const int M=2010;
const int N=20010;
const int MOD=1000000007;
const double eps=1e-6;
const double pi=acos(-1.0);
using namespace std;
int n,m,u,v;
int index,scc_cnt,dfn[N],low[N],instk[N],sum_low[N],rev[N],ok;
vector<int> edge[N];
stack<int> stk;
void init()
{
    ok=1;
    scc_cnt=index=0;
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(edge,0,sizeof(edge));
    memset(rev,0,sizeof(rev));
    memset(sum_low,0,sizeof(sum_low));
}
void tarjan(int u)
{
    dfn[u]=low[u]=++index;
    stk.push(u);instk[u]=1;
    for(int i=0;i<edge[u].size();i++)
    {
        int v=edge[u][i];
        if(dfn[v]>dfn[u]) ok=0;
        if(!dfn[v])
        {
            tarjan(v);
            if(low[u]>low[v]) sum_low[u]++,low[u]=low[v];
        }
        else if(instk[v]) rev[u]++,low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u])
    {
        scc_cnt++;
        while(1)
        {
            int now=stk.top();
            stk.pop();
            instk[now]=0;
            if(now==u) break;
        }
    }
}
int main()
{
    int i,j,T;
    scanf("%d",&T);
    while(T--)
    {
        init();
        scanf("%d",&n);
        while(~scanf("%d%d",&u,&v)&&(u||v)) edge[u].push_back(v);
        for(i=0;i<n;i++) if(!dfn[i]) tarjan(i);
        for(i=0;i<n;i++) if(sum_low[i]+rev[i]>1) {ok=0;break;}
        if(scc_cnt==1&&ok==1) puts("YES");
        else puts("NO");
    }
}