Uva 派 (Pie,NWERC 2006,LA 3635)

依然是一道二分查找

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

const double PI=acos(-1.0);
int N,F;
double r[10001];

bool ok(double area)
{
    int sum=0;
    for(int i=0;i<N;i++)
        sum+=floor(r[i]/area);
    return sum>=F+1;
}

int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        cin>>N>>F;
        double maxn=-1;
        for(int i=0;i<N;i++)
        {
            int a;
            cin>>a;
            r[i]=PI*a*a;
            maxn=max(maxn,r[i]);
        }
        double L=0,R=maxn;
        while(R-L>1e-4)
        {
            double M=(L+R)/2;
            if(ok(M)) L=M;else R=M;
        }
        printf("%.4lf\n",L);
    }
    return 0;
}