HDU 5144 NPY and shot(三分法)

当时做这道题时一直想退出物理公式来，但是后来推到导数那一部分，由于数学不好，没有推出来那个关于Θ的最值，后来直接暴力了，很明显超时了，忘了三分法的应用，这道题又是典型的三分求最值，是个单峰曲线，下面是代码

#include <stdio.h>
#include <math.h>
#define PI 3.1415926
int v, h;
double f(double i)//推倒物理公式 
{
    return v*v*1.0*sin(2*i)/9.8+(sqrt(v*v*1.0*sin(i)*sin(i)+2.0*9.8*h)-v*sin(i))/9.8 * v*cos(i);
}
int main()
{
    int T;
    scanf("%d", &T);
    while (T--)
    {
        scanf("%d %d", &h, &v);
        double ans = 0.0;
        double l = 0.0, r = PI/2;
        double midl, midr;
        while (r - l >= 1e-6)//三分求最值 
        {
            midl=(2*l+r)/3;//是由midl = (r - l) / 3 + l化简得来的 
            midr=(l+2*r)/3;//是由midr = (r - l) * 2 / 3 + l化简得来的 
            if(f(midl)>f(midr))//画图像比较好看这个关系 
                r=midr;
            else
                l=midl;
        }
        printf("%.2lf\n", f((r+l)/2));
    }
    
    return 0;
}