001:派

我的生日要到了!根据习俗,我需要将一些派分给大家。我有N个不同口味、不同大小的派。有F个朋友会来参加我的派对,每个人会拿到一块派(必须一个派的一块,不能由几个派

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

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

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

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

3 3
4 3 3

样例输出

25.133

对最大派的体积二分就可以了。
应当注意一点,二分的下界应该是0,而不是最小派的体积。因为最小派的体积可能比答案大,就像样例一样。

#include
#include
#define PI acos(-1)
using namespace std;
int n,f;
double pie[10005];
double EPS = 1e-6;
inline double Capacity(int r)
{return 1.0*r*r*PI;
}
//判断二分的体积是否可以将所有蛋糕分成f+1块
bool Check(double v)
{int cnt = 0;for(int i = 0; i<n; ++i){cnt += (int)(pie[i]/v);}if(cnt>=f+1)return true;return false;
}
//二分蛋糕的最大体积
double BinarySearch(double l,double r)
{double mid=0;while(r-l>EPS){mid = l + (r-l)/2.0;if(Check(mid))l = mid;elser = mid;}return mid;
}
int main()
{int t;scanf("%d%d",&n,&f);double r = 0;for(int i = 0; i<n; ++i){scanf("%d",&t);pie[i] = Capacity(t);if(pie[i]>r)r = pie[i];}printf("%.3f\n",BinarySearch(0.0,r));return 0;
}