标签:upload font ntc tco 递增 讲解 class article 序列
简单点说二分是查找区间,相当于一次函数,三分就是二次函数了,求它的极值,怎么做,数学常用的是求导,计算机就用查找咯,那么请看下面的简单概述吧!
mid = (left + right) / 2;
midmid = (mid + right) / 2;
(3)、我们mid比midmid更靠近最值,我们就舍弃右区间,否则我们舍弃左区间?。
if (cal(mid) > cal(midmid)) right = midmid; else left = mid;
1 double three_devide(double low,double up) 2 { 3 double m1,m2; 4 while(up-low>=eps) 5 { 6 m1=low+(up-low)/3; 7 m2=up-(up-low)/3; 8 if(f(m1)<=f(m2)) 9 low=m1; 10 else 11 up=m2; 12 } 13 return (m1+m2)/2; 14 }
1 const double EPS = 1e-10; 2 double calc(double x) 3 { 4 // f(x) = -(x-3)^2 + 2; 5 return -(x-3.0)*(x-3.0) + 2; 6 } 7 8 double ternarySearch(double low, double high) 9 { 10 double mid, midmid; 11 while (low + EPS < high) 12 { 13 mid = (low + high) / 2; 14 midmid = (mid + high) / 2; 15 double mid_value = calc(mid); 16 double midmid_value = calc(midmid); 17 if (mid_value > midmid_value) 18 high = midmid; 19 else 20 low = mid; 21 } 22 return low; 23 }
调用ternarySearch(0, 6),返回的结果为3.0000
我们都知道 二分查找 适用于单调函数中逼近求解某点的值。
如果遇到凸性或凹形函数时,可以用三分查找求那个凸点或凹点。
下面的方法应该是三分查找的一个变形。
如图所示,已知左右端点L、R,要求找到白点的位置。
思路:通过不断缩小 [L,R] 的范围,无限逼近白点。
做法:先取 [L,R] 的中点 mid,再取 [mid,R] 的中点 mmid,通过比较 f(mid) 与 f(mmid) 的大小来缩小范围。
当最后 L=R-1 时,再比较下这两个点的值,我们就找到了答案。
1、当 f(mid) > f(mmid) 的时候,我们可以断定 mmid 一定在白点的右边。
反证法:假设 mmid 在白点的左边,则 mid 也一定在白点的左边,又由 f(mid) > f(mmid) 可推出 mmid < mid,与已知矛盾,故假设不成立。
所以,此时可以将 R = mmid 来缩小范围。
2、当 f(mid) < f(mmid) 的时候,我们可以断定 mid 一定在白点的左边。
反证法:假设 mid 在白点的右边,则 mmid 也一定在白点的右边,又由 f(mid) < f(mmid) 可推出 mid > mmid,与已知矛盾,故假设不成立。
同理,此时可以将 L = mid 来缩小范围。
1 int SanFen(int l,int r) //找凸点 2 { 3 while(l < r-1) 4 { 5 int mid = (l+r)/2; 6 int mmid = (mid+r)/2; 7 if( f(mid) > f(mmid) ) 8 r = mmid; 9 else 10 l = mid; 11 } 12 return f(l) > f(r) ? l : r; 13 }
标签:upload font ntc tco 递增 讲解 class article 序列
原文地址:http://www.cnblogs.com/ECJTUACM-873284962/p/6536414.html