OpenJudge 008:派

008:派

总时间限制: 1000ms 内存限制: 65536kB

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

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

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

输入
第一行包含两个正整数N和F，1 ≤ N, F ≤ 10 000，表示派的数量和朋友的数量。

第二行包含N个1到10000之间的整数，表示每个派的半径。

输出
输出每个人能得到的最大的派的体积，精确到小数点后三位。

样例输入

3 3
4 3 3

样例输出

25.133

思想：
本题同样是二分法， right就是将输入数据排序后的最大值（即最大半径），left设为0；
然后就在于验证mid所对应的体积是否能够分够F+1个人，对于算能够分的人数可以用——每块蛋糕的体积除以mid对应体积然后int强制转化后累加——用得到的结果和F+1比较；
如果小于不够分则说明mid太大，设置right=mid-MIN（很小的数），否则就是够分mid较小，记录下mid对应的体积，然后设置left=mid+MIN。

本题的槽点是小数点后的几位总是有那么一点点的差别，我也是在不断的修改PIE的值的精度以及MIN的精度，才最终AC的。

代码：

#include<iostream>
#include<algorithm>
#include<iomanip>
#include<cmath>
using namespace std;

//精度很重要！！！
const double PI=3.141592653589793;
const double MIN=1e-9;

//计算体积
double get_V(double r){
    return PI * r * r;
}

//计算mid对应体积能分多少人
int cut_count(double V, const int* ar, int N){
    int cnt=0;
    for(int i=0;i<N;++i)
        cnt+=int(get_V(ar[i])/V);
    return cnt;
}

//二分
double process(const int N,const int F){
    int *ar=new int[N];
    for(int i=0;i<N;++i)
        cin >> ar[i];

    sort(ar,ar+N);

    double maxv=0;
    double left=0.0,right=get_V(ar[N-1]);
    while(left<=right){
        double mid=left+(right-left)/2;
        if(cut_count(get_V(mid), ar, N) < (F + 1))
            right=mid-MIN;
        else {
            maxv=get_V(mid);
            left = mid + MIN;
        }
    }

    delete[] ar;
    return maxv;
}

int main(){
    int N,F;
    //freopen("in-4.txt","r",stdin);
    cin>>N>>F;
    cout<<setiosflags(ios::fixed)
        <<setprecision(3)<<process(N,F);
    return 0;
}