强联通分量算法入门

强联通分量，听着名字挺高大上的，一直听别人说但是自己也不知道是什么东西，昨天晚上，看见别人在学习这个东西，闲来无事的无就也学了学。

强连通分量在图论问题中得到广泛的应用，往往可以将有向图缩点，得到一个 DAG，于是避免了原图中可能有环造成后效性，可以在上面进行动态规划求解。 
 

强联通分量究竟是干什么用的呢？我感觉吧他就是为了节约时间，因为有的题图给的会比较大，如果我们按正常的方法去写，一般都会超时的，既然图这么大，我们能不能把图给简化一下呢，简单的说就是缩点，把一些可以放到因为一起的点给他们变成一个点，这样这个图就会小一点了，那么什么样的点才算是相同的点呢，假如一个子图的里面的点的集合，任意两个点都是可达的，那么我们就称这为强联通分量，就相当于是一个环。

强连通分量分解可以通过两次简单的DFS3实现。第一次DFS时，选取任意顶项点作为起点，遍历所看清来访过的顶点，井在回谢前给顶点标号 (pot oder,后序遍历)o对剩余的未访问过的顶
点，不断重复上述过程。
完成标号后，越接近图的尾部(搜索树的叶子)，顶点的标号越小。第二次DFS时，先将所有边反向，然后以标号最大的顶点为起点进行DFS。这样DFS所遍历的顶点集合就构成了一个强连通分量。之后，只要还有尚未访问的顶点，就从中选取标号最大的顶点不断重复上述过程。

下面给出模板，

#include<bits/stdc++.h>
using namespace std;
#define met(Q,QQ) memset(Q,QQ,sizeof(Q))
const int maxn=1e4+7;
vector<int> G[maxn];//存图
vector<int> rG[maxn];//把边反向后的图
vector<int> vec;//后序遍历顺序的图
bool vis[maxn];//访问标记
int cmp[maxn];//所属强联通分量的拓扑序列
int n,m;

void add_edge(int from,int to)//存图
{
    G[from].push_back(to);
    rG[to].push_back(from);
}

void dfs(int k)     //对图进行dfs，并访问到的点依次存入vec中,
{
    vis[k]=true;
    int kk=G[k].size();
    for(int i=0;i<kk;i++)
    {
        int kkk=G[k][i];
        if(!vis[kkk]) dfs(kkk);
    }
    vec.push_back(k);
}

void rdfs(int k,int cnt)//构建拓扑序
{
    vis[k]=true;
    cmp[k]=cnt;
    int kk=rG[k].size();
    for(int i=0;i<kk;i++)
    {
        int kkk=rG[k][i];
        if(!vis[kkk]) rdfs(kkk,cnt);
    }
}
int main()
{
    scanf("%d %d",&n,&m);
    int from,to;
    for(int i=1;i<=m;i++)
    {
        scanf("%d %d",&from,&to);
        add_edge(from,to);
    }
    met(vis,0);
    for(int i=1;i<=n;i++)
    {
        if(!vis[i]) dfs(i);
    }

    met(vis,0);
    int kk=vec.size();
    int cnt=0;
    for(int i=kk-1;i>=0;i--)
    {
        if(!vis[vec[i]])
        {
            cnt++;
            int qq=vec[i];
            rdfs(qq,cnt);
        }
    }


    return 0;
}

 

感觉这个模板的最终的结果都存到了cmp数组中了，也就是图的拓扑序，遇到题的时候只要灵活的运用cmp就行了