几何。。。。。。
好讨厌。。。。。感觉自己用的公式不对。。。。
#include <stdio.h> #include <math.h> #define eps 1e-9 #define pai acos(-1) double r,R; int main() { while(scanf("%lf%lf",&r,&R)!=EOF) { if(2.0*r>R) {printf("NO Solution!\n");continue;} double di=R*sqrt(3)*0.5; double low=0,hig=di,mid=di/2.0; while(hig-low>eps) { // double tmp=R*R-di*di; double t1=asin(mid/R); double len=2*R*sin((pai-t1)/2.0); double rr=2*len*len*sin(mid/R)/(len*2+2*mid); printf("bian=%lf\n",len); if(rr-r>eps) hig=mid; else low=mid; mid=(low+hig)/2.0; } // printf("mid=%lf\n",mid); double t1=asin(mid/R); double len=2*R*sin((pai-t1)/2); // printf("mid/r=%lf asin=%lf len=%lf\n",mid/R,asin(mid/R),len); printf("%.18lf %.18lf %.18lf\n",mid*2,len,len); } return 0; }