Tarjan——强连通分量

是的你没有看错,又是Tarjan(说过他发明了很多算法),那么什么是Tarjan

首先看看度娘的解释:

有向图强连通分量:在有向图G中,如果两个顶点vi,vj间(vi>vj)有一条从vi到vj的有向路径,同时还有一条从vj到vi的有向路径,
则称两个顶点强连通(strongly connected)。如果有向图G的每两个顶点都强连通,称G是一个强连通图。有向图的极大强连通子图,
称为强连通分量(strongly connected components)。

是不是没有看懂,我也是

首先了解几个概念:强连通,强连通图,强连通分量

强连通:在一个有向图G中,两个点a,b,a可以走到b,b可以走到a,我们就说(a,b)强连通

强连通图:在一个有向图G中,任意两个点都是强连通

强连通分量:在一个有向图G中,有一个子图,它任意两个点都是强连通,我们就说这个子图为强连通分量,特别的,一个点也是一个强连通分量

如图所示:

显然可得:1,2,3,5 构成了一个强连通分量(一个点也是)

下面进入正题,厉害的Tarjan发明的厉害的Tarjan

首先引进一个概念(以前应该都学过),时间戳,用数组dfn表示(应该不用讲吧),也就是搜索这个图的顺序(最后还是讲了)

每个节点的时间戳不同

再是一个low数组,low[x]表示以x为根的子树中,每个节点中连接的点的时间戳的最小值(可能有点难懂,但这个非常重要,是核心思想,等一下的模拟过程会详细讲述)

low的初值:low[x]=dfn[x]

那么如何储存强连通分量呢,可以用栈(学c++的有福利啦,我们有STL)

可惜我不会用,(手写栈挺好的),每次遍历到一个新节点,就把它放进栈,如果这个点有出度,就继续往下找,直到不能再找

每一次回来都要更新low值,当然是取小的那个,如果发现low[x]=dfn[x]那么它的子节点中肯定有一个连上来,既然可以过去又可以回来,很明显是一个强连通分量

那么这个x就是这个强连通分量的根节点,那么栈中间,比这个x晚进来的点就是x的子节点,那么这些点全部出栈,就组成了一个强连通分量

到这来就完了,但是好像还是没有理解透彻(反正我是这样)

那么就模拟一下

还是这张图5——>4

low[1]=dfn[1]=1,1入栈

low[2]=dfn[2]=2,2入栈

low[3]=dfn[3]=3,3入栈

low[5]=dfn[5]=4,4入栈

然后发现5连接着1,已经寻找过得了,那么就看看,谁才是真正的祖先

low[1]=1,low[5]=4,好吧,5输了,所以1是5的根节点,low[5]=min(low[5],low[1])=1

继续发现还有4,low[4]=dfn[4]=5,4入栈

但是4已经没有了出度,往回退

发现low[4]=dfn[4]那么4就是一个强连通分量的根节点(其实也就它一个),4退栈

继续往回退:low[5]=min(low[5],low[4])=1;

继续:一直到1,low[3]=min(low[3],low[5])=1;low[2]=min(low[2],low[3])=1;

low[1]=min(low[1],low[2])=1;

发现此时low[1]=dfn[1],所以1也是一个强连通分量的根,此时发现栈里还有1,2,3,5

所以这个强连通分量为1,2,3,5

1还有一个出度:4

寻找4,low[4]=dfn[4]=6,发现没有出度

low[4]=dfn[4],所以4是一个强连通分量的根节点(还是只有他一个),退栈

往回退,low[1]=min(low[1],dfn[4])=1;

这样就完了吗?

万一还有图没有遍历到呢

所以要加一个语句:

for(int i=1;i<=number;i++)
       if(!dfn[i])Tarjan(i);

以上的模拟过程用程序表示就是:

void Tarjan(int x){
    dfn[x]=low[x]=++deep;
    q[++hd]=x;res[x]=1;
    for(int i=head[x];i;i=next[i]){
        int y=ver[i];
        if(!dfn[y]){
            Tarjan(y);
            low[x]=min(low[x],low[y]);
        }
        else if(res[y])low[x]=min(low[x],dfn[y]);
    }
    if(low[x]==dfn[x]){
        int sum=0;
        do{
            res[q[hd]]=0;hd--;sum++;
        }while(x!=q[hd+1]);
        if(sum>1)ans++;
    }
    
}

来一个模板题:洛谷

P2863 [USACO06JAN]牛的舞会The Cow Prom

大意是求强连通分量数量,不过注意的是,每一个点不能算,所以要加一个特判

以下是AC代码:

#include<bits/stdc++.h>
using namespace std;

const int N=50002;
int number,m,ans=0;
int tot,next[N],ver[N],head[N];
int dfn[N],low[N],deep,res[N],q[N],hd;

int read(){
    int s=0,w=1;char ch=getchar();
    while(ch<'0'||ch>'9')w=(ch=='-')?-1:1,ch=getchar();
    while(ch>='0'&&ch<='9')s=s*10+ch-'0',ch=getchar();
    return s*w;
}

void add(int x,int y){
    ver[++tot]=y;next[tot]=head[x];head[x]=tot;
}

void Tarjan(int x){
    dfn[x]=low[x]=++deep;
    q[++hd]=x;res[x]=1;
    for(int i=head[x];i;i=next[i]){
        int y=ver[i];
        if(!dfn[y]){
            Tarjan(y);
            low[x]=min(low[x],low[y]);
        }
        else if(res[y])low[x]=min(low[x],dfn[y]);
    }
    if(low[x]==dfn[x]){
        int sum=0;
        do{
            res[q[hd]]=0;hd--;sum++;
        }while(x!=q[hd+1]);
        if(sum>1)ans++;
    }
    
}

int main(){
    number=read();m=read();
    for(int i=1;i<=m;i++){
        int x=read(),y=read();
        add(x,y);
    }
    for(int i=1;i<=number;i++)
       if(!dfn[i])Tarjan(i);
    cout<<ans;
    return 0;
}
原文地址:https://www.cnblogs.com/GMSD/p/11320154.html