#define rep(i,n) for (int i=0;i < (int)(n);i++) #define ALL(v) v.begin(),v.end() typedef long long ll; #include using namespace std; double r,h,d; //被積分関数 (acos(x)でx>1 や sqrt(x)でx<0 だと nan になるので端っこ注意) double f(double x) { double p=min(1.0,d*h/(2*r*x)); double q=max(0.0,4*r*r*x*x-d*d*h*h); return 2*r*r/(h*h)*x*x*acos(p)-d/h*sqrt(q)/2; } int main() { cin>>r>>h>>d; double A=d*h/(2*r); //下端 double B=h; //上端 int N=100000; const double delta=(B-A)/N; double value=0.0; //i=1からでも大して影響ない for(int i=0;i