洛谷P1726 上白泽慧音（Tarjan强连通分量）

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强联通分量
处理的时候两个数组记录答案，一个是当前的答案，一个是以前的答案
更新的时候只要当前节点多或者节点数相同当前字典序小就可以更新
每次更新前数组初始别忘了赋0
输出答案的时候也不要把0输出来 
*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define maxn 50001
#define inf 0x3f3f3f3f

using namespace std;
int dis[maxn],head[maxn],vis[maxn];
int dfn[maxn],low[maxn],stack[maxn],Ans1[maxn],Ans2[maxn];
bool in_stack[maxn];
int n,m,s,t,x,y,z,num,ans,cnt,top;
struct node
{
    int u,v,w,next;
}e[maxn<<1];

inline int init()
{
    int x=0,f=1;char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}

inline void add(int u,int v)
{
    e[++num].u=u;
    e[num].v=v;
    e[num].next=head[u];
    head[u]=num;
}

void Tarjan(int u)
{
    dfn[u]=low[u]=++cnt;
    stack[++top]=u;in_stack[u]=true;
    for(int i=head[u];i;i=e[i].next)
    {
        int v=e[i].v;
        if(!dfn[v]) 
        {
            Tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(in_stack[v])
        {
            low[u]=min(low[u],dfn[v]);
        }
    }
    if(dfn[u]==low[u])
    {
        int cnt1=0;
        while(stack[top]!=u)
        {
            cnt1++;
            in_stack[stack[top]]=false;
            Ans1[cnt1]=stack[top];//当前答案数组 
            top--;
        }
        cnt1++;
        in_stack[stack[top]]=false;
        Ans1[cnt1]=stack[top];
        top--;        
        if(cnt1>ans)//节点数多 
        {
            memset(Ans2,0,sizeof Ans2);//以前答案数组 
            ans=cnt1;
              for(int i=1;i<=cnt1;i++)
                Ans2[i]=Ans1[i];
            sort(Ans2+1,Ans2+n+1);
        }
        else if(cnt1==ans)//节点数相同 
        {
            sort(Ans1+1,Ans1+n+1);
            if(Ans2[1]>Ans1[1])//字典序小 
            memset(Ans2,0,sizeof Ans2);
              for(int i=1;i<=cnt1;i++)
                Ans2[i]=Ans1[i];
        }
    }
}

int main()
{
    n=init();m=init();
    for(int i=1;i<=m;i++)
    {
        x=init();y=init(),z=init();
        if(z==1)add(x,y);
        else add(x,y),add(y,x);
    }
    for(int i=1;i<=n;i++)
     if(!dfn[i]) Tarjan(i);
    printf("%d
",ans);
    for(int i=1;i<=n;i++)
      if(Ans2[i]!=0) printf("%d ",Ans2[i]);//不要输出零 
    return 0;
}