2019_GDUT_新生专题I选集 三分

K
题意：函数F（x）定义如下：F（x）= max（Si（x）），i = 1 … n。 x的定义域是[0，1000]。 Si（x）是二次函数。求F（x）的最小值。
t组测试数据 t(t<=100) n (n ≤ 10000)
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3714

显然每个二次函数都是一个凸函数，在取最大值后拼接起来的F显然也是一个凸函数，此时可以用三分法来查找函数的最小值。
三分法说简单点就是在二分的基础上再多分一个点出来比较。
取区间[l,r]的两个三等分点
lmid=(left*2+right)/3 rmid=(left+2*right)/3；
显然
如果 lmid 和 rmid 在最大（小）值的同一侧： 那么由于单调性，一定是二者中较大（小）的那个离最值近一些，较远的那个点对应的区间不可能包含最值，所以可以舍弃。
如果在两侧： 由于最值在二者中间，我们舍弃两侧的一个区间后，也不会影响最值，所以可以舍弃。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<iostream>
using namespace std;
int n;
double a[10005],b[10005],c[10005];
double func(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;
	scanf("%d",&t);
	while (t--)
	{
		scanf("%d",&n);
		for (int i=1;i<=n;i++)
			scanf("%lf %lf %lf",&a[i],&b[i],&c[i]);
		double l=0,r=1000;
		while (r-l>=1e-9)
		{
			double p=(r+2*l)/3,q=(2*r+l)/3;
			if (func(p)>=func(q)) l=p;
			else r=q;
		}
		printf("%.4lf\n",func(l));
	}
	return 0;
}