【洛谷 p3382】模板-三分法(算法效率)

题目:给出一个N次函数,保证在范围[l,r]内存在一点x,使得[l,x]上单调增,[x,r]上单调减。试求出x的值。

解法:与二分法枚举中点使区间分成2份不一样,三分法是枚举三分点,再根据题目的情况修改 l , r 。

P.S.嘻嘻,其实我就是为了凑够我这3个月博客刚好150篇才补了这个我NOIP比赛临出发打的题。

 1 #include<cstdio>
 2 #include<cstdlib>
 3 #include<cstring>
 4 #include<iostream>
 5 using namespace std;
 6 #define eps 1e-8
 7 
 8 int n;
 9 double a[15];
10 
11 double mabs(double x) {return x>0?x:-x;}
12 double f(double x)
13 {
14     double h=0.0,tmp=1.0;
15     for (int i=0;i<=n;i++)
16     {
17       h+=tmp*a[i];
18       tmp*=x;
19     }
20     return h;
21 }
22 int main()
23 {
24     double l,r;
25     scanf("%d%lf%lf",&n,&l,&r);
26     for (int i=0;i<=n;i++)
27       scanf("%lf",&a[n-i]);
28     while (mabs(r-l)>eps)
29     {
30       double ll=l+(r-l)/3,rr=l+(r-l)/3*2;
31       if (f(ll)>f(rr)) r=rr;//
32       else l=ll;
33     }
34     printf("%.5lf",(double)l);
35     return 0;
36 }

 

posted @ 2016-11-21 18:31  konjac蒟蒻  阅读(366)  评论(0编辑  收藏  举报