【20171005】Openjudge 1.11-09 膨胀的木棍

不得不说，写这道题的时候，我也快膨胀了~~~

总时间限制:

1000ms

内存限制:

65536kB

描述

当长度为L的一根细木棍的温度升高n度，它会膨胀到新的长度L'=(1+n*C)*L，其中C是热膨胀系数。

当一根细木棍被嵌在两堵墙之间被加热，它将膨胀形成弓形的弧，而这个弓形的弦恰好是未加热前木棍的原始位置。

你的任务是计算木棍中心的偏移距离。

输入

三个非负实数：木棍初始长度（单位：毫米），温度变化（单位：度），以及材料的热膨胀系数。
保证木棍不会膨胀到超过原始长度的1.5倍。

输出

木棍中心的偏移距离（单位：毫米），保留到小数点后第三位。

样例输入

1000 100 0.0001

样例输出

61.329

来源

Waterloo local 2004.06.12

生活总是充满了欢乐，这根小木棍感受到了夏日浓浓的热情，也禁不住膨胀起来了~~~

那么，是谁那么无聊让我们去算它膨胀了多少呢？

这个神奇的“偏移距离”

#include<stdio.h>
#include<math.h>
double searchN(double l,double ln,double head,double end)
{
    double mid,r,ln_temp,temp;
    mid=(head+end)/2.0;//mid是猜测的弧度 
    while(fabs(ln_temp-ln)>=1e-14)
    {
    mid=(head+end)/2;
        r=(l/2.0)/sin(mid/2.0);
        ln_temp=mid*r;
        if(ln_temp-ln>=1e-14)
        {
            temp=mid;
            end=temp;
            continue;
        }
        else if(ln_temp-ln<-1e-14)
        {
            temp=mid;
            head=temp;
            continue;
        }
    }
    return mid;
}
int main()
{
    double c,n,l;
    double ln;//L_New
    scanf("%lf%lf%lf",&l,&n,&c);
    ln=(1+n*c)*l;
    
    double k=searchN(l,ln,0,M_PI);
    double h=(l/2.0)/sin(k/2)-(l/2.0)/tan(k/2);
    
    printf("%.3f
",h);
    return 0;
}