HDU 5144 NPY and shot(三分法)

 1 #include <stdio.h>
 2 #include <math.h>
 3 #define PI 3.1415926
 4 int v, h;
 5 double f(double i)//推倒物理公式 
 6 {
 7     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);
 8 }
 9 int main()
10 {
11     int T;
12     scanf("%d", &T);
13     while (T--)
14     {
15         scanf("%d %d", &h, &v);
16         double ans = 0.0;
17         double l = 0.0, r = PI/2;
18         double midl, midr;
19         while (r - l >= 1e-6)//三分求最值 
20         {
21             midl=(2*l+r)/3;//是由midl = (r - l) / 3 + l化简得来的 
22             midr=(l+2*r)/3;//是由midr = (r - l) * 2 / 3 + l化简得来的 
23             if(f(midl)>f(midr))//画图像比较好看这个关系 
24                 r=midr;
25             else
26                 l=midl;
27         }
28         printf("%.2lf\n", f((r+l)/2));
29     }
30     
31     return 0;
32 }