题目描述
给你N个二次函数,记第i个为:fi(x)=aix^2+bix+ci,
(1≤N≤100000,0 ≤ ai ≤ 100, |bi| ≤ 5000, |ci| ≤ 5000)。
设函数F(x)=max{f1(x), f2(x) … fn(x)}。
请你求出F(x)的在区间[0, 1000]上的最小值,保留3位有效数字。
题解
因为F(X)是个凸函数所以可以三分
三分
考虑分情况考虑,如果lmid和rmid在极值的一侧,让远离极值的一侧往上跳肯定没问题,否则哪边往上跳都可以
代码
#include <bits/stdc++.h>
#define maxn 1000005
#define LL long long
#define INF 0x3f3f3f3f
#define I inline
#define re register
#define eps 1e-10
using namespace std;
int read(){
int res; bool f=1; char c;
while(!isdigit(c=getchar())) if(c=='-') f=0; res=c^48;
while(isdigit(c=getchar())) res=(res<<3)+(res<<1)+(c^48);
return f?res:-res;
}
int n;
double a[maxn],b[maxn],c[maxn];
int main(){
n=read();
for(int i=1;i<=n;i++){
scanf("%lf%lf%lf",&a[i],&b[i],&c[i]);
}
double l=0,r=1000.0,sl,sr;
while(l<r-eps){
double lm=(r-l)/3.0+l,rm=2.0*(r-l)/3.0+l;
sl=-1e100,sr=-1e100;
for(int i=1;i<=n;i++){
sl=max(sl,a[i]*lm*lm+b[i]*lm+c[i]);
sr=max(sr,a[i]*rm*rm+b[i]*rm+c[i]);
}
if(sl>sr) l=lm;
else r=rm;
}
printf("%.3lf\n",sl);
return 0;
}