【noip2010模拟赛】classroom

【noip2010模拟赛】classroom
2017年7月29日
强连通分量

---

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

Input Format
第1行：两个正整数N,M
第2..M+1行：每行三个正整数a,b,t, t=1表示存在从村庄a到b的单向道路，t=2表示村庄a,b之间存在双向通行的道路。保证每条道路只出现一次。
Output Format
第1行： 1个整数，表示最大的绝对连通区域包含的村庄个数。
第2行：若干个整数，依次输出最大的绝对连通区域所包含的村庄编号。

---

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

int N, M;
struct Edge{
    int x, y, next;
}E[51000];    int ans = 0;
int V[11000];
//
int dfn[11000], low[11000], ind = 0;//DFS
int stack[11000], tails = 0;//栈
bool Jud[11000];
int numSCC[11000], sum = 0, countsofSCC[11000];

void AssignEdge(int xi, int yi)
{
    E[++ans].next = V[xi];    E[ans].x = xi;    E[ans].y = yi;
    V[xi] = ans;
}

void Putin()
{
    memset(V, -1, sizeof(V));
    memset(Jud, false, sizeof(Jud));
    memset(dfn, 0, sizeof(dfn));
    memset(low, 0, sizeof(low));
    memset(countsofSCC, 0, sizeof(countsofSCC));
    //
    cin >> N >> M;
    for(int i = 1, a, b, k; i <= M; i++){
        cin >> a >> b >> k;
        if(k == 1)  AssignEdge(a, b);
        else    AssignEdge(a, b), AssignEdge(b, a);
    }
}

void SCC(int x)
{
    dfn[x] = low[x] = ++ind;
    stack[++tails] = x;
    Jud[x] = true;
    for(int i = V[x]; i > 0; i = E[i].next)
    {
        int y = E[i].y;
        if(!dfn[y]){
            SCC(y);
            low[x] = min(low[x], low[y]);
        }
        else
            if(Jud[y] == true)    low[x] = min(low[x], dfn[y]);
    }
    if(low[x] >= dfn[x])
    {
        int k;    sum++;
        do{
            k = stack[tails--];
            Jud[k] = false;
            numSCC[k] = sum;
            countsofSCC[sum]++;
         } while(k != x);
    }
}

void Printout()
{
    int maxx = -1, outteam;
    for(int i = 1; i <= sum; i++)
        maxx = max(maxx, countsofSCC[i]);
    cout << maxx << endl;
    for(int i = 1; i <= N; i++)
        if(countsofSCC[numSCC[i]] == maxx)
            outteam = numSCC[i];
    for(int i = 1; i <= N; i++)
        if(numSCC[i] == outteam)
            cout << i << ' ';
    cout << endl;
}

int main()
{
    Putin();
    for(int i = 1; i <= N; i++)
        if(dfn[i] == 0)    SCC(i);
    Printout();
    return 0;
}