题目链接：http://noi.openjudge.cn/ch0111/05/

总时间限制: 1000ms

    内存限制: 65536kB

描述

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

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

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

输入

第一行包含两个正整数N和F，1 ≤ N, F ≤ 10 000，表示派的数量和朋友的数量。

第二行包含N个1到10000之间的整数，表示每个派的半径。

输出

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

样例输入

3 34 3 3

样例输出

25.133

思路：在min和max之间进行二分。min是0，max是最大的那一个派的体积。

这个题要注意精度的控制，我这里二分继续的条件是max-min还大于0.0000001。

另外，关于∏的取值，应该用数学库函数里面的反余弦函数acos()。就是PI=acos(-1.0)。   （解析：余弦cos(∏)=-1,所以反余弦acos(-1)=∏）

还有一点，看下面的代码：

while( max-min > cmpPrecision )

{

    mid=min+(max-min)/2;

    count=compute(v,N,mid,F);

    if(count>=0) min=mid;

    else max=mid;

}

红色部分，当等号成立时为何不结束二分呢？这个问题主要是因为这里是对浮点数区间进行二分，而且分割出来的块要尽量大。试想，比如有一个大块总体积11，分割为每块体积3.1可以分3块，但是分割为每一块体积3.5同样也可以分成3块。所以，为了满足分割出来的块尽量大这一条件，必须用if(count>=0) min=mid。

1 #include 
     
       2 #include
      
        3 #include
       
         4 #include
        
          5  6 #define cmpPrecision 0.0000001 7  8 int compute(double *v,long N,double mid,long F);//假定分发的派的体积是mid，统计能够分割的个数ans。根据ans与F的大小关系返回1,0，-1. 9 int main()10 {11     long N,F;12     long r,h=1;13     double *v=NULL,min=0,max=0,mid=0;14     double PI=acos(-1.0);//获取PI的值15     long i;16     long count;17 18     freopen("5.in","r",stdin);19     scanf("%ld%ld",&N,&F);20     F++;21     v=(double *)malloc(N*sizeof(double));22     memset(v,0,N);23     for(i=0;i
         
          max) max=v[i];28     }29 30     while( max-min > cmpPrecision )31     {32         mid=min+(max-min)/2;33         count=compute(v,N,mid,F);34         if(count>=0) min=mid;//注意，当count==0的时候，并不一定就可以结束二分。可能可以分割更大的块却得到相同的数量。（主要是因为块的体积是浮点数。比如有一块总体积11，假如每一块体积3.1可以分3块。但是每一块3.5也同样只能分3块。题目要求块大小尽量大，所以当分割结果刚好合适时不一定就应该结束二分。）35         else max=mid;36     }37     printf("%.3lf\n",min);38     return 0;39 }40 int compute(double *v,long N,double mid,long F)//假定分发的派的体积是mid，统计能够分割的个数ans。根据ans与F的大小关系返回1,0，-1.41 {42     long ans=0,i;43     for(i=0;i
          
           F) return 1;47 }48 if(ans==F) return 0;49 else return -1;50 }
          
         
        
       
      
     

额这段代码还要注意一些库函数，比如malloc，memset，acos，floor等。