【noiOJ】p8211 (PS:二分浮点数的精度问题)

时间:2023-03-09 01:56:08
【noiOJ】p8211 (PS:二分浮点数的精度问题)

05:派

总时间限制: 
1000ms

内存限制: 
65536kB
描述

我的生日要到了!根据习俗,我需要将一些派分给大家。我有N个不同口味、不同大小的派。有F个朋友会来参加我的派对,每个人会拿到一块派(必须一个派的一块,不能由几个派的小块拼成;可以是一整个派)。

我的朋友们都特别小气,如果有人拿到更大的一块,就会开始抱怨。因此所有人拿到的派是同样大小的(但不需要是同样形状的),虽然这样有些派会被浪费,但总比搞砸整个派对好。当然,我也要给自己留一块,而这一块也要和其他人的同样大小。

请问我们每个人拿到的派最大是多少?每个派都是一个高为1,半径不等的圆柱体。

输入
第一行包含两个正整数N和F,1 ≤ N, F ≤ 10 000,表示派的数量和朋友的数量。
第二行包含N个1到10000之间的整数,表示每个派的半径。
输出
输出每个人能得到的最大的派的体积,精确到小数点后三位。
样例输入
3 3
4 3 3
样例输出
25.133
惨痛的教训:浮点数最好转化成整数再算,不然精度有问题。
#include <iostream>
#include <cstdio>
using namespace std;
double r,l,mid;
double a[];
int n,k,sum;
int judge(double x)
{
int ans,i;
ans=;
for (i=;i<=n;i++)
ans+=a[i]/x;
return ans;
}
int main()
{
int i;
scanf("%d%d",&n,&k);
l=; r=;
for (i=;i<=n;i++)
{
scanf("%lf",&a[i]);
a[i]=a[i]*a[i]*3.141592653589793*;
r=max(r,a[i]);
}
//r=3.14159263589793*100000*100000;
sum=;
while (l<r)
{
mid=(r+l)/;
sum=judge(mid);
if (sum>=k+)
l=mid;
else
r=mid;
if (r-l<)
break;
}
if (judge(r)==k+)
printf("%.3lf",r/);
else
printf("%.3lf",l/);
return ;
}