最优化方法:斐波那契搜索

适用于无约束、单变量或一个自由度问题

#include<iostream>
#include<vector>
#include<limits>
#include<cmath>
#include<algorithm>
using namespace std;

vector<int> vect(46);
int feibo(int n){
    if(vect[n] != 0) return vect[n];
    int ans;
    if(n <= 1) ans = 1;
    else ans = feibo(n-1) + feibo(n-2);
    vect[n] = ans;
    return ans;
}
double func(double x){
    return pow(x,2) - 6.0*x + 2.0;
}
int main(){
    feibo(45);
    // for(auto it : vect) cout<<it<<" ";
    double a,b,error;
    double pre_x=INT_MIN,pre_y=INT_MIN;
    cin>>a>>b>>error;
    int n = lower_bound(vect.begin(),vect.end(),ceil(1/error)) - vect.begin();
    cout<<n<<endl;
    while(n>0){
        double l = (pre_x==INT_MIN)?(vect[n-1]*1.0*(b-a)/vect[n]):(pre_x - a);
        double x1 = a+l,x2 = b-l;
        double y1 = (x1=
  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值