搜索—二分三分搜索—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;
}