Pie POJ 3122(二分找最值)

题目大意：给定n,f分别为Pie的数量和朋友的数量，接下来是Pie的半径，求在每个朋友只能拿一块蛋糕的条件下，能分得的最大体积的Pie是多少。

思路：直接二分体积，（PS：注意二分之后枚举区间的移动l=mid,r=mid!不要写成没有精度问题时的l=mid+1，PI=acos(-1.0)）. 判断枚举区间的位置依据条件当前的体积可以够多少人分的。

#include<map>
#include<queue>
#include<cmath>
#include<iostream>
#include<cstdio>
#include<stack>
#include<cstring>
#include<algorithm>
#define inf 0x3f3f3f3f
#define eps 1e-5
using namespace std;
const double PI=acos(-1.0);
int r[100010],n,m;

int op(int a,int b){
    return a>b;
}
int fi(double t){
    int tmp,sum=0;
    for(int i=0;i<n;i++){
        sum+=(r[i]*r[i]*PI)/t;
    }

    if(sum>=m)
        return 1;
    return 0;
}
int main(){
    int i,j,k,cla;
    scanf("%d",&cla);
    while(cla--){
        scanf("%d%d",&n,&m);
        ++m;
        int mi=0x3f3f3f3f;
        double o=0;
        for(i=0;i<n;i++){
            scanf("%d",&r[i]);
            o+=r[i]*r[i]*PI;
        }
        sort(r,r+n,op);
        if(n>m){
            printf("%.4f\n",r[m-1]*r[m-1]*PI);
            continue;
        }
        if(n==1){
            printf("%.4f\n",(r[0]*r[0]*PI)/m);continue;
        }
        double ls,rs,ans;
        ls=0;
        rs=o/m;
        while(rs-ls>eps){
            double mid=(ls+rs)/2;
            if(fi(mid)){
                ls=mid;
                ans=mid;
            }
            else
                rs=mid;
        }
        printf("%.4f\n",ans);
    }
    return 0;
}