求斐波那契数列的三种方法
什么是斐波那契数列,1,1,2,3,5,8,13...这样一个数列就是斐波那契数列,求第n项的值。
一、经典求法
观察数列可得,除了第一项和第二项,所有的数列的值都是前一项和前一项的前一项的加和,转换成函数也就是f(n) = f(n-1) + f(n-2)
public static int f1(int n) {if(n < 1) {return 0;}else if(n == 1 || n == 2) {return 1;}return f1(n-1) + f1(n-2);
}
显然,递归n次,时间复杂度O(2^n),太恐怖,所以,必须优化。
二、顺序求法
因为斐波那契数列可以从左到右顺序的求出每一项的值,因此只需要顺序计算到n项即可,时间复杂度为O(n)的,我们可以把它看成在单链表的最后插入一个右最后一个和倒数第二个指针指向的值来决定的。
public static int f2(int n) {if(n < 1) {return 0;}else if(n == 1 || n == 2) {return 1;}int res = 1;int pre = 1;int temp = 0;for(int i = 3; i < n; i++) {temp = res;res = pre + res;pre = temp;}return res;
}
三、状态矩阵相乘求法
这是一个时间复杂度为O(log n)的算法。因为斐波那契数列在n大于等于三的时候严格遵守递推数列f(n) = f(n-1) + f(n-2),而对于一个二阶的递推数列来说,我们可以用矩阵乘法来表示,且状态矩阵是2阶的方阵:
= A
(F(N),F(N-1))= (F(N-1),F(N-2))*A
这个矩阵是怎么得到的呢,其实回顾一下我们所学的遇到一个递推公式,如何求通项公式的问题,是不是就是构造一个等比数列啊,然后换元累加求解,举个例子,温习一下:
---------------------------------------------------------------------------------------------------------------------
使用特征方程计算:
可以看到,这两种方法计算的结果相同,且都是正确的,但使用特征方程求解,十分方便.使用特征方程的一个问题是,如何计算得到递推公式的特征方程,
如上图,计算一个递推公式的通项公式,只要将c1和c2的值带入r^2 = c1 * r + c2即可,然后用上述的方法求解通项公式即可.注意,上述方法仅针对
这种形式的递推公式.
原理解决,下面继续看具体求法:
----------------------------------------------------------------------------------------------------------------------
把斐波那契数列的前四项F(1) = 1、F(2) = 1、F(3) = 2、F(4) = 3带进去可得状态矩阵的为:
当求出状态矩阵之后,当n >= 2时,原来的公式可简化为:
(F(3),F(2))= (F(2),F(1))*A = (1,1)*A
(F(4),F(3))= (F(3),F(2))*A = (1,1)*A^2
.....
(F(N),F(N-1))= (F(N-1),F(N-2))*A = (1,1) *A^(n-2)
所以其斐波那契数列n项和的问题就转化为如何用一个最快的方法求一个矩阵的N次方的问题,而求矩阵N方的问题显然能够在O(log n)时间内解决的问题,因为一个矩阵的实质其实就是一个值,假设这个数为a,如何快速的求a^75次幂的值,因为75次幂的二进制形式为1001011,那么a^75就可以写成a^64*a^8*a^2*a^1,在这过程中,我们先求出a的一次幂,再求2次幂,在根据2次幂去求4次幂以此类推,最后根据32次幂求出64次幂,即75的二进制形式共有多少位我们就使用了几次乘法,而且只乘相位上为1的位置即可。
代码实现:
//特征矩阵法//求矩阵m的p次幂的值public static int[][] matrixPower(int[][] m, int p){int [][] res = new int[m.length][m[0].length];//先把res设为单位矩阵,相当于整数中的1for (int i = 0; i < res.length; i++) {res[i][i] = 1;}//临时矩阵int [][] tmp = m;for (; p != 0; p >>= 1) {if ((p & 1) != 0) {//按位与操作,其实也就是位置为1的乘把这个矩阵乘一下res = muliMatrix(res, tmp);}//等于0 的时候我临时矩阵自己乘一下tmp = muliMatrix(tmp, tmp);}return res;}//两矩阵相乘public static int[][] muliMatrix(int[][] m1, int[][] m2) {//构造出一个m1行和m2列的矩阵(矩阵相乘规则)int [][] res = new int[m1.length][m2[0].length];//给我们要得到的矩阵的每一个元素的值进行赋值操作for(int i = 0; i < m1.length; i++) {for(int j = 0; j < m2[0].length; j++) {for(int k = 0; k < m2.length; k++) {//i,j位置元素的值就是i行元素和j行元素乘积的加和res[i][j] += m1[i][k] * m2[k][j];}}}return res;}//利用矩阵乘法求解斐波那契数列第n项的值public int f3(int n){if(n < 1) {return 0;}else if(n == 1 || n == 2) {return 1;}//状态矩阵int[][] base = {{1,1},{1,0}};//最后的结果需要求矩阵的n-2次int [][] res = matrixPower(base, n-2);return res[0][0] + res[1][0];}
求斐波那契数列的三种方法相关推荐
- C语言--斐波那契数列(三种方法)
文章目录 一·介绍 二·代码实现 1·递归实现 2·迭代实现 3·数组实现 一·介绍 斐波那契数列,就是前两个数是1,之后从第三个数开始等于前面两个数的和,请用代码方式求出第n个斐波那契数列的大小. ...
- 斐波拉契数列 Java三种实现
对于很多Java初学者来说,求解斐波拉契数列,基本上所使用的方式都是递归.对于求解斐波拉契数列来说,递归的实现方式的效率是极其低下的. 在此,我贴出三种实现(递归,迭代,动态规划): 我们先看代码的执 ...
- 剑指offer——面试题9:求斐波那切数列的四种方法
剑指offer--面试题9:求斐波那切数列的四种方法 另一个相关的链接:https://blog.csdn.net/Allenlzcoder/article/details/80297333 总结下求 ...
- 斐波拉契数列的三种解法
斐波那契数列: f(n)=f(n-1)+f(n-2); n>=2 f(0)=0; f(1)=1; 即有名的兔子繁衍问题. 斐波那契数列共有三种解法,因而写这篇文章总结一下. 1. 递归求解 递归 ...
- 斐波拉契数列的三种实现方法
百度解释 斐波那契数列(Fibonacci sequence),又称黄金分割数列.因数学家列昂纳多·斐波那契(Leonardoda Fibonacci)以兔子繁殖为例子而引入,故又称为"兔子 ...
- C语言实现求解斐波那契数列的四种方法及优化处理(递归,迭代,特殊性质公式,矩阵快速幂)
众所周知,斐波那契数列是非常经典的一个数列,它的数学公式如下 为了便于观察,我们列出它的几项:0 1 1 2 3 5 8 13 21...... 下面我们将介绍四种方法来用C语言计算机 ...
- C++输出斐波那契数列的几种方法
定义: 斐波那契数列指的是这样一个数列:0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, ... 这个数列从第三项开始,每一项都等于前两项之和. 以输出斐波那 ...
- C语言实现斐波那契数列的几种方法
斐波那契数列指的是这样一个数列:1 .1.2.3.5.8.13.21.34.55.89--是从第三项开始每一项都是前两项之和. 用递推的方法定义为: f(0)=1;f(1)=1; f(n)=f(n-1 ...
- C语言斐波那契数列求和两种方法
递归调用函数: #define _CRT_SECURE_NO_WARNINGS int fun(int n); #include<stdio.h> int main() {//斐波那契数列 ...
- 求斐波那契数列前n项的值
Description 输入n,求斐波那契数列前n项的值.斐波那契数列规律如下:1, 1, 2, 3, 5, 8, 13,21, 34,55-,从第三项开始,每一项都是前面两项的和. Input 输入 ...
最新文章
- telegraf input的配置
- 成功解决pyinstaller生成exe缺少各种包的问题
- 关于setTimeout
- Entity Framework查询,EF执行SQl
- java设计模式迭代器模式_迭代器模式和Java
- 剖析如何自己写框架Struts和Mybatis,学其原理而非重复造轮子
- Linux开发标准LSB简介:Linux Standard Base
- c语言无线网络抓包程序,c语言实现抓包
- 我的台式机的几张照片
- 单元测试框架TestableMock快速入门(三):校验Mock调用
- 浅析机关单位人力资源内部控制
- PCL:RANSAC 空间直线拟合
- 计算机内码汉字,gbk内码(汉字机内码在线查询)
- php生成excel带图片格式,php导出excel单元格带图片显示方法
- Win2008建立kms服务器vlmcsd
- 基于element-ui的Vue计算工作日组件
- 同步011序列检测电路
- 汉语属于哪个语系_汉语语族的所属语系
- 什么是凸多边形和凹多边形
- Extremely hard RSA 低加密指数攻击
热门文章
- 【转帖】mysql锁机制
- deficit记忆_背单词:不要让“捷径”欺骗你的记忆
- private的使用
- java游戏繁体字名字_繁体字游戏名(精选500个)_繁体字游戏名字大全_繁体字游戏网名...
- 用户空间和内核空间的区别
- html浏览xsd,那点你不知道的XHtml(Xml+Html)语法知识(DTD、XSD)
- Python提取CSV数据统计四分位数
- discuz 模板php代码,自定义HTML模板DIY支持PHP代码解析
- python 抓百度招聘信息,将搜索结果解析下来,用来跟踪招聘趋势数据,免费数据变为有用的东东
- 【线性代数】n阶行列式