浮点数二分答案 HDU1969

#include <iostream>
#include <cstdio>
#include <queue>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

double PI=acos(-1.0);
double arr[10010];
int n,m;

bool test(double mid)
{
int ans=0;
for(int i=0;i<n;i++)
{
ans+=(int)(arr[i]/mid);
}
if(ans>=m+1)
return true;
else
return false;
}

int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
double sum=0;
double rr;
for(int i=0;i<n;i++)
{
scanf("%lf",&rr);
arr[i]=PI*rr*rr;
sum+=arr[i];
}
double r=sum/(m+1);
double mid=0.0;
double l=0.0;
while(r-l>1e-6)
{
mid=(r+l)/2;
if(test(mid))
{
l=mid;
}
else
{
r=mid;
}
}
printf("%.4lf\n",mid);
}
return 0;
}