强连通分量scc

  • Post author:
  • Post category:其他


给出一个有向图 , 要求这个图中的所有强连通分量 (一个单独的点也是强连通分量 , 这点和无向图的双连通分量不同 , 他是一条边就是双连通分量)

这里延伸出一个题目 , 给出一个有向图 , 请问最小需要连接多少边才能让这个图变成强连通图 。

解法:先把图中的强连通分量都求出来 , 然后再把每个强连通分量缩成一个点 , 就能得到一个DAG图 , 然后最小值 = max(a , b) , a表示入度为0的点 , b表示进度为0的点。

要注意的一点:当整个图就是一个强连通图时 , 则需要增加的边为0;

代码:

#include

#include

#include

#include

#include

#include

using namespace std;

const int MAXN = 1000;

vectorgrap[MAXN] ;

int pre[MAXN] , lowlink[MAXN] , sccno[MAXN] , dfs_clock , scc_cnt;

stacks;

int n , m;

void init()

{




for(int i = 0 ; i <= n; i++)








grap[i].clear();




memset(pre , 0 , sizeof(pre));




memset(sccno , 0 , sizeof(sccno));




dfs_clock = scc_cnt = 0;

}

void dfs(int u)

{




pre[u] = lowlink[u] = ++dfs_clock;




s.push(u);




for(int i = 0 ; i < grap[u].size() ; i++)




{








int v = grap[u][i];








if(!pre[v])








{












dfs(v);












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








}








else if(!sccno[v])








{












lowlink[u] = min(lowlink[u] , pre[v]) ;








}




}




if(lowlink[u] == pre[u])




{








scc_cnt++;








for(; 😉








{












int x = s.top() ; s.pop();












sccno[x] = scc_cnt;












if(x == u)

break;








}




}

}

int main()

{




while(scanf(“%d %d” , &n , &m) != EOF)




{








init();








int i , x , y;








for(i = 0 ; i < m; i++)








{












scanf(“%d %d” , &x , &y);












grap[x].push_back(y);








}








dfs(0);








int in0[MAXN] , out0[MAXN] ;








memset(in0 , 0 , sizeof(in0));








memset(out0 , 0 , sizeof(out0));








for(int u = 0 ; u < n ; u++)












for(int i = 0 ; i < grap[u].size() ; i++)












{
















int v = grap[u][i] ;
















if(sccno[u] != sccno[v])

in0[sccno[v]] = out0[sccno[u]] = 1;












}








int a = 0 , b = 0;








for(int i = 1; i <= scc_cnt ; i++)








{












if(!in0[i])

a++;












if(out0[i])

b++;








}








int ans = max(a , b);








if(scc_cnt == 1)

ans = 0;








printf(“%d\n” , ans);




}




return 0;

}



版权声明:本文为zengchenacmer原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。