POJ Expanding Rods

几何问题二分法求解

最新推荐文章于 2017-12-11 13:28:12 发布

原创最新推荐文章于 2017-12-11 13:28:12 发布 · 400 阅读

0 ·

CC 4.0 BY-SA版权

数学同时被 2 个专栏收录

9 篇文章

订阅专栏

分治

3 篇文章

订阅专栏

本文介绍了一种使用二分法解决特定几何问题的方法，通过给定的红色线段L和绿色弧线L(1+n*C)，利用正弦定理和二分查找算法求解两线中点之间的距离。

点击打开题目

题目大意

这里写图片描述

给定L，n，C，L为红色线段，L(1+n*C)为绿色弧，求两者中点的距离

二分圆心角度数，接下来就是几何的能力了

根据正弦定理，可得：
$\frac{L}{sin\theta}=\frac{r}{sin(90°-\theta)}$

则弧长：
$a=\pi r·\frac{\theta}{180}$

将a与nL作比较来二分

精度满天飞 QWQ

代码如下：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
double pi=3.1415926535897932384626433832795,eps=0.000000001;
double L,C,n,nL;
using namespace std;
double work(double mid)
{
    double a1=sin(2*pi*mid/360),a2=sin(2*pi*(90-mid/2)/360);
    double r=L/a1*a2;
    return 2*r*pi*mid/360;
}
int main() 
{
    while(scanf("%lf%lf%lf",&L,&n,&C))
    {
        if(L==-1&&n==-1&&C==-1)return 0;
        nL=(1+n*C)*L;
        if(nL==L){printf("0.000\n");continue;}
        double l=0,r=180,mid;
        while(r-l>eps)
        {
            mid=(l+r)/2;
            if(work(mid)<nL)l=mid;
            else r=mid;
        }
        printf("%.3lf\n",L/sin(2*pi*l/360)*sin(2*pi*(90-mid/2)/360)-sqrt(pow(L/sin(2*pi*l/360)*sin(2*pi*(90-mid/2)/360),2)-pow(L/2,2)));
    }
}