【NOIP2011模拟9.3】旅行 (Standard IO)

Description

　　Z 小镇是一个景色宜人的地方，吸引来自各地的观光客来此旅游观光。Z 小镇附近共有N 个景点（编号为1,2,3,…,N），这些景点被M 条道路连接着，所有道路都是双向的，两个景点之间可能有多条道路。也许是为了保护该地的旅游资源，Z 小镇有个奇怪的规定，就是对于一条给定的公路Ri，任何在该公路上行驶的车辆速度必须为Vi。
　　速度变化太快使得游客们很不舒服，因此从一个景点前往另一个景点的时候，大家都希望选择行使过程中最大速度和最小速度的比尽可能小的路线，也就是所谓最舒适的路线。

Input

　　第一行包含两个正整数，N 和M。
　　接下来的M 行每行包含三个正整数：x，y 和v。表示景点x 到景点y 之间有一条双向公路，车辆必须以速度v 在该公路上行驶。
最后一行包含两个正整数s，t，表示想知道从景点s 到景点t 最大最小速度比最小的路径。s 和t 不可能相同。

Output

　　如果景点s 到景点t 没有路径，输出“IMPOSSIBLE”。否则输出一个数，表示最小的速度比。如果需要，输出一个既约分数。

Sample Input

【样例输入1】

4 2

1 2 1

3 4 2

1 4



【样例输入2】

3 3

1 2 10

1 2 5

2 3 8

1 3



【样例输入3】

3 2

1 2 2

2 3 4

1 3

Sample Output

【样例输出1】

IMPOSSIBLE



【样例输出2】

5/4



【样例输出3】

2

Data Constraint

Hint

【数据范围】

　　对于100%的数据，1 < N≤500,1≤x,y≤N，0

题解：

用最小生成树（Kruskal），先用每条边的权值排个序，然后就枚举加边，只要s和t之间是联通的就可以更新答案退出。

如果s 到t 没有路径，就输出“IMPOSSIBLE”。如果它们的比是整数就直接输出，否则就要用gcd约分一下输出分数。

代码：

var
  n,m,i,s,t,xi,yi,flag:longint;
  fa:array[0..5000]of longint;
  x,y,v:array[0..5000]of longint;
procedure sort(l,r:longint);
var
  i,j,t,z:longint;
begin
  randomize;
  i:=l;j:=r;
  z:=v[random(r-l)+l];
  repeat
    while v[i]>z do inc(i);
    while v[j]<z do dec(j);
    if i<=j then
    begin
      t:=x[i];x[i]:=x[j];x[j]:=t;
      t:=y[i];y[i]:=y[j];y[j]:=t;
      t:=v[i];v[i]:=v[j];v[j]:=t;
      inc(i);dec(j);
    end;
  until i>j;
  if i<r then sort(i,r);
  if j>l then sort(l,j);
end;
function find(son:longint):longint;
begin
  if fa[son]=son then exit(son);
  fa[son]:=find(fa[son]);
  exit(fa[son]);
end;
procedure try(step:longint);
var
  i:longint;
begin
  for i:=1 to n do
    fa[i]:=i;
  for i:=step to m do
    begin
      fa[find(x[i])]:=find(y[i]);
      if find(s)=find(t) then
      begin
        inc(flag);
        if flag=1 then
        begin
          xi:=v[step];
          yi:=v[i];
        end;
        if (v[step]/v[i])<(xi/yi) then
        begin
          xi:=v[step];
          yi:=v[i];
        end;
      end;
    end;
end;
function gcd(a,b:longint):longint;
begin
  if b=0 then exit(a);
  exit(gcd(b,a mod b));
end;
procedure print;
begin
  if xi mod yi=0 then
  begin
    writeln(xi div yi);
    exit;
  end;
  writeln(xi div gcd(xi,yi),'/',yi div gcd(xi,yi));
end;
begin
  randomize;
  readln(n,m);
  for i:=1 to m do
    readln(x[i],y[i],v[i]);
  readln(s,t);
  sort(1,m);
  for i:=1 to m do
    try(i);
  if flag=0 then writeln('IMPOSSIBLE')
  else print;
end.