【洛谷 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 }