MOOC程序设计与算法二第4周作业2 派 二分

#include<bits/stdc++.h>
#define PI acos(-1)
#define eps 1e-5 
using namespace std;
int N,F;
int a[10010];
double size(int r);
int good(double d);
int main()
{
	double L,R;//面积 
	double D,result,sumsize;
	while(scanf("%d%d",&N,&F)!=EOF)
	{
		F=F+1;//加上自己 
		sumsize=0;
		for(int i=0;i<N;i++)
		{
			scanf("%d",&a[i]);
			sumsize+=size(a[i]);
		}
		sort(a,a+N);
		L=0;
		R=sumsize/F;
		D=(L+R)/2;
		while(R-L>eps)
		{
			if(good(D))
			{
				L=D;
				result=D;
			}
			else
			{
				R=D;
			}
			D=(L+R)/2;
		}
		printf("%.3lf",result); 
	}
}
int good(double d)
{
	int people=0;
	for(int i=0;i<N;i++)
	{
		people+=(int)(size(a[i])/d);
	}
	if(people>=F)
		return 1;
	return 0;
}
double size(int r)
{
	return r*r*PI;
}