hdu 2298(推公式、解方程)

/*
    hdu 2298
    题目大意就是
    你站在(0,0)点射箭要射到(x,y)点
    告诉你箭的速度v
    问通过调整角度能不能射到(x,y)点
    如果可以问最小角度是多少
    如果不行输出-1
    这就是一道高中物理题
    当箭刚好到(x,y)点那么所需的角度最小
    设与x方向的夹角为α
    则
    y方向有y=Vy*t-g*t^2*1/2
    x方向有x=Vx*t
    将速度正交分解得
    Vx=V*cos(α)
    Vy=V*sin(α)
    然后化简成为tan(α)的一元二次方程,得
    g*x^2*(tan(α))^2-2*x*v^2+(2*v^2*y+g*x^2)=0
    (这题好像不能用反向的平抛运动来求,不知道是我求错了还是出题者根本没考虑)


    最后解方程求反正切，得到两个角
    输出
    范围在atan(y/x)~pi/2(至于为什么是这个范围，应该挺好想的，角度要符合事实)之间且小的那个角就行
*/
#include <iostream>
#include <algorithm>
#include <stdio.h>
#include <math.h>
#include <string.h>
#define mod 1000000007
#define MAX 100005
#define ll long long
#define PI acos(-1)
#define g 9.8
using namespace std;


int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        double x,y,v;
        cin>>x>>y>>v;
        double a,b,c;
        a=g*x*x;
        b=-2*v*v*x;
        c=2*v*v*y+g*x*x;
        double d=b*b-4*a*c;
        if(d<0)
        {
            printf("-1\n");
            continue;
        }
        double res1=(-b+sqrt(d))/(2*a);
        double res2=(-b-sqrt(d))/(2*a);
        double ans1=atan(res1);
        double ans2=atan(res2);
        double angle=atan(y/x);
        if(ans1>angle&&ans1<PI/2&&ans2>angle&&ans2<PI/2)
        {
            if(ans1>ans2)
                printf("%.6lf\n",ans2);
            else
                printf("%.6lf\n",ans1);
        }
        else if(ans1>angle&&ans1<PI/2&&ans2<=angle&&ans2>=PI/2)
        {
            printf("%.6lf\n",ans1);
        }
        else
        {
            printf("%.6lf\n",ans2);
        }


    }
    return 0;
}