POJ Expanding Rods <简单二分>

题目链接

题目大意：两墙之前的细杆在加热之后会膨胀，若增加的温度为$n$，膨胀系数为$C$，那么原先长为$L$的细杆膨胀后长度为$L'=(1+n*C)*L$，于是细杆会被挤压成弧形，求细杆前后中心点的高度差是多少。

分析：可以对高度$h$使用二分，由于题目中有一句

Input data guarantee that no rod expands by more than one half of its
original length

所以可以将$h$上限放大到$L/2$，先求圆弧对应的半径$r=\frac{h^2+L^2/4}{2h}$，然后得到半个圆心角$\theta=arcsin(\frac{L}{2r})$，这样就可以求得当前$h$下对应的弧长$L_{now}=2r\theta$，与实际弧长相比就知道 h <script type="math/tex" id="MathJax-Element-1716">h</script>是大了还是小了。

代码：

#include <iostream>
#include <cmath>
#include <algorithm>
#include <iomanip>
using namespace std;

const double eps=1e-5;

double l,n,c,L;

bool judge(double h){
    double r=(h*h+l*l/4.0)/(2.0*h);
    return asin(l/(2.0*r))*r<L/2.0;
}

int main(){
    while(cin>>l>>n>>c,l!=-1&&n!=-1&&c!=-1){
        L=(1+n*c)*l;
        double lb=0,rb=l/2,mid;
        while(rb-lb>eps){
            mid=(lb+rb)/2;
            if(judge(mid)) lb=mid;
            else rb=mid;
        }
        cout<<fixed<<setprecision(3)<<mid<<endl;
    }
}