[三分]求凹凸函数最值

https://www.luogu.org/problemnew/show/P3382
二分用于单调函数, 三分用于凹凸函数.
三分就是把区间分成三份, 设左右端点为L和R, 有两种分法。
第一种是区间内的第一个点lmid = (L + R) / 2, 第二个点rmid = (lmid + R) / 2。
第二种是lmid = l + (r – l) / 3, rmid = l + 2 * (r – l) / 3。
然后判断函数在这两个位置的值，根据实际情况缩小答案的区间。
比如求凸函数最大值，如果f(lmid) < f(rmid), 则L = lmid，否则R = rmid。然后循环一定次数保证精度，最后的R就是最大值点。如果是整数三分，那就是L = lmid + 1，否则R = rmid - 1。注意1e-5的精度至少64次，1e-6建议100次。实数三分：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;const int maxn = 1e3 + 5, mod = 1e9 + 7;int n;
double a[15];
double F(double x)
{double sum = 0;for (int i = 0; i <= n; ++i){double prod = 1;for (int j = 0; j < i; ++j){prod = prod * x;}sum += a[i] * prod;}return sum;
}int main()
{ios::sync_with_stdio(0); cin.tie(0);double l, r;cin >> n >> l >> r;for (int i = n; i >= 0; --i) cin >> a[i];for (int i = 0; i < 32; ++i){//注意观察区别double lmid = l + (r - l) / 3, rmid = l + 2 * (r - l) / 3;//double lmid = (l + r) / 2, rmid = (lmid + r) / 2;if (F(lmid) < F(rmid)) l = lmid;else r = rmid;}cout << fixed << setprecision(5) << r << endl;return 0;
}

整数三分：
这里有个小技巧，为什么要r – l > 3呢，是为了防止l和r离的太近使得l或r直接跳过了极值点，二分不会出现这样的问题，最后求出来的[l, r]是一个区间，极值点一定在这个区间里，挨个暴力就好了，可以根据不同的题适当改变区间大小，这样check函数不需要那么严格。

int l = mi, r = ma;
while (r - l > 3){int t = (r - l) / 3;int lmid = l + t, rmid = l + 2 * t;if (check(lmid, rmid)) l = lmid + 1;else r = rmid - 1;
}

[三分]求凹凸函数最值相关推荐

7：C语言之求分段函数的值
#include <stdio.h>//定义main函数 int main() {//请在此添加'求分段函数的值'的代码/*****************Begin*********** ...
matlab求分段函数的值.,如何用MATLAB求分段函数的最小值和最大值？
7.1.1 分段线性插值所谓分段线性插值就是通过插值点用折线段连接起来逼近原曲线,这也是计算机绘制图形的基本原理.实现分段线性插值不需编制函数程序,MATLAB自身提供了内部函数interp1其主要 ...
求Ackerman函数的值
原题:已知Akerman函数的定义如下: { m+1 n=0 Akm(n,m)= { akm(n-1,1) n>0,m=0 {ak ...
遗传算法求取函数最值问题
目录 1. 关于遗传算法 2. 遗传算法的步骤 3. 代码实现 3.1 工具函数 3.1.1 目标函数 3.1.2 解码 3.1.3 交叉 3.1.4 变异 3.2 主函数部分 3.3 代码 4. 其 ...
差分进化算法求解函数最值问题
差分进化算法求解函数最值问题声明: 1.本文源代码来自书目<智能优化算法及其MATLAB实例(第3版)>,目的在于为MATLAB初学者提供更简明的代码解析,方便读者了解算法及MATLAB ...
Java 求分段函数g(x)的值
题目描述输入一个实数x,按下列方式求分段函数g(x)的值,输出时保留2位小数. x的取值范围 g(x)的值X<-40 0.75X-40≤x≤20 x^2(x的平方)x>20 x-6 输入 ...
求带便函数权值的最短路径
function [distance,path]=zuiduanlujing(D,start1,terminal); %求带便函数权值的最短路径 %嵌套调用函数 [B,F]=datainfun; ...
Python实现遗传算法求函数最值
Python实现遗传算法求函数最值详细源代码:GA.py 1.算法过程图解 2.详细过程举例说明 (1)待求解方程 (2)确定编码方案主要是确定编码长度: def segment_length(s ...
蚁群算法求最值c语言实现,蚁群算法代码(求函数最值)
<蚁群算法代码(求函数最值)>由会员分享,可在线阅读,更多相关<蚁群算法代码(求函数最值)(4页珍藏版)>请在人人文库网上搜索. 1.function F=F(x1,x2) % ...

[三分]求凹凸函数最值

[三分]求凹凸函数最值相关推荐

最新文章

热门文章