Fibonacci数列求余数和尾递归的一次实践

Fibonacci数列求余数和尾递归的一次实践

作者：LOYISHEN
时间：2020年01月19日
关键词：Fibonacci数列、求余数、递归、尾递归

背景

在练习蓝桥杯的试题时，有一道题是求Fibonacci数列中第n个数对10007取余数的题目（轻声说一句，这是入门题）：
图片来源于蓝桥杯网页，侵删

实际操作

像我这种对递归算法只有初步了解的初等程序员，第一时间想到的肯定是用递归算法来计算Fibonacci数列的第n个数是啥，然后再对10007取余，得到想要的数。因为我并没有想象到其中会有一个巨大的坑。

可是看到练习题中的提示，说到的“不需要先计算出Fn的准确值，再将计算的结果除以10007取余数”，我一开始并不理解到底是什么意思，因为在我当时看来，求余数就是直接取余，怎么能不需要先计算出具体值来取余呢？当时我就一直很迷惑这件事，加上准备期末考试和自己有一些事，我就没有再继续想下去。

期间有一天我睡觉前又想了一下，觉得这件事肯定有什么数学基础在里面。比如，求余数可以把一个数拆开多个来分开求？（意思是能够把7拆开成3和4）

思考和验证算法

到现在，我才拿出来又看了一下，想着能不能证明一下我那个想法。经过证明（并不严格），发现真的可以：

n	1	2	3	4	5	6	7	8	9	10
Fn	1	1	2	3	5	8	13	21	34	55
mod(Fn, 10)	1	1	2	3	5	8	3	1	4	5
mod(mod(F(n-1), 10)+F(n-2), 10)	-	-	2	3	5	8	3	1	4	5

这里只是找了一个小的除数来验证，因为这又不是去严格证明，只需要把结果算出来就行了，所以就只是对10取余，代码是会改回来的。

总结出的公式就是：
$$F_n=mod(mod(F_{n-1},10)+F_{n-2}))$$
然后我就写出了下面这个代码：

#include <stdio.h>
#include <stdlib.h>

int fibonacci_mod(int n)
{
	if (1 == n || 2 == n)
	{
		return 1;
	}
	return fibonacci_mod(n-1) % 10007 + fibonacci_mod(n-2);
}

int main(int argc, char *argv[])
{
	int n = 0;
	scanf("%d", &n);
	printf("%d", fibonacci_mod(n));
	return 0;
}

运行的时候，输入题目的示例，有正确的输出。但是在尝试输入100的时候，我等了很久都没有获得输出，这里就看到了这个递归函数的弊端了：当递归的层数达到一定的数量时，即使电脑计算速度非常快，但也不能在短时间计算出这个结果的，即使有足够的时间来计算，但也会因为内存不足而出错。

其实看到了上面的代码，就可以知道其主要原因在于上面设计的递归算法。因为在上面的递归算法中，每一层的调用和计算都需要一定的时间和内存，程序运行后当输入的n足够大时，计算机的计算力就不够了，再加上题目要求有时间限制，那么这一个递归算法就不是一个好的递归算法了。

优化代码

其实还有一个更好的递归算法，叫做尾递归。很幸运，我以前在看某乎的时候了解过尾递归，知道尾递归能够很好地解决上面递归算法遇到的问题。关于尾递归，其实就是把所有的计算工作都留到最后一步来计算，而递归本身只传递每一个需要运算的数值到最后一步。

然后我写出了下面的代码：

#include <stdio.h>
#include <stdlib.h>

int fibonacci_mod(int n, int ret1, int ret2)
{
	if (0 == n)
	{
		return ret1;
	}

	return fibonacci_mod(n-1, ret2 % 10007, ret1 + ret2);
}


int main(int argc, char *argv[]) {
	
	int n = 0;
	scanf("%d", &n);
	printf("%d", fibonacci_mod(n, 0, 1));
	
	return 0;
}

这个就是这个程序用到的尾递归算法。当我输入100时，并没有出现像前面那样的需要等很长时间的情况，而是直接出来结果。

这就可以看出尾递归的好处了，同样的运算结果，运算速度却有巨大的差别，幸好之前我了解过尾递归 :）

但是！当我输入题目给出的最大的数值时，程序给出的是一个错误！

哇塞，这个不是最好的递归算法了吗？为什么还会出错，而且还是直接退出？

其实这里应该就是堆栈溢出了吧！当递归层数太多的时候，程序把所有的堆栈空间都占用了都还不够用，可是程序还没有运行完，还需要获取空间，于是程序就会退出有报错了。即使尾递归是最后的时候才运算，但它依然是递归，在调用函数的时候依然需要用到堆栈的空间，即使尾递归的一个函数调用用到的空间很少，但一旦调用得多了，累积起来那就肯定是不够用了。

再一次优化代码

然后我仔细去看了一下尾递归函数，这好像就是一个条件循环吧？我直接写成循环不就可以不去调用函数了吗！

于是又有了下面这个代码：

#include <stdio.h>
#include <stdlib.h>

int main(int argc, char *argv[]) {
	
	int n = 0;
	scanf("%d", &n);
	
	int ret1 = 0;
	int ret2 = 1;
	int tmp = 0;
	while (n > 0)
	{
		n--;
		tmp = ret1 + ret2;
		ret1 = ret2 % 10007;
		ret2 = tmp;
	}
	printf("%d", ret1);
	
	return 0;
}

运行一下，输入100，跟前面的结果一样，没问题；输入最大的数1000000，也并没有报错，计算结果也没有本质上的出错（因为第100个Fibonacci数实在太大了，具体对10007取余求得的结果也不知道是啥，智能验证它并没有比10007大），于是就提交这段代码到平台上了。

其实利用循环能够减少不必要的堆栈操作，言外之意就是既能减少空间的使用，又能加快运算速度。

总结

其实在一开始的时候，我们都能够想到最直接最暴力的递归计算，因为这个是我们一开始学习编程语言的时候接触到的一个计算方式（好像说是计算方式也不太对）。

但是在使用递归的过程中，我们需要知道使用递归调用时不能有太多的递归层数，因为一旦递归得太深，那么程序运算的过程中就会需要用到更多的堆栈空间，并且也会因为需要用到堆栈空间而减慢运算的速度。

到最后使用循环的方法把这两个问题给解决了，这也说明在对待同一个问题的时候，并不会只有一个算法，而且最容易想到的算法也不一定会好。

所以，在解决问题的时候我们需要不断地去发现问题，然后根据这些问题来改进我们的算法，直到我们能够顺利完成我们的任务。

欢迎来访问我的个人博客