OpenJudge 派

目录

 

派

要求：

 描述：

 输入：

输出：

样例输入： 

样例输出： 

问题分析： 

细节处理：

最终代码：

总结：

---

派

要求：

总时间限制: 1000ms

内存限制: 65536kB

 描述：

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

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

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

 输入：

第一行包含两个正整数N和F，1 ≤ N, F ≤ 10 000，表示派的数量和朋友的数量。
第二行包含N个1到10000之间的整数，表示每个派的半径。

输出：

 输出每个人能得到的最大的派的体积，精确到小数点后三位。

样例输入： 

3 3
4 3 3

样例输出： 

25.133

问题分析： 

我们还是先来分析题意，对于N个给定的圆柱形派，我们怎么分割使得将这些派平均且最大分给含“我”在内的F+1个人？

那么我们首先面对的问题是如何处理题目中的派（圆柱体）？我们知道圆柱的体积公式为

，而对于所有的派，高均为1，所以决定V的大小的关键因素就是，所以在存储这些派的相关信息时，我们不要像给出的数据那样存不同派的半径，我们要存的是，在后续操作中可以直接对进行划分。

所以对于样例，我们存了这样的一个数组，我们的目标就是尽可能大并且平均地将这三个数分成四份，由于派的数量N最大可以到10000，所以一般的搜索枚举一定是会超时的。（此题精度是三位小数）所以我们在没有好的贪心算法的情况下，利用二分法来快速搜索最优的情况。

细节处理：

这一部分十分重要，关乎是否能够AC

1）Pi 不要用自己#define的3.1415926……，这个的结果是Wrong Answer，可以采用acos(-1.0)的方法来得到Pi。

2）二分时终止条件的精度需要调整到1e-5，1e-4以及1e-3都是错误的。

3）小优化：在判断当前情况下派的划分是否过多时，如果超过了F+1，可以直接跳出循环，不用在进行下去；在定义二分上限时，可以定义成SumPie/(F+1)，使上限尽可能小又不失准确性。

最终代码：

#include<iostream>
#include<cmath>
using namespace std;
const double Pi=acos(-1.0);//小细节
int main(){
    int N;//派的数量
    int F;//朋友数量
    scanf("%d %d",&N,&F);
    int* Pie=new int[N];
    double SumPie=0;
    for(int i=0;i<N;i++){
        scanf("%d",&Pie[i]);
        Pie[i]*=Pie[i];
        SumPie+=Pie[i];
    }
    double Left=0;double Right=SumPie/(F+1);
    double Middle=Left+(Right-Left)/2;
    while(Right-Left>1e-5){//小细节
        int Division=0;
        for(int i=0;i<N;i++){
            Division+=int(Pie[i]/Middle);
            if(Division>=F+1){
                Left=Middle;
                break;
            }
        }
        if(Division<F+1) {
            Right = Middle;
        }
        Middle = Left + (Right - Left) / 2;
    }
    printf("%.3lf\n",Middle*Pi);
    delete []Pie;
}

总结：

这个题目的二分比上一个题目“月度开销”要简单明了得多。如果对二分法还不了解，可以参考我的上一篇文章。此题最关键的地方就是需要将，而非存起来进行划分。希望这篇文章可以为您带来帮助！