Luogu 3382(三分)

传送门

模板题。好想有个什么叫黄金分割比例优化三分,洛谷的题解里有一位大佬提到过orz。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=15;
const double eps=1e-10;
double a[N],L,R;
int n;
inline double fpow(double a,int b) {
	double ret=1.0;
	while (b) {
		if (b&1) ret=ret*a;
		a*=a,b>>=1;
	}
	return ret;
}
inline double cal(double x) {
	double ret=0.0;
	for (int i=n;~i;--i)
		ret+=a[i]*fpow(x,i);
	return ret;
}
int main() {
//	freopen("P3382.in","r",stdin);
	scanf("%d%lf%lf",&n,&L,&R);
	for (int i=n;~i;--i) scanf("%lf",&a[i]);
	while (R-L>eps) {
		double lm=L+(R-L)/3,rm=R-(R-L)/3;
		if (cal(lm)<cal(rm)) L=lm;
		else R=rm;
	}
	printf("%.5lf\n",L);
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值