POJ2181强连通分量

POJ 2181 强连通分量

  (2009-07-11 00:00:00)

转载

标签：  强连通分量   杂谈

分类： POJ解题报告

题目大意 ：

有一些牛，然后他们互相仰慕，为有向图，仰慕可传递， a仰慕b， b仰慕c，a就仰慕c，让你求
解被所有牛仰慕的牛的个数。

思路：

先做出强连通分量，然后将一个分量作为一个节点（分量中的每一个节点都具有相同的性质），重新构图。找出是否有出度为0的点。如果出度为0的点不是一，则不存在这样的cow。

如果存在设其id为k

    则 剩下的强连通中 至少存在一个强连通分量i指向k   剩下的num-2个强连通分量中有指向k或者i的。。。（这里即可转化为 判断剩下的图是否是连通

 

代码

#include <iostream>
#include<vector>
#include<map>
using namespace std;
const int MAXN = 10001;//表示节点的最大值
//typedef int AdjTable[MAXN]; //邻接表类型
int      n;//n表示有多少个点
bool     flag[MAXN]; //访问标志数组
int      belg[MAXN]; //存储强连通分量,其中belg[i]表示顶点i属于第belg[i]个强连通分量
int      numb[MAXN]; //结束时间标记,其中numb[i]表示离开时间为i的顶点
vector<int> adj[MAXN],radj[MAXN];//邻接表,逆邻接表  
//用于第一次深搜,求得numb[1..n]的值  （注意顶点id从1开始）
int m;
int out[MAXN];
int f[MAXN];
int num;
void VisitOne(int cur, int &sig)
{
flag[cur] = true;
for ( int i=0; i<adj[cur].size(); ++i )
{
     if ( false==flag[adj[cur][i]] )
     {
         VisitOne(adj[cur][i],sig);
     }
}
numb[++sig] = cur;
}
//用于第二次深搜,求得belg[1..n]的值
void VisitTwo(int cur, int sig)
{
flag[cur] = true;
belg[cur] = sig;
for ( int i=0; i<radj[cur].size(); ++i )
{
     if ( false==flag[radj[cur][i]] )
     {
         VisitTwo(radj[cur][i],sig);
     }
}
}
//Kosaraju算法,返回为强连通分量个数
int Kosaraju_StronglyConnectedComponent()
{
int i, sig;
//第一次深搜
memset(flag,0,sizeof(flag));
for ( sig=0,i=1; i<=n; ++i )
{
     if ( false==flag[i] )
     {
         VisitOne(i,sig);
     }
}
//第二次深搜
memset(flag,0,sizeof(flag));
for ( sig=0,i=n; i>=1; --i )
{
     if ( false==flag[numb[i]] )
     {
         VisitTwo(numb[i],++sig);
     }
}
return sig;  
}
int getf(int k)
{
    int root=k;
    while(root!=f[root])root=f[root];
    return root;
}
void re_graph(void)
{
     memset(out,0,sizeof(out));
     int i,j,b;
     for(i=1;i<=n;i++)
     {
         out[i]=0;
         f[i]=i;
     }
     for(i=1;i<=n;i++)
     {
         int t=adj[i].size();
         for(j=0;j<t;j++)
         {
             b=adj[i][j];
             if(belg[i]!=belg[b])
             {  
                 out[belg[i]]++;
                 f[belg[i]]=getf(belg[b]);
             }
         }
     }
}
int pd()
{
    int i;
    int kk=0;
    int pp=0;
    int qq=getf(1);
    for(i=1;i<=num;i++)
    {
        if(out[i]==0)
        {  
            kk++;
            pp=i;
        }
        if(getf(i)!=qq)
           return 0;
    }
    if(kk!=1)
       return 0;
    int result=0;
    for(i=1;i<=n;i++)
        if(belg[i]==pp)
            result++;
    return result;
}
int main()
{
    scanf("%d%d",&n,&m);
    int i,j,a,b;
    for(i=0;i<=n;i++)
    {
        adj[i].clear();
        radj[i].clear();            
    }   
    for(i=0;i<m;i++)
    {  
        scanf("%d%d",&a,&b);
        adj[a].push_back(b);
        radj[b].push_back(a);
    }
    num=Kosaraju_StronglyConnectedComponent();  //用此表示出有几个强连通分量  强连通分量的id从1开始
    re_graph();
    printf("%d/n",pd());
    return 1;
}