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：复制

输出样例#1：复制

3

1 3 5

说明

对于60%的数据：N <= 200且M <= 10,000

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

tarjan塔尖算法强联通

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

#define inf 2147483647

const ll INF = 0x3f3f3f3f3f3f3f3fll;

#define ri register int

template <class T> inline T min(T a, T b, T c)

{

    return min(min(a, b), c);

}

template <class T> inline T max(T a, T b, T c)

{

    return max(max(a, b), c);

}

template <class T> inline T min(T a, T b, T c, T d)

{

    return min(min(a, b), min(c, d));

}

template <class T> inline T max(T a, T b, T c, T d)

{

    return max(max(a, b), max(c, d));

}

#define scanf1(x) scanf("%d", &x)

#define scanf2(x, y) scanf("%d%d", &x, &y)

#define scanf3(x, y, z) scanf("%d%d%d", &x, &y, &z)

#define scanf4(x, y, z, X) scanf("%d%d%d%d", &x, &y, &z, &X)

#define pi acos(-1)

#define me(x, y) memset(x, y, sizeof(x));

#define For(i, a, b) for (int i = a; i <= b; i++)

#define FFor(i, a, b) for (int i = a; i >= b; i--)

#define bug printf("***********\n");

#define mp make_pair

#define pb push_back

const int N = ;

const int M=;

// name*******************************

int dfn[N],low[N];

int ins[N];

int link[N],gett[N];

int fst[N];

struct edge

{

    int to,nxt;

} e[M];

int cnt=,tot=;

stack<int>S;

int n,m,num=;

// function******************************

void add(int u,int v)

{

    e[++tot].to=v;

    e[tot].nxt=fst[u];

    fst[u]=tot;

}

void tarjan(int u)

{

    dfn[u]=low[u]=++cnt;

    ins[u]=;

    S.push(u);

    for(int p=fst[u]; p; p=e[p].nxt)

    {

        int v=e[p].to;

        if(dfn[v])

        {

            if(ins[v])

                low[u]=min(low[u],dfn[v]);

        }

        else

        {

            tarjan(v);

            low[u]=min(low[u],low[v]);

        }

    }

    if(dfn[u]==low[u])

    {

        num++;

        while(S.top()!=u)

        {

            int t=S.top();

            S.pop();

            link[t]=num;

            gett[num]++;

            ins[t]=;

        }

        int t=S.top();

        S.pop();

        link[t]=num;

        gett[num]++;

        ins[t]=;

    }

}

//***************************************

int main()

{

//    ios::sync_with_stdio(0);

//    cin.tie(0);

    // freopen("test.txt", "r", stdin);

    //  freopen("outout.txt","w",stdout);

    cin>>n>>m;

    For(i,,m)

    {

        int u,v,t;

        scanf("%d%d%d",&u,&v,&t);

        if(t==)

            add(u,v);

        else

        {

            add(u,v);

            add(v,u);

        }

    }

    For(i,,n)

    if(link[i]==)

        tarjan(i);

    int pos,maxx=;

    For(i,,num)

    if(maxx<gett[i])

    {

        maxx=gett[i];

        pos=i;

    }

    cout<<maxx<<endl;

    For(i,,n)

    {

        if(link[i]==pos)

            cout<<i<<" ";

    }

    return ;

}