002：棍的膨胀

Messj

已知一根长为L的细棍被加热了n摄氏度，那么其新的长度为L'=(1+n*C)*L。中间的C是热膨胀系数。当一根细棍被夹在两面墙中间然后被加热，它会膨胀，其形状会变成一个弧，而原来的细棍(加热前的细棍)就是这个弧所对的弦。你的任务是计算出弧的中点与弦的中点的距离。

                               
登录/注册后可看大图

输入
包含多组数据(每组占一行)。每一行包含三个非负数：细棍的长度，温度的变化值和细棍材料的热膨胀系数C。输入数据保证细棍不会膨胀超过自己的一半。输入数据以三个连续的-1结尾。

输出
对于每一组数据，输出计算的距离，答案保留三位小数。

样例输入
1000 100 0.0001
15000 10 0.00006
10 0 0.001
-1 -1 -1

样例输出
61.329
225.020
0.000

Messj

解析：https://www.cnblogs.com/SinGuLaRiTy2001/p/7171242.html#_label0

之所以可以使用分治策略中的二分渐进法，是由于偏移的距离与角度成正比，有单调性。

1. #include<cstdio>
   
2. #include<cmath>
   
3. 
   
4. using namespace std;
   
5. #define eps 1e-4
   
6. 
   
7. double l,ll,n,c,r,sita;//l-棍长 n-温度 c-膨胀常量 ll-弧长 r-半径 sita-对应角
   
8. 
   
9. double solve(double a,double b) //二分渐进法
   
10. {
    
11.         double left,right,mid;
    
12.         left=a;
    
13.         right=b;
    
14.         while(left+eps<right)
    
15.         {
    
16.                 mid=(left+right)/2;
    
17.                 r = l*l/8/mid+mid/2;
    
18.                 sita = 2*asin(l/2/r);
    
19.                 if(sita*r>ll)  // ll=sita*r  l=sin(sita)*r
    
20.                         right=mid;
    
21.                 else
    
22.                         left=mid;
    
23.         }
    
24.         return left;
    
25. }
    
26. int main()
    
27. {
    
28.         while(scanf("%lf%lf%lf",&l,&n,&c))
    
29.         {
    
30.                 if(l<0||n<0||c<0)break;
    
31.                 ll=l*(1+n*c);
    
32.                 printf("%.3f",solve(0,l/2));
    
33.         }
    
34. }

复制代码