[复习][poj3737]三分 UmBasketella

题目背景
poj 3737

题目描述
给定圆锥的表面积S ，求这个圆锥的最大体积V ，以及此时它的高h与底面半径r。

输入格式
输入包含几个测试数据。每个案例包含只有一个实数，代表圆锥的表面积。保证1≤S≤10000。

输出格式
对于每一个测试数据，输出应包含三行。
第一行一个实数代表最大锥体体积。
在第二、第三行输出圆锥高度和圆锥底面半径。
所有实数都四舍五入到0.01。

样例数据
输入
30
输出
10.93
4.37
1.55

分析：一道三分题，证明如下：
$S_表=\pi*r*l+\pi*r^2$
$\Rightarrow$$l=\frac{S_表}{π∗r}−r$ （l是母线长）
$\Rightarrow$$h=\sqrt{l^2-r^2}=\sqrt{\frac{{S_表}^2}{π^2∗r^2}−\frac{2∗{S_表}}{π}}$
$\Rightarrow$$V=\frac{1}{3}∗π∗r^2∗h=\frac{1}{3}∗π∗r^2∗\sqrt{\frac{{S_表}^2}{π^2∗r^2}−\frac{2∗{S_表}}{π}}=\frac{1}{3}∗r∗\sqrt{{S_表}^2−2∗{S_表}∗π∗r^2}$
不难看出最后推出的体积表达式是一个V关于r的单峰函数，所以我们可以三分r的取值从而得到答案。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<cctype>
#include<iomanip>
#include<queue>
#include<set>
using namespace std;

int getint()
{
    int f=1,sum=0;
    char ch;
    for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar());
    if(ch=='-')
    {
        f=-1;
        ch=getchar();
    }
    for(;ch>='0'&&ch<='9';ch=getchar())
        sum=(sum<<3)+(sum<<1)+ch-48;
    return sum*f;
}

const double pi=acos(-1.0);//反函数求pi更精确
const double eps=1e-6;
double s,l,r,midl,midr;
double ansr,ansl,ansh,ansv;

double find(double r)
{
    double l=(s-r*r)/r;
    double h=sqrt(l*l-r*r);
    double v=pi*r*r*h/3;
    return v;
}

int main()
{
    freopen("UmBasketella.in","r",stdin);
    freopen("UmBasketella.out","w",stdout);

    while(scanf("%lf",&s)!=EOF)
    {
        s=s/pi;
        l=0,r=sqrt(s);
        while(r-l>eps)
        {
            midl=l+(r-l)/3,midr=r-(r-l)/3;
            if(find(midl)<find(midr))
            {
                l=midl;
                ansr=midr;
            }
            else
            {
                r=midr;
                ansr=midl;
            }
        }

        ansl=(s-ansr*ansr)/ansr;
        ansh=sqrt(ansl*ansl-ansr*ansr);
        ansv=pi*ansr*ansr*ansh/3;
        printf("%0.2f\n%0.2f\n%0.2f\n",ansv,ansh,ansr);
    }

    return 0;
}

本题结。