洛谷1726 上白泽慧音

洛谷1726 上白泽慧音
本题地址： http://www.luogu.org/problem/show?pid=1726 题目描述
在幻想乡，上白泽慧音是以知识渊博闻名的老师。春雪异变导致人间之里的很多道路都被大雪堵塞，使有的学生不能顺利地到达慧音所在的村庄。因此慧音决定换一个能够聚集最多人数的村庄作为新的教学地点。人间之里由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算法求出scc（强连通分量），判断比较最大值即可。需要注意的是rank的比较，这里以scc中的最小标号作为scc的rank。

【代码】

1 #include<cstdio>
2 #include<cstring>
3 #include<stack>
4 using namespace std;
5
6 const int maxn = 5000+10 , maxm=50000+10;
7 struct Edge{
8     int v,w,next;
9 }e[2*maxm];
10 int en=-1,front[maxn];
11
12 int n,m;
13
14 inline void AddEdge(int u,int v) {
15     en++; e[en].v=v; e[en].next=front[u]; front[u]=en;
16 }
17 stack<int> S;
18 int scc_cnt,dfs_clock;
19 int sccno[maxn],pre[maxn],lowlink[maxn],scc_N[maxn],scc_rank[maxn];
20
21 void dfs(int u) {
22     pre[u]=lowlink[u]=dfs_clock++;
23     S.push(u);
24     for(int i=front[u];i>=0;i=e[i].next) {
25         int v=e[i].v;
26         if(!pre[v]) {
27             dfs(v);
28             lowlink[u]=min(lowlink[u],lowlink[v]);
29         }
30         else if(!sccno[v]) {
31             lowlink[u]=min(pre[v],lowlink[u]);
32         }
33     }
34     if(pre[u]==lowlink[u]) {
35         scc_cnt++;
36         for(;;) {
37             int v=S.top(); S.pop();
38             sccno[v]=scc_cnt;
39             if(v==u) break;
40         }
41     }
42 }
43
44 void find_scc(int n) {
45     dfs_clock=scc_cnt=0;
46     memset(pre,0,sizeof(pre));
47     memset(sccno,0,sizeof(sccno));
48     for(int i=1;i<=n;i++) if(!pre[i])
49        dfs(i);
50 }
51
52 int main() {
53     memset(front,-1,sizeof(front));
54     scanf("%d%d",&n,&m);
55     int u,v,w;
56     for(int i=0;i<m;i++) {
57         scanf("%d%d%d",&u,&v,&w);
58         if(w==1) AddEdge(u,v);
59         else {
60             AddEdge(u,v);
61             AddEdge(v,u);
62         }
63     }
64     find_scc(n);
65     for(int i=1;i<=n;i++) {
66         scc_N[sccno[i]]++;
67         if(!scc_rank[sccno[i]]) scc_rank[sccno[i]]=i;
68     }
69     int _max=0,k;
70     for(int i=1;i<=scc_cnt;i++)
71        if(scc_N[i]>_max || scc_N[i]==_max&&scc_rank[i]<scc_rank[k]) {
72               _max=scc_N[k=i];
73        }
74     printf("%d\n",_max);
75     for(int i=1;i<=n;i++){
76         if(sccno[i]==k) printf("%d ",i);
77     }
78     putchar('\n');
79     return 0;
80 }