欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

洛谷 P3382 【模板】三分法

程序员文章站 2022-07-13 13:55:16
...

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

输入格式
第一行一次包含一个正整数N和两个实数l、r,含义如题目描述所示。

第二行包含N+1个实数,从高到低依次表示该N次函数各项的系数。

输出格式
输出为一行,包含一个实数,即为x的值。四舍五入保留5位小数。

输入输出样例
输入 #1复制
3 -0.9981 0.5
1 -3 -3 1
输出 #1复制
-0.41421
说明/提示
时空限制:50ms,128M

数据规模:

对于100%的数据:7<=N<=13

样例说明:

洛谷 P3382 【模板】三分法

如图所示,红色段即为该函数f(x)=x3-3x2-3x+1在区间[-0.9981,0.5]上的图像。

当x=-0.41421时图像位于最高点,故此时函数在[l,x]上单调增,[x,r]上单调减,故x=-0.41421,输出-0.41421。

(Tip.l&r的范围并不是非常大ww不会超过一位数)

解法:三分

就是个模板,也没什么可以说的

AC代码

#include<cstdio>
#define re register int
using namespace std;
int n; double L,R,a[15];
inline double check(double x) {
	double sum=0;
	for(re i=n;i>=0;i--) {
		sum=sum*x+a[i];
	}
	return sum;
}
int main() {
	scanf("%d%lf%lf",&n,&L,&R);
	for(re i=n;i>=0;i--) {
		scanf("%lf",&a[i]);
	}
	while(R-L>=1e-7) {
		double mid=(L+R)/2;
		if(check(mid+1e-7)>check(mid-1e-7)) L=mid;
		else R=mid;
	}
	printf("%.5lf",L);
	return 0;
}