POJ 3737 UmBasketella G++ 三分 背

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;
//英语 看博友分析 抄博友程序  不均匀三分未实现  均匀三分  背 
// 圆锥侧面积公式 s=ip*r*l  圆锥体积公式 v=(1/3.0)*h*ip*r*r
double s;//圆锥表面积 
double ip;
double getH(double r)
{
	double l=(s-ip*r*r)/(r*ip);
	return sqrt(l*l-r*r);
}
double getV(double r)
{
	double h=getH(r);
	return (1.0/3.0)*h*ip*r*r;
}
int main()
{
	ip=3.141592653589793284626;
	while(cin>>s)
	{
		double lb=0;
		double rb=sqrt(s/ip);
		double jg;
		for(int i=0;i<40;i++)
		{
			double ml=lb+(rb-lb)/3.0;
			double mr=rb-(rb-lb)/3.0;
			if(getV(ml)<getV(mr))
			{
				lb=ml;
				jg=mr;
			}else
			{
				rb=mr;
				jg=ml;
			}
		}
		cout<<fixed<<setprecision(2)<<getV(jg)<<endl<<getH(jg)<<endl<<jg<<endl;
	}
	return 0;
}