二分法解题报告

题目：派

描述：

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

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

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

输入：

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

输出：

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

样例输入：

3 3

4 3 3

样例输出：

25.133

 

 

 

代码：

#include<iostream>

#include<algorithm>

#define_USE_MATH_DEFINES

#include<math.h>

#include<iomanip>

usingnamespace std;

#defineEPS 0.00000000001

 

bool judge(int *a, doubler, intN, intF)

{

       int num = 0;

       for (int i = 0; i < N; i++)

       {

              num +=(a[i] * a[i]) / (r*r);

              if (num >= F)

                     returntrue;

       }

       returnfalse;

}

 

double twiceserch(int* a, intmaxr, intN, intF)

{

       double MAX = maxr, MIN = 0;

       while (MAX - MIN > EPS)

       {

              double mid = (MAX + MIN)/ 2;

              if (judge(a, mid, N, F))

                     MIN= mid;

              else

                     MAX= mid;

       }

       return MIN;

}

 

int main()

{

       int N, F, maxr = 0;

       cin >> N >> F;

       F += 1;

       int* a = newint[N];

       for (int i = 0; i < N;i++)

       {

              cin >> a[i];

              if (a[i] > maxr)

                     maxr= a[i];

       }

       double r = twiceserch(a,maxr, N, F);

       cout << fixed << setprecision(3) <<M_PI*r*r << endl;

       return 0;

}