洛谷 P1726 上白泽慧音 tarjan求强连通分量

https://www.luogu.org/problem/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行：若干个整数，依次输出最大的绝对连通区域所包含的村庄编号。

思路：很明显题目要求最大的强连通分量，那么求出所有的强连通分量再取最大的那个就可以了，所以要用到$tarjan$算法，说一下字典序的问题，因为任意两个强连通分量都没有交集，所以只要知道每个强连通分量的最小的节点编号就可以确定答案了。

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define eps 1e-8
#define pr pair<int,int>
using namespace std;
typedef long long ll;

const int maxn=5005;

int s[maxn][maxn];
int dfn[maxn],low[maxn],Stack[maxn],ins[maxn],id[maxn];
vector<int> scc[maxn];
int n,m,num,top,cnt,ans;

void tarjan(int x)
{
    dfn[x]=low[x]=++num;
    Stack[++top]=x,ins[x]=1;
    for(int i=1;i<=n;i++)
    {
        if(!s[x][i])
            continue ;
        if(!dfn[i])
        {
            tarjan(i);
            low[x]=min(low[x],low[i]);
        }
        else if(ins[i])
            low[x]=min(low[x],dfn[i]);
    }
    if(dfn[x]==low[x])
    {
        ++cnt;
        int y;
        do
        {
            y=Stack[top--],ins[y]=0;
            id[y]=cnt,scc[cnt].push_back(y);
        }while(x!=y);
        if(scc[cnt].size()>ans)
            ans=scc[cnt].size();
    }
}

int main()
{
    scanf("%d %d",&n,&m);
    int u,v,t;
    for(int i=0;i<m;i++)
    {
        scanf("%d %d %d",&u,&v,&t);
        if(t==1)
            s[u][v]=1;
        else
            s[u][v]=s[v][u]=1;
    }
    for(int i=1;i<=n;i++)
        if(!dfn[i])
            tarjan(i);
    int len=0,MIN=INF,re=INF,idx=0;
    for(int i=1;i<=cnt;i++)
    {
        if(scc[i].size()==ans)
        {
            MIN=INF;
            for(int j=0;j<ans;j++)
                MIN=min(MIN,scc[i][j]);
            if(MIN<re)
                re=MIN,idx=i;
        }
    }
    sort(scc[idx].begin(),scc[idx].end());
    printf("%d\n",scc[idx].size());
    for(int i=0;i<ans;i++)
        printf("%d ",scc[idx][i]);
    return 0;
}