POJ 3737 UmBasketella(三分）

题意：给出圆锥的表面积（包含底面）。求其最大体积，以及此时的底面半径及高

#include <cmath>
#include <iostream>
using namespace std;

#define eps 0.00000005
#define PI acos(-1.0)  // PI用反三角函数比较准确

int main()
{
	double s;
	while ( scanf("%lf",&s) != EOF )
	{
		double v = 0.0;
		double low = 0, high = sqrt ( s / PI );
		while ( high - low > eps )
		{
			double mid1 = low + ( high - low ) / 3;
			double mid2 = high - ( high - low ) / 3;
			double v1 = PI * mid1 * mid1 * sqrt ( pow((s-PI*mid1*mid1)/(PI*mid1),2) - mid1*mid1 ) / 3;
			double v2 = PI * mid2 * mid2 * sqrt ( pow((s-PI*mid2*mid2)/(PI*mid2),2) - mid2*mid2 ) / 3;
			if ( v1 > v2 )
			{
				v = v1;
				high = mid2;
			}
			else
			{
				v = v2;
				low = mid1;
			}
		}
		double r = high;
		double h = sqrt ( pow((s-PI*r*r)/(PI*r),2) - r*r );
		printf("%.2lf\n%.2lf\n%.2lf\n", v, h, r);
	}
	return 0;
}