POJ 3122 Pie【二分答案】

题意：给出n个圆，分别的半径，现在需要把它们分成f+1块，问每一块的面积

和上面一题一样，二分面积，上限是这里面最大的圆的面积

另外可以先只二分r*r,到最后再乘上pi,不过这题先乘后乘都能过

不过不懂的是,输出是lf就一直wa,后来看了discuss,改成f就过了

#include<iostream>  
#include<cstdio>  
#include<cstring> 
#include <cmath> 
#include<stack>
#include<vector>
#include<map> 
#include<set>
#include<queue> 
#include<algorithm>  
using namespace std;

typedef long long LL;
const int INF = (1<<30)-1;
const int mod=1000000007;
const int maxn=1000005;
const double pi=3.14159265359;

double a[maxn],r[maxn];
int n,f;

int ok(double x){
    int ans=0;
    for(int i=1;i<=n;i++){
        ans+=(int)(a[i]/x);
    }
    if(ans<(f+1)) return 0;
    return 1;
}

int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d %d",&n,&f);
        double smax=-1;
        for(int i=1;i<=n;i++){
            scanf("%lf",&r[i]);
            a[i]=r[i]*r[i];
            
        //    printf("a[%d]=%lf
",i,a[i]);
            smax=max(smax,a[i]);
        }
        
        double l=0.000,r=smax,mid;
        while(r-l > 1e-6){
            mid=(l+r)/2;
            if(ok(mid)) l=mid;
            else r=mid;
    //        printf("l=%lf
",l);
    //        printf("r=%lf
",r);
    //        printf("mid=%lf
",mid);
        }
        printf("%.4f
",l*pi);
        
    }
    return 0;
}