【openjudge】派 二分查找

05:派

总时间限制: 1000ms 内存限制: 65536kB
描述
我的生日要到了！根据习俗，我需要将一些派分给大家。我有N个不同口味、不同大小的派。有F个朋友会来参加我的派对，每个人会拿到一块派（必须一个派的一块，不能由几个派的小块拼成；可以是一整个派）。

我的朋友们都特别小气，如果有人拿到更大的一块，就会开始抱怨。因此所有人拿到的派是同样大小的（但不需要是同样形状的），虽然这样有些派会被浪费，但总比搞砸整个派对好。当然，我也要给自己留一块，而这一块也要和其他人的同样大小。

请问我们每个人拿到的派最大是多少？每个派都是一个高为1，半径不等的圆柱体。

输入
第一行包含两个正整数N和F，1 ≤ N, F ≤ 10 000，表示派的数量和朋友的数量。
第二行包含N个1到10000之间的整数，表示每个派的半径。
输出
输出每个人能得到的最大的派的体积，精确到小数点后三位。
样例输入
3 3
4 3 3
样例输出
25.133

思路：枚举每个人可能分得的蛋糕体积。注意卡精度。

代码：

#include<iostream>
#include<cstdio>
using namespace std;
double l,r,mid;
long long cnt,i,n,f;
double fg[10001];
double x;
long long love(double mid)
{
    long long sum=0,i;
    for (i=1;i<=n;i++)
        sum+=fg[i]/(mid*1.0);
    return sum;
}
int main()
{
    scanf("%d%d",&n,&f);
    f++;
    for (i=1;i<=n;i++)
    {
        scanf("%lf",&x);
        fg[i]=x*x*3.141592653589;
        r=max(r,fg[i]);
    }
    l=0;
    while (l+0.00001<r)
    {
        mid=(l+r)/2.0;
        cnt=love(mid);
        if (cnt<f) r=mid-0.00001;
        else l=mid;
    }
    printf("%.3lf",l);
}