这一看不就直接tarjan求强连通分量，然后最大的强连通分量不就是答案了吗，然后就开始码，然后就过了。。。
在幻想乡，上白泽慧音是以知识渊博闻名的老师。春雪异变导致人间之里的很多道路都被大雪堵塞，使有的学生不能顺利地到达慧音所在的村庄。因此慧音决定换一个能够聚集最多人数的村庄作为新的教学地点。人间之里由N个村庄（编号为1..N）和M条道路组成，道路分为两种一种为单向通行的，一种为双向通行的，分别用1和2来标记。如果存在由村庄A到达村庄B的通路，那么我们认为可以从村庄A到达村庄B，记为(A,B)。当(A,B)和(B,A)同时满足时，我们认为A,B是绝对连通的，记为<A,B>。绝对连通区域是指一个村庄的集合，在这个集合中任意两个村庄X,Y都满足<X,Y>。现在你的任务是，找出最大的绝对连通区域，并将这个绝对连通区域的村庄按编号依次输出。若存在两个最大的，输出字典序最小的，比如当存在1,3,4和2,5,6这两个最大连通区域时，输出的是1,3,4。
输入格式
第1行：两个正整数N,M
第2..M+1行：每行三个正整数a,b,t, t = 1表示存在从村庄a到b的单向道路，t = 2表示村庄a,b之间存在双向通行的道路。保证每条道路只出现一次。
输出格式
第1行： 1个整数，表示最大的绝对连通区域包含的村庄个数。
第2行：若干个整数，依次输出最大的绝对连通区域所包含的村庄编号。
输入输出样例
输入 #1
5 5 1 2 1 1 3 2 2 4 2 5 1 2 3 5 1
输出 #1
3 1 3 5
说明/提示
对于60%的数据：N <= 200且M <= 10,000
对于100%的数据：N <= 5,000且M <= 50,000
#include <cstdio>
#include<iostream>
#include <cstring>
#include <algorithm>
#include <stack>
#include<set>
#include<vector>
using namespace std;
#define ll long long
const int maxn=3e5;
stack<int>q;
set<int>ans;
int cnt,n,m,head[maxn],dfn[maxn],low[maxn],co[maxn],tot,id,out[maxn],root;
vector<int>a[maxn];
struct Edge
{
int to,next;
}e[maxn];
void add(int u,int v)
{
e[cnt].to=v;
e[cnt].next=head[u];
head[u]=cnt++;
}
void init()
{
ans.clear();
memset(out,0,sizeof(out));
memset(co,0,sizeof(co));
memset(head,-1,sizeof(head));
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
while(!q.empty()) q.pop();
cnt=0;id=0,tot=0;
}
void tarjan(int u)
{
q.push(u);
dfn[u]=low[u]=++tot;int ch=0;
for(int i=head[u];i!=-1;i=e[i].next)
{
int v=e[i].to;
if(!dfn[v])
{
ch++;
tarjan(v);
low[u]=min(low[u],low[v]);
if(low[v]>=dfn[u]&&u!=root) ans.insert(u);
}
else if(!co[v])
low[u]=min(low[u],dfn[v]);
}
if(root==u&&ch>=2) ans.insert(u);
if(low[u]==dfn[u])
{
id++;
while(q.top()!=u)
{
int x=q.top();
q.pop();
co[x]=id;
}
co[u]=id;
q.pop();
}
}
bool cmp(vector<int> a,vector<int> b)
{
return a.size()>b.size();
}
int main()
{
scanf("%d %d",&n,&m);
init();
for(int i=1;i<=m;i++)
{
int u,v;int op;scanf("%d %d",&u,&v);
scanf("%d",&op);
add(u,v); if(op==2)add(v,u);
}
for(int i=1;i<=n;i++)
{
if(dfn[i]==0) {
root=i;
tarjan(i);
}
}
for(int i=1;i<=n;i++) a[co[i]].push_back(i);
sort(a+1,a+1+id,cmp);
sort(a[1].begin(),a[1].end());
printf("%d\n",(int)a[1].size());
for(int i=0;i<a[1].size();i++)
{
printf("%d",a[1][i]);
printf("%c",i==a[1].size()-1?'\n':' ');
}
}