迷宫城堡（强连通分量）

为了训练小希的方向感，Gardon建立了一座大城堡，里面有N个房间(N<=10000)和M条通道(M<=100000)，每个通道都是单向的，就是说若称某通道连通了A房间和B房间，只说明可以通过这个通道由A房间到达B房间，但并不说明通过它可以由B房间到达A房间。Gardon需要请你写个程序确认一下是否任意两个房间都是相互连通的，即：对于任意的i和j，至少存在一条路径可以从房间i到房间j，也存在一条路径可以从房间j到房间i。
Input
输入包含多组数据，输入的第一行有两个数：N和M，接下来的M行每行有两个数a和b，表示了一条通道可以从A房间来到B房间。文件最后以两个0结束。
Output
对于输入的每组数据，如果任意两个房间都是相互连接的，输出"Yes"，否则输出"No"。
Sample Input

3 3
1 2
2 3
3 1
3 3
1 2
2 3
3 2
0 0

Sample Output

Yes
No

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<vector>
#include<stack> 
using namespace std;
const int maxn=1e+5;
int cnt,num,n,m;
int low[maxn],dfn[maxn];
int book[maxn];
vector<int>a[maxn];
stack<int>s; 
void tarjan(int u)
{
	dfn[u]=low[u]=++num;
	s.push(u);
	int i;
	book[u]=1;
	for(i=0;i<a[u].size();i++)
	{
		int v=a[u][i];
		if(dfn[v]==0)
		{
			tarjan(v);
			low[u]=min(low[u],low[v]);
		}
		else if(book[v]==1)
		low[u]=min(low[u],low[v]);
	}
		if(dfn[u]==low[u])
		{
			cnt++;
			int vt;
			do
			{
				vt=s.top();
				s.pop();
				book[vt]=0;
			} while(u!=vt);
	 	} 
}
int main()
{
	int i,c,b;
	while(~scanf("%d%d",&n,&m))
	{
	if(n==0&&m==0)
	break;
	//memset(book,0,sizeof(book));
	memset(dfn,0,sizeof(dfn));
	memset(low,0,sizeof(low));
	for(i=1;i<=n;i++)
	a[i].clear();
	for(i=1;i<=m;i++)
	{
		scanf("%d%d",&c,&b);
		a[c].push_back(b); 
	}
	cnt=0,num=0;
	for(i=1;i<=n;i++)
	if(dfn[i]==0)
	tarjan(i);
	if(cnt==1)
	printf("Yes\n");
	else
	printf("No\n",cnt);
	}
}

kosaraju

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<vector>
using namespace std;
int N,V,M;
//vector<int> v[n]:表示n个vector v ，二维数组
vector<int>g[90000]; 
vector<int>rg[90000];
vector<int>vs;//后序遍历顺序的顶点列表 
int book[90000];
int a[90000];
int x[900000],y[900000];
void add(int u,int v)
{
	g[u].push_back(v);//在数组的最后添加一个数据
	rg[v].push_back(u);
}
void dfs(int v)
{
	book[v]=1;
	int i;
	for(i=0;i<g[v].size();i++)
	if(!book[g[v][i]])
	dfs(g[v][i]);
	vs.push_back(v); 
}
void rdfs(int v,int k)
{
	book[v]=1;
	a[v]=k;
	int i;
	for(i=0;i<rg[v].size();i++)
	{
		if(!book[rg[v][i]])
		rdfs(rg[v][i],k);
	}
}
int scc()
{
	memset(book,0,sizeof(book));
	vs.clear();
	int i;
	for(i=0;i<V;i++)
	if(!book[i])
	dfs(i);
	memset(book,0,sizeof(book));
	int k=0;
	for(i=vs.size()-1;i>=0;i--)
	if(!book[vs[i]])
	rdfs(vs[i],k++);
	return k;
}
int main()
{
	while(~scanf("%d%d",&N,&M))
	{
	//memset(a,-1,sizeof(a));
	//memset(y,-1,sizeof(y));
	//memset(x,-1,sizeof(x));
	if(N==0&&M==0)
	break;
	int i;
	for(i=0;i<N;i++)
	{
        g[i].clear();
        rg[i].clear();
    }
	for(i=0;i<M;i++)
	scanf("%d%d",&x[i],&y[i]);
	V=N;
	for( i=0;i<M;i++)
	add(x[i]-1,y[i]-1);
	int n=scc();
	if(n==1)
	printf("Yes\n");
	else
	printf("No\n");
	}
}

//