hdu2298 Toxophily

题目来源:Problem - 2298 (hdu.edu.cn)

题目大意:从圆心以速度v抛出某物,求通过(x0,y0)的最小角度。

解题思路:

如图所示,根据题意,不妨设x0=2,则从原点(0,0)抛出物体时,

当θ=0时,是平抛运动(黑线),则抛物线与x=2的交点的y值小于0;

当θ=Π/2时,是上抛运动(蓝线),则抛物线与x=2的交点的y值也必然小于0;

当0<θ<Π/2时,是抛物运动(红线),与x=2的交点的y值可能大于0,也就是存在解的可能。

当θ从0逐渐变化到Π/2时,交点y的值是先变大再变小,如图所示:

不妨设最大值ymax,对应的角度为θmax。

如果ymax>=y0,则存在解的可能(若不等于,则可能存在两个解。分别在[0,θmax] 和[θmax,Π/2]之间,显然答案在[0,θmax]。

因此,算法可以分为两步:

1)先在[0,Π/2]内找到θmax。这个是从某个范围找最大值,且线具有凸性,可以使用三分法。

2)再在[0,θmax]内找到θ,使得经过(x0,y0)点。这个单调线上查找某个根,因此可以使用二分法。

参考代码如下:

#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;

#define G 9.8
#define PI 3.1415926
#define E 1e-8

double x,y,v;
//返回:以xita角度抛出时,与x=x0垂直线的交点的y值
double f(double xita){
    return x*tan(xita)-1/2.0*G*x*x/v/v/cos(xita)/cos(xita);
}
double triPartition(double (*pf)(double xita),double l,double r,double e);
double biPartition(double (*pf)(double x),double l,double r,double e);
int main()
{
    int count;
    cin>>count;
    while(count--){
       cin>>x>>y>>v;
       double argMax = triPartition(f,0,PI/2,E);
       if (f(argMax)<y){
           cout<<-1<<endl;
       }
       else{
          double arg = biPartition(f,0,argMax,E);
          cout<<fixed<<setprecision(6)<<arg<<endl;
       }
    }
    return 0;
}

//三分法
//pf 函数曲线,凸性,向上凸。...
//l表示左边界,r表示有边界,e表示精度
double triPartition(double (*pf)(double xita),double l,double r,double e){
    while(abs(r-l)>e){
        double x1= l/3*2+r/3;
        double x2 = l/3 +r/3*2;
        if (pf(x1)<pf(x2))
            l= x1;
        else //if (pf(x1)>pf(x2))
            r = x2;
    }
    return r;
}

//二分法:
//pf 函数曲线,凸性,向上凸。...
//l表示左边界,r表示有边界,e表示精度
double biPartition(double (*pf)(double x),double l,double r,double e){
    while(abs(r-l)>=e){
        double m= (l+r)/2;
        if (pf(m)>y)
            r = m;
        else
            l = m;
    }
    return r;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值