LA 3635 Pie --- 二分答案

ICPC Live Archive 题目链接 https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1636

题意：F+1个人分N个圆形蛋糕，每个人需要分到一片蛋糕，这片是从一块蛋糕中分出来而不是多个蛋糕中的部分拼出来的，要求每个人分到的蛋糕面积相同，求每人分到的蛋糕最大面积是多少？

题解：二分答案即可，答案x，min(x) = 0, max(x) = max(蛋糕面积)，然后对x进行二分，每次判断x是否满足要求，就是遍历蛋糕，用蛋糕面积除以x再取整，算出蛋糕最多能否分给x个人。 注意因为面积是浮点数，所以二分循环条件要写成r-l>tiny，tiny是一个很小的浮点数，不过这个浮点数别设置得太小了，不然会TLE 

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#define MAXN 10010
#define PI acos(-1)
using namespace std;
int F,N;
double res = 0.0;
int r[MAXN];

// res:每个人分到的面积大小
bool isok(double res) {
    int sum = 0;
    for(int i = 0;i < N;i++) {
        double s = PI * (double)r[i] * r[i];
        int t = s / res;
        sum += t;
    }
    return sum >= F+1;
}

int main() {
    int T;
    scanf("%d",&T);
    while(T--) {
        F = N = 0;
        res =0.0;
        memset(r,0,sizeof(r));
        scanf("%d%d",&N,&F);
        for(int i = 0;i < N;i++) {
            scanf("%d",&r[i]);
        }
        sort(r,r+N);
        double L = 0;
        double R = PI * (double)r[N-1] * r[N-1];
        while(R-L>1e-5) {
            double m = (R+L)/2.0;
            if(isok(m)) {
                L = m;
            } else {
                R = m;
            }
        }
        printf("%.4f\n",R);
    }

    return 0;
}