第一次接触三分算法,我认为是比较好理解的。三分算法用于求单峰函数的最值,与二分不同,二分用于求单调函数中趋近某个值的值。
三分搜索的实现主要是判断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)就等于了极值。(摘自博主huzujun)
来看这道题,有两条线段AB,CD,在AB上行走的速度为p,CD上行走的速度为q,平面里行走的速度为r(r>p,q)。求从A走到D的最短距离。
无论如何,路线都是AE—EF—FD,不妨先把E点定下来,现在要求EF/r+FD/q的最小值,根据初中数学知识,这个东西用勾股定理之类的表示出来发现他确实是一个单峰函数(或者可以用几何画板画一画),是可以用三分法求出来的。
不妨记f(E)为min(EF/r+FD/q),也就是说,现在我们又要求AE/p+f(E)的最小值,这个呢,感性理解下,它似乎也是单峰函数,这个我不知道如何证明,如果有同学会的话欢迎评论。
所以,综合起来,就是一个三分套三分。
#include<bits/stdc++.h>
using namespace std;
const double eps=1e-8;
double ax,ay,bx,