三分查找

算法分析

---

二分查找是利用了函数的单调性，通过不停二分状态直到找到合适的值，但是如果遇到凸性或者凹性函数这种方法就不使用了，这时候就需要三分法了。
*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;
}