B - 曲线

B - 曲线

明明做作业的时候遇到了 nn 个二次函数 S_i(x)= ax^2 + bx + cSi​(x)=ax2+bx+c，他突发奇想设计了一个新的函数 F(x) = \max{S_i(x)}, i = 1\ldots nF(x)=max{Si​(x)},i=1…n。

明明现在想求这个函数在 [0,1000][0,1000] 的最小值，要求精确到小数点后四位，四舍五入。

输入格式

输入包含 TT 组数据，每组第一行一个整数 nn；

接下来 nn 行，每行 33 个整数 aa, bb, cc ，用来表示每个二次函数的 33 个系数。注意：二次函数有可能退化成一次。

输出格式

每组数据输出一行，表示新函数 F(x)F(x) 的在区间 [0,1000][0,1000] 上的最小值。精确到小数点后四位，四舍五入。

样例

Input	Output
2 1 2 0 0 2 2 0 0 2 -4 2	0.0000 0.5000

数据范围与提示

对于 50\%50% 的数据，1 \leq n \leq 1001≤n≤100；

对于 100\%100% 的数据，1 \leq T \leq 101≤T≤10, 1 \leq n \leq 10^51≤n≤105, 0 \leq a \leq 1000≤a≤100, 0 \leq |b| \leq 50000≤∣b∣≤5000, 0 \leq |c| \leq 50000≤∣c∣≤5000。

Sponsor

 

#include<iostream>
#include<algorithm>
#include<iomanip>
using namespace std;
typedef long long ll;
double* a, * b, * c;//用inline可以减少时间和空间的浪费
inline double value(double x, int i) {//求S（x）的函数值
	return a[i] * x * x + b[i] * x + c[i];
}
double fvalue(double mid, int n) {//求F（x）的函数值
	double low = value(mid, 0);
	for (int i = 1; i < n; i++) {
		low = max(low, value(mid, i));
	}
	return low;
}
int main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	int t;
	cin >> t;
	double pre = 1e-9;
	while (t--) {
		int n;
		cin >> n;
	    double left = 0, right = 1000, mid1, mid2;
		a = new double[n];
		b = new double[n];
		c = new double[n];
		for (int i = 0; i < n; i++) {
			cin >> a[i] >> b[i] >> c[i];
		}
		while (left + pre < right) {//三分,这里必须补充一个值弥补误差否则会造成死循环
			mid1 = left + (right - left) / 3.0;//这个精度至少是1E-9（我试出来的），比这个小就行（也别太小）
			mid2 = right - (right - left) / 3.0;
			if (fvalue(mid1, n) > fvalue(mid2, n)) left = mid1;
			else right = mid2;
		}
		cout << setiosflags(ios::fixed) << setprecision(4) << fvalue(left, n) << endl;
	}//注意输出的是函数值
}