codevs-1332 上白澤慧音
阿新 • • 發佈:2018-12-08
題意不用多說,有向圖,求強聯通分量。找到強連通分量中定點最多的分量,並輸出該最大分量中包含的頂點。
是一道求強連通分量的模板題,用的是Kosaraju演算法。
Kosaraju演算法步驟:
1.先根據題意建原圖與原圖的反圖
2.根據原圖進行第一次DFS,(DFS1)得到一個ord[]陣列,存放DFS1的遍歷順序
3.再根據ord[]陣列,按照第一次的遍歷順序去DFS(DFS2)反圖,染色。染色數就是強連通分量的數量,每種顏色的使用數量就是這種顏色(該強連通分量)包含的頂點。
/* *looooop * Do not go gentle into that good night * -Dylan Thomas */ #include <iostream> #include <stdio.h> #include <string.h> #include <stack> #include <queue> #include <map> #include <set> #include <vector> #include <math.h> #include <bitset> #include <algorithm> #include <climits> using namespace std; #define lson 2*i #define rson 2*i+1 #define LS l,mid,lson #define RS mid+1,r,rson #define UP(i,x,y) for(i=x;i<=y;i++) #define DOWN(i,x,y) for(i=x;i>=y;i--) #define MEM(a,x) memset(a,x,sizeof(a)) #define W(a) while(a) #define gcd(a,b) __gcd(a,b) #define LL long long #define N 1000005 #define MOD 1000000007 #define INF 0x3f3f3f3f #define EXP 1e-8 #define lowbit(x) (x&-x) #define MAX 10005 vector<int>G[MAX],re_G[MAX]; int ord[MAX]; //正向搜尋的dfs次序 int num[MAX]; // int vis[MAX]; int belong[MAX]; //當前頂點屬於哪個集合,相當於染色,當前頂點被染成了什麼顏色 int color;int n,m;int out[MAX]; //out是縮點轉換成DAG後每個縮點的出度 int ans[MAX]; //每種顏色包含多少個頂點,即強連通數量 int no; //正向搜尋的編號 void dfs1(int st){ //dfs1用來遍歷原圖得到dfs序 vis[st] = 1; for(int i = 0; i < G[st].size(); i++){ int v = G[st][i]; if(!vis[v]) dfs1(v); } ord[no++] = st; } void dfs2(int st){ vis[st] = 1; belong[st] = color; for(int i = 0; i < re_G[st].size(); i++){ int v = re_G[st][i]; if(!vis[v]){ ans[color]++; dfs2(v); //color++; } } } void Kosaraju(){ MEM(vis,0); MEM(ord,0); MEM(belong,0); MEM(out,0); no = 1;color = 1; for(int i = 1; i<=MAX; i++) ans[i] = 1; for(int i = 1; i <= n; i++){ if(!vis[i]) vis[i] = 1,dfs1(i); } MEM(vis,0); for(int i = no-1; i>= 1; i--){ int v = ord[i]; if(!vis[v]){ dfs2(v); color++; } } int temp=1; for(int i = 1; i <= color; i++){ if(ans[i] >= ans[temp]) temp = i; } printf("%d\n",ans[temp]); for(int i = 1; i <= n; i++){ if(belong[i] == temp) printf("%d ",i); } printf("\n"); } int main(int argc,char *argv[]){ scanf("%d%d",&n,&m); for(int i = 0 ;i < m; i++){ int x,y,t; scanf("%d%d%d",&x,&y,&t); G[x].push_back(y); re_G[y].push_back(x); if(t==2){ G[y].push_back(x); re_G[x].push_back(y); } } Kosaraju(); return 0; }