09:膨胀的木棍

09:膨胀的木棍


总时间限制: 
1000ms 
内存限制: 
65536kB
描述

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

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

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

09:膨胀的木棍_第1张图片

输入
三个非负实数:木棍初始长度(单位:毫米),温度变化(单位:度),以及材料的热膨胀系数。
保证木棍不会膨胀到超过原始长度的1.5倍。
输出
木棍中心的偏移距离(单位:毫米),保留到小数点后第三位。
样例输入
1000 100 0.0001
样例输出
61.329
来源
Waterloo local 2004.06.12

点评:二分求角的极限法 参考 http://www.cnblogs.com/huashanqingzhu/p/5662002.html


#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
using namespace std;
typedef long long LL;
#define MAX 100001
int a[MAX];
int b[MAX];

int main() {
    //freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
    double L_before; //原始长度
    double L_After; //形变后的长度
    double rig, lef, mid; //角的极大极小值
    double n, c;
    scanf("%lf%lf%lf", &L_before, &n, &c);
    if(L_before < 1e-14){
        printf("0.000\n");
        return 0;
    }

    L_After = L_before * (1 + n * c);

    lef = 0.0; //角的极小值
    rig = asin(1.0); //角的极大值

   while(rig - lef>1e-14){
        mid = (rig + lef) / 2;
        if(L_After * sin(mid) / mid <= L_before){

            rig = mid;
        }else {
            lef = mid;
        }
    }

    printf("%.3lf\n", L_before / 2 * tan(lef / 2));

    return 0;

}





你可能感兴趣的:(POJ水题)