上白泽慧音（luogu P1726

                        题目传送门

题目描述

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

                 

　　tarjan模板题。（学习tarjan推荐byvoid的博客:传送门）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<stack>
#include<bitset>
#define LL long long
#define RI register int
using namespace std;
const int INF = 0x7ffffff ;
const int N = 5000 + 10 ;
const int M = 50000 + 10 ;

inline int read() {
    int k = 0 , f = 1 ; char c = getchar() ;
    for( ; !isdigit(c) ; c = getchar())
      if(c == '-') f = -1 ;
    for( ; isdigit(c) ; c = getchar())
      k = k*10 + c-'0' ;
    return k*f ;
}
struct Edge {
    int to, next ;
}e[M<<1] ;
int n, m, dclock, tot ; int head[N], dfn[N], low[N], num[N], bl[N] ;
stack<int>s ;  bitset<N>ins ;
inline void add_edge(int x,int y) {
    static int cnt = 1 ;
    e[++cnt].to = y, e[cnt].next = head[x], head[x] = cnt ;
}

void tarjan(int x) {
    dfn[x] = low[x] = ++dclock ; s.push(x) ; ins[x] = 1 ;
    for(int i=head[x];i;i=e[i].next) {
        int y = e[i].to ;
        if(!dfn[y]) {
            tarjan(y) ; low[x] = min(low[x],low[y]) ;
        } else if(ins[y]) low[x] = min(low[x],dfn[y]) ;
    }
    if(dfn[x] == low[x]) {
        tot++ ;
        while(1) {
            int xx = s.top() ; s.pop() ; ins[xx] = 0 ; 
            bl[xx] = tot, num[tot]++ ;
            if(xx == x) break ;
        }
    }
}

int main() {
    n = read(), m = read() ;
    for(int i=1;i<=m;i++) {
        int x = read(), y = read(), t = read() ;
        add_edge(x,y) ; if(t == 2) add_edge(y,x) ;
    }
    for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i) ;
    int maxx = 0 ;
    for(int i=1;i<=tot;i++) maxx = max(maxx,num[i]) ;
    int k ;
    for(int i=1;i<=n;i++) if(num[bl[i]] == maxx) { k = bl[i] ; break ; }
    printf("%d
",maxx) ;
    for(int i=1;i<=n;i++) if(bl[i] == k) printf("%d ",i) ;
    return 0 ;
}