3402: [Usaco2009 Open]Hide and Seek 捉迷藏

Time Limit: 3 Sec Memory Limit: 128 MB

Submit: 78 Solved: 64

Description

贝茜在和约翰玩一个“捉迷藏”的游戏．

她正要找出所有适合她躲藏的安全牛棚．一共有N(2≤N≤20000)个牛棚，被编为1到N号．她知道约翰（捉牛者）从牛棚1出发．所有的牛棚由M(1≤M≤50000)条双向路连接，每条双向路连接两个不同的牛棚．所有的牛棚都是相通的．贝茜认为同牛棚1距离最远的的牛棚是安全的．两个牛棚间的距离是指，从一个牛棚到另一个牛棚最少需要通过的道路数量．请帮贝茜找出所有的安全牛棚．

Input

第1行输入两个整数N和M，之后M行每行输入两个整数，表示一条路的两个端点．

Output

仅一行，输出三个整数．第1个表示安全牛棚（如果有多个，输出编号最小的）；第2个表示牛棚1和安全牛棚的距离；第3个表示有多少个安全的牛棚．

Sample Input

6 7 3 6 4 3 3 2 1 3 1 2 2 4 5 2

Sample Output

4 2 3

HINT

Source

Silver

题解：不用多说，明显的单源最短路，但是出现了一个很神奇的小插曲——一开始用dijkstra写，结果TLE得很惨；于是换成spfa，然后192ms就AC了？！？！

于是，本人打算就此展开下一次实验研究，同时将可能考虑引入用堆优化的dijkstra，以及在进行全源最短路时引入Floyd算法进行对比，敬请期待

 1 /**************************************************************
 2     Problem: 3402
 3     User: HansBug
 4     Language: Pascal
 5     Result: Accepted
 6     Time:192 ms
 7     Memory:13904 kb
 8 ****************************************************************/
 9  
10 type
11     point=^node;
12     node=record
13                g,w:longint;
14                next:point;
15     end;
16 var
17    i,j,k,l,m,n,f,r:longint;
18    a:array[0..100000] of point;
19    b,c,d:array[0..1000000] of longint;
20    p:point;
21 procedure add(x,y,z:longint);
22           var p:point;
23           begin
24                new(p);p^.g:=y;p^.w:=z;
25                p^.next:=a[x];a[x]:=p;
26           end;
27 function min(x,y:longint):longint;
28          begin
29               if x<y then min:=x else min:=y;
30          end;
31 begin
32      readln(n,m);
33      for i:=1 to n do a[i]:=nil;
34      for i:=1 to m do
35          begin
36               readln(j,k);
37               add(j,k,1);add(k,j,1);
38          end;
39      f:=1;r:=2;
40      fillchar(c,sizeof(c),0);
41      b[1]:=1;c[1]:=1;
42  
43      while f<r do
44            begin
45                 p:=a[b[f]];
46                 while p<>nil do
47                       begin
48                            if (c[p^.g]=0) or (c[p^.g]>(c[b[f]]+p^.w)) then
49                               begin
50                                    b[r]:=p^.g;
51                                    c[p^.g]:=c[b[f]]+p^.w;
52                                    inc(r);
53                               end;
54                            p:=p^.next;
55                       end;
56                 inc(f);
57            end;
58      j:=0;k:=0;l:=0;
59      for i:=1 to n do
60          if c[i]>j then
61             begin
62                  j:=c[i];
63                  k:=1;
64                  l:=i;
65             end
66          else if c[i]=j then inc(k);
67      writeln(l,' ',j-1,' ',k);
68 end.