B:派

B:派

总时间限制: 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 <string.h>
#include <sstream>
#include <vector>
#include <algorithm>
#define PI 3.141592653589//没想到这个PI的精度居然很重要，少几位的话会WA
#define exp 0.0000001
using namespace std;

int n = 0, f = 0;
double CirS[10001] = {};

bool check(double s){
    int cnt = f + 1;
    for(int i = 0; i < n; ++i){
        int tmp = CirS[i]/s;
        cnt -= tmp;
        if(cnt <= 0)
            return 1;
    }
   return 0;//cnt > 0
}
int main(){
    cin >> n >> f;
    double s = 0, e = 0;
    double sum  = 0;
    double ans = 0;
    double r = 0;
    for(int i = 0; i < n; ++i){
        cin >> r;
        CirS[i] = r*r;
        sum += CirS[i];
    }
    e = sum/(f + 1);//加上我
    while(e - s > exp){
        double mid = s+(e-s)/2;
        if(check(mid)){
            ans = mid;
            s = mid;
        }
        else{
            e = mid;
        }
    }
    ans *= PI;
    printf("%.3f\n",ans);
    return 0;
}