洛谷 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自动忽略大环中的小环（yeah）

---

【代码】

//洛谷 上白泽慧音
#include<stack>
#include<queue> 
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=100005;
stack <int> s;
struct edge {int to,next;} f[mxn];
int head[mxn],dfn[mxn],low[mxn],siz[mxn],cun[mxn],be[mxn];
bool in[mxn];
int n,m,a,b,t;
int cnt,ans,tim,num;
inline void add(int u,int v)
{
    f[++cnt].to=v,f[cnt].next=head[u],head[u]=cnt;
}
inline void tarjan(int u)
{
    int i,v;
    dfn[u]=low[u]=++tim;
    in[u]=1;
    s.push(u);
    for(i=head[u];i;i=f[i].next)
    {
        v=f[i].to;
        if(!dfn[v]) tarjan(v),low[u]=min(low[u],low[v]);
        else if(in[v]) low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u])
    {
        num++;
        do
        {
            v=s.top();
            s.pop();
            in[v]=0;
            be[v]=num;
            siz[num]++;
        }while(u!=v);
    }
}
int main()
{
    int i,j,k,u,v,ans=0;
    scanf("%d%d",&n,&m);
    while(m--)
    {
        scanf("%d%d%d",&a,&b,&t);
        if(t==1) add(a,b);
        else add(a,b),add(b,a);
    }
    fo(i,1,n)
      if(!dfn[i])
        tarjan(i);
    fo(i,1,num) ans=max(ans,siz[i]);
    memset(cun,0x3f,sizeof cun);
    fo(i,1,num) if(ans==siz[i])
    {
        int tmp=0;
        fo(j,1,n) if(be[j]==i) break;
        if(j<cun[1]) fo(k,j,n) if(be[k]==i) cun[++tmp]=k;
    }
    printf("%d\n",ans);
    fo(i,1,ans) printf("%d ",cun[i]);
    printf("\n");
    return 0;
}