二分法 Pie HDU - 1969

#include< cstdio>
#include< cmath>
#include< iostream>
#include< algorithm>
using namespace std;
#define eps 1e-7
const double PI=acos(-1.0); //精确计算π的值，减小误差
int n,f;
double cake[10005];
bool check(double h)
{
int sum=0;
for(int i=0;i<n;i++)
{
sum=sum+cake[i]/h; //可以分给几个人
}
if(sum>=(f+1))
{
return true;
}
else
{
return false;
}
}
int main ()
{
int T;
scanf("%d",&T);
while(T–)
{
scanf("%d%d",&n,&f);
for(int i=0;i<n;i++)
{
scanf("%lf",&cake[i]);
cake[i]=cake[i]*cake[i]PI1.0; //每个派的体积
}
sort(cake,cake+n); //找出体积最大个派，目的为下面取有极限做准备
double left=0, right=cake[n-1]; //利用二分法（浮点型）
while(left + eps <right)
{
double mid=(left+right)/2.0;
if(check(mid))
{
left=mid;
}
else
{
right=mid;
}
}
printf("%.4lf\n",left);
}
}