pku2728 最优比率生成树

额，这个东西我是最近才听说的，然后想了半天没想出来，看了各种大牛的题解发现原来是二分答案然后用prim，真的很神奇。

题目大意：

  给你一堆点，每个点有x，y，d三个参数，分别是xy坐标和高度？额，反正就这意思，两个点距离为欧几里得距离，然后高度距离也是欧几里得距离。求sum点距离/sum高度距离最小。

简略题解：

  其实就是把除法乘过来然后转化成判定性问题，用二分答案的方法然后prim解决。注意精度。

附程序：

program pku_2728; var a,b,w,d:array[0..1000] of double; c,dis,data:array[0..1000,0..1000] of double; f:array[0..1000] of boolean; n:longint; sum:double; {-------------------------------------} procedure init; var i,j:longint; begin sum:=0; for i:=1 to n do readln(a[i],b[i],w[i]); for i:=1 to n do sum:=sum+w[i]; for i:=1 to n do for j:=1 to n do begin dis[i,j]:=sqrt(sqr(a[i]-a[j])+sqr(b[i]-b[j])); data[i,j]:=abs(w[i]-w[j]); end; end; {-------------------------------------} function mst:double; var i,j,k:longint; min:double; begin for i:=1 to n do f[i]:=false; for i:=1 to n do d[i]:=c[1,i]; f[1]:=true;mst:=0; for i:=1 to n-1 do begin min:=-maxlongint; for j:=1 to n do if not f[j] and (d[j]>min) then begin min:=d[j];k:=j; end; f[k]:=true; mst:=mst+d[k]; for j:=1 to n do if not f[j] and (c[k,j]>d[j]) then d[j]:=c[k,j]; end; end; {-------------------------------------} procedure main; var l,r,mid,ans:double; i,j:longint; begin l:=0;r:=sum/(n-1); while r-l>0.00001 do begin mid:=(l+r)/2; for i:=1 to n do for j:=1 to n do c[i,j]:=mid*dis[i,j]-data[i,j]; ans:=mst; if ans>0.00001 then r:=mid else l:=mid end; writeln(l:0:3); end; {-------------------------------------} begin readln(n); while n<>0 do begin init; main; readln(n); end; end.