斐波那契数列与卡特兰数

       斐波那契数列

       斐波那契数列（Fibonacci sequence），又称黄金分割数列、因数学家列昂纳多·斐波那契（Leonardoda Fibonacci）以兔子繁殖为例子而引入，故又称为“兔子数列”，指的是这样一个数列：0、1、1、2、3、5、8、13、21、34、……在数学上，斐波纳契数列以如下被以递归的方法定义：F（0）=0，F（1）=1，F（n）=F(n-1)+F(n-2)（n≥2，n∈N*）在现代物理、准晶体结构、化学等领域，斐波纳契数列都有直接的应用，为此，美国数学会从1963起出版了以《斐波纳契数列季刊》为名的一份数学杂志，用于专门刊载这方面的研究成果。

斐波那契数列指的是这样一个数列 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233，377，610，987，1597，2584，4181，6765，10946，17711，28657，46368..... 这个数列从第3项开始，每一项都等于前两项之和。

推广

编辑

斐波那契—卢卡斯数列

卢卡斯数列1、3、4、7、11、18…，也具有斐波那契数列同样的性质。（我们可称之为斐波那契—卢卡斯递推：从第三项开始，每一项都等于前两项之和f(n) = f(n-1)+ f(n-2）。

卢卡斯数列的通项公式为 f(n)=[(1+√5）/2]^n+[(1-√5)/2]^n

这两个数列还有一种特殊的联系（如下表所示），F（n）*L（n）=F(2n），及L（n）=F（n-1)+F(n+1)

n	1	2	3	4	5	6	7	8	9	10	…
斐波那契数列F(n)	1	1	2	3	5	8	13	21	34	55	…
卢卡斯数列L(n)	1	3	4	7	11	18	29	47	76	123	…
F(n)*L(n)	1	3	8	21	55	144	377	987	2584	6765	…

类似的数列还有无限多个，我们称之为 斐波那契—卢卡斯数列。

如1，4，5，9，14，23…，因为1，4开头，可记作F[1，4]，斐波那契数列就是F[1，1]，卢卡斯数列就是F[1，3]，斐波那契—卢卡斯数列就是F[a，b]。

斐波那契—卢卡斯数列之间的广泛联系

①任意两个或两个以上斐波那契—卢卡斯数列之和或差仍然是斐波那契—卢卡斯数列。

如：F[1,4]n+F[1,3]n=F[2,7]n，F[1,4]n-F[1,3]n=F[0,1]n=F[1,1](n-1），

n	1	2	3	4	5	6	7	8	9	10	…
F[1,4]n	1	4	5	9	14	23	37	60	97	157	…
F[1,3]n	1	3	4	7	11	18	29	47	76	123	…
F[1,4]n-F[1,3]n	0	1	1	2	3	5	8	13	21	34	…
F[1,4]n+F[1,3]n	2	7	9	16	25	41	66	107	173	280	…

②任何一个斐波那契—卢卡斯数列都可以由斐波那契数列的有限项之和获得，如

n	1	2	3	4	5	6	7	8	9	10	…
F[1,1](n)	1	1	2	3	5	8	13	21	34	55	…
F[1,1](n-1)	0	1	1	2	3	5	8	13	21	34	…
F[1,1](n-1)	0	1	1	2	3	5	8	13	21	34	…
F[1,3]n	1	3	4	7	11	18	29	47	76	123	…

黄金特征与孪生斐波那契—卢卡斯数列

斐波那契—卢卡斯数列的另一个共同性质：中间项的平方数与前后两项之积的差的绝对值是一个恒值，

斐波那契数列：|1*1-1*2|=|2*2-1*3|=|3*3-2*5|=|5*5-3*8|=|8*8-5*13|=…=1

卢卡斯数列：|3*3-1*4|=|4*4-3*7|=…=5

F[1，4]数列：|4*4-1*5|=11

F[2，5]数列：|5*5-2*7|=11

F[2，7]数列：|7*7-2*9|=31

斐波那契数列这个值是1最小，也就是前后项之比接近 黄金比例最快，我们称为黄金特征，黄金特征1的数列只有斐波那契数列，是独生数列。卢卡斯数列的黄金特征是5，也是独生数列。前两项 互质的独生数列只有斐波那契数列和卢卡斯数列这两个数列。

而F[1，4]与F[2，5]的黄金特征都是11，是孪生数列。F[2，7]也有孪生数列：F[3，8]。其他前两项互质的斐波那契—卢卡斯数列都是孪生数列，称为孪生斐波那契—卢卡斯数列。

广义斐波那契数列

斐波那契数列的黄金特征1，还让我们联想到佩尔数列：1，2，5，12，29，…，也有|2*2-1*5|=|5*5-2*12|=…=1（该类数列的这种 特征值称为勾股特征）。

佩尔数列Pn的递推规则：P1=1，P2=2，Pn=P(n-2)+2P(n-1).

据此类推到所有根据前两项导出第三项的通用规则：f(n) = f(n-1) * p + f(n-2) * q，称为广义斐波那契数列。

当p=1，q=1时，我们得到斐波那契—卢卡斯数列。

当p=1，q=2时，我们得到佩尔—勾股弦数（跟边长为整数的 直角三角形有关的数列集合）。

当p=2，q=-1时，我们得到等差数列。其中f1=1，f2=2时，我们得到自然数列1，2，3，4…。自然数列的特征就是每个数的平方与前后两数之积的差为1（等差数列的这种差值称为 自然特征）。

具有类似黄金特征、勾股特征、自然特征的广义——斐波那契数列p=±1。

当f1=1，f2=2，p=2，q=0时，我们得到 等比数列1，2，4，8，16……

相关数学

编辑

排列组合

有一段楼梯有10级台阶，规定每一步只能跨一级或两级，要登上第10级台阶有几种不同的走法?

这就是一个斐波那契数列：登上第一级台阶有一种登法；登上两级台阶，有两种登法；登上三级台阶，有三种登法；登上四级台阶，有五种登法……

1，2，3，5，8，13……所以，登上十级，有89种走法。

类似的，一枚均匀的硬币掷10次，问不连续出现正面的可能情形有多少种？

答案是（1/√5)*{[(1+√5)/2]^(10+2) - [(1-√5)/2]^(10+2)}=144种。

求递推数列a⑴=1，a(n+1)=1+1/a(n）的 通项公式

由 数学归纳法可以得到：a(n)=F(n+1)/F(n），将斐波那契数列的通项式代入，化简就得结果。

兔子繁殖问题

斐波那契数列又因数学家 列昂纳多·斐波那契以兔子繁殖为例子而引入，故又称为“ 兔子数列”。

一般而言，兔子在出生两个月后，就有繁殖能力，一对兔子每个月能生出一对小兔子来。如果所有兔子都不死，那么一年以后可以繁殖多少对兔子？

我们不妨拿新出生的一对小兔子分析一下：

第一个月小兔子没有繁殖能力，所以还是一对

两个月后，生下一对小兔对数共有两对

三个月以后，老兔子又生下一对，因为小兔子还没有繁殖能力，所以一共是三对

－－－－－－

依次类推可以列出下表：

经过月数	0	1	2	3	4	5	6	7	8	9	10	11	12
幼仔对数	1	0	1	1	2	3	5	8	13	21	34	55	89
成兔对数	0	1	1	2	3	5	8	13	21	34	55	89	144
总体对数	1	1	2	3	5	8	13	21	34	55	89	144	233

幼仔对数=前月成兔对数

成兔对数=前月成兔对数+前月幼仔对数

总体对数=本月成兔对数+本月幼仔对数

可以看出幼仔对数、成兔对数、总体对数都构成了一个数列。这个数列有关十分明显的特点，那是：前面相邻两项之和，构成了后一项。

这个数列是意大利中世纪数学家斐波那契在<算盘全书>中提出的，这个 级数的通项公式，除了具有a(n+2)=an+a(n+1）的性质外，还可以证明通项公式为：an=(1/√5)*{[(1+√5)/2]^n-[(1-√5)/2]^n}（n=1,2,3.....）

数列与矩阵

对于斐波那契数列1、1、2、3、5、8、13、……。有如下定义

F(n)=F(n-1)+F(n-2)

F(1)=1

F(2)=1

对于以下 矩阵乘法

F(n+1) = 11 F(n)

F(n) 10 F(n-1)

它的运算就是右边的矩阵 11乘以矩阵 F(n) 得到：

10 F(n-1)

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

F(n)=F(n)

可见该矩阵的乘法完全符合斐波那契数列的定义

设矩阵A=1 1 迭代n次可以得到：F(n+1) =A^(n) * F(1)= A^(n)*1

1 0 F(n) F(0) 0

这就是斐波那契数列的矩阵乘法定义。

另矩阵乘法的一个运算法则A^n(n为偶数) = A^(n/2)* A^(n/2)，这样我们通过二分的思想，可以实现对数复杂度的矩阵相乘。

因此可以用递归的方法求得答案。

数列值的另一种求法：

F(n) = [ (( sqrt ( 5 ) + 1 ) / 2) ^ n ]

其中[ x ]表示取距离 x 最近的整数。

卡特兰数

卡特兰数又称卡塔兰数，英文名Catalan number，是组合数学中一个常出现在各种计数问题中出现的数列。以比利时的数学家欧仁·查理·卡塔兰 (1814–1894)的名字来命名，其前几项为 : 1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900, 2674440, 9694845, 35357670, 129644790, 477638700, 1767263190, 6564120420, 24466267020, 91482563640, 343059613650, 1289904147324, 4861946401452, ...

令h(0)=1,h(1)=1，catalan数满足递推式 [1]   ：

h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)h(0) (n>=2)

例如：h(2)=h(0)*h(1)+h(1)*h(0)=1*1+1*1=2

h(3)=h(0)*h(2)+h(1)*h(1)+h(2)*h(0)=1*2+1*1+2*1=5

另类递推式 [2]   ：

h(n)=h(n-1)*(4*n-2)/(n+1);

应用

编辑

实质上都是递推等式的应用

括号化

矩阵连乘： P=a1×a2×a3×……×an，依据乘法结合律，不改变其顺序，只用括号表示成对的乘积，试问有几种括号化的方案？(h(n-1)种) [3]  

出栈次序

一个栈(无穷大)的 进栈序列为1，2，3，…，n，有多少个不同的 出栈序列? [4-5]  

常规分析

首先，我们设f（n）=序列个数为n的出栈序列种数。（我们假定，最后出栈的元素为k，显然，k取不同值时的情况是相互独立的，也就是求出每种k最后出栈的情况数后可用加法原则，由于k最后出栈，因此，在k入栈之前，比k小的值均出栈，此处情况有f(k-1)种，而之后比k大的值入栈，且都在k之前出栈，因此有f(n-k)种方式，由于比k小和比k大的值入栈出栈情况是相互独立的，此处可用乘法原则，f(n-k)*f(k-1)种，求和便是Catalan递归式。ps.author.陶百百）

首次出空之前第一个出栈的序数k将1~n的序列分成两个序列，其中一个是1~k-1，序列个数为k-1，另外一个是k+1~n，序列个数是n-k。

此时，我们若把k视为确定一个序数，那么根据 乘法原理，f（n）的问题就等价于——序列个数为k-1的出栈序列种数乘以序列个数为n - k的出栈序列种数，即选择k这个序数的f（n）=f（k-1）×f（n-k）。而k可以选1到n，所以再根据 加法原理，将k取不同值的序列种数相加，得到的总序列种数为：f（n）=f（0）f（n-1）+f（1）f（n-2）+……+f（n-1）f（0）。

看到此处，再看看卡特兰数的递推式，答案不言而喻，即为f（n）=h（n）= C（2n,n）/（n+1）= c（2n,n）-c（2n,n-1）（n=0，1，2，……）。

最后，令f（0）=1，f（1）=1。

非常规分析

对于每一个数来说，必须进栈一次、出栈一次。我们把 进栈设为状态‘1’，出栈设为状态‘0’。n个数的所有状态对应n个1和n个0组成的2n位 二进制数。由于等待入栈的操作数按照1‥n的顺序排列、入栈的操作数b大于等于 出栈的操作数a(a≤b)，因此输出序列的总数目=由左而右扫描由n个1和n个0组成的2n位二进制数，1的累计数不小于0的累计数的方案种数。

在2n位二进制数中填入n个1的方案数为c(2n,n),不填1的其余n位自动填0。从中减去不符合要求（由左而右扫描，0的累计数大于1的累计数）的方案数即为所求。

不符合要求的数的特征是由左而右扫描时，必然在某一奇数位2m+1位上首先出现m+1个0的累计数和m个1的累计数，此后的2(n-m)-1位上有n-m个 1和n-m-1个0。如若把后面这2(n-m)-1位上的0和1互换，使之成为n-m个0和n-m-1个1，结果得1个由n+1个0和n-1个1组成的2n位数，即一个不合要求的数对应于一个由n+1个0和n-1个1组成的排列。

反过来，任何一个由n+1个0和n-1个1组成的2n位 二进制数，由于0的个数多2个，2n为 偶数，故必在某一个奇数位上出现0的累计数超过1的累计数。同样在后面部分0和1互换，使之成为由n个0和n个1组成的2n位数，即n+1个0和n-1个1组成的2n位数必对应一个不符合要求的数。

因而不合要求的2n位数与n+1个0，n－1个1组成的排列一一对应。

显然，不符合要求的方案数为c(2n,n+1)。由此得出输出序列的总数目=c(2n,n)-c(2n,n+1)=c(2n,n)/(n+1)=h(n)。

类似问题 买票找零

有2n个人排成一行进入剧场。入场费5元。其中只有n个人有一张5元钞票，另外n人只有10元钞票，剧院无其它钞票，问有多少中方法使得只要有10元的人买票，售票处就有5元的钞票找零？(将持5元者到达视作将5元入栈，持10元者到达视作使栈中某5元出栈)

凸多边形三角划分

在一个 凸多边形中，通过若干条互不相交的对角线，把这个多边形划分成了若干个三角形。任务是键盘上输入凸多边形的边数n，求不同划分的方案数f（n）。比如当n=6时，f（6）=14。

分析

如果纯粹从f（4）=2，f（5）=5，f（6）=14，……，f（n）=n慢慢去归纳，恐怕很难找到问题的递推式，我们必须从一般情况出发去找规律。

因为凸多边形的任意一条边必定属于某一个三角形，所以我们以某一条边为基准，以这条边的两个顶点为起点P1和终点Pn（P即Point），将该凸多边形的顶点依序标记为P1、P2、……、Pn，再在该凸多边形中找任意一个不属于这两个点的顶点Pk（2<=k<=n-1），来构成一个三角形，用这个三角形把一个凸多边形划分成两个凸多边形，其中一个凸多边形，是由P1，P2，……，Pk构成的凸k边形（顶点数即是边数），另一个凸多边形，是由Pk，Pk+1，……，Pn构成的凸n-k+1边形。

此时，我们若把Pk视为确定一点，那么根据 乘法原理，f（n）的问题就等价于——凸k多边形的划分方案数乘以凸n-k+1多边形的划分方案数，即选择Pk这个顶点的f（n）=f（k）×f（n-k+1）。而k可以选2到n-1，所以再根据加法原理，将k取不同值的划分方案相加，得到的总方案数为：f（n）=f（2）f（n-2+1）+f（3）f（n-3+1）+……+f（n-1）f（2）。看到此处，再看看卡特兰数的递推式，答案不言而喻，即为f（n）=h（n-2） （n=2，3，4，……）。

最后，令f（2）=1，f（3）=1。

此处f（2）=1和f（3）=1的具体缘由须参考详尽的“卡特兰数”，也许可从 凸四边形f（4）=f（2）f（3）+ f（3）f（2）=2×f（2）f（3）倒推，四边形的划分方案不用规律推导都可以知道是2，那么2×f（2）f（3）=2，则f（2）f（3）=1，又f（2）和f（3）若存在的话一定是整数，则f（2）=1，f（3）=1。（因为我没研究过卡特兰数的由来，此处仅作刘抟羽的臆测）。

类似问题

一位大城市的律师在她住所以北n个街区和以东n个街区处工作。每天她走2n个街区去上班。如果她从不穿越（但可以碰到）从家到办公室的对角线，那么有多少条可能的道路？

在 圆上选择2n个点,将这些点成对连接起来使得所得到的n条线段不相交的方法数？

给定节点组成二叉搜索树

给定N个 节点，能构成多少种不同的 二叉搜索树？

（能构成h（N）个）

（这个公式的下标是从h(0)=1开始的）