码迷,mamicode.com
首页 > 其他好文 > 详细

三分法

时间:2018-08-25 11:30:52      阅读:155      评论:0      收藏:0      [点我收藏+]

标签:href   图片   nal   png   amp   show   int   img   ima   

三分法用来求解函数的极值,极值左右区间满足单调性。

我们使用类似二分的思想来求解 极大值(极小值同理):

我们定义 \(\displaystyle mid = \frac{l+r}{2}\)\(\displaystyle mmid = \frac{mid+r}{2}\),分类讨论:

  • 如果 \(f(mid) > f(mmid)\)
    • 如果二者在同侧,\(f(mid)\) 更靠近 极大值,舍去 \((mmid, r]\) 区间;
    • 如果二者在极大值两侧,舍去其中一边都可以,不影响答案区间。
  • 如果 \(f(mid) < f(mmid)?\)
    • 如果二者在同侧,\(f(mmid)\) 更靠近 极大值,舍去 \([l, mid)\) 区间;
    • 如果二者在极大值两侧,舍去其中一边都可以,不影响答案区间。
  • 如果 \(f(mid) = f(mmid)\)
    • 二者一定在极大值两侧,舍去其中一边都可以。

容易实现。

以下是代码:(题面:LG3382 技术分享图片

#include <cstdio>
#define eps 1e-6
int n; double l, r, a[14], mid, mmid;

inline double calc(double x) {
    double res = 0;
    for (int i=0; i<=n; i++) res = res * x + a[i];
    return res;
}

int main() {
    scanf("%d%lf%lf", &n, &l, &r);
    for (int i=0; i<=n; i++) scanf("%lf", &a[i]);
    while (r > l + eps) {
        mid = (l + r) / 2.0;
        mmid = (mid + r) / 2.0;
        if (calc(mid) > calc(mmid)) r = mmid;
        else l = mid; 
    }
    printf("%.5f\n", l);
    return 0;
}

三分法

标签:href   图片   nal   png   amp   show   int   img   ima   

原文地址:https://www.cnblogs.com/greyqz/p/9532893.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!