[LuoguP1883]函数 三分法

Posted 2021-01-13 zhinv

三分法的模板，一直提交，一直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
",F(l));
    }
    return 0;
}