[图论多解法]hdoj 3342：Legal or Not

大致题意：

    很简单，就是给出一个有向图，判断这个图是否含有环。

 

大致思路：

    正解应该是拓扑排序判断环，不过Tarjan求强联通分量，floyd也可以解。

实验室快关门了，先贴上Tarjan的～～

 

/*
    Tarjan算法求强连通分量
*/
#include<iostream>
#include<cstdio>
#include <algorithm>
#include<cstring>
using namespace std;
const int inf=1<<30;
const int nMax=3015;
const int mMax=5001000;
class edge{
public:Legal or Not
    int v,nex;
};edge e[mMax];
int k,head[nMax];//head[i]是以点i为起点的链表头部

void addedge(int a,int b){//向图中加边的算法，注意加上的是有向边//b为a的后续节点既是a---->b
    e[k].v=b;
    e[k].nex=head[a];
    head[a]=k;k++;
}

int dfn[nMax],low[nMax],sta[nMax],top,atype,belon[nMax],dep;   //atype 强连通分量的个数
bool insta[nMax];

void Tarjan(int u){                 //我的Tarjan模版
    int i,j;
    dfn[u]=low[u]=++dep;
    sta[++top]=u;
    insta[u]=1;
    for(i=head[u];i;i=e[i].nex){
        int v=e[i].v;
        if(!dfn[v]){
            Tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else{
            if(insta[v])low[u]=min(low[u],dfn[v]);
        }
    }
    if(dfn[u]==low[u]){
        atype++;              //强连通分量个数
        do{
            j=sta[top--];
            belon[j]=atype;   //第j个点属于第type个连通块
            insta[j]=0;
        }while(u!=j);
    }
}

void init(){
    k=1;
    dep=1;
    top=atype=0;
    memset(insta,0,sizeof(insta)); //是否在栈中
    memset(head,0,sizeof(head));   //静态链表头指针
    memset(low,0,sizeof(low));     //Tarjan的low数组
    memset(dfn,0,sizeof(dfn));     //Tarjan的dfn数组
    memset(belon,0,sizeof(belon)); //记录每个点属于哪一个强连通分量
}

int main(){
    int n,m,a,b,i,j,flag;
    while(cin>>n>>m&&n)
    {
        flag=1;
        init();
        for(i=0;i<m;i++)
        {
            cin>>a>>b;
            if(a==b)flag=0;
            addedge(a,b);
        }
        for(i=0;i<n;i++)
        {
            if(!dfn[i])Tarjan(i);
        }
        if(atype!=n)flag=0;
        if(flag)printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}