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
*/