【codevs 1332】上白泽慧音

题目描述 Description
在幻想乡，上白泽慧音是以知识渊博闻名的老师。春雪异变导致人间之里的很多道路都被大雪堵塞，使有的学生不能顺利地到达慧音所在的村庄。因此慧音决定换一个能够聚集最多人数的村庄作为新的教学地点。人间之里由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。

输入描述 Input Description
第1行：两个正整数N,M

第2..M+1行：每行三个正整数a,b,t, t = 1表示存在从村庄a到b的单向道路，t = 2表示村庄a,b之间存在双向通行的道路。保证每条道路只出现一次。

输出描述 Output Description
第1行： 1个整数，表示最大的绝对连通区域包含的村庄个数。

第2行：若干个整数，依次输出最大的绝对连通区域所包含的村庄编号。

样例输入 Sample Input
5 5

1 2 1

1 3 2

2 4 2

5 1 2

3 5 1

样例输出 Sample Output
3

1 3 5

数据范围及提示 Data Size & Hint
对于60%的数据：N <= 200且M <= 10,000

对于100%的数据：N <= 5,000且M <= 50,000

找出最大的强连通分量，输出。

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<stack>using namespace std;const int maxn=50000+5;int first[maxn],nxt[maxn],tot=0;int dfn[maxn],low[maxn],dfs_clock,scc_num;int scc[maxn],size[maxn];int n,m;stack<int>s;struct edge{    int f,t;}l[maxn];void build(int f,int t){    l[++tot]=(edge){f,t};    nxt[tot]=first[f];    first[f]=tot;    return;}void dfs(int u){    low[u]=dfn[u]=++dfs_clock;    s.push(u);    for(int i=first[u];i;i=nxt[i])    {        int w=l[i].t;        if(!dfn[w])        {            dfs(w);            low[u]=min(low[u],low[w]);        }        else if(!scc[w])        {            low[u]=min(low[u],dfn[w]);        }    }    if(low[u]==dfn[u])    {        scc_num++;        while(!s.empty())        {            int x=s.top();            s.pop();            scc[x]=scc_num;            size[scc_num]++;            if(x==u) break;        }    }    return;}int main(){    int f,t,k;    scanf("%d%d",&n,&m);    for(int i=1;i<=m;i++)    {        scanf("%d%d%d",&f,&t,&k);        if(k==1) build(f,t);        else if(k==2)        {            build(f,t);            build(t,f);        }    }    for(int i=1;i<=n;i++)    {        if(!dfn[i]) dfs(i);    }    int sz=0;    for(int i=1;i<=scc_num;i++)    {        sz=max(sz,size[i]);    }    cout<<sz<<'\n';    int g;    for(int i=1;i<=scc_num;i++)    {        if(size[i]==sz)        {            g=i;            break;        }    }    for(int i=1;i<=n;i++)    {        if(scc[i]==g) printf("%d ",i);    }    return 0;}