POJ 1905 Expanding Rods【二分答案】

题意：已知木棍原长度，温度，系数，加热后弯曲，求弯曲前后木棍中点的距离

如图，蓝色为杆弯曲前，长度为L

红色为杆弯曲后，长度为s

h是所求

依题意知

S=(1+n*C)*L

 

又从图中得到三条关系式;

（1）       角度→弧度公式  θr = 1/2*s

（2）       三角函数公式  sinθ= 1/2*L/r

（3）       勾股定理  r^2 – ( r – h)^2 = (1/2*L)^2

 

把四条关系式化简可以得到

逆向思维解二元方程组：

要求（1）式的h，唯有先求r

但是由于（2）式是三角函数式，直接求r比较困难

 

因此要用顺向思维解方程组：

在h的值的范围内枚举h的值，计算出对应的r，判断这个r得到的(2)式的右边  与 左边的值S的大小关系  （ S= (1+n*C)*L ）

 

很显然的二分查找了。。。。。

那么问题只剩下 h 的范围是多少了

下界是0 (不弯曲)

确定上界：题意已告知不能延伸超过其初始长度的一半

就是说 S max = 3/2 L

理论上把上式代入（1）(2)方程组就能求到h的最小上界，但是实际操作很困难

因此这里可以做一个范围扩展，把h的上界扩展到 1/2L  ，不难证明这个值必定大于h的最小上界，那么h的范围就为  0<=h<1/2L

这样每次利用下界low和上界high就能得到中间值mid，寻找最优的mid使得(2)式左右两边差值在精度范围之内，那么这个mid就是h

以上转自 http://blog.csdn.net/lyy289065406/article/details/6648562

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

int main()
{
   // freopen("in.txt","r",stdin);
    double newL,c;
    double r,l,n;
    while(scanf("%lf%lf%lf",&l,&n,&c)!=EOF)
    {
        double low,high,mid;
        if(l<0&&n<0&&c<0) break;

        newL = (1+n*c)*l;
        low = 0.0,high = 0.5*l;

        while( high - low > 0.0001)
     // while( low - high < 0.0001) //脑残这样写了，明明是非正数肯定小于它的
        {
            mid = (low+high)/2.0;
            r = (4*mid*mid+l*l)/(8*mid);
            if( newL > asin(l/(2*r))*2*r  )  low = mid ;//范围要求高，+1的话加的太大
            else   high = mid;
        }
        printf("%.3lf\n",high);
    }
     return 0;
}