1004-搜索初步

将射箭的方向正交分解为x轴和y轴。

 则x=v*t*cosa----->(1)  y=v*t*sina-(g*t^2)/2----->（2）

消去t得到：y=x*tana-(g/2)*[(x^2)/v^2*cosa*cosa]----->(3)

1/(cosa)^2=(sina)^2+(cosa)^2/(cosa)^2=1+(tana)^2-----(4)

将（4）带入（3）得出公式。

g*x*x*(tana)^2-2*v*v*x*tana+g*x*x+2*v*v*y=0

一元二次方程，得出a，b，c；

#include<iostream>
#include<cmath>
#include<iomanip>
#include<cstdio>
using namespace std;
#define PI acos(-1)
#define g 9.8
int main()
{

    int T;
    cin>>T;
    while(T--)
    {
        double x,y,v;//x，y描述物品的坐标。v表示箭的速度。
        cin>>x>>y>>v;
        //在推导出的方程当中，a=g*x^2,b=-2*x*v^2,c=g*x*x+2*v*v*y;
        double a,b,c,dyx,wyx,temp;
        int wrong=-1;
        a=g*x*x;
        b=(-2)*x*v*v;
        c=g*x*x+2*v*v*y;
        temp=b*b-4*a*c;
        if(temp<0)
        cout<<wrong<<endl;
        else
        {
            dyx=atan((-1*b+sqrt(temp))/(2*a));
            wyx=atan((-1*b-sqrt(temp))/(2*a));
            if((dyx<0||dyx>PI/2.0)&&(wyx<0||wyx>PI/2.0))
            cout<<wrong<<endl;
            else if(dyx<0||dyx>PI/2.0)
            cout<<fixed<<setprecision(6)<<wyx<<endl;
            else if(wyx<0||wyx>PI/2.0)
            cout<<fixed<<setprecision(6)<<dyx<<endl;
            else
            cout<<fixed<<setprecision(6)<<(wyx>dyx?dyx:wyx)<<endl;
        }
    }
    return 0;
}