题目背景
poj 3737
题目描述
给定圆锥的表面积S ,求这个圆锥的最大体积V ,以及此时它的高h与底面半径r。
输入格式
输入包含几个测试数据。每个案例包含只有一个实数,代表圆锥的表面积。保证1≤S≤10000。
输出格式
对于每一个测试数据,输出应包含三行。
第一行一个实数代表最大锥体体积。
在第二、第三行输出圆锥高度和圆锥底面半径。
所有实数都四舍五入到0.01。
样例数据
输入
30
输出
10.93
4.37
1.55
分析:一道三分题,证明如下:
S表=π∗r∗l+π∗r2
⇒
l=S表π∗r−r
(l是母线长)
⇒
h=l2−r2−−−−−√=S表2π2∗r2−2∗S表π−−−−−−−−−−−√
⇒
V=13∗π∗r2∗h=13∗π∗r2∗S表2π2∗r2−2∗S表π−−−−−−−−−−−√=13∗r∗S表2−2∗S表∗π∗r2−−−−−−−−−−−−−−−−−√
不难看出最后推出的体积表达式是一个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;
}
本题结。