我的生日要到了!根据习俗,我需要将一些派分给大家。我有N个不同口味、不同大小的派。有F个朋友会来参加我的派对,每个人会拿到一块派(必须一个派的一块,不能由几个派的小块拼成;可以是一整个派)。
我的朋友们都特别小气,如果有人拿到更大的一块,就会开始抱怨。因此所有人拿到的派是同样大小的(但不需要是同样形状的),虽然这样有些派会被浪费,但总比搞砸整个派对好。当然,我也要给自己留一块,而这一块也要和其他人的同样大小。
请问我们每个人拿到的派最大是多少?每个派都是一个高为1,半径不等的圆柱体。
输入
第一行包含两个正整数N和F,1 ≤ N, F ≤ 10 000,表示派的数量和朋友的数量。
第二行包含N个1到10000之间的整数,表示每个派的半径。
输出
输出每个人能得到的最大的派的体积,精确到小数点后三位。
样例输入
3 3
4 3 3
样例输出
25.133
对最大派的体积二分就可以了。
应当注意一点,二分的下界应该是0,而不是最小派的体积。因为最小派的体积可能比答案大,就像样例一样。
#include<cstdio>
#include<cmath>
#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;
else
r = 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;
}