#include using namespace std; #define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i> R >> H >> D; constexpr int Z = 1 << 9; const double T = D / (2 * R); double ret = 0.0; const double lo = T; const double dz = 1.0 / Z; REP(t, Z) { double u = 1.0 - (dz * dz * (t + 0.5) * (t + 0.5)) * (1 - lo); double dd = 2 * dz * dz * (t + 0.5) * (1 - lo); double theta = acos(T / u); ret += u * u * (theta * 0.5 - sin(theta * 2) * 0.25) * dd; } cout << fixed << setprecision(20) << ret * 4 * R * R * H << '\n'; }