1. 首页
2. 模板
3. 本页

【模板】有向图的强连通分量

const int MAXN = 10010;  //点数
const int MAXM = 100100;  //边数

struct edge
{
    int to,next;

} edge[MAXM];
int head[MAXN],tot;
int Low[MAXN],DFN[MAXN],Stack[MAXN];
int Belong[MAXN]; //标记当前点属于第几个强连通分量

int index,top;
int scc;         //强连通分量的个数
bool instack[MAXN];
int num[MAXN];   //各个强连通分量包含点的个数

void addedge(int u,int v)
{
    edge[tot].to=v;
    edge[tot].next=head[u];
    head[u]=tot++;
}

void Tarjan(int u)
{
    int v;
    Low[u]=DFN[u]=++index; //访问时间
    Stack[top++]=u;
    instack[u]=true;
    for(int i=head[u]; ~i; i=edge[i].next)
    {
        v=edge[i].to; 
        if(!DFN[v])   //未访问过
        {
            Tarjan(v);
            if(Low[u]>Low[v])
                Low[u]=Low[v];
        }
        else if(instack[v]&&Low[u]>DFN[v]) //已经在栈  判断是否要更新low
            Low[u]=DFN[v];                  

    }
    if(Low[u]==DFN[u])
    {
        scc++;
        do
        {
            v=Stack[--top];
            instack[v]=false; //出栈
            Belong[v]=scc;
            num[scc]++;

        }
        while(v!=u);
    }

}

void solve(int n)
{
    memset(DFN,0,sizeof(DFN));
    memset(instack,0,sizeof(instack));
    memset(num,0,sizeof(num));
    index=scc=top=0;
    for(int i=1; i<=n; ++i)
    {
        if(!DFN[i])
            Tarjan(i);
    }

}

void init()
{
    tot=0;
    memset(head,-1,sizeof(head));

}