草庐IT

C语言实现求解斐波那契数列的四种方法及优化处理(递归,迭代,特殊性质公式,矩阵快速幂)

Ty_Man 2024-04-21 原文

        众所周知,斐波那契数列是非常经典的一个数列,它的数学公式如下

        为了便于观察,我们列出它的几项:0  1  1  2  3  5  8  13  21......

        下面我们将介绍四种方法来用C语言计算机代码实现对斐波那契数列的求解,分别是:递归法,迭代法,矩阵求解法以及特殊性质公式。

一、递归法

        (PS:没有递归基础的建议先学习递归的基础概念,在此我仅简要介绍一下递归的思想和求解代码)

        在递归的实现中,我们知道,递归有两个要求:(1)进行递归这一操作所需要满足的条件 (2)此条件需要最终不被满足,使得函数的嵌套调用能够返回。在斐波那契数列中,我们知道当x=0时,返回值为0;当x=1时,返回值为1。所以,要求(1)中所需要满足的条件即x>=2;而要求(2)中要使该条件最终不被满足,即x不断减小,且最终为0或1。

        于是我们给出代码的实现:

int Fibonacci(int x)
{
    if (x >= 2)
        return Fibonacci(x - 1) + Fibonacci(x - 2);    //当x>=2时,继续按照公式f(x)=f(x-1)+f(x-2)嵌套调用
    else if (x == 1)
        return 1;    //当x=1时,返回1,同时嵌套调用开始返回
    else
        return 0;    //当x=0时,返回0,同时嵌套调用开始返回
}

        用递归的方法求解斐波那契数列是不是非常简单呢?

        然而! 用递归的方法也有一个很大的缺点:栈溢出!

        栈:函数及其形参的空间由栈分配,在函数调用未完成时,为其分配的空间并不会销毁!于是当递归中函数的调用次数过多时,会导致栈溢出的问题。在本题的递归中,我们能发现,当我们求

   F(5)= F(4)+F(3)

= F(3)+F(2)+F(2)+F(1)

= F(2)+F(1)+F(1)+F(0)+F(1)+F(0)+F(1)时,不仅对F(2)进行了多次调用,甚至还进行了多次计算!这就导致计算效率的大幅降低以及资源占用的大幅增加。

        于是我们给出第二种方法。

二、迭代法

                                                什么是迭代法?

        迭代法也称辗转法,是一种不断用变量的旧值递推新值的过程。

        也就是说,这是一种各个变量直接相互作用的过程。

        分析思路:设有a、b、c三个变量,以及我们输入的第x项的x,令a=0(第一个斐波那契数)b=1(第二个斐波那契数),当x>=2时,c=a+b,同时再令a、b分别往右移一个数,如此往复x-2次,c的值就变成了第x个数的值。

        下面给出代码:

int Fibonacci(int x)
{
    int a = 0;
    int b = 1;
    int c = 0;
    if (x == 1)
        return 1;    //当x=1,返回1
    if (x == 0)
        return 0;    //当x=0,返回0
    while (x >= 2)   //输入x>=2时,进行迭代
    {
        c = a + b;   //每次迭代令c=a+b,即进行f(x)=f(x-1)+f(x-2)
        a = b;       //使得a,b往后移一个数字
        b = c;
        x--;
    }
    return c;
}

三、特殊性质公式法

        同递归法,只不过需要用到公式:

                F(n+m)=F(m)*F(n+1)+F(m-1)*F(n)

        在这个公式下,F(x)可以转换为F((x/2)   +   x-(x/2)),即一种二分思想,可以将时间复杂度为O(n)的算法转化为时间复杂度为O(log n)的算法。

        下面给出代码:

int Fibonacci(int x)
{
    if (x == 1 || x == 2)    //当x=2时也要返回1,否则将进入死循环。
        return 1;
    if (x == 0)
        return 0;
    int a = x / 2;
    int b = x - a;
    return Fibonacci(a + 1) * Fibonacci(b) + Fibonacci(a) * Fibonacci(b - 1);
}

         然而虽然但是事实上这个方法还是不如迭代法来得实在的QoQ

  四、矩阵法

                        

         于是,只要先算出系数矩阵的n-1次幂,再用与之相乘即可求出F(n)

         下面是用矩阵快速幂的方法求出所需矩阵n-1次幂的代码:

(PS:小编只学了快速幂,这个矩阵快速幂的算法是自己琢磨出来的,可能不标准,如果有心可以去学一下矩阵快速幂的算法!!下面会给出一般快速幂的算法给各位小可爱们参考一下~)

int Fibonacci(int x)
{
    if (x == 0)
        return 0;
    if (x == 1 || x == 2)
        return 1;
    int a[2][2] = { {1,1},{1,0} };
    int b[2][2] = { {1,1},{1,0} };
    int c[2][2] = { {1,0},{0,1} };
    int d[2][2] = { {1,0},{0,1} };
    int t = x - 1;
    while (t)
    {
        if (t%2==1)
        {
            d[0][0] = c[0][0] * b[0][0] + c[0][1] * b[1][0];
            d[0][1] = c[0][0] * b[0][1] + c[0][1] * b[1][1];
            d[1][0] = c[1][0] * b[0][0] + c[1][1] * b[1][0];
            d[1][1] = c[1][0] * b[0][1] + c[1][1] * b[1][1];
            c[0][0] = d[0][0];
            c[0][1] = d[0][1];
            c[1][0] = d[1][0];
            c[1][1] = d[1][1];
        }
        b[0][0] = a[0][0] * a[0][0] + a[0][1] * a[1][0];
        b[0][1] = a[0][0] * a[0][1] + a[0][1] * a[1][1];
        b[1][0] = a[1][0] * a[0][0] + a[1][1] * a[1][0];
        b[1][1] = a[1][0] * a[0][1] + a[1][1] * a[1][1];  
        a[0][0] = b[0][0];
        a[0][1] = b[0][1];
        a[1][0] = b[1][0];
        a[1][1] = b[1][1];
        t >>= 1;
    }
    return c[0][0];
}

         其中F(x) = c[0][0]*F(1)+c[1][0]*F(0) = c[0][0]

快速幂代码:

int f(int a,int n)
{
    int ans = 1;
    while (n)
    {
        if (n & 1 == 1)    //判断n的最后一位是否为1
        {
            ans *= a;
        }
        a *= a;            //a每次自乘a,使得a在每次循环中分别为a,a^2,a^3...
        n >>= 1;           //n右移一,使得n的二进制表示中1所在的位置所代表的位全与a当前的次幂相同
    }
    return ans;
}

好啦这就是全部内容啦!!!

 

有关C语言实现求解斐波那契数列的四种方法及优化处理(递归,迭代,特殊性质公式,矩阵快速幂)的更多相关文章

  1. ruby - 如何将脚本文件的末尾读取为数据文件(Perl 或任何其他语言) - 2

    我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚

  2. ruby - 如何根据特征实现 FactoryGirl 的条件行为 - 2

    我有一个用户工厂。我希望默认情况下确认用户。但是鉴于unconfirmed特征,我不希望它们被确认。虽然我有一个基于实现细节而不是抽象的工作实现,但我想知道如何正确地做到这一点。factory:userdoafter(:create)do|user,evaluator|#unwantedimplementationdetailshereunlessFactoryGirl.factories[:user].defined_traits.map(&:name).include?(:unconfirmed)user.confirm!endendtrait:unconfirmeddoenden

  3. ruby - 为什么 Ruby 的 each 迭代器先执行? - 2

    我在用Ruby执行简单任务时遇到了一件奇怪的事情。我只想用每个方法迭代字母表,但迭代在执行中先进行:alfawit=("a".."z")puts"That'sanalphabet:\n\n#{alfawit.each{|litera|putslitera}}"这段代码的结果是:(缩写)abc⋮xyzThat'sanalphabet:a..z知道为什么它会这样工作或者我做错了什么吗?提前致谢。 最佳答案 因为您的each调用被插入到在固定字符串之前执行的字符串文字中。此外,each返回一个Enumerable,实际上您甚至打印它。试试

  4. ruby - 寻找通过阅读代码确定编程语言的ruby gem? - 2

    几个月前,我读了一篇关于ruby​​gem的博客文章,它可以通过阅读代码本身来确定编程语言。对于我的生活,我不记得博客或gem的名称。谷歌搜索“ruby编程语言猜测”及其变体也无济于事。有人碰巧知道相关gem的名称吗? 最佳答案 是这个吗:http://github.com/chrislo/sourceclassifier/tree/master 关于ruby-寻找通过阅读代码确定编程语言的rubygem?,我们在StackOverflow上找到一个类似的问题:

  5. 华为OD机试用Python实现 -【明明的随机数】 2023Q1A - 2

    华为OD机试题本篇题目:明明的随机数题目输入描述输出描述:示例1输入输出说明代码编写思路最近更新的博客华为od2023|什么是华为od,od薪资待遇,od机试题清单华为OD机试真题大全,用Python解华为机试题|机试宝典【华为OD机试】全流程解析+经验分享,题型分享,防作弊指南华为o

  6. Unity 热更新技术 | (三) Lua语言基本介绍及下载安装 - 2

    ?博客主页:https://xiaoy.blog.csdn.net?本文由呆呆敲代码的小Y原创,首发于CSDN??学习专栏推荐:Unity系统学习专栏?游戏制作专栏推荐:游戏制作?Unity实战100例专栏推荐:Unity实战100例教程?欢迎点赞?收藏⭐留言?如有错误敬请指正!?未来很长,值得我们全力奔赴更美好的生活✨------------------❤️分割线❤️-------------------------

  7. 7个大一C语言必学的程序 / C语言经典代码大全 - 2

    嗨~大家好,这里是可莉!今天给大家带来的是7个C语言的经典基础代码~那一起往下看下去把【程序一】打印100到200之间的素数#includeintmain(){ inti; for(i=100;i 【程序二】输出乘法口诀表#includeintmain(){inti;for(i=1;i 【程序三】判断1000年---2000年之间的闰年#includeintmain(){intyear;for(year=1000;year 【程序四】给定两个整形变量的值,将两个值的内容进行交换。这里提供两种方法来进行交换,第一种为创建临时变量来进行交换,第二种是不创建临时变量而直接进行交换。1.创建临时变量来

  8. 基于C#实现简易绘图工具【100010177】 - 2

    C#实现简易绘图工具一.引言实验目的:通过制作窗体应用程序(C#画图软件),熟悉基本的窗体设计过程以及控件设计,事件处理等,熟悉使用C#的winform窗体进行绘图的基本步骤,对于面向对象编程有更加深刻的体会.Tutorial任务设计一个具有基本功能的画图软件**·包括简单的新建文件,保存,重新绘图等功能**·实现一些基本图形的绘制,包括铅笔和基本形状等,学习橡皮工具的创建**·设计一个合理舒适的UI界面**注明:你可能需要先了解一些关于winform窗体应用程序绘图的基本知识,以及关于GDI+类和结构的知识二.实验环境Windows系统下的visualstudio2017C#窗体应用程序三.

  9. MIMO-OFDM无线通信技术及MATLAB实现(1)无线信道:传播和衰落 - 2

     MIMO技术的优缺点优点通过下面三个增益来总体概括:阵列增益。阵列增益是指由于接收机通过对接收信号的相干合并而活得的平均SNR的提高。在发射机不知道信道信息的情况下,MIMO系统可以获得的阵列增益与接收天线数成正比复用增益。在采用空间复用方案的MIMO系统中,可以获得复用增益,即信道容量成倍增加。信道容量的增加与min(Nt,Nr)成正比分集增益。在采用空间分集方案的MIMO系统中,可以获得分集增益,即可靠性性能的改善。分集增益用独立衰落支路数来描述,即分集指数。在使用了空时编码的MIMO系统中,由于接收天线或发射天线之间的间距较远,可认为它们各自的大尺度衰落是相互独立的,因此分布式MIMO

  10. 【Java入门】使用Java实现文件夹的遍历 - 2

    遍历文件夹我们通常是使用递归进行操作,这种方式比较简单,也比较容易理解。本文为大家介绍另一种不使用递归的方式,由于没有使用递归,只用到了循环和集合,所以效率更高一些!一、使用递归遍历文件夹整体思路1、使用File封装初始目录,2、打印这个目录3、获取这个目录下所有的子文件和子目录的数组。4、遍历这个数组,取出每个File对象4-1、如果File是否是一个文件,打印4-2、否则就是一个目录,递归调用代码实现publicclassSearchFile{publicstaticvoidmain(String[]args){//初始目录Filedir=newFile("d:/Dev");Datebeg

随机推荐