旅行

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

Hint

【数据范围】
1<N<=500
1<=x,y<=N,0<v<30000,x≠y
0<M<=5000

分析
按速度由大到小排序,把依次加入边。

程序:

type node=record
x,y,v:longint;
end;
var
e:array[1..6000]of node;
r:array[1..6000]of longint;
n,m,s,t,minx,miny:longint;
procedure init;
var
i:longint;
begin
    readln(n,m);
    for i:=1 to m do
    readln(e[i].x,e[i].y,e[i].v);
    readln(s,t);

end;
procedure swap(x,y:longint);
var
t:node;
begin
    t:=e[x];e[x]:=e[y];e[y]:=t;
end;
procedure sort(s,t:longint);
var
i,j:longint;
x:node;
begin
    i:=s;j:=t;x:=e[(s+t)shr 1];
    repeat
         while (i<=j)and(x.v>e[i].v) do inc(i);
         while (i<=j)and(x.v<e[j].v) do dec(j);
         if i<=j then
         begin
             swap(i,j);
             inc(i);dec(j);
         end;
    until i>j;
    if s<j then sort(s,j);
    if i<t then sort(i,t);
end;
function find(x:longint):longint;
begin
    if r[x]<0 then find:=x else
    begin
        r[x]:=find(r[x]);
        find:=r[x];
    end;
end;
procedure add(x,y:longint);
begin
    x:=find(x);
    y:=find(y);
    if x<>y then
    begin
        r[x]:=r[x]+r[y];
        r[y]:=x;
    end;
end;
function gcd(a,b:longint):longint;
begin
    if b=0 then gcd:=a else gcd:=gcd(b,a mod b);
end;
procedure print;
var
d:longint;
begin
    if minx*e[1].v>miny*e[m].v then writeln('IMPOSSIBLE') else
    begin
        d:=gcd(minx,miny);
        minx:=minx div d;
        miny:=miny div d;
        if minx mod miny=0 then writeln(minx div miny) else writeln(minx,'/',miny);
    end;
end;
procedure main;
var
i,j,k:longint;
begin
    randomize;
    sort(1,m);
    minx:=e[m].v+1;
    miny:=e[1].v;
    for i:=1 to m do
    begin
        fillchar(r,sizeof(r),$ff);
        j:=i-1;
        while (j<m)and(find(s)<>find(t)) do
        begin
            inc(j);
            add(e[j].x,e[j].y);
        end;
        if find(s)<>find(t) then break;
        if e[j].v*miny<e[i].v*minx then
        begin
            minx:=e[j].v;
            miny:=e[i].v;
        end;
    end;
end;
begin
    init;
    main;
    print;
end.
posted @ 2017-12-14 16:48  银叶草  阅读(103)  评论(0编辑  收藏  举报
Live2D