搜索—二分三分搜索—1004_三分搜索二分-优快云博客

本文介绍了一种结合三分法和二分法求解抛物线弹道问题的方法。通过三分法找到使抛物线在特定水平距离下达到最大高度的角度，再利用二分法精确确定能够覆盖指定垂直距离的发射角度。

这个题主要是一个公式的推导还有二分法和三分法的综合运用。。

思路：先三分出一个角度，让抛物线在x=x0出的y取得最大值。然后再到这个y一下去二分一个合适的值。

代码：

#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<math.h>
#define PI acos(-1.0)
#define G 9.8
#define eps 1e-12

using namespace std;

double X,Y,v,sita;

double cal(double x)
{
double t;
t=X/(v*cos(x));
return v*sin(x)*t-0.5*G*t*t;
}

double triplediv()
{
int i;
double mid1,mid2,left,right,h1,h2;
left=0,right=0.5*PI;
for(i=1;i<=100;i++)
{
mid1=(2*left+right)/3;
mid2=(left+2*right)/3;
h1=cal(mid1);
h2=cal(mid2);
if(h1>h2)
{
right=mid2;
}
else
{
left=mid1;
}
}
sita=left;
return cal(left);
}

void doublediv(double maxy)
{
int i;
double left,right,mid,h;
left=0,right=sita;
for(i=1;i<=100;i++)
{
mid=(left+right)/2;
h=cal(mid);
if(h>Y)
right=mid;
else
left=mid;
}
printf("%.6lf\n",left);
}

int main()
{
int t;
double maxy;
scanf("%d",&t);
while(t--)
{
scanf("%lf%lf%lf",&X,&Y,&v);
maxy=triplediv();
if(maxy<Y)
printf("-1\n");
else
doublediv(maxy);
}
return 0;
}