算法分析
二分查找是利用了函数的单调性,通过不停二分状态直到找到合适的值,但是如果遇到凸性或者凹性函数这种方法就不使用了,这时候就需要三分法了。
*hint:三分函数不仅适用于凸函数,还适用于所有单峰函数。所谓单峰函数就是先严格递增后严格递减(此时存在唯一的最大值),或者先严格递减后严格递增(此时存在唯一的最小值)。具体理解还是自己画一个抛物线理解一下,其实很简单,看了一会就理解了。
三分函数
int SanFen(int l,int r){
while(l < r-1){
int mid = (l+r)/2;
int mmid = (mid+r)/2;
if(f(mid) > f(mmid)) r = mid;
else l = mid;
}
return f(l) > f(r)?l:r;
}
LA Error Curves
虽然有n条直线,但是组成的图形依然是单峰函数,三分即可。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 1e4+10;
int n,a[maxn], b[maxn], c[maxn];
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 T;
scanf("%d", &T);
while(T--){
scanf("%d", &n);
for(int i = 0; i < n; i++) scanf("%d%d%d", &a[i], &b[i], &c[i]);
double L = 0.0, R = 1000.0;
for(int i = 0; i < 100; i++){
double m1 = L+(R-L)/3;
double m2 = R-(R-L)/3;
if(F(m1) < F(m2)) R = m2;
else L = m1;
}
printf("%.4lf\n", F(L));
}
return 0;
}