poj 2186(强连通分量)
题意：
n头奶牛，给出若干个欢迎关系a b，表示a欢迎b，欢迎关系是单向的，但是是可以传递的,如：a欢迎b，b欢迎c，那么a欢迎c 。另外每个奶牛都是欢迎他自己的。求出被所有的奶牛欢迎的奶牛的数目.
方法：对有向图求强连通分量，然后找出所有独立的强连通分量（所谓独立，就是该连通分量里面的点到外面的点没有通路，当然，连通分量外的点是可以有路到强连通分量内的点的），如果独立的强连通分量的数目只有一个，那么，就输出这个强连通分量内解的个数，否则输出无解。
kosaraju:
http://www.nocow.cn/index.php/Kosaraju%E7%AE%97%E6%B3%95
View Code
1 // File Name: 2186.cpp 2 // Author: Missa 3 // Created Time: 2013/2/5 星期二 20:47:56 4 5 #include<iostream> 6 #include<cstdio> 7 #include<cstring> 8 #include<algorithm> 9 #include<cmath> 10 #include<queue> 11 #include<stack> 12 #include<string> 13 #include<vector> 14 #include<cstdlib> 15 #include<map> 16 using namespace std; 17 18 //***********************kosaraju************** 19 const int maxn = 1e4+5; 20 const int maxm = 5e4+5; 21 int n,m; 22 vector <int> adj[maxn];//正向邻接表 23 vector <int> radj[maxn];//反向邻接表 24 vector <int> ord;//后续访问顺序 25 bool vis[maxn]; 26 int ma[maxn];//表示点v属于哪个连通分量 27 int cnt;//有向图强连通分量的个数 28 void init() 29 { 30 for(int i=0;i<maxn;i++) 31 { 32 adj[i].clear(); 33 radj[i].clear(); 34 } 35 ord.clear(); 36 cnt=0; 37 } 38 void dfs1(int v) 39 { 40 vis[v]=1; 41 for(int i=0;i<adj[v].size();i++) 42 if(!vis[adj[v][i]]) 43 dfs1(adj[v][i]); 44 ord.push_back(v); 45 } 46 void dfs2(int v) 47 { 48 vis[v]=1; 49 ma[v]=cnt; 50 for(int i=0;i<radj[v].size();i++) 51 if(!vis[radj[v][i]]) 52 dfs2(radj[v][i]); 53 } 54 void kosaraju() 55 { 56 memset(vis,0,sizeof(vis)); 57 ord.clear(); 58 for(int i=1;i<=n;i++)//顶点从1--n 59 if(!vis[i]) 60 dfs1(i); 61 memset(vis,0,sizeof(vis)); 62 cnt=0; 63 for(int i=ord.size()-1;i>=0;i--) 64 if(!vis[ord[i]]) 65 { 66 cnt++; 67 dfs2(ord[i]); 68 } 69 } 70 //********************************************* 71 void solve() 72 { 73 kosaraju(); 74 int ans[maxn]; 75 memset(vis,0,sizeof(vis)); 76 memset(ans,0,sizeof(ans)); 77 for(int i=1;i<=n;i++) 78 { 79 ans[ma[i]]++;//统计每个连通分量的大小 80 for(int j=0;j<adj[i].size();j++) 81 { 82 int v=adj[i][j]; 83 if(ma[i] != ma[v])//有向边(i,v),若i与v不在一个连通分量 84 vis[ma[i]]=1;//i所在的强连通分量出度不为0. 85 } 86 } 87 int flag=0,cow; 88 for(int i=1;i<=cnt;i++) 89 { 90 if(!vis[i])//出度为0 91 { 92 cow=i; 93 flag++; 94 } 95 } 96 if(flag==1) 97 printf("%d\n",ans[cow]); 98 else 99 printf("0\n");//存在多个出度为0。表明至少两组是不互相羡慕的 100 } 101 int main() 102 { 103 while(~scanf("%d%d",&n,&m)) 104 { 105 init(); 106 for(int i=0;i<m;i++) 107 { 108 int x,y; 109 scanf("%d%d",&x,&y); 110 adj[x].push_back(y); 111 radj[y].push_back(x);//反向边 112 } 113 solve(); 114 } 115 return 0; 116 } 117 118 /* 119 3 3 120 1 2 121 2 1 122 2 3 123 * /
