http://acm.hdu.edu.cn/showproblem.php?pid=2899
三分法
我的代码
1 #include <stdio.h>
2 #include <math.h>
3 const double eps=1e-8;
4 double y;
5 double f(double x)
6 {
7 return 6*pow(x,7)+8*pow(x,6)+7*x*x*x+5*x*x-y*x;
8 }
9 int main()
10 {
11 int T;
12 double l,r,ml,mr;
13 scanf("%d",&T);
14 while (T--)
15 {
16 scanf("%lf",&y);
17 l=0; r=100;
18 while (fabs(f(l)-f(r))>eps)
19 {
20 ml=l+(r-l)/3;
21 mr=l+2*(r-l)/3;
22 if (f(ml)<f(mr)) r=mr;
23 else l=ml;
24 }
25 printf("%.4lf\n",f(l));
26 }
27 return 0;
28 }
