【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[10010];
int n,k,sum;
int judge(double x)
{
    int ans,i;
    ans=0;
    for (i=1;i<=n;i++)
        ans+=a[i]/x;
    return ans;
}
int main()
{
    int i;
    scanf("%d%d",&n,&k);
    l=0; r=0;
    for (i=1;i<=n;i++)
    {
        scanf("%lf",&a[i]);
        a[i]=a[i]*a[i]*3.141592653589793*100000;
        r=max(r,a[i]);
    }
    //r=3.14159263589793*100000*100000;
    sum=0;
    while (l<r)
    {
        mid=(r+l)/2;
        sum=judge(mid);
        if (sum>=k+1)
            l=mid;
        else
            r=mid;    
        if (r-l<1)
            break;
    }
    if (judge(r)==k+1)
        printf("%.3lf",r/100000);
    else
        printf("%.3lf",l/100000);
    return 0;
}