【codevs 1332】上白泽慧音

点击就送屠龙宝刀

题目描述

在幻想乡，上白泽慧音是以知识渊博闻名的老师。春雪异变导致人间之里的很多道路都被大雪堵塞，使有的学生不能顺利地到达慧音所在的村庄。因此慧音决定换一个能够聚集最多人数的村庄作为新的教学地点。人间之里由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行：若干个整数，依次输出最大的绝对连通区域所包含的村庄编号。

样例输入

5 5
1 2 1
1 3 2
2 4 2
5 1 2
3 5 1

样例输出

3
1 3 5

数据范围

对于60%的数据：N <= 200且M <= 10,000
对于100%的数据：N <= 5,000且M <= 50,000
scc模板题

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <vector>
#include <stack>
using namespace std;
const int MAXN = 50000 + 5;
struct edge
{
    int f,t;
}l[MAXN << 1];
int first[MAXN],nxt[MAXN << 1],tot;
void init()
{
    memset(first,0xfff,sizeof(first));
    tot = 0;
    return;
}
void build(int f,int t)
{
    l[++tot] = (edge){f,t};
    nxt[tot] = first[f];
    first[f] = tot;
    return;
}
int low[MAXN],dfn[MAXN];
int dfs_clock=0,scc_num[MAXN],scc_cnt=0;
stack<int> s;
void dfs(int u)
{
    dfn[u]=low[u]=++dfs_clock;
    s.push(u);
    for(int i=first[u];~i;i=nxt[i])
    {
        int v=l[i].t;
        if(!dfn[v])
        {
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!scc_num[v])
        {
            low[u]=min(low[u],dfn[v]);
        }
    }
    if(low[u]==dfn[u])
    {
        scc_cnt++;
        while(true)
        {
            int x=s.top();
            s.pop();
            scc_num[x]=scc_cnt;
            if(x==u) break;
        }
    }
    return;
}
bool cmp(vector <int> a,vector <int> b)
{
    return a.size() > b.size() || (a.size() == b.size() && a[0] < b[0]);
}
int n,m;
int q,f,t;
vector <int> scc[5000];
int main()
{
    init();
    scanf("%d %d",&n,&m);
    for(int i = 1;i <= m;i ++)
    {
        scanf("%d %d %d",&f,&t,&q);
        build(f,t);
        if(q == 2)
            build(t,f);
    }
    for(int i=1;i<=n;i++)
        if(!dfn[i])
            dfs(i);
    for(int i = 1;i <= n;i ++)
        scc[scc_num[i]].push_back(i);
    sort(scc + 1,scc + scc_cnt + 1,cmp);
    printf("%d\n",scc[1].size());
    for(int i = 0;i < scc[1].size();i ++)
        printf("%d ",scc[1][i]);
    puts("");
    return 0;
}