POJ 2762 弱联通

<pre name="code" class="cpp">//晕死，随便敲也能过，这是真心数据弱吧。。。。。

//poj 2762
//判断有向图的弱连通
//【题目大意】
//给出T个山洞，对于每个山洞，如果任意选择两点s,e，都满足s可以到达e或者e可以到达s，则输出Yes,否则输出No。
//
//【解题思路】
//这题实际上是判断是否是弱连通的，所以首先强连通，然后缩点，对缩点形成的图最多只能有一个入度为0的点，
//如果有多个入度为0的点，则这两个连通分量肯定是不连通的。
//缩点后形成的图形是一棵树，入度为0的点是这颗树的根，这棵树只能是单链，不能有分叉，如果有分叉，则这些分叉之间是不可达的，所以就对这棵树进行DFS，如果是单链则是YES。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<queue>
#define maxn 1050
using namespace std;

int n, m;//n m 为点数和边数

struct Edge
{
    int from, to, nex;
    bool sign;//是否为桥
}edge[10000];


int head[maxn], edgenum;
bool visit[maxn];
int num;
bool judge;

void addedge(int u, int v)
{//边的起点和终点
    Edge E={u, v, head[u], false};
    edge[edgenum] = E;
    head[u] = edgenum++;
}

int DFN[maxn], Low[maxn], Stack[maxn], top, Time; //Low[u]是点集{u点及以u点为根的子树} 中(所有反向弧)能指向的(离根最近的祖先v) 的DFN[v]值(即v点时间戳)
int taj;//连通分支标号，从1开始
int Belong[maxn];//Belong[i] 表示i点属于的连通分支
bool Instack[maxn];
vector<int> bcc[maxn]; //标号从1开始

void tarjan(int u ,int fa)//求出图的强联通分量
{
    DFN[u] = Low[u] = ++ Time ;
    Stack[top ++ ] = u ;
    Instack[u] = 1 ;

    for (int i = head[u] ; ~i ; i = edge[i].nex )
	{
        int v = edge[i].to ;
        if(DFN[v] == -1)
        {
            tarjan(v , u) ;
            Low[u] = min(Low[u] ,Low[v]) ;
            if(DFN[u] < Low[v])
            {
                edge[i].sign = 1;//为割桥
            }
        }
        else if(Instack[v]) Low[u] = min(Low[u] ,DFN[v]) ;
    }
    if(Low[u] == DFN[u])
	{
        int now;
        taj ++ ;
		bcc[taj].clear();
        do
		{
            now = Stack[-- top] ;
            Instack[now] = 0 ;
            Belong [now] = taj ;
            bcc[taj].push_back(now);
        }while(now != u) ;
    }
}


void tarjan_init(int all)
{
    memset(DFN, -1, sizeof(DFN));
    memset(Instack, 0, sizeof(Instack));
    top = Time = taj = 0;
    for(int i=1;i<=all;i++)if(DFN[i]==-1 )tarjan(i, i); //注意开始点标！！！
}


vector<int>G[maxn];
int du[maxn];

void  dfslian(int x)
{
   // cout<<x<<" mi "<<G[x].size()<<" "<<num<<" "<<taj<<endl;
    visit[x]=true;
    if(num==taj){ judge=true;return ;}
    if(G[x].size()>1)
        return ;
    for(int i=0;i<G[x].size();i++)
    {
        int temp=G[x][i];
        if(!visit[temp]){visit[temp]=true;num++;}
        dfslian(temp);
    }
}

void suodian()
{
    memset(du, 0, sizeof(du));
    for(int i = 1; i <= taj; i++)G[i].clear();
    for(int i = 0; i < edgenum; i++)
	{
        int u = Belong[edge[i].from], v = Belong[edge[i].to];
        if(u!=v)G[u].push_back(v), du[v]++;///g[u]存的是出度，du[v]存的是入度
    }
}

void init(){memset(head, -1, sizeof(head)); edgenum=0;}


int main()
{
    int cas;
    scanf("%d",&cas);
    for(int ca=1;ca<=cas;ca++)
    {
        bool flag=false;
        scanf("%d%d",&n,&m);
        int tempa,tempb;
        init();
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&tempa,&tempb);
            addedge(tempa,tempb);
        //    cout<<tempa<<" "<<tempb<<endl;
        }
        tarjan_init(n);
       // cout<<"taj  "<<taj<<endl;
        suodian();
        int in=0;
        for(int i=1;i<=taj;i++)
            if(du[i]==0)
            in++;
      //  cout<<in<<" "<<in<<endl;
        num=1;int index=0;
     //   cout<<in<<endl;
      //  for(int i=1;i<=taj;i++)
        //    cout<<du[i]<<" "<<i<<" "<<Belong[i]<<endl;
        memset(visit,false,sizeof(visit));
        if(in>1) flag=false;
        else if(in==0) flag=true;
        else
        {
          //  cout<<"miao"<<endl;
            for(int i=1;i<=taj;i++)
                if(du[i]==0)
                index=i;
    //        cout<<"index"<<index<<endl;
            judge=false;
            dfslian(index);
            if(judge)
                flag=true;
            else
                flag=false;
        }
         if(flag )
            cout<<"Yes"<<endl;
         else
            cout<<"No"<<endl;
    }

}