HDU 2298 Toxophily

题目地址： http://acm.hdu.edu.cn/showproblem.php?pid=2298

 

x=v*cosα*t;

y=v*sinα*t-1/2*g*t^2;

化简消去t得 ：y=x*tanα-x^2*g/(2*v^2*cosα*cosα).

三分+二分

可知该方程先递增后递减，所以先三分求出最大值，如果y小于最大值，则输出-1，否则二分求出α值。

 

代码如下：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
using namespace std;

/*
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
*/

double x,y,v;
const double PI=acos(-1.0),g=9.8;

double f(double a)
{
    return x*tan(a)-x*x*g/(2*v*v*cos(a)*cos(a));
}

double sanfen(double l,double r)//三分查找
{
    double mid,midmid;
    while((r-l)>1e-9)//这里注意精度，1e-8就WA了
    {
        mid=(l+r)/2;
        midmid=(mid+r)/2;
        if(f(mid)>f(midmid))
            r=midmid;
        else
            l=mid;
    }
    return (mid+midmid)/2;
}

double erfen(double l,double r)//二分查找
{
    double mid;
    while((r-l)>1e-9)
    {
        mid=(l+r)/2;
        if(f(mid)>y)
            r=mid;
        else
            l=mid;
    }
    return mid;
}

int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        scanf("%lf%lf%lf",&x,&y,&v);
        double max=sanfen(0,PI/2);
        if(f(max)<y)
        {
            printf("-1\n");
        }
        else
        {
            printf("%.6f\n",erfen(0,max));
        }
    }
    return 520;
}

 

 另外一种直接用数学方法求解。

y=x*tanα-x^2*g/(2*v^2*cosα*cosα)

三角函数换元得

.x^2*g/(2*v^2)*tanα-x*tanα+x^2*g/(2*v^2)+y=0;

a=x^2*g/(2*v^2);

b=-x;

c=x^2*g/(2*v^2)+y;

通过判断deta和α范围确定答案。

 

代码如下：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
using namespace std;

/*
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
*/

double x,y,v;
const double PI=acos(-1.0),g=9.8;

double f(double a)
{
    return x*tan(a)-x*x*g/(2*v*v*cos(a)*cos(a));
}

int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        scanf("%lf%lf%lf",&x,&y,&v);
        double a,b,c,deta,a1,a2;
        a=x*x*g/(2*v*v);
        b=-x;
        c=a+y;
        deta=b*b-4*a*c;
        if(deta<0)
            printf("-1\n");
        else
        {
            a1=atan((-b+sqrt(deta))/2/a);
            a2=atan((-b-sqrt(deta))/2/a);
            if((a1>PI/2||a1<0)&&(a2>PI/2||a2<0))
                printf("-1\n");
            else if((a1>PI/2||a1<0))
                printf("%.6f\n",a2);
            else if(a2>PI/2||a2<0)
                printf("%.6f\n",a1);
            else
                printf("%.6f\n",a1<a2?a1:a2);
        }
    }
    return 520;
}