Tarjan的强联通分量

　　求强联通分量有很多种。 《C++信息学奥赛一本通》  中讲过一个dfs求强联通分量的算法Kosdaraju，为了骗字数我就待会简单的说说。然而我们这篇文章的主体是Tarjan，所以我肯定说完之后再赞扬一下Tarjan大法好是不是

　　首先我们讲一下强联通分量

　　强联通分量指的是图的一个子图。在这个子图中，任意两个节点都可以互相到达。从定义上我们就可以看出是一个有向图来，因为任意一个无向图都符合该定义。

　　而它的标准定义是：有向图中任意两点都联通的最大子图。

 　　　　　　　　　　　　　　　　　　　　

　　咳咳，首先庆祝一下哈——本人博客的第一张图。绘图历时3分钟。

　　在咱们举的例子中，可以看出1 、2 、3 、5 通过边可以相互到达，它们算一个强联通分量，但4却被它们隔绝在外。从图中可以看出，从4点出发不能到达任意一个点。所以它单个节点也算一个强联通分量。所以图中的强联通分量有两个：一个是1-2-3-5，一个是4。

　　ok看完了强联通分量是什么我们就讲一下Kosaraju。

　　这个算法的思路是，对图进行DFS并记录每个点的退出顺序。再构造反图(就是有向边的方向全都反过来),按照退出顺序的逆序DFS反图，对得到的点进行染色即为强联通分量。

　　讲完思路开始模拟。以起点1为起点遍历顺序如下：

　　[ 1 2 3 5 4  5 3 2 4 4 1 ]

　　加粗斜体外带下划线的部分就是本图的退出顺序。

　　于是我们得到这样一个数组：[ 5 3 2 4 1 ] 。按照这个数组的逆序对反图遍历得到：

　　[ 5 3 2 1 退出 4 退出 ]

　　即得到要求的两个强联通分量。

　　还要两遍DFS，麻烦的一比。看我大Tarjan一遍DFS就能求出强联通分量

　　首先我们要明确Tarjan要用到的两个数组：dfn[] 和 low[] 

　　dfn指的是在DFS过程中访问到该点的顺序。从1开始DFS全图，那么1的dfn值就是1，2的dfn值是2,5的dfn值是4，4的dfn值是5。剩下的一个类推

　　那么low呢？low指的是如果逆着DFS序往前回溯，该节点最早是由哪个节点走过来的。

　　比如在上图中2 、3 、5 、4 最早都是由1走过来的，所以它们的low值都是1

　　下面贴出dfn和low的算法

每次dfs(点u){

　　dfn[u] = 进入 dfs() 函数的次数  （自己定义一个时间戳记录 如 time）

       枚举与其相邻的点v{

   　　     如果 没有 访问过点v {   ( 就是dfs树上的树边 )

　　　　　　　　dfs(v);

　　　　　　　　如果 v 能追溯 到 比“u 追溯到的最早的点” 更早的点；

　　　　　　　　那么 u 就能 通过 v 来追溯到 那个点；

　　　　　　　　low[u]=min(low[u],low[v]);

　　　　  }

 　　　　  如果 访问过点v && v在栈中

　　　　　　　low[u]=min(low[u],dfn[v]); 

     }

　　缩点

}

　　上面那些伪代码是从伟大的GeneralLiu那里带过来的，在此先%%%

　　然后  假设我们走到一个节点i，发现这个i不能继续扩展了，也就是dfn[i]==low[i]

　　于是我们开始往回走。往回走的过程中，我们就把和它一个分量的节点进行染色，给它们统一的标记。　　最后统计有多少种不同的标记即是强联通分量个数

　　luogu的一道题刻录光盘非常好，可以用于练手。

　　放代码

#include<iostream>
#include<cstring>
using namespace std;
int head[10000],num;
struct Edge{
    int next,to;
}edge[100000];
int stack[10000],top;
int color[10000],cnt;
int dfn[10000],low[10000];
int ID;
bool jd[10000];
int vis[10000];
inline void add(int from,int to){
    edge[++num]={head[from],to};
    head[from]=num;
}

void tarjan(int x){
    dfn[x]=++ID;
    low[x]=ID;
    jd[x]=1;
    stack[++top]=x;
    for(int i=head[x];i;i=edge[i].next){
        int to=edge[i].to;
        if(!dfn[to]){
            tarjan(to);
            low[x]=min(low[x],low[to]);
        }
        else if(jd[to])    low[x]=min(low[x],dfn[to]);
    }
    if(dfn[x]==low[x]){
        jd[x]=0;
        color[x]=++cnt;
        while(stack[top]!=x){
            color[stack[top--]]=cnt;
            jd[stack[top+1]]=0;
            color[stack[top+1]]=cnt;
        }
        top--;
    }
    
}

int main(){
    int n;
    cin>>n;
    int x;
    for(int i=1;i<=n;++i){
        while(cin>>x&&x!=0){
            add(i,x);
        }
    }
    for(int i=1;i<=n;++i){
        if(!dfn[i])    tarjan(i);
    }
    memset(jd,0,sizeof(jd));
    for(int i=1;i<=n;++i){
        for(int j=head[i];j;j=edge[j].next){
            if(color[i]!=color[edge[j].to]){
                jd[color[edge[j].to]]=1;
            }
        }
    }
    int ans=0;
    for(int i=1;i<=cnt;++i)    if(!jd[i])    ans++;
    cout<<ans<<endl;
    return 0;
}