openjudge8211_派

openjudge8211_派

时空限制    1000ms/64MB

描述

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

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

输入

第一行包含两个正整数N和F，1 ≤ N, F ≤ 10 000，表示派的数量和朋友的数量。
第二行包含N个1到10000之间的整数，表示每个派的半径。

输出

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

样例输入

3 3
4 3 3

样例输出

25.133

代码

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
const double PI = acos(-1), eps = 1e-7;
int n,f;
double v[10005];

bool can(double x){		//能否切成f+1块
	int sum=0;
	for (int i=1; i<=n; i++) sum+=int(v[i]/x);
	return sum>=f+1;
}

double solve(double low,double high){
	double mid;
	while (low+eps<high){
		mid=(low+high)/2;
		if (can(mid)) low=mid;
		else high=mid;
	}
	return high;
}

int main(){
	cin>>n>>f;
	for (int i=1,r; i<=n; i++){
		cin>>r;
		v[i]=PI*r*r;	//存派体积
	}
	printf("%.3lf\n",solve(0.0,400000000.0));
	return 0;
}