【模板】 强连通分量分解

#include<vector>
#include<cstring>
#define MAX_V 10050

int V,E;                  //V表示顶点数，E表示边数（顶点标号为0~V-1)
vector<int> G[MAX_V];     //图的邻接表表示
vector<int> rG[MAX_V];    //把边反向后的图
vector<int> vs;           //后序遍历顺序的顶点列表（即用一个栈模拟了回溯时的标号）
bool used[MAX_V];         //DFS中用到的访问标记
int cmp[MAX_V];           //顶点所属强连通分量的拓扑序（0~k-1）

//向图中添加一条从_from到_to的边
void add_edge(int _from,int _to)
{
    G[_from].push_back(_to);
    rG[_to].push_back(_from);
}

//第一次在原图上的DFS
void dfs(int v)
{
    used[v]=true;
    for(int i=0;i<G[v].size();i++)
        if(!used[G[v][i]])
            dfs(G[v][i]);
    vs.push_back(v);
}

//第二次在反向图上的DFS
void rdfs(int v,int k)
{
    used[v]=true;
    cmp[v]=k;
    for(int i=0;i<rG[v].size();i++)
        if(!used[rG[v][i]])
            rdfs(rG[v][i],k);
}

//求解强连通分量（SCC:Strongly Connected Component），返回值为图中强连通分量的个数
int scc()
{
    memset(used,false,sizeof(used));
    vs.clear();
    for(int v=0;v<V;v++)
        if(!used[v])
            dfs(v);

    memset(used,false,sizeof(used));
    int k=0;
    for(int i=vs.size()-1;i>=0;i--)
        if(!used[vs[i]])
        {
            rdfs(vs[i],k);
            k++;
        }

    return k;
}