【三分+精度问题】G. Toxophily

https://www.bnuoj.com/v3/contest_show.php?cid=9154#problem/G

【题意】

已知人的坐标在(0,0),靶的位置在(x,y),人以速度v射箭并且射中靶，求v与x的夹角。

【思路】

• 经典的三分+二分题
• 先三分求出纵坐标最大时的sita，如果这时的纵坐标都小于y,就输出-1
• 三分确定sita后，在二分求出正解
• 这道题我的代码eps=1e-9才能过

【Accepted】

#include<iostream>
#include<cstdio> 
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>

using namespace std;
const double pi=acos(-1.0);
double x,y,v,sita;
const double eps=1e-8;
const double G=9.8;
double calc(double ang)
{
    double t=x/(v*cos(ang));
    return v*sin(ang)*t-0.5*G*t*t;
}
double triplediv()
{
    double l=0.0;
    double r=0.5*pi;
    while(r-l>eps)
    {
        double mid1=(2*l+r)/3.0;
        double mid2=(l+2*r)/3.0;
        if(calc(mid1)<calc(mid2))
        {
            l=mid1+eps; 
        } 
        else
        {
            r=mid1-eps;
        }
    }
    sita=l;
    return calc(l);
}

double doublediv()
{
    double l=0;double r=sita;
    while(r-l>eps)
    {
        double mid=(l+r)/2.0;
        if(calc(mid)<y)
        {
            l=mid+eps;
        }
        else
        {
            r=mid-eps;
        }
    }
    return l;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lf%lf%lf",&x,&y,&v);
        if(triplediv()<y)
        {
            puts("-1");
        }
        else
        {
            double ans=doublediv();
            printf("%.6f
",ans);
        }
    }
    return 0;
}