[洛谷 1883]函数 三分法

三分法的模板,一直提交,一直WA

结果,今天把精度判断从l-r>=t改成 fabs(F(l)-F(r))>=t ,过了三个点,再将 \(t\) 改成 \(10^{-7}\)\(AC\) 了。。。

这给我一个十分重要的启示:精度问题十分重要!!!

#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;
const int MAXN = 10000+5;
int a[MAXN],b[MAXN],c[MAXN],n;

inline double F(double x) {
	double ans=a[1]*x*x+b[1]*x+c[1];
	for(int i=2; i<=n; i++) ans = max(ans ,a[i]*x*x+b[i]*x+c[i]);
	return ans;
}

int main() {
	int T;
	cin >> T;
	while(T--) {

		cin >> n;
		for(int i=1; i<=n; i++) cin >> a[i] >> b[i] >> c[i];

		double l=0,r=1000,t=1e-7;
		while(fabs(F(l)-F(r))>t) {
			double p=(2*l+r)/3,q=(l+2*r)/3;
			if(F(p) > F(q)) {
				l=p;
			} else r=q;
		}

		printf("%.4lf\n",F(l));
	}
	return 0;
}
posted @ 2018-10-19 15:27  昤昽  阅读(201)  评论(0编辑  收藏  举报