#include #include #include #include using namespace std; const double eps= 1e-5; double H,h,D; double f(double d) { return D-d+H-(H-h)*D/d; // return D-d+h-(H-h)*(D-d)/d; } double sanfen(double l,double r) { double mid,midr,ans; while(fabs(r-l)>eps) { mid=(l+r)/2; midr=(mid+r)/2; if(f(mid)