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=;
   eps=1e-8;
 var
   u,v,w:array[..]of longint;
   father:array[..]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)>>];
    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<> 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:= to m do readln(u[i],v[i],w[i]);
    readln(s,t);
    sort(,m);
    min:=inf;
    for i:= to m do
       begin
          for j:= 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= then write(q) else write(q,'/',p); end;
    close(input);close(output);
 end.