【NOIP2016练习】T1 挖金矿（二分答案）

题意：

思路：二分答案A

合法的答案 sigma(s[i][xi])/sigma(xi)>=a i<=m

sigma(s[i][xi]-a*xi)>=0

对于每个i找到xi使s[i][xi]-a*xi最大累加验证答案即可

 var a,s:array[..]of int64;
     n,m,i,j:longint;
     l,r,mid,last:real;

 function isok(a:real):boolean;
 var i,j:longint;
     t,tmp:real;
 begin
  tmp:=;
  for i:= to n do
  begin
   t:=s[(i-)*m+]-a;
   for j:= to m do
    if s[(i-)*m+j]-a*j>t then t:=s[(i-)*m+j]-a*j;
   tmp:=tmp+t;
  end;
  if tmp>= then exit(true)
   else exit(false);
 end;

 begin
  assign(input,'gold.in'); reset(input);
  assign(output,'gold.out'); rewrite(output);
  readln(n,m);
  for i:= to n*m do read(a[i]);
  for i:= to n do
  begin
   s[(i-)*m+]:=a[(i-)*m+];
   for j:= to m do s[(i-)*m+j]:=s[(i-)*m+j-]+a[(i-)*m+j];
  end;
  l:=; r:=; last:=;
  while r-l>1e-5 do
  begin
   mid:=(l+r)/;
   if isok(mid) then begin last:=mid; l:=mid; end
    else r:=mid;
  end;
  writeln(last::);
  close(input);
  close(output);
 end.