题目来源: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;
}