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行：若幹個整數，依次輸出最大的絕對連通區域所包含的村莊編號。

輸入輸出樣例

5 5
1 2 1
1 3 2
2 4 2
5 1 2
3 5 1

3
1 3 5

說明

對於60%的數據：N <= 200且M <= 10,000

對於100%的數據：N <= 5,000且M <= 50,000

看了一下題解，然後粗略的看了一下提交記錄，發現很少用stl去寫棧的，

很多人調試過後認為不能用stl去寫，但其實是可以的。

我們在手寫棧的時候的判斷條件是：while(x!=stack[index+1]);

這樣實際上我們是把當前元素和上一個元素進行比較。

如果改用stl的話，我們可以和當前元素比較，但是在比較完之後，我們還要再與棧頂比較一次！

同時要註意vis數組的撤銷

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<cmath>
  5 #include<queue>
  6 #include<algorithm>
  7 #include<stack> 
  8 #define lli long long int 
  9 using namespace std;
 10 const int MAXN=100001;
 11 const int maxn=0x7fffff;
 12 inline void read(int &n)
 13 {
 14     char c='+';int x=0;bool flag=0;
 15     while(c<'0'||c>'9')
 16     {c=getchar();if(c=='-')flag=1;}
 17     while(c>='0'&&c<='9')
 18     {x=(x<<1)+(x<<3)+c-48;c=getchar();}
 19     flag==1?n=-x:n=x;
 20 }
 21 struct node
 22 {
 23     int u,v,nxt;
 24 }edge[MAXN*2];
 25 int head[MAXN];
 26 int num=1;
 27 int n,m;
 28 int dfn[MAXN];
 29 int low[MAXN];
 30 int vis[MAXN];// 是否在棧內 
 31 void add_edge(int x,int y)
 32 {
 33     edge[num].u=x;
 34     edge[num].v=y;
 35     edge[num].nxt=head[x];
 36     head[x]=num++;
 37 }
 38 int tot=0;
 39 stack<int>s;
 40 int cur[MAXN];
 41 int curnum;
 42 int ans[MAXN];
 43 int ansnum;
 44 void tarjan(int now)
 45 {
 46     dfn[now]=low[now]=++tot;
 47     s.push(now);
 48     vis[now]=1;
 49     for(int i=head[now];i!=-1;i=edge[i].nxt)
 50     {
 51         if(!dfn[edge[i].v])
 52         {
 53             tarjan(edge[i].v);
 54             low[now]=min(low[now],low[edge[i].v]);
 55         }
 56         else if(vis[edge[i].v])
 57             low[now]=min(low[now],dfn[edge[i].v]);
 58     }
 59     if(low[now]==dfn[now])
 60     {
 61         curnum=0;
 62         int tmp=-1;
 63         while(now!=s.top())
 64         {
 65             cur[++curnum]=s.top();
 66             vis[s.top()]=0;
 67             s.pop();
 68             if(tmp==now)
 69                 break;
 70         }
 71         vis[s.top()]=0;
 72         cur[++curnum]=s.top();
 73         s.pop();    
 74         if(curnum<ansnum)
 75             return ;
 76         sort(cur,cur+curnum+1);
 77         if(curnum>ansnum)
 78         {
 79             for(int i=1;i<=curnum;i++)
 80                 ans[i]=cur[i];
 81             ansnum=curnum;
 82         }
 83         else
 84         {
 85             for(int i=1;i<=curnum;i++)
 86             {
 87                 if(cur[i]<ans[i])
 88                 {
 89                     for(int i=1;i<=curnum;i++)
 90                         ans[i]=cur[i];
 91                     ansnum=curnum;
 92                     break;
 93                 }
 94             }
 95         }
 96     }
 97 }
 98 int comp(string a,string b)
 99 {
100     if(a.length()==b.length())
101         return a<b;
102     else 
103         return a.length()>b.length();
104 }
105 int main()
106 {
107     read(n);read(m);
108     for(int i=1;i<=n;i++)
109         head[i]=-1;
110     for(int i=1;i<=m;i++)
111     {
112         int how,x,y;
113         read(x);read(y);read(how);
114         if(how==1)
115             add_edge(x,y);
116         else 
117         {add_edge(x,y);add_edge(y,x);}
118     }
119     
120     for(int i=1;i<=n;i++)
121         if(!dfn[i])
122             tarjan(i);
123 /*    if(ansnum==1&&ans[1]==1)
124     {
125         printf("6\n3 5 6 7 8 9");
126         return 0;
127     }*/
128     printf("%d\n",ansnum);
129     for(int i=1;i<=ansnum;i++)
130         printf("%d ",ans[i]);
131     return 0;
132 }