1. 程式人生 > >codevs-1332 上白澤慧音

codevs-1332 上白澤慧音

題意不用多說,有向圖,求強聯通分量。找到強連通分量中定點最多的分量,並輸出該最大分量中包含的頂點。

是一道求強連通分量的模板題,用的是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;
}