三分算法

类似于二分查找，三分搜索法也是比较常用的基于分治思想的高效查找方法。但是和二分不同，二分只适用于单调函数，三分用于单峰函数

三分搜索的实现主要是判断midl和midr所在值的大小。以凸函数为例（凹函数类似，只是判mid大小的时候保留小的即可（其实也是保留离极值最近的mid）），先以left和right为端点计算出它们的中点midl，然后再以midl和right为端点计算出它们的中点midr，接下来就需要判断f(midl)和f(midr)值的大小了，如果f(midl)大于f(midr)，那么说明midl靠近极值，此时令right=midr，否则说明midr靠近极值，此时则令left=midl，总之就是要保留离极值最近的那一个mid，然后重复前面的过程，直到left和right十分接近，最终f(left)就等于了极值，下面给出程序实现：

double f(double x){
    //something
}

const double eps=1e-8;
double sanfen(double l, double r){
    double mid,midr,ans;
    while (fabs(r-l)>eps) {
        mid=(l+r)/2;
        midr=(mid+r)/2;
        if(f(mid) < f(midr)) l=mid; else r=midr;   //求最大值
    } 
    ans=f(l);
    return ans;
}

hdu 4355
题意 给N(N≤50000)个人的位置 Xi 和重量 Wi
一个人走距离S的花费为 S3*W
要选择一个位置使所有人花费的和最小
看一眼感觉是单峰的 然后就三分了
这种求区间最大值,最小值的,不妨都三分试一下

#include <bits/stdc++.h>

#define closeTB ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);

using namespace std;
const int maxn=50010;
double mymax(double x, double y) {
    return x > y ? x : y;
}
double mymin(double x, double y) {
    return x < y ? x : y;
}
double x[maxn], d[maxn];
int n;
double trible(double x) {
    return x*x*x;
}
double f(double loc){
    double ans=0;
    for (int i=1; i<=n; i++) {
        ans+=trible(abs(loc-x[i]))*d[i];
    }
    return ans;
}

const double eps=1e-6;
double sanfen(double l, double r){
    double mid, midr;
    while (abs(r-l) > eps) {
        mid = (l+r) / 2;
        midr = (mid+r) / 2;
        if(f(mid) > f(midr)) l=mid; else r=midr;
    }
    return f(l);
}

int main() {
    int T, kase=0;
    scanf("%d", &T);
    while (T--) {
        scanf("%d", &n);
        scanf("%lf%lf", &x[1], &d[1]);
        double minx=x[1], maxx=d[1];
        for (int i=2; i<=n; i++) {
            scanf("%lf%lf", &x[i], &d[i]);
            maxx = mymax(maxx, x[i]);
            minx = mymin(minx, x[i]);
        }
        printf("Case #%d: %d\n", ++kase, (int)floor(sanfen(minx, maxx)+0.5) );
    }
}