Fibonacci数列的递归与非递归实现

原创 2012年03月23日 11:24:29

Fibonacci(斐波那契)数列定义:





(1)递归实现

思想简单,无需赘述,实现如下:

// 递归方法复杂度O(2^n),计算到n=41已达到计算机的极限
int F( int n )
{
	if ( n <= 0 ) return 1;
	else if ( n == 1 ) return 1;
	else return F(n-1)+F(n-2);
} 

int* Fibonacci1( int len )
{
	int* arr = (int*)malloc( sizeof(int) * len);
	if ( arr == NULL ) return arr;

	for ( int i = 0; i < len; ++i ) arr[i] = F(i);

	return arr;
}

该方法存在很多重复计算,虽然是典型的递归,但效率很低,资源消耗大。经分析,改算法时间复杂读为指数阶(O(2^n)),远高于多项式阶;占用计算机空间资源大,经测试在n=41时已达到PC的极限。


(2)采用通项公式计算

递推公式为x(n+2)=x(n+1)+x(n),是一个二阶常系数齐次线性差分方程,利用该差分方程的特征方程t^2=t+1,以及初始条件x(0)=1, x(1)=1易解出x(n)的通项表达式。

f(n)的通项表达式:





C实现:

struct data
{
	double f1;
	double f2;
};

// 利用通项公式计算,时间复杂度O(n)
int* Fibonacci2( int len )
{
	const double sqrt_5 = sqrt(5.0);
	const double p1 = (1+sqrt_5)*0.5;
	const double p2 = (1-sqrt_5)*0.5;
	data d = {(1+1.0/sqrt_5)/2.0, (1.0/sqrt_5-1)/2.0};

	int* arr = (int*)malloc( sizeof(int) * len);
	if ( arr == NULL ) return arr;

	for ( int i = 0; i < len; ++i )
	{
		arr[i] = (int)(d.f1-d.f2+0.5);
		d.f1 *= p1;
		d.f2 *= p2;
	}

	return arr;
}

虽然采用通项公式的算法能以O(1)的时间复杂度计算Fibonacci数列,但引入了浮点数,计算精度无法保证。


(3)采用二阶递推的方法

注意到,存在矩阵A=(1,1;1,0)使得下式成立



据此得到Fibonacci数列的实现代码:

// 用于迭代操作的2x2矩阵
class matrix
{
public:
	matrix():a(1),b(1),c(1),d(0){}
	matrix(int a1, int b1, int c1, int d1):a(a1),b(b1),c(c1),d(d1){}
	matrix( const matrix& B ):a(B.a),b(B.b),c(B.c),d(B.d){}
	const matrix& operator*=( const matrix& B)
	{
		int a1 = a, b1 = b, c1 = c, d1 = d;
		int a2=B.a, b2=B.b, c2=B.c, d2=B.d;
		a = a1 * a2 + b1 * c2;
		b = a1 * b2 + b1 * d2;
		c = c1 * a2 + d1 * c2;
		d = c1 * b2 + d1 * d2;
		return *this;
	}
	
	unsigned int a;
	unsigned int b;
	unsigned int c;
	unsigned int d;
};

const matrix operator*( const matrix& A, const matrix& B)
{
	matrix C(A);
	C *= B;
	return C;
}

void printMatrix( const matrix& A )
{
	printf("%d  %d\n", A.a, A.b);
	printf("%d  %d\n", A.c, A.d);
}

matrix matrixPower( const matrix& A, int n )
{
	matrix m(1,0,0,1); // 单位矩阵
	matrix temp(A);

	for ( ; n > 0; n >>= 1 )
	{
		if ( n&1 ) m *= temp;
		temp *= temp;
	}

	return m;
}

int F3(int n)
{
	matrix B(1,1,1,0);
	matrix A = matrixPower( B, n-1 );
	return A.a + A.c;
}

int* Fibonacci3( int len )
{
	int* arr = (int*)malloc( sizeof(int) * len );
	if ( arr == NULL ) return arr;

	for ( int i = 0; i < len; ++i )
	{
		arr[i] = F3(i);
	}

	return arr;
}
该算法时间复杂度为O(logN).

相关文章推荐

非递归实现fibonacci数列

  • 2012年10月20日 18:10
  • 542KB
  • 下载

Fibonacci递归与非递归实现

  • 2013年02月06日 15:13
  • 766B
  • 下载

斐波那契数列递归与非递归实现(JAVA语言描述)

斐波那契数列(Fibonacci sequence),又称黄金分割数列、因数学家列昂纳多·斐波那契(Leonardoda Fibonacci[1] )以兔子繁殖为例子而引入,故又称为“兔子数列”,指...

裴波那契数列的递归实现与非递归实现

斐波那契数列是数学家列昂纳多·斐波那契(Leonardoda Fibonacci[1]  )以兔子繁殖为例子而引入,也称为“兔子数列”。 指的是这样一个数列:0、1、1、2、3、5、8、13、21、3...

斐波那契数列的非递归实现(Java)

package chapter_4; public class Fibonacci { int i = 1; int j = 1; public void fibonacci(int n) ...

由斐波那契数列(Fibonacci)谈递归(Recursion)内部实现

一般说到递归,很多人第一时间可能都会想到斐波那契数列(Fibonacci),(注:斐波纳契,又称黄金分割数列,指的是这样一个数列:1、1、2、3、5、8、13、21、……数学上,斐波纳契数列以如下被以...

迭代与递归(及异同点的简单分析)以及Fibonacci数列的两种实现

迭代 迭代算法是用计算机解决问题的一种基本方法。它利用计算机运算速度快、适合做重复性操作的特点,让计算机对一组指令(或一定步骤)进行重复执行,在每次执行这组指令(或这些步骤)时,都从变量的原值推出它的...

[学习笔记]Fibonaci数列的递归与非递归算法实现

前一阵一直在忙实验室的事情,博客就暂时放下了,这次和大家分享一下Fibonaci数列的递归与非递归两种算法实现。        想必Fibonaci数列大家肯定都特熟,当输入正整数n为0或1时,f(0...

非递归实现广度遍历生成二叉树

  • 2014年11月05日 19:16
  • 1021B
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Fibonacci数列的递归与非递归实现
举报原因:
原因补充:

(最多只允许输入30个字)