【三分】51nod 1629 B君的圆锥

51nod 1629 B君的圆锥（三分）

/*
time:2016/12/2
write:Swords

题意：给你一个圆锥的表面积大小，让你求这样表面积的圆锥体积最大是多少
题解：三分法根据选半径确定最大体积，r的范围不大于10^5（因为S不大于10^9）。
*/
#include
#include
#include
#define PI 3.1415926535898
#define eps 1e-7
using namespace std;
int s;
double val(double r)
{
    double l=s*1.0/(PI*r)-r;
    if(l<=r) return -1;
    double h=sqrt(l*l-r*r);
    double v=PI*r*r*h/3.0;
    return v;
}
double solve(double l,double r)
{
    while(l + eps < r)
    {
        double lmid = l + (r-l)/3,rmid = r - (r-l)/3;
        if(val(lmid) >= val(rmid))
            r = rmid;
        else
            l = lmid;
    }
    return val(l);
}
int main()
{
    while(~scanf("%d",&s))
        printf("%lf\n",solve(0,100000));
    return 0;
}