1 /* 2 二分 + 几何 3 弧长L, 圆半径R, 弧度 q, L=R*q; 4 二分: 弧度(0~PI) 或者 高度(L/2~L) 5 */ 6 #include<cstdio> 7 #include<iostream> 8 #include<cmath> 9 using namespace std; 10 const double PI = acos(-1.0); 11 double L, L1, T, C, R, Q; 12 13 int main(){ 14 15 while(scanf("%lf%lf%lf", &L, &T, &C) && (L!=-1 || T!=-1 || C!=-1)){ 16 L1 = (1+T*C)*L; 17 double ld=0.0, rd=PI, tmp; 18 Q = (ld+rd)/2; 19 R = L/2/sin(Q); 20 while(fabs(tmp = R*2*Q-L1)>1e-8){ 21 if(tmp>0) rd=Q; 22 else ld=Q; 23 Q=(ld+rd)/2; 24 R = L/2/sin(Q); 25 } 26 printf("%.3lf ", R-sqrt(R*R - (L/2)*(L/2))); 27 } 28 return 0; 29 }