还是用到二分法,和上一个组装电脑的非常像,去二分每个人得到的那一块派x,
看是否满足》=f+1块,满足继续二分,不满足,继续二分。
这里用到的是浮点数的二分方法:
double L=0,R=maxa;
while(R-L > 1e-5){
double M=L+(R-L)/2;
if(ok(M)) L=M; else R=M;
}
#include<bits/stdc++.h>
using namespace std;
const double PI = acos(-1.0);
const int maxn = 10000 + 5;
int n,f;
double A[maxn];
bool ok(double area){
int sum=0;
for(int i=0;i<n;i++) sum+=floor(A[i]/area);
return sum >= f+1;
}
int main()
{
int T;
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&f);
double maxa = -1;
for(int i=0;i<n;i++){
int r;
scanf("%d",&r);
A[i]=PI*r*r;
maxa=max(maxa,A[i]);
}
double L=0,R=maxa;
while(R-L > 1e-5){
double M=L+(R-L)/2;
if(ok(M)) L=M; else R=M;
}
printf("%.5lf\n",L);
}
return 0;
}