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(*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/;
         double midl, midr;
         while (r - l >= 1e-)//三分求最值
         {
             midl=(*l+r)/;//是由midl = (r - l) / 3 + l化简得来的
             midr=(l+*r)/;//是由midr = (r - l) * 2 / 3 + l化简得来的
             if(f(midl)>f(midr))//画图像比较好看这个关系
                 r=midr;
             else
                 l=midl;
         }
         printf("%.2lf\n", f((r+l)/));
     }

     return ;
 }