【暖*墟】#三分法# 求单峰函数极值

最新推荐文章于 2022-09-11 17:22:08 发布

flora715_hss

最新推荐文章于 2022-09-11 17:22:08 发布

阅读量1.4k

点赞数 1

分类专栏： C++ 知识点

本文链接：https://blog.csdn.net/flora715/article/details/80983382

版权

本文介绍了利用三分法寻找单峰函数极值的思路，通过BZOJ1857 传送带问题作为实例，阐述了如何在上凸函数中找到极值点。文章首先解释了三分法的基本原理，然后详细分析了如何确定最优路径，最后讨论了如何应用三分法解决实际问题。

摘要由CSDN通过智能技术生成

【三分法】

1. 对于任意一个上凸函数，选取函数上任意两个点A,B（xA<xB），

若满足yA<yB,那么该函数的极值点必然在[xA,+∞)中，

若满足yA>yB，那么该函数极值点必然在(-∞,xB]中，

若满足yA=yB，那么该函数的极值点必然在[xA,xB]中。

2. 对于任意一个下凸函数，选取函数上任意两个点A,B（xA<xB），

若满足yA<yB,那么该函数的极值点必然在(-∞,xB]中，

若满足yA>yB，那么该函数极值点必然在[xA,+∞)中，

若满足yA=yB，那么该函数的极值点必然在[xA,xB]中。

三分法的思路与二分法很类似，不过其用途没有那么广泛，主要用于求单峰函数的极值。

这里写图片描述

此处，任选的点为m1，m2。与l，r把区间分成三段。

void Solve(){
    double left, right, m1, m2, m1_value, m2_value;
    left = MIN;
    right = MAX;
    while (left + EPS < right){
        m1 = left + (right - left)/3;
        m2 = right - (right - left)/3;
        m1_value = f(m1);
        m2_value = f(m2);
        if (m1_value >= m2_value)
            right = m2;  //假设求解极大值
        else  le