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塔尖算法 强联通

#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("***********
");
#define mp make_pair
#define pb push_back
const int N = 5010;
const int M=200005;
// 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=0,tot=0;
stack<int>S;
int n,m,num=0;
// 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]=1;
    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]=0;
        }
        int t=S.top();
        S.pop();
        link[t]=num;
        gett[num]++;
        ins[t]=0;
    }
}

//***************************************
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,1,m)
    {
        int u,v,t;
        scanf("%d%d%d",&u,&v,&t);
        if(t==1)
            add(u,v);
        else
        {
            add(u,v);
            add(v,u);
        }
    }
    For(i,1,n)
    if(link[i]==0)
        tarjan(i);
    int pos,maxx=0;
    For(i,1,num)
    if(maxx<gett[i])
    {
        maxx=gett[i];
        pos=i;
    }
    cout<<maxx<<endl;
    For(i,1,n)
    {
        if(link[i]==pos)
            cout<<i<<" ";
    }

    return 0;
}