《二分》hdoj acm 4.1.3

#include<stdio.h>
#include<math.h>
double pi=acos(-1.0);
int main()
{
int t,n,f,i,ri,num;
double m[10005];
double maxn,r,l,mid;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&f);
maxn=0;
f++;
for(i=1;i<=n;i++)
{
scanf("%d",&ri);
m[i]=ri*ri*pi;
if(m[i]>maxn)
maxn=m[i];
}
l=0;
r=maxn;
while(l+0.00001<r)
{
mid=(r+l)/2;
num=0;
for(i=1;i<=n;i++)
{
num+=(int)floor(m[i]/mid);
}
if(num<f)
r=mid;
else
l=mid;
}
printf("%.4lf\n",l);
}
return 0;
}