题意:
有n块圆形馅饼,每块的馅饼的半径为r,可以任意切分,剩下的都扔掉。
问分给n+1个人,使得每个人获得馅饼的体积最大
解析:
先让每个数变成r*r,然后二分答案x
O(n)计算当前答案为x,可以分得多少块 a[i]/x(1<=i<=n) 假设可以获得cnt块
如果cnt>=n+1 说明x过小,使得cnt变大了,所以我们要调整下界
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define PI 3.141592653589793238462643383279502884197169399
using namespace std;
const int N=1e5+1000;
const double eps=1e-8;
double a[N];
int t,n,m;
bool check(double x)
{
int cnt=0;
for(int i=1;i<=n;i++)
{
cnt=cnt+(int)(a[i]/x);
}
if(cnt>=m+1) return false;
return true;
}
int main()
{
scanf("%d",&t);
while(t--)
{
scanf("%d %d",&n,&m);
double l=0,r=0;
for(int i=1;i<=n;i++) scanf("%lf",&a[i]),a[i]=a[i]*a[i],r=max(r,a[i]);
while (r - l >eps)
{
double mid = (l + r) / 2;
if (check(mid)) r = mid;
else l = mid;
}
printf("%.4f\n",l*PI);
}
}