LA3635-二分法

14 篇文章 0 订阅
13 篇文章 0 订阅

解释都在大白书上,我就不再造车轮了,上代码:

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

const int maxn=1e4+5;
const double pi=acos(-1);
double a[maxn];
int n,f;

bool ok(double s)
{
    int sum=0;
    for(int i=0;i<n;i++)sum+=floor(a[i]/s);
    return sum>=f+1;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&f);
        double r;
        double m=-1;
        for(int i=0;i<n;i++){
            scanf("%lf",&r);
            a[i]=pi*r*r;
            m=max(m,a[i]);
        }
        double R=0,L=m;
        while(L-R>1e-5){
            double mid=(L+R)/2;
            if(ok(mid))R=mid;
            else L=mid;
        }
        printf("%.4f\n",L);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值
>