LA5009三分法

/*LA5009:
定义 fi(x)=a[i]*x^2+b[i]*x+c[i](a[i]>=0),F(x)=max(fi(x)),0<=x<=1000
求F(x)在区间上的最小值
三分法模板题：总结一下：
1、三分法适用类型，a:凹函数的最大值的最小值，b:凸函数的最小值的最大值。
必须是单峰函数啊
2、特点：升维，记住三分法不是用来解零点的
3、解法:a类型，取三分点m1，m2，m1<m2,若F(m1)<F(m2) 则 解在(l,m2)上 ，否则在（m1，r)上
        b类型，解在靠近大的值得区间上
*/
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
#define eps 1e-10
#define LL long long
using namespace std;
int n,T;
double a[10000+5],b[10000+5],c[10000+5];
double f(int i,double x)
{
    return a[i]*x*x+b[i]*x+c[i];
}

double F(double x)
{
    double m=f(0,x);
    for(int i=1;i<n;i++) m=max(m,f(i,x));
    return m;
}
int main()
{
    cin>>T;
    for(int cas=1;cas<=T;cas++)
    {
        cin>>n;
        for(int i=0;i<n;i++) cin>>a[i]>>b[i]>>c[i];
        double l=0,r=1000;
        while(r-l>eps)
        {
            double m1=l+(r-l)/3,m2=m1+(r-l)/3;
            if (F(m1)<F(m2)) r=m2;else l=m1;
        }
        printf("%.4lf
",F(l));//注意最后输出的是函数值
    }
    return 0;
}