HOJ 2615 Pie ---------------二分查找的运用

//题意：有f+1个人分n块披萨，每个人要求分得的面积一样，
//      且披萨只能被切开而不能重新组合，求每个人能分到的最大面积v。
//思路：对于每个确定的v，可以计算出最多能满足的人数p。
//      因此得到一个单调递减的函数关系，并且v的范围也可以确定为0~max(size(i)),i=1...n。（图中）
//      一般要求什么酒对什么进行二分，对某一个面积进行判断，如果满足条件，那么找比这个面积大的满足条件的
//调试：编译的时候出现 undefined reference to `WinMain@16'
//     由于把main写成 mian啦
//     这个题貌似很恶心 Pi  还是用反三角函数比较好
//      人数是输入的F+1个儿不是F个他本人也要分
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#define  Pi acos(-1)
#define maxlen 10010
using namespace std;
long long high,mid,low,ans,s[maxlen];
int T,N,i,F;
double r;
bool judge(long long mid)
{
int i;
bool  flag;
long long p = 0;
for (i = 0; i < N; i++)
p += s[i] / mid;
if(p>=F+1)flag=true;//注意这里是f+1
else flag=false;
return flag;
}//判断分成这个面积可不可以
int main()
{
cin >> T;
while(T--)
{
low=0,high=0,ans=0;
cin >> N >> F;
memset(s,0,sizeof(s));
for(i=0; i<N; i++)
{
cin >> r;
s[i]=r*Pi*r*1000000;//先乘以一个大的树，再除回来，怕精度问题
high+=s[i];
}
while(low<=high)
{
mid=(low+high)/2;
if (judge(mid))
{
low = mid + 1;
ans = mid;
}
else
high = mid - 1;
}
printf("%.4lf\n", (double)ans/1000000);
}
return 0;
}