UVALIVE 3635 Pie 二分

有m个派， 要分给n + 1个人， 要求每个人拿到的大小相同， 并且每个人的派必须是一整块， 不能说分别从两个派切出一块凑成。 现在给出m个派的半径(派均为圆柱，高为1），输出每人可以拿到的最大派的体积。

注意double精度等问题直接二分就OK

#include <map>
#include <set>
#include <list>
#include <cmath>
#include<cctype>
#include <ctime>
#include <deque>
#include <stack>
#include <queue>
#include <cstdio>
#include <string>
#include <vector>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
#define PI 3.1415926535897932626
using namespace std;
int gcd(int a, int b)
{
	return a % b == 0 ? b : gcd(b, a % b);
}
#define MAXN 10005
double src[MAXN],limit,Maxv;
int N,F;
int cmp(const int a,const int b) {return a>b;}
bool judge(double v)
{
   int sum=0;
   for (int i=0;i<N;i++)
   {
       double tmp=PI*src[i]*src[i];
       sum+=floor(tmp/v);
       if (sum>=F+1) return true;
   }
   return false;
}
int main()
{
    //freopen("sample.txt","r",stdin);
    int T;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d",&N,&F);
        limit=0.0;
        for (int i=0;i<N;i++)
        {
            scanf("%lf",&src[i]);
            if (src[i]>limit) limit=src[i];
        }
        sort(src,src+N);
        Maxv=PI*limit*limit;
        //printf("%.4lf\n",limit);
        double left=0.0,right=Maxv,mid;
        while (right-left>=1e-5)
        {
            mid=left+(right-left)/2;
            if (judge(mid)) left=mid;
            else right=mid;
        }
       // printf("%.4lf %.4lf %.4lf\n",left,mid,right);
       printf("%.4lf\n",mid);
    }
    return 0;
}