poj Pie 二分

#include<stdio.h> #include<string.h> #include<algorithm> #define PI 3.1415926535897932 using namespace std; double a[20000]; int T,sum; int main() { scanf("%d",&T); while(T--) { int n,f; scanf("%d%d",&n,&f); f=f+1; sum=0; double Max=0.0; for(int i=0;i<n;i++) { scanf("%lf",&a[i]); a[i]=1.0*a[i]*a[i]; sum+=a[i]; } double low,high,mid; int num=0; low=0.0;//每个人能分到的下界 high=sum/f;//上界 while(high-low>0.000001) { mid=(low+high)/2; num=0; for(int i=0;i<n;i++) { num+=(int)(a[i]/mid); } if(num<f)//进行判断 改变上界与下界的值 { high=mid; } else { low=mid; } } printf("%.4f\n",mid*PI); } return 0; }