poj 1905 Expanding Rods

最新推荐文章于 2020-05-02 19:02:25 发布

G_rass

最新推荐文章于 2020-05-02 19:02:25 发布

阅读量343

点赞数

分类专栏：分治&二分

本文链接：https://blog.csdn.net/G_rass/article/details/46971617

版权

分治&二分专栏收录该内容

8 篇文章 0 订阅

订阅专栏

Description

When a thin rod of length L is heated n degrees, it expands to a new length L'=(1+n*C)*L, where C is the coefficient of heat expansion.
When a thin rod is mounted on two solid walls and then heated, it expands and takes the shape of a circular segment, the original rod being the chord of the segment.

Your task is to compute the distance by which the center of the rod is displaced.

Input

The input contains multiple lines. Each line of input contains three non-negative numbers: the initial lenth of the rod in millimeters, the temperature change in degrees and the coefficient of heat expansion of the material. Input data guarantee that no rod expands by more than one half of its original length. The last line of input contains three negative numbers and it should not be processed.

Output

For each line of input, output one line with the displacement of the center of the rod in millimeters with 3 digits of precision.

Sample Input

1000 100 0.0001
15000 10 0.00006
10 0 0.001
-1 -1 -1

Sample Output

61.329
225.020
0.000

思路转自：優YoU http://user.qzone.qq.com/289065406/blog/1301845324
大致题意：
一根两端固定在两面墙上的杆 受热弯曲后变弯曲
求前后两个状态的杆的中点位置的距离
 
解题思路：
几何和二分的混合体
 
 
 
 
如图，蓝色为杆弯曲前，长度为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 (不弯曲)
关键确定上界
题中提及到
Input data guarantee that no rod expands by more than one half of its original length.
意即输入的数据要保证没有一条杆能够延伸超过其初始长度的一半
就是说 S max = 3/2 L
理论上把上式代入（1）(2)方程组就能求到h的最小上界，但是实际操作很困难
因此这里可以做一个范围扩展，把h的上界扩展到 1/2L  ，不难证明这个值必定大于h的最小上界，那么h的范围就为  0<=h<1/2L
这样每次利用下界low和上界high就能得到中间值mid，寻找最优的mid使得(2)式左右两边差值在精度范围之内，那么这个mid就是h
 
精度问题是必须注意的
由于数据都是double，当low无限接近high时， 若二分查找的条件为while(low<high)，会很容易陷入死循环，或者在得到要求的精度前就输出了不理想的“最优mid”
精度的处理方法参考程序
#include<iostream>
#include<stdio.h>
#include<math.h>
using namespace std;
main()
{
    double pp = 1e-5;//精度限制
	double l,n,c,x;
	while(~scanf("%lf%lf%lf",&l,&n,&c))//注意下double的输入。。。
	{
		if(l<0||n<0||c<0)
			break;
		x=(1+n*c)*l;
        double high=l/2.0,low=0.0;
        double r,mid,tmp;
        while(high-low>pp)
       {
    	     mid = (heigh+low)/2;
            r = (l*l/4.0+mid*mid)/(2.0*mid);
            tmp = 2.0*r*asin(l/(2.0*r));
            if(tmp>x)
           {
        	    heigh = mid;
           }
            else
        	    low = mid;
        }
     printf("%.3f\n",mid);

	}
}