浮点数二分答案 HDU1969

#include <iostream>
#include <cstdio>
#include <queue>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

double PI=acos(-1.0);
double arr[10010];
int n,m;

bool test(double mid)
{
    int ans=0;
    for(int i=0;i<n;i++)
    {
        ans+=(int)(arr[i]/mid);
    }
    if(ans>=m+1)
        return true;
    else
        return false;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        double sum=0;
        double rr;
        for(int i=0;i<n;i++)
        {
            scanf("%lf",&rr);
            arr[i]=PI*rr*rr;
            sum+=arr[i];
        }
        double r=sum/(m+1);
        double mid=0.0;
        double l=0.0;
        while(r-l>1e-6)
        {
            mid=(r+l)/2;
            if(test(mid))
            {
                l=mid;
            }
            else
            {
                r=mid;
            }
        }
        printf("%.4lf
",mid);
    }
    return 0;
}