算法——三分法

A2274

于 2023-12-05 18:41:02 发布

阅读量1.5k

点赞数 30

分类专栏：算法与数据结构文章标签：算法 c++

本文链接：https://blog.csdn.net/qq_61350148/article/details/134813768

版权

算法与数据结构专栏收录该内容

16 篇文章 0 订阅

订阅专栏

三分法

简介

三分法求单峰(或单谷)的极值是二分法的简单拓展，单峰(或单谷)可以简单的理解为二次函数含最值的区间，左右两边严格单调，对于这样的二次函数：y = ax² + bx + c(a≠0)，最值在 x = - b/2a 处取得，则给定一个区间，该区间包含这个 x 值，就可以使用三分法求解取得最值的 x 值。

原理

拿上面的二次函数举例，若二次项系数a大于0，则该二次函数的图像是一个单谷，给定区间包含最值点 x = - b/2a ，则我们可以将该区间三等分，左三等分点为 $mi d L$ ，右三等分点为 $mi d R$ 。
若 $mi d L$ 对应的函数值大于 $mi d R$ 对应的函数值，则有两种情况：第一种是 $mi d L$ 和 $mi d R$ 都在最值点的左边，第二种是 $mi d L$ 在最值点的左边，而 $mi d R$ 在最值点的右边。可见，两种情况都是 $mi d L$ 在最值点左边。
若 $mi d L$ 对应的函数值小于 $mi d R$ 对应的函数值，则有两种情况：第一种是 $mi d L$ 和 $mi d R$ 都在最值点的右边，第二种是 $mi d L$ 在最值点的左边，而 $mi d R$ 在最值点的右边。可见，两种情况都是 $mi d R$ 在最值点右边。
结合上面两个推断，可以得出一个策略：对于二次项系数a大于0的二次函数，当 $mi d L$ 对应的函数值大于 $mi d R$ 对应的函数值时，让子区间的左端点取 $mi d L$ ；当 $mi d L$ 对应的函数值小于 $mi d R$ 对应的函数值时，让子区间的右端点取 $mi d R$ 。这样可以保证在缩小区间的时候让最值点一直在区间内。
对于二次项系数a小于0的二次函数，结论刚好是相反的。

复杂度

每次将长度为 $n$ 的区间缩小为 1/3 ，算法的时间复杂度为 O(log₃n) 。

建模

对于三分法能解决的问题，最难的是建立一个单峰(或单谷)函数，其它的有如下的模版：

// type_name是基本数据类型中 整型 或 浮点型 的一种，根据精度和题目而改变
type_name func_x(type_name x); // 根据传入的x取其对应的函数值
type_name left, right, midL, midR;
// 接收left和right的值
while ( /* 条件，实数三分和整数三分不同 */ )
{
	midL = left + (right - left) / 3;
	midR = right - (right - left) / 3;
	if (func_x(midL) > func_x(midR))
		// left = midL; 或 right  = midR; 根据函数的单调性而改变
	else
		// right  = midR; 或 left = midL; 根据函数的单调性而改变
}
// 输出结果