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行:若干个整数,依次输出最大的绝对连通区域所包含的村庄编号。
输入输出样例
说明
对于60%的数据:N <= 200且M <= 10,000
对于100%的数据:N <= 5,000且M <= 50,000
tarjan塔尖算法 强联通
#include <bits/stdc++.h> using namespace std; typedef long long ll; #define inf 2147483647 const ll INF = 0x3f3f3f3f3f3f3f3fll; #define ri register int template <class T> inline T min(T a, T b, T c) { return min(min(a, b), c); } template <class T> inline T max(T a, T b, T c) { return max(max(a, b), c); } template <class T> inline T min(T a, T b, T c, T d) { return min(min(a, b), min(c, d)); } template <class T> inline T max(T a, T b, T c, T d) { return max(max(a, b), max(c, d)); } #define scanf1(x) scanf("%d", &x) #define scanf2(x, y) scanf("%d%d", &x, &y) #define scanf3(x, y, z) scanf("%d%d%d", &x, &y, &z) #define scanf4(x, y, z, X) scanf("%d%d%d%d", &x, &y, &z, &X) #define pi acos(-1) #define me(x, y) memset(x, y, sizeof(x)); #define For(i, a, b) for (int i = a; i <= b; i++) #define FFor(i, a, b) for (int i = a; i >= b; i--) #define bug printf("***********\n"); #define mp make_pair #define pb push_back const int N = 5010; const int M=200005; // name******************************* int dfn[N],low[N]; int ins[N]; int link[N],gett[N]; int fst[N]; struct edge { int to,nxt; } e[M]; int cnt=0,tot=0; stack<int>S; int n,m,num=0; // function****************************** void add(int u,int v) { e[++tot].to=v; e[tot].nxt=fst[u]; fst[u]=tot; } void tarjan(int u) { dfn[u]=low[u]=++cnt; ins[u]=1; S.push(u); for(int p=fst[u]; p; p=e[p].nxt) { int v=e[p].to; if(dfn[v]) { if(ins[v]) low[u]=min(low[u],dfn[v]); } else { tarjan(v); low[u]=min(low[u],low[v]); } } if(dfn[u]==low[u]) { num++; while(S.top()!=u) { int t=S.top(); S.pop(); link[t]=num; gett[num]++; ins[t]=0; } int t=S.top(); S.pop(); link[t]=num; gett[num]++; ins[t]=0; } } //*************************************** int main() { // ios::sync_with_stdio(0); // cin.tie(0); // freopen("test.txt", "r", stdin); // freopen("outout.txt","w",stdout); cin>>n>>m; For(i,1,m) { int u,v,t; scanf("%d%d%d",&u,&v,&t); if(t==1) add(u,v); else { add(u,v); add(v,u); } } For(i,1,n) if(link[i]==0) tarjan(i); int pos,maxx=0; For(i,1,num) if(maxx<gett[i]) { maxx=gett[i]; pos=i; } cout<<maxx<<endl; For(i,1,n) { if(link[i]==pos) cout<<i<<" "; } return 0; }