读者也许听说过“循环定义”,即在定义概念 A 的时候直接或间接地用到了 A 自身。例 如将“逻辑学”定义成“研究逻辑的科学”,这实际上是同语反复,并未揭示任何新的内涵; 又如将“美丽”定义成“漂亮”,再将“漂亮”定义成“美丽”,这种循环定义实际上也是同 语反复。循环定义是一种常见的逻辑错误,应尽量避免使用,但在数学和程序设计中,我们 经常在一个函数的定义中直接或间接地用到该函数自身,这种函数称为递归(recursive②) 函数。通过下面的讨论我们会看到,这种递归定义不同于循环定义,它能够明确地定义出函 数的意义。
递归是一种强大的算法设计思想和方法,利用递归可以轻松解决很多难题。下面我们通 过例子来介绍这种方法。
阶乘
数学中的阶乘运算通常用下式定义:
注意,当 n 为 0 时,其阶乘被定义为 1。
如果要编程计算 n 的阶乘,可以采用以前介绍过的累积算法模式来实现,累积算法的关 键部分是一个循环语句。下面是此方法的实现代码及执行实例:
24
>>> fac(40)
815915283247897734345611269596115894272000000000L
下面我们用另一种方式来观察阶乘的定义。在阶乘定义式中,等号右边的第一个 n 之后 的部分是什么?稍加思考即可看出就是(n-1)的阶乘,即阶乘定义式可写成:
n! = n x (n - 1)!
这个式子的含义是:n 的阶乘定义为 n 乘(n-1)的阶乘。我们看到,“阶乘”的定义中用到了 “阶乘”本身,这就是递归定义。
现代编程语言都支持递归函数,Python 也不例外。读者也许会将上面的递归定义式直 接翻译成如下 Python 函数:
def fac(n):
return n * fac(n-1)
但这个定义是错误的。如果执行这个函数,将会形成如下调用链条:
fac(n) => fac(n-1) => ... => fac(1) => fac(0) => fac(-1) => fac(-2) => ...
显然,递归将无穷无尽地延续下去①。
有效递归定义的关键是具有终止条件,使得到达某一点后不再递归,否则会导致像无穷循环一样的后果。对阶乘来说,当 n=0 时,n!的值定义为 1,此时无需递归。在上面的 fac 函数中添加这个终止条件,即可得正确的递归版阶乘函数:
为了理解递归函数的执行过程,需要回顾第 4 章中介绍的函数调用与返回的知识。图10.1 展示了 fac(2)的计算过程。
图 10.1 fac(2)的计算过程
计算 fac(n)时,由于每次递归都导致计算更小的数的阶乘,因此这个递归过程迟早会到 达计算 fac(0)的情形。而根据 fac()的定义,fac(0)直接返回 1,无需递归计算。我们称这种情 形为递归定义的奠基情形。对于奠基情形,递归函数能够直接计算结果。
要说明的是,上面的阶乘函数定义其实仍然有 bug:当 n 的初始值小于 1 时,调用 fac(n) 会导致无穷递归!解决这个问题很容易,只需在程序开始处检查 n 是否为负数即可,并且仅 当 n 为非负自然数时才能计算阶乘。编写递归程序时很容易在终止条件上面犯错误,作为好 的编程习惯,我们应当围绕递归奠基情形测试各种情形。
还要说明一点,每次递归调用 fac()都相当于调用一个新的函数,系统将为该函数的局 部变量和参数分配新的空间,与其他 fac()调用的局部变量和参数完全没有关系。初学者在 这一点上也会经常犯错误,以为各递归调用中使用的变量是全局共享的。在图 10.1 中有三 次对 fac(n)的调用,这三次调用应当视为独立的三个函数,其中用到的参数 n 应当视为三个 相互独立的局部变量。
列表处理
递归对于处理列表是非常有用的,因为列表本身就是“递归结构”——任一列表都可看 作是由第一个数据成员与剩余数据列表组成的,即:[a1,a2,…,an]可视为由 a1 和[a2,…,an]组成。 编程处理这个列表时,只需要单独考虑如何处理 a1,而对[a2,…,an]的处理可以通过递 归调用来解决。显然,每次递归都导致处理一个更短的列表,如此递归下去终将到达空列表 情形,这正可作为奠基情形。在 Python 中通过索引很容易取得列表 list 的第一个数据和剩 余数据列表,它们分别是 list[0]和 list[1:]。
作为例子,下面我们写一个递归函数来逆向显示列表的数据,即将[a1,a2,…,an]显 示为
an,...,a2,a1
根据列表的“递归结构”性质,不难形成这样的递归思考:为了逆向显示 list,只需先 逆向显示 list[1:],然后显示 list[0];当剩余数据列表为空时停止递归。这个递归思考可以直 接翻译成如下 Python 代码:
>>> def revList(list):
if list != []:
revList(list[1:])
print list[0],
else:
return
>>> revList([1,2,3,4,5])
5 4 3 2 1
对于简单列表的处理任务,用 for 循环语句也很容易实现;但当列表很复杂(例如列表 中的元素本身可能是列表),用循环语句就很难编程,而用递归则可以很容易地解决问题。 作为练习,读者不妨思考一下如何逆向显示如下形状的列表:
[1,[2,3],4,[5,6,[7,8],9]]
二分搜索
10.1 节中介绍了线性搜索算法,读者已经知道线性搜索的优点是适合无序的数据列表,缺点是不适合大量数据。当列表中的数据有序时,存在更好的搜索策略,这个策略的基本思 想可以通过一个猜数游戏展现出来。
假设某甲心中想好了一个 100 以内的自然数,让某乙来猜这个数。某乙每猜一次,某甲 都会告诉他猜对了、猜大了或猜小了三种情形之一。某乙该采用什么策略来玩这个游戏呢? 某乙可以每次都随机猜一个数,也可以系统化地按 1、2、3、……的顺序猜(此即线性搜索), 但这两种策略平均需要猜很多次才能猜中。最好的策略是先猜 1~100 的中间数 50,如果猜 中自不必说,如果猜大了则接下去猜 1~49 的中间数 25,如果猜小了则接下去猜 51~100 的中间数 75。依此类推,每次都猜可能值范围的中间值,直至猜中。这个策略就是我们要 介绍的二分搜索(binary search)算法。
下面我们利用二分搜索来解决在一个有序数据列表 list 中查找指定数据 x 的问题。先看 如何利用循环来实现二分搜索。算法的核心是一个循环,每一次循环都检查当前搜索范围的 中间数据是否等于 x;不等的话,根据大小信息重新设定搜索范围;如果找到了 x,或者没 有剩余数据了,则循环终止。为了便于设定搜索范围,可以用两个变量 low 和 high 分别记 录搜索范围的两端,每次循环后根据比较结果调整这两个变量即可重新设定搜索范围。代码 如下:
def binary(list,x): low = 0
high = len(list) - 1
while low <= high:
mid = (low + high) / 2
if list[mid] == x:
high = mid - 1
else:
low = mid + 1
return -1
再看二分搜索的递归实现。二分搜索算法可以这样表达:检查当前搜索范围的中间数据, 如果该数据就是目标数据,则算法结束;如果不是,则选择某一半范围重新进行二分搜索。 这段话可以翻译成如下伪代码:
这个算法中有三处(见划线部分)涉及几乎相同的操作,这正是二分搜索的递归性质的 体现。奠基情形是找到了目标值或者检查完所有数据都未找到目标值,这时将不再递归。由 于每次递归调用都将搜索空间减小了一半,因此迟早会到达奠基情形。下面给出递归版本二 分搜索的 Python 代码实现。注意与循环版本不同的是,每次递归都需要指明搜索范围,因 此我们将搜索范围的两个端点 low 和 high 也作为函数参数。
>>> def recBinSearch(list,x,low,high):
if low > high:
return -1
mid = (low + high) / 2
m = list[mid]
if m == x:
return mid
elif x < m:
return recBinSearch(list,x,low,mid-1)
else:
return recBinSearch(list,x,mid+1,high)
>>> recBinSearch([1,3,5,7,9],5,0,4) 2
>>> recBinSearch([1,3,5,7,9],6,0,4)
-1
阶乘和二分搜索这两个例子说明,许多问题既可用循环(或称迭代)来实现,也可用递 归来实现。很多情况下,两种方法在设计上都很容易;但对有些问题,迭代算法很难设计, 而递归算法则非常容易得到,例如下面的 Hanoi 塔问题。
Hanoi 塔问题
Hanoi 塔问题是体现递归方法强大能力的经典例子,该问题涉及如下故事:在某个寺庙 里有三根柱子(不妨称为 A、B、C 柱),A 柱上有 n 个同心圆盘,圆盘尺寸各不相同,并且 小盘叠在大盘之上,B 柱和 C 柱为空(参见图 10.2)。寺庙的僧侣们有一项任务:将 n 个圆 盘从 A 柱移到 C 柱,移动过程中可以利用 B 柱作为临时存放柱。具体的移动圆盘的规则是:
圆盘只能放在柱子上;
大圆盘永远不能置于小圆盘之上。
图 10.2 Hanoi 塔问题
下面我们来设计解决此问题的算法,该算法能够给出搬运步骤。例如对于 n = 3 的情形,算法将显示如下移动过程(其中 A -> C 表示将 A 柱顶部圆盘移至 C 柱顶部,余类推):
A -> C
A -> B
C -> B
B -> A
B -> C
A -> C
Hanoi 塔问题看上去有点难度,但如果采用递归方法,算法是非常简单的。稍加思考即 可明白,为了将 A 柱上的所有圆盘移到 C 柱,必然需要有一步将底部的最大圆盘从 A 柱移 到 C 柱,而为此又必须先将最大圆盘上面的 n - 1 个圆盘从 A 柱移到 B 柱,从而形成最大 圆盘之上没有其他圆盘、同时 C 柱上也没有圆盘的局面(参见图 10.3)。
图 10.3 为移动最大圆盘必须形成的格局
至此,可以将最大圆盘从 A 柱移到 C 柱,显然接下去再也不需要移动这个圆盘了,因此可视为不存在。这时形成的局面(图 10.4)与初始局面非常相似,即:B 柱上有 n - 1 个 圆盘,A 柱和 C 柱为空(无视最大圆盘)。于是任务变成了:将 n - 1 个圆盘从 B 柱移动到 C 柱,移动过程中可以利用 A 柱作为临时存放柱。将这里的黑体部分文字与前面的初始问 题文字相比较,即可看出 Hanoi 塔问题的递归性质:一旦最大圆盘到达目的地,剩下来的问 题恰好又是初始 Hanoi 塔问题,只不过问题规模变成了 n - 1,并且 A 柱和 B 柱的角色相互 交换。
图 10.4 最大圆盘就位之后的格局
根据以上分析,容易得到解决 Hanoi 塔问题的算法。下面是算法的伪代码:
算法:将 n 个圆盘从 A 柱移到 C 柱,以 B 柱作为临时存放柱。
将 n - 1 个圆盘从 A 柱移到 B 柱,以 C 柱作为临时存放柱;
将最大圆盘从 A 柱移到 C 柱;
将 n - 1 个圆盘从 B 柱移到 C 柱,以 A 柱作为临时存放柱。
从算法中可见,通过递归,我们将规模为 n 的问题转化成了两个规模为 n – 1 的问题。 如此递归下去,最终将转化成规模为 1 的问题。而 n = 1 的 Hanoi 塔问题是平凡的,直接移 动一步即可,不再需要递归,这就是奠基情形。有了奠基情形,每次递归又导致问题规模变 小,可见上述递归算法能正确终止。下面我们给出对上述算法的 Python 实现,并对 n = 3 的 情形进行验证。代码中 hanoi 函数的参数分别表示圆盘个数和三根柱子(源、目的地、临时 存放)。
>>> def hanoi(n,source,dest,temp):
if n == 1:
print source,"->",dest
else:
hanoi(n-1,source,temp,dest)
hanoi(1,source,dest,temp)
hanoi(n-1,temp,dest,source)
>>> hanoi(3,"A","C","B")
A -> C
A -> B
C -> B
A -> C
B -> A
B -> C
A -> C
至此,一个看上去挺难的问题通过递归就轻松地解决了。读者有兴趣的话不妨试试如何 利用循环(迭代)来解决 Hanoi 塔问题。
最后对递归方法做个小结。
递归是非常重要的算法设计方法,在解决很多具有递归性质的问题、结构的时候,设计
递归算法往往是直接而简单的。递归定义必须满足以下条件才是良定义的:
有一个或多个无需递归的奠基情形;
递归总是针对规模更小的问题。
虽然递归算法容易设计、实现,也容易理解,但递归是有代价的。由于递归涉及大量的 函数调用,因此需要耗费较多的内存和较长的执行时间,即递归算法的效率较差。而迭代算 法不涉及函数调用,故速度更快,更节省内存。