Toxophily HDU

代码+解析：

//题意：
//有一个大炮在(0,0)位置，为你可不可以把炮弹射到(x,y)这个位置
//题目给你炮弹初始速度，让你求能不能找出来一个炮弹射出时角度满足题意
//题解：
//由物理公式分析可知：
//Vx=v*cos(a)
//Vy=v*sin(a)
//t=x/Vx=x/(v*cos(a))
//y=-(1/2)*g*t*t+Vy*t=-(1/2)*g*t*t+v*sin(a)*t
//一看是一个一元二次函数那他的图像不是先减后增就是先增后减。那就用三分找出来它的最高位置对应度数是多少(假设是ans)
//那么(0,ans)与(ans,90)这两个区间所对应曲线肯定一致，那我们就挑一个区间进行二分查找答案就完了
#include<stdio.h>
#include<stdlib.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<math.h>
using namespace std;
const int maxn=500005;
const int INF=0x3f3f3f3f;
#define eps 1e-11
#define PI acos(-1.0)  //180度的弧度值
typedef long long ll;
double x,y,v;
double f(double a)  //求大炮角度为a时，大炮能射多高
{
    double t = x/(v*cos(a));
    return v*sin(a)*t - 9.8/2.0*t*t;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lf%lf%lf",&x,&y,&v);
        double l=0,r=PI/2,m1,m2;   //先找出来大炮的能射最高所对应的角度，只需要在(0,90)度范围内找就行了
        while(r-l>eps)
        {
            m1=(l+r)/2;
            m2=(m1+r)/2;
            if(f(m1)<f(m2))
                l=m1;
            else r=m2;
        }
        if(f(l)<y)  //如果大炮最高也不到y，那就根本无法打到那个位置
        {
            printf("-1
");
            continue;
        }
        double mid,ll=0,rr=l;
        while(rr-ll>eps)  //找到之后在二分出来那个角度能刚好射到高度为y
        {
            mid=(ll+rr)/2;
            if(f(mid)<y)
                ll=mid;
            else rr=mid;
        }
        printf("%.6lf
",ll);
    }
    return 0;
}