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;
}