hdu2298（二分+三分）

hdu-2298

 

题意：

一个人站在（0,0）向（x,y）射箭，x代表水平距离，y代表竖直距离，g=9.8。

求出射箭的角度，如果到达不了，就输出-1.

 

思路：

（1）总体思路

先求出箭头可以到达的最大高度，然后求出最大高度时所需的角度，

如果最大的高度无法到达，就输出-1.

然后在这个角度范围内查找合适的高度。

（2）计算公式

竖直方向：H = Vx*t + 1/2*g*t*t;  Vx =  V * sin(du);

水平方向：X = Vy*t; Vy = V * cos(du);

 

复习：

利用y=Vy*t - 1/2 *g *t*t; x = vx*t;消去t，得到公式：y = x*tanu - 1/2*g*x*x*(1+tanu*tan*u)/(v*v)。（u表示发射方向与x的角度）因为y对tanu是二次函数，不能直接用二分法求解，所以可以先用三分法求出y的最大值时对应的角度u`,然后再[0,u']区间内用

二分法求解对应的u``。注意：考虑90度的情况。

 

参考文章

 

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
const double eps = 1e-11;
double x,y,sp;
double f(double du) //返回高度的函数 
{
	double t=x/(sp*cos(du)); //求出时间t，是关于t的二次函数。 
	return sp*sin(du)*t-4.9*t*t;
}
int main(void)
{
	int T,i;
	scanf("%d",&T);
	while(T--){
		scanf("%lf%lf%lf",&x,&y,&sp);
		double l=0,r=acos(-1)/2-eps; //l,r分别表示第一象限的角度范围 ，不能是90度，开区间 
		double m1,m2;
		if(x==0){ //90度的情况 
			if(sp*sp/2/9.8>y) printf("%.6lf\n",r); //如果最大距离超过y 
			else printf("-1\n");
		}
		else{
			for(i=1;i<=100;i++){ //三分100次 ，求出最大高度时所用的角度 
				m1=(l+r)/2;
				m2=(r+m1)/2;
				if(f(m1)>f(m2)) r=m2;
				else l=m1;
			}
			if(f(l)<y){
				printf("-1\n");continue; //最大高度无法到达 
			}
			r=l;l=0;
			for(i=1;i<=100;i++){ //在所需的高度中查找所需的角度，100次二分，尽量精确 
				m1=(l+r)/2;
				if(f(m1)<y) l=m1;
				else r=m1;
			}
			printf("%.6lf\n",l);
		}
	}
	return 0;
}