学步园

tarjan这个人很神奇呀，发表的算法非常的多，这个求有向图的强连通量，也是人家发明的。这个算法可以求出图中的强连通分量数目，和分别每个顶点所对应的强连通分量的编号，十分的好用呀，结合上其他的算法，就可以将图论的连通性展现的淋漓尽致呀。

直接贴出代码吧：

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <vector>
#include <stack>
#include <string>

/*
Õâ¾ÍÊǷdz£¾­µäµÄtarjanËã·¨£¬Õâ¸öËã·¨µÄÖ÷Òª×÷ÓþÍÊÇÇó³öÓÐÏòͼµÄÇ¿Á¬Í¨ÊýÁ¿£¬
²¢ÇÒ½«Ã¿¸öÇ¿Á¬Í¨·ÖÁ¿·Ö±ð±êºÅ£¬ÄÇôֻҪ±êºÃºÅÁËÖ®ºó£¬ÄǾͿÉÒ԰ɱêºÅÊÇÒ»ÑùµÄµãÈ«²¿
Ëù³ÆÒ»¸öµãÁË£¬ÒòΪһ¸öÇ¿Á¬Í¨Í¼ÉϵÄËùÓе㶼ÊÇ¿ÉÒÔ»¥Ïൽ´ïµÄ£¬¿ÉÒÔËõ³ÉÒ»¸öµã£¬
Èç¹ûÒªÇóËùÓеÄÇ¿Á¬Í¨·ÖÁ¿£¬Á¬³ÉÒ»¸ö£¬ÕûÌåµÄÇ¿Á¬Í¨·ÖÁ¿£¬¾ÍҪͳ¼Æÿ¸öËõµãºóµÄµãµÄÈë¶ÈΪ0µÄµã£¬
³ö¶ÈΪ0µÄµã¡£È»ºó±È½ÏÁ½ÕߵĴóС£¬´óµÄ¾ÍÊÇÒªÌí¼ÓµÄ±ß¡£ÕâÊÇÒòΪÿ¸öµãÈôÏëÁªÍ¨ËùÓеĵ㣬
±ØÐëÊÇÒªÓиö³ö¶ÈºÍÈë¶ÈµÄ¡£ 
*/

using namespace std;

const int MAXN = 100 + 11;

int N, M;

vector <int> G[MAXN];

int pre[MAXN], lowlink[MAXN], sccno[MAXN], dfs_clock, scc_cnt;

stack <int> S;

void init()
{
	for (int i = 0; i <= N; i++)
	{
		G[i].clear();
	}
}

void DFS(int u)
{
	pre[u] = lowlink[u] = ++dfs_clock;
	S.push(u);
	for (int i = 0; i < G[u].size(); i++)
	{
		int v = G[u][i];
		if (!pre[v])
		{
			DFS(v);
			lowlink[u] = min(lowlink[u], lowlink[v]); //Çó³ö´Ë½ÚµãÒÔϵÄ×Ó½ÚµãµÄ×îСʱ¼ä´Á¡£ 
		}
		else if (!sccno[v])//Èç¹ûÕâ¸öµã²»ÊôÓÚÈκÎÁ¬Í¨·ÖÁ¿¡£ 
		{
			lowlink[u] = min(lowlink[u], pre[v]);
		}
	}
	if (lowlink[u] == pre[u])
	{
		scc_cnt++;
		for (;  
		{
			int x = S.top();
			S.pop();
			sccno[x] = scc_cnt;
			if (x == u)
			{
				break;
			}
		}
	}
}

void find_scc()
{
	dfs_clock = scc_cnt = 0;
	memset(sccno, 0, sizeof(sccno)); //sccnoÊÇÓÃÀ´´æÿ¸öµã¶ÔÓ¦µÄÁ¬Í¨·ÖÁ¿µÄ±êºÅµÄ¡£ 
	memset(pre, 0, sizeof(pre)); //pre»¹ÊÇÀϹæ¾Ø£¬ÓÃÀ´±ê¼ÇÊÇ·ñ·ÃÎÊ£¬ºÍµ±×öʱ¼ä´Á¡£ 
	for (int i = 1; i <= N; i++)
	{
		if (!pre[i])
		{
			DFS(i);
		}
	}
}

int main()
{
	while (scanf("%d%d", &N, &M) != EOF)
	{
		init();
		int a, b;
		for (int i = 0; i < M; i++)
		{
			scanf("%d%d", &a, &b);
			G[a].push_back(b);
		}
		find_scc();
		printf("the scc_cnt = %d\n", scc_cnt);
	} 
	return 0;
} 

/*
6 8
1 2
2 3
2 4
2 5
5 2
5 6
3 6
6 3

*/