Error Curves UVALive - 5009 （三分法求单峰函数的极值）

传送门

#include<stdio.h>
#include<iostream>
#include<sstream>
#include<cmath>
#include<math.h>
#include<algorithm>
#include<string.h>
#include<string>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include <ctype.h>
#define INF 0x3f3f3f3f
#define NINF 0xc0c0c0c0
using namespace std;
const int N = 10010;
int T,n,a[N],b[N],c[N];

double f(double x){
    double res = a[0]*x*x+b[0]*x+c[0];
    for(int i=1;i<n;i++){
        res = max(res,a[i]*x*x+b[i]*x+c[i]);
    }
    return res;
}

int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i=0;i<n;i++){
            scanf("%d%d%d",&a[i],&b[i],&c[i]);
        }
        double l = 0, r = 1000;
        for(int i=0;i<100;i++){
            double m1 = l+(r-l)/3;
            double m2 = r-(r-l)/3;
            if(f(m1)<f(m2)){
                r = m2;
            } else{
                l = m1;
            }
        }
        printf("%.4f\n",f(l));
    }
    return 0;
}