POJ 3122 Pie 二分查找 +高精度

 点击打开链接

题意：作者要开一个生日party，他现在拥有n块高度都为1的圆柱形奶酪，

 已知每块奶酪的底面半径为r不等，作者邀请了f个朋友参加了他的party，
 他要把这些奶酪平均分给所有的朋友和他自己（f+1人），
 每个人分得奶酪的体积必须相等（这个值是确定的），形状就没有要求。
 现在要你求出所有人都能够得到的最大块奶酪的体积是多少？
 
一个人只能从其中一个蛋糕中 取一部分
蛋糕可以有剩余
 思路:
     二分枚举出蛋糕的面积
     注意
     1 精度问题   ///高精度加在较小的数上
     2 mid =high 或mid= low 不能 +1
     3 %f

     4 acos(-1.0)  为

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<queue>
#include<set>
#define PI acos(-1)
using namespace std;
double a[20000];
int n,T,m;
double solve(double low,double high)
{
    double sum,zmax=0.0;
    int num=0;
  //  while(high-low>)
    while(low+0.000001<high)         ///高精度加在较小的数上
    {
        double mid=(low+high)/2;
        num=0;
        for(int i=1;i<=n;i++)
        {
            num+=(a[i]/mid);
        }
        if(num<m)
            high=mid;
        else
        {
            low=mid;
            zmax=max(zmax,mid);
        }
    }
    return zmax;
}
int main()
{
        scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        m=m+1;
        double Max=0.0;
        for(int i=1; i<=n; i++)
        {
            scanf("%lf",&a[i]);
            a[i]=1.0*a[i]*a[i];
            Max=max(Max,a[i]);
        }
        printf("%.4f\n",solve(0,Max)*PI);
    }
    return 0;
}