C++二分学习派

总时间限制: 

1000ms

内存限制: 

65536kB

描述

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

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

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

输入

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

输出

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

样例输入

3 3
4 3 3

  题目的意思是总共有N个派，N个派的半径是不同的，分给K个人，每个人获得的体积是一样的，要使得获得的体积最大，因为 

   每个派的半径是不同的，依然用二分的想法，直接枚举这个最大体积，显然最大体积的极值必然没最大的派的体积，因为最大派都被分完了，就不可能再分了。最小体积为0。开始二分，如果每个派的体积/（每个人得到的最大体积）》总人数，说明每个人得到的最大体积还不够大，否则说明体积过大。                                                         

#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
const double PI=acos(-1.0); 
const double eps = 1e-6;
int N;
int F;
double r[10001];
bool Valid(double V)
{
	int flag=0;
	for(int i=0;i<N;i++)
	{
		flag+=r[i]*r[i]/V;
	}
	if(flag>=F)
	{
		return true;
	}
	return false;
}
int main()
{
	cin>>N>>F;
	F++;
	double R=0;
	for(int i=0;i<N;i++)
	{
		cin>>r[i];
		R=max(R,r[i]);
	}
	double L=0;
	R=R*R;
	double last=0;
	while(R-L>eps)
	{
		double mid=L+(R-L)/2;
		if(Valid(mid))
		{
			last=mid;
			L=mid;
		}
		else
			R=mid;
	}
	printf("%.3f\n",last*PI);
	return 0;
}