P1726 上白澤慧音(0分)

題目描述

在幻想鄉，上白澤慧音是以知識淵博聞名的老師。春雪異變導致人間之里的很多道路都被大雪堵塞，使有的學生不能順利地到達慧音所在的村莊。因此慧音決定換一個能夠聚集最多人數的村莊作為新的教學地點。人間之里由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

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<cmath>
  5 #include<queue>
  6 #include<algorithm>
  7 using namespace std;
  8 const int MAXN=5001;
  9 void read(int &n)
 10 {
 11     char c='+';int x=0;bool flag=0;
 12     while(c<'0'||c>'9')
 13     {c=getchar();if(c=='-')flag=1;}
 14     while(c>='0'&&c<='9')
 15     {x=x*10+c-48;c=getchar();}
 16     flag==1?n=-x:n=x;
 17 }
 18 struct node
 19 {
 20     int u,v,w,nxt;
 21 }edge[MAXN*4];
 22 int head[MAXN];
 23 int num=1;
 24 void add_edge(int x,int y,int z)
 25 {
 26     edge[num].u=x;
 27     edge[num].v=y;
 28     edge[num].w=z;
 29     edge[num].nxt=head[x];
 30     head[x]=num++;
 31 }
 32 bool have[MAXN][MAXN];
 33 int pre[MAXN];
 34 int nex[MAXN];
 35 int out[MAXN];
 36 void cler()
 37 {
 38     memset(have,0,sizeof(have));
 39     memset(pre,-1,sizeof(pre));
 40     memset(nex,-1,sizeof(nex));
 41 }
 42 int ans=0;
 43 int ansbg=0;
 44 int n,m;
 45 void calc_ans(int bg,int ed)
 46 {
 47     int num=0;
 48     int hh=ed;
 49     while(pre[ed]!=-1)
 50     {
 51         num++;
 52         ed=pre[ed];
 53     }
 54     num+=1;
 55     if(num>ans)
 56     {
 57         ans=num;
 58         ansbg=bg;
 59         for(int i=0;i<=n;i++)
 60             out[i]=nex[i];
 61     }
 62 }
 63 void dfs(int bg,int now)
 64 {
 65     for(int i=head[now];i!=-1;i=edge[i].nxt)
 66     {
 67         if(!have[now][edge[i].v])
 68         {
 69             if(edge[i].v==bg)
 70             {
 71                 calc_ans(bg,now);
 72                 return ;
 73             }
 74             have[now][edge[i].v]=1;
 75             have[edge[i].v][now]=1;
 76             pre[edge[i].v]=now;
 77             nex[now]=edge[i].v;
 78             dfs(bg,edge[i].v);
 79             have[now][edge[i].v]=0;
 80             have[edge[i].v][now]=0;
 81             pre[edge[i].v]=-1;
 82             nex[now]=-1;
 83         }
 84     }
 85 }
 86 int main()
 87 {
 88     
 89     read(n);read(m);
 90     for(int i=1;i<=n;i++)
 91         head[i]=-1;
 92     for(int i=1;i<=m;i++)
 93     {
 94         int x,y,how;
 95         read(x);read(y);read(how);
 96         if(how==1) add_edge(x,y,1);
 97         else {add_edge(x,y,1);add_edge(y,x,1);}
 98     }
 99     for(int i=1;i<=n;i++)
100     {
101         cler();
102         dfs(i,i);
103     }
104     printf("%d\n",ans);
105     if(ansbg)
106         printf("%d ",ansbg);
107     while(out[ansbg]!=-1)
108     {
109         ansbg=out[ansbg];
110         printf("%d ",ansbg);
111     }
112     return 0;
113 }