1. 程式人生 > >[洛谷P3382]【模板】三分法

[洛谷P3382]【模板】三分法

esp ace print 註意 兩個 分法 ont define nbsp

題目大意:給出一個N次函數,保證在範圍[l,r]內存在一點x,使得[l,x]上單調增,[x,r]上單調減。試求出x的值。

解題思路:三分法。像我這種什麽函數都不知道的,只知道要三分。取兩個“mid”,然後分別求出答案。如果ansl<ansr,則l=midl,否則r=midr。

基本和二分一樣。

註意eps選得小一點,我用$10^{-14}$。

C++ Code:

#include<cstdio>
using namespace std;
#define eps 1e-14
int n;
double l,r,a[12222];
inline double calc(double x){
	double p=0,now=1;
	for(int i=0;i<=n;++i){
		p+=now*a[i];
		now*=x;
	}
	return p;
}
int main(){
	scanf("%d%lf%lf",&n,&l,&r);
	for(int i=n;i>=0;--i)scanf("%lf",&a[i]);
	while(l+eps<r){
		double midl=(r-l)/3+l,midr=(r-l)/3*2+l;
		double ansl=calc(midl),ansr=calc(midr);
		if(ansl+eps<ansr)l=midl;else
		r=midr;
	}
	printf("%.5f\n",l);
}

[洛谷P3382]【模板】三分法