POJ 3737
题目大意
暂无
题目解法
暂无
RTFC
#include <cmath>
#include <cstdio>
const double pi = acos(-1.0), eps = 1e-6;
double S;
inline double calc(double r)
{
double l = (S - r * r) / r;
double h = sqrt(l * l - r * r);
return pi * r * r * h / 3;
}
int main()
{
while (~scanf("%lf", &S))
{
S /= pi;
double l = 0, r = sqrt(S), m1, m2, d;
while ((d = r - l) > eps)
calc(m1 = l + d / 3) < calc(m2 = r - d / 3) ? l = m1 : r = m2;
l = (S - r * r) / r;
double h = sqrt(l * l - r * r);
double V = pi * r * r * h / 3;
printf("%.2f\n%.2f\n%.2f\n", V, h, r);
}
return 0;
}