Tarjan

一、Tarjan算法求LCA

二、Tarjan算法求强连通分量【栈里存点---------每个点都属于一个强连通分量】

(着重理解一下第11行的else if语句:dfn[ i ]!=0 且 instk[i]==0的点一定已经是另一个强连通分量里面的点了,所以就不用考虑了,所以用else if(instk[ i ])来判断是不是已经是别的强连通分量里的点了,如果不是,则是当前点的祖先,可以用来更新low[ u ],这也是强连通分量的更新low[ u ]的精髓。)

 1 void Tarjan(int u) {
 2     dfn[u] = low[u] = ++tim;
 3     in_stk[u] = 1;
 4     stk[tp++] = u;
 5     for (int i = head[u]; i != -1; i = e[i].next) {
 6         int v = e[i].to;
 7         if (!dfn[v]) {
 8             Tarjan(v);
 9             if (low[v] < low[u])low[u] = low[v];
10         }
11         else if (in_stk[v]) {
12             if (low[v] < low[u])low[u] = low[v];//这里使用dfn[v]和low[v]都可以
13         }
14     }
15     if (dfn[u] == low[u]) {
16         scc_cnt++;
17         do {
18             tp--;
19             int tt = stk[tp];
20             col[tt] = scc_cnt;
21             in_stk[tt] = 0;
22             vc[scc_cnt].push_back(tt);
23         } while (stk[tp] != u);
24     }
25 }

 

三、Tarjan算法求双连通分量【栈里存边------每条边都属于一个点双连通分量】

(由于是无向图,这里需要注意的是不能重复访问边,用else if(dfn[ u ]>dfn[ v ]&&v!=fa)来防止重复访问)

 1 void Tarjan(int u, int fa) {
 2     dfn[u] = low[u] = ++tim;
 3     int child = 0;
 4     for (int i = head[u]; i != -1; i = e[i].next) {
 5         int v = e[i].to;
 6         struct edge e1;
 7         e1.fr = u; e1.to = v;
 8         if (!dfn[v]) {
 9             stk.push(e1);
10             child++;
11             Tarjan(v);
12             if (low[v] >= dfn[u]) {
13                 iscut[u] = true;
14                 bcc_cnt++; bcc[bcc_cnt].clear();
15                 do {
16                     struct edge e2 = stk.top(); stk.pop();
17                     if (col[e2.fr]!=bcc_cnt) {
18                         bcc[bcc_cnt].push_back(e2.fr);
19                         col[e2.fr] = bcc_cnt;
20                     }
21                     if (col[e2.to]!=bcc_cnt) {
22                         bcc[bcc_cnt].push_back(e2.to);
23                         col[e2.to] = bcc_cnt;
24                     }
25                 } while (e2.fr != u || e2.to != v);
26             }
27         }
28         else if (dfn[v] < dfn[u] && v != fa) {
29             stk.push(e1);
30             low[u] = min(low[u], dfn[v]);//注意这里必须是dfn[v],而不能是low[v],因为low[v]很可能在另一个双连通分量已经被更新,这个时候用会出错
31         }
32         if (fa < 0 && chile == 1)iscut[u] = false;
33     }
34 }

 

posted @ 2018-09-13 20:22  MekakuCityActor  阅读(1454)  评论(0编辑  收藏  举报