luoguP1726 上白泽慧音

上白泽慧音

题目描述

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

坑点:

　　1)更新的时候只要当前节点多或者节点数相同并且当前字典序小就可以更新答案
　　2)每次更新前的时候数组初始别忘了赋0!
　　3)输出答案的时候千万不要把0输出来!

上代码:

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#define INF 99999
using namespace std;

const int N = 5050;
const int M = 100050;///双向边 
int n,m,ans1[N],ans2[N];
///ans1数组存放现在的答案,ans2数组存放以前最优答案 
int answer;
int top,head[N];
int DFN[M],LOW[M];
int index0,tot,stack[M];
bool vis[M];
/*
注意: 
1)更新的时候只要当前节点多
或者节点数相同并且当前字典序小就可以更新答案 
2)每次更新前数组初始别忘了赋0
3)输出答案的时候也不要把0输出来 
*/
struct A {
    int next,v;
}t[M];

void add(int u,int v)
{///链表 
    ++top;
    t[top].v=v;
    t[top].next=head[u];
    head[u]=top;
}

void tarjan(int x)
{///模板 
    DFN[x]=LOW[x]=++tot;
    stack[++index0]=x;
    vis[x]=true;
    for(int i=head[x];i!=-1;i=t[i].next)
    {
        int v=t[i].v;
        if(!DFN[v])
        {
            tarjan(v);
            LOW[x]=min(LOW[x],LOW[v]);
        }
        else if(vis[v])
            LOW[x]=min(LOW[x],DFN[v]);
    }
    int y;
    if(LOW[x]==DFN[x])
    {
        int cnt=0;
        do
        {
            cnt++;
            y=stack[index0--];
            vis[y]=false;
            ans1[cnt]=y;
        } while(y!=x);
        if(cnt>answer)
        {
            memset(ans2,0,sizeof(ans2));
            answer=cnt;
            for(int i=1;i<=cnt;++i)
                ans2[i]=ans1[i];///更新 
            sort(ans2+1,ans2+1+cnt);///+cnt
        }
        else if(cnt==answer)
        {
            sort(ans1+1,ans1+1+cnt);
            if(ans2[1]>ans1[1])///选取字典序较小的进行更新 
                memset(ans2,0,sizeof(ans2));
            for(int i=1;i<=cnt;++i)
                ans2[i]=ans1[i];
        }
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    memset(head,-1,sizeof(head));
    for(int i=1,a,b,t;i<=m;++i)
    {
        scanf("%d%d%d",&a,&b,&t);
        add(a,b);
        if(t==2)
            add(b,a);
    }
    for(int i=1;i<=n;++i)
        if(!DFN[i])
            tarjan(i);
    printf("%d\n",answer);
    for(int i=1;i<=n;++i)
    {
        if(ans2[i]==0)///不能输出0 
            break;
        printf("%d ",ans2[i]);
    }
    return 0;
}