ZOJ 4806

几何。。。。。。

好讨厌。。。。。感觉自己用的公式不对。。。。

#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;
}