【算法笔记】三分法

前言

三分法通常用于求单峰函数的极值，就像这样：

那么，三分法具体是怎样一个算法呢？

大致思路

我们拿求单峰函数中最小值举例。
首先确定一个区间，然后将它平均分成三份：

这时我们可以发现，当lmid对应的值<rmid对应的值时，min一定不会在rmid的右边，因此r到rmid中间可以被排除，我们可以把rmid赋值给r；反之，我们把lmid的值赋值给l。
重复进行上述操作，直到找到min值。

代码实现

我们以洛谷上P3382为例。

#include <iostream>
using namespace std;

const int N=20;
double l,r,a[N];
int n;
double f(double x)
{
	double t=0;
	for(int i=0;i<=n;i++)
	{
		t=t*x+a[i];//秦九韶算法，可以自行推理一下，当然按照题意模拟也可以
	}
	return t;
}
int main()
{
	cin>>n>>l>>r;
	for(int i=0;i<=n;i++) cin>>a[i];
	while(r-l>=1e-5)//根据题目要求的相对或绝对误差决定停止条件
	{
		double t=(r-l)/3.0;//求区间的小的1/3是多少
		double lmid=l+t,rmid=r-t;//lmid是在l的基础上在向右t，rmid同理
		if(f(lmid)>f(rmid)) r=rmid;
		else l=lmid;
	}
	cout<<l;
	return 0;
}

一些补充

此外，任何三分问题都可以转换为二分问题解决，并且二者运行时间接近。所以只要掌握好二分也就能解决此类问题了，当然博主本人认为三分解决求单峰极值一类问题时更加直观易懂，故作本文。