求有向图的强连通分量个数（kosaraju算法）

有向图的连通分量的求解思路

 

kosaraju算法

逛了很多博客，感觉都很难懂，终于找到一篇能看懂的，摘要记录一下

原博客https://www.cnblogs.com/nullzx/p/6437926.html

关于连通分量是什么自行百度，这里主要说明连通分量的求解方法

基本思路：第一次DFS得出顶点的顺序，根据顶点顺序进行第二次DFS，也就是逆后序遍历（手动模拟一下堆栈就知道第二次DFS的过程就能得出答案）。

为什么要两次DFS?

如果从连通分量A中任意一个定点DFS，得不到正确结果。应该按照被指向的强连通分量的定点排在前面的顺序进行DFS。上图按照B3，B4，B5，A0，A1，A2的顺序DFS。实际中我们只要保证被指向的强连通分量的至少一个顶点排在指向这个连通分量的所有顶点前面即可，比如B3，A0，A1，A2，B4，B5；B3排在强连通分量A所有定点的前面。

如何得到满足要求的顶点顺序：对原图取反，从反向图的任意节点开始进行DFS的逆后序遍历

DFS的逆后序遍历指：如果当前顶点没被访问，先遍历完与当前顶点相连的且未被访问的所有其他顶点，然后将当前顶点加入栈，最后从栈顶到栈底的顺序是我们需要的顶点顺序。

其实它是利用了有向图的方向性：例如在上图中，强连通分量A和B不管正反图都能自己跑一圈，但是从A到B就只能从A2跑到B3，不可能从B3跑到A2，所以将图取反（注意图取反，不能从A2跑到B3，顶点顺序被记录，按顶点顺序一个个弹出遍历，之前遍历过的点就不用再遍历），做成反向图，再逆后序遍历，A0在栈顶，遍历一圈，不能从A2跑到B3，就得到一个连通分量（如下图）

接下来原博主的代码看不懂，还是百度百科kosaraju算法里面的代码好

#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <stack>
#include <algorithm>
using namespace std;
#define INF 999999999
#define MAXN 551
#define MOD 1000000009
int n;
int mp[MAXN][MAXN], nmp[MAXN][MAXN], vis[MAXN];
stack<int> s;
void dfs_1(int v)
{
//    cout << v <<endl;
    vis[v] = 1;
    for(int i = 1; i <= n; ++i)
        if(!vis[i] && mp[v][i])
        dfs_1(i);
    s.push(v);
}
void dfs_2(int v)
{
    vis[v] = 1;
    for(int i = 1; i <= n; ++i)
        if(!vis[i] && nmp[v][i])
        dfs_2(i);
}
int kosaraju()
{
    while(!s.empty())
        s.pop();
    memset(vis, 0, sizeof(vis));
    
    for(int i = 1; i <= n; ++i)
        if(!vis[i])
        dfs_1(i);
        
    int ans=0;
    memset(vis, 0, sizeof(vis));
    
    while(!s.empty())
    {
        int v = s.top();
        s.pop();
        if(!vis[v])
        {
//            cout << v <<endl;
            ans++;
            dfs_2(v);
        }
    }
    return ans;
}
int main()
{
    int m, a, b;
    //n个点，m条边
    cin >> n >>m;
    memset(mp, 0, sizeof(mp));
    memset(nmp, 0, sizeof(nmp));
    for(int i = 0; i < m; ++i)
    {
        cin >> a >>b;
        mp[a][b] = 1;
        nmp[b][a] = 1;
    }
    cout << kosaraju() <<endl;
    return 0;
}
/*
5 5
1 2
2 1
2 3
3 4
4 1
样例输出：
2
*/