推公式hdu2298

(0,0)点到(x,y)以v射箭 问 最小的角度

不能射到 -1

x=v*cos(a)*t;

y=v*sin(a)*t-0.5*g*t*t;

消去t;

y=x*tan(a)-0.5*g*x*x/(v*v*cos(a)*cos(a));

v*v*y=v*v*x*tan(a)-0.5*g*x*x*sec(a)*sec(a);

根据 1+tan(a)*tan(a)=sec(a)*sec(a);

化简得到 g*x*x*tan(a)*tan(a)-2*x*v*v*tan(a)+2*v*v*y+g*x*x=0;

然后二元一次方程

几个特判一下

a=g*x*x;   a!=0     所以x要判

然后b*b-4*a*c;

然后解中找小的但是>=0

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>

using namespace std;
typedef long long ll;
#define MAXN 100010
#define g   9.8
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        double x,y,v,a,b,c,d,x1,x2;
        scanf("%lf%lf%lf",&x,&y,&v);
        if(x==)
        {
            if(y==)
                printf("%.6lf\n",);
            else
            {
                double s1=v*v/(*g);
                if(s1>=y)
                    printf("%.6lf\n",);
                else
                    printf("-1\n");
            }
            continue;
        }
        a=g*x*x;
        b=(-)*x*v*v;
        c=*v*v*y+g*x*x;
        if(b*b-*a*c<)
        {
            printf("-1\n");
            continue;
        }
        d=sqrt(b*b-*a*c);
        x1=(-b-d)/(*a);
        x2=(-b+d)/(*a);
        if(x1>=)
            printf("%.6lf\n",atan(x1));
        else if(x2>=)
            printf("%.6lf\n",atan(x2));
        else
            printf("-1\n");
    }
    return  ;
}