la 5009 - Error Curves(三分)

la 5009 - Error Curves

单峰函数用三分

#include<stdio.h>
#include<string.h>
#include<math.h>
#define MAXN 10005
const double eps=1e-9;
int a[MAXN],b[MAXN],c[MAXN];
int n;
double max(double a,double b) {return a>b?a:b;}
double f(double x)
{
	double ans=a[0]*x*x+b[0]*x+c[0];
	for(int i=1;i<n;i++)
	{
		ans=max(ans,a[i]*x*x+b[i]*x+c[i]);
	}
	return ans;
}
int main()
{
	int cas;
	scanf("%d",&cas);
	while(cas--)
	{
		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;
		while(l+eps<r)
		{
			double m1=(l+r)/2;
			double m2=(m1+r)/2;
			if(f(m1)+eps<f(m2)) r=m2;
			else l=m1;
		}
		printf("%.4lf\n",f(l));
	}
	return 0;
}