hdu 5144

三分法加物理公式……

推出来了公式，但是发现不会求解……QAQ ，以为要求导……然后发现不会求导……

原来是传说中的三分法！

然后……就出来了……

三分法就是在二分的情况下再分一次。

#include <stdio.h>
#include <math.h>
#define eps 1e-12
double v,g,h;
double pi;

int sig(double x)
{
    return (x>eps)-(x<-eps);
}
double f(double ang)
{
    double jump=2*v*sin(ang)*v*cos(ang)/g;
    double fall=(sqrt(2*g*h+v*v*sin(ang)*sin(ang))-v*sin(ang))*v*cos(ang)/g;
    return jump+fall;
}
int main()
{
    int T;
    g=9.8;
    pi=3.14145926;
    scanf("%d",&T);
    while(T--)
    {
        int i,j,k;
        scanf("%lf%lf",&h,&v);
        double l,r;
        l=0.0,r=pi/2;
        while(sig(r-l)>0)
        {
            double mid=(l+r)/2;
            double mmid=(mid+r)/2;
            double t1=f(mid);
            double t2=f(mmid);
            if(t1>t2) r=mmid;
            else l=mid;
        }
        printf("%.2lf\n",f(l));
    }
    return 0;
}