bzoj1050[HAOI2006]旅行comf（枚举+贪心+并查集）

Description

　　给你一个无向图，N(N<=500)个顶点, M(M<=5000)条边，每条边有一个权值Vi(Vi<30000)。给你两个顶点S和T，求一条路径，使得路径上最大边和最小边的比值最小。如果S和T之间没有路径，输出”IMPOSSIBLE”，否则输出这个比值，如果需要，表示成一个既约分数。 备注： 两个顶点之间可能有多条路径。

Input

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

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

 

 

 

题目中要求的最小比值涉及两个值，因此枚举其中一个，再贪心求另一个。比如先枚举最小的，然后将比它大的边按从小到大的顺序加入并查集，直到s和t连通，最后加入的边就是最好的最大边。

program rrr(input,output);
const
  inf=123456789;
  eps=1e-8;
var
  u,v,w:array[0..5050]of longint;
  father:array[0..505]of longint;
  n,m,s,t,i,j,c,d,p,q:longint;
  min:double;
procedure sort(q,h:longint);
var
  i,j,x,t:longint;
begin
   i:=q;j:=h;x:=w[(i+j)>>1];
   repeat
     while w[i]<x do inc(i);
     while x<w[j] do dec(j);
     if i<=j then
        begin
           t:=u[i];u[i]:=u[j];u[j]:=t;
           t:=v[i];v[i]:=v[j];v[j]:=t;
           t:=w[i];w[i]:=w[j];w[j]:=t;
           inc(i);dec(j);
        end;
   until i>j;
   if j>q then sort(q,j);
   if i<h then sort(i,h);
end;
function find(k:longint):longint;
begin
   if father[k]=k then exit(k);
   father[k]:=find(father[k]);
   exit(father[k]);
end;
function gcd(a,b:longint):longint;
var
  r:longint;
begin
   r:=a mod b;
   while r<>0 do begin a:=b;b:=r;r:=a mod b; end;
   exit(b);
end;
begin
   assign(input,'r.in');assign(output,'r.out');reset(input);rewrite(output);
   readln(n,m);
   for i:=1 to m do readln(u[i],v[i],w[i]);
   readln(s,t);
   sort(1,m);
   min:=inf;
   for i:=1 to m do
      begin
         for j:=1 to n do father[j]:=j;
         for j:=i to m do
            begin
               c:=find(u[j]);d:=find(v[j]);
               if c<>d then father[c]:=d;
               c:=find(s);d:=find(t);
               if c=d then break;
            end;
         if c=d then if w[j]/w[i]<min then begin min:=w[j]/w[i];p:=w[i];q:=w[j]; end;
      end;
   if abs(min-inf)<eps then write('IMPOSSIBLE')
   else begin c:=gcd(p,q);p:=p div c;q:=q div c;if p=1 then write(q) else write(q,'/',p); end;
   close(input);close(output);
end.