poj-3737(UmBasketella)

三分法入门题。

题目大意：给出圆锥的表面积，求最大体积的圆锥的高和底面半径。附链接：http://poj.org/problem?id=3737

大体思路：这道题由圆锥体积可知，当高h不变时，体积和底面半径 r 有二次关系，直接用三分法求出体积最大时的底面半径。或者数学好的话可以直接求出表面积和圆锥高、体积、底面半径的公式，在poj的 Disscus 看到的，见识了，还可以这样求，厉害厉害。

以下是ac代码：

#include<iostream>
#include<cmath>
#include<cstring>
#include<stdio.h>
using namespace std;

const double error=1e-4;
const double pi=acos(-1.0);

int main(){
    double s;
    while(~scanf("%lf",&s)){
        double r1,r2,h1,h2;
        double left=0;
        double right=sqrt(s/pi);
        double v1,v2;
        while(right-left>error){
            r1=left+(right-left)/3.0;
            h1=sqrt((s/pi/r1-r1)*(s/pi/r1-r1)-r1*r1);
            r2=right-(right-left)/3.0;
            h2=sqrt((s/pi/r2-r2)*(s/pi/r2-r2)-r2*r2);
            v1=pi*r1*r1*h1/3.0;
            v2=pi*r2*r2*h2/3.0;
            if(v1<v2)
                left=r1;
            else
                right=r2;

        }
        printf("%.2lf\n%.2lf\n%.2lf\n",v1,h1,r1);
    }
    return 0;
}

这代码在 oj 上用C++提交，如果用G++会wa，原因是编译器自身的问题，这里就不解释了，自行百度吧（我也是踩坑过来的）。