POJ 1905 Expanding Rods

http://poj.org/problem?id=1905

---------------------------------------------------------------------------------

思路：这题感觉难点在几何上，数学真心无语啊；

r=(4*h*h+L*L)/(8*h);

s=2*r*arc sin(L/(2*r));

求h得先求r，由s比较难以求得r的值，于是，直接用二分来做，先将中点带入第一个式子中，求得r的值，再带入2式中，求s的值，与s的原值比较，接下来就是二分没啥好说的了；

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<stdlib.h>
const double esp=1e-5;//最小精度
int main()
{
    double L,n,c,s;
    double h;
    double r;
    double mid;
    while(scanf("%lf%lf%lf",&L,&n,&c)!=EOF)
    {
        if(L==-1&&n==-1&&c==-1)
        break;
        double low=0;
        double high=0.5*L;
        double mid;
        s=(1+n*c)*L;
        while(high-low>esp)//这里不能直接用hign>low,否则当low与high很接近时会陷入死循环
        {
            mid=(high+low)/2;
            r=(4*mid*mid+L*L)/(8*mid);
            if(2*r*asin(L/(2*r))<s)
            low=mid;
            else high=mid;
        }
        h=mid;
        printf("%.3lf\n",h);
    }
    return 0;
}