[POJ 3737]UmBasketella(三分)

题目链接：http://poj.org/problem?id=3737
题目大意：已知圆锥表面积$n$，求圆锥体积$V$最大时的圆锥高度$h$,底面半径$r$,以及圆锥体积$V$。
思路：圆锥体积$V$与底面半径$r$成单调上凸关系(换句话说即函数$V=f(r)$的导函数为单调减并有零点)。那么我们可以采用三分来求$V$的极值。下图是三分的过程，其中$lowerBound、upperBound$是当前三分的上下界，而$mid1、mid2$则各为当前区间$[lowerBound,upperBound]$的$\frac{1}{3}$点以及$\frac{2}{3}$点。$val1、val2$各为$mid1、mid2$的对应值。

然后三分的过程比较简单，类似二分，若$val1<val2$，则$lowerBound=mid1$，否则$upperBound=mid1$

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
#include <cmath>

#define PI (3.1415926535897384626)
#define EPS (1e-7)

using namespace std;

double S; //圆锥的表面积

double calc(double r) //已知圆锥底面半径r，计算当前圆锥的体积
{
    double R=(S-PI*r*r)/PI/r; //侧面展开的扇形半径
    double h=sqrt(R*R-r*r); //
    return PI*r*r*h/3;
}

int main()
{
    while(scanf("%lf",&S)!=EOF)
    {
        double lowerBound=0,upperBound=sqrt(S/PI);
        double mid1,mid2,val1,val2;
        while(lowerBound+EPS<upperBound)
        {
            mid1=lowerBound+(upperBound-lowerBound)/3;
            mid2=upperBound-(upperBound-lowerBound)/3;
            val1=calc(mid1);
            val2=calc(mid2);
            if(val1<val2) lowerBound=mid1;
            else upperBound=mid2;
        }
        double R=(S-PI*upperBound*upperBound)/PI/upperBound; //圆锥的侧面展开后的扇形半径
        double h=sqrt(R*R-upperBound*upperBound); //最终的圆锥高度
        double V=PI*upperBound*upperBound*h/3; //最终的圆锥最大体积
        printf("%.2f\n%.2f\n%.2f\n",V,h,upperBound);
    }
    return 0;
}