【模板】强连通分量分解

const int MAX_V=50005;  
vector<int> G[MAX_V];  
vector<int> rG[MAX_V];  
int V;  
vector<int> vs;  
int comp[MAX_V];  
bool used[MAX_V];  
void add_edge(int from,int to)  
{  
    G[from].push_back(to);  
    rG[to].push_back(from);  
}  
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);  
}  
void rdfs(int v,int k)  
{  
    comp[v]=k;  
    used[v]=true;  
    for(int i=0;i<rG[v].size();i++)  
    {  
        if(!used[rG[v][i]]) rdfs(rG[v][i],k);  
    }  
}  
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++);  
    }  
    return k;  
}