题意不用多说,有向图,求强联通分量。找到强连通分量中定点最多的分量,并输出该最大分量中包含的顶点。
是一道求强连通分量的模板题,用的是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;
}