浅谈递推与生成函数的应用

最新推荐文章于 2022-07-25 23:24:50 发布

zhouchengyunew

最新推荐文章于 2022-07-25 23:24:50 发布

阅读量4.3k

点赞数 2

文章标签：扩展 c

本文链接：https://blog.csdn.net/zhouchengyunew/article/details/6249478

版权

这几天看了点组合数学的东西，以前一直以为排列计数之类的东西挺简单的，研究了几天感觉真是相当的难，那本书没坚持看完，看过的也没记住多少，感觉递推和生成函数那章有点意思，就写点这方面的，其他的有机会再好好钻研一下。

递推关系，相信大家都接触过，经常用于定义一个数列，比如著名的婓波那契序列：f₀ =0,f₁ =1,......f_n =f_n-1 +f_n-2 (n>=3).给定数列的初始元素值，其它元素根据前面的元素值可以推算出来。但是我现在希望能找到一个f_n 的通项公式，使它不依赖于前面元素的值，这样我们可以直接写出这个数列的任意元素的值，而不用一个一个的推算得来。如果有一个公式可以直接把婓波那契序列的第n项写出来，会不会觉得很神奇呀，下面我们就来把这个公式给找出来。

中心思想就是模仿微分方程的解题思路，先假设，再带入，设f_n =c * Xⁿ ,然后再将其带入递推关系式中，就以婓波那契序列为例，将f_n =c * Xⁿ 带入f_n =f_n-1 +f_n-2，可得 c * Xⁿ =c * X^n-1 +c * X^n-2 ,化简得，X² -X-1=0，这就是这个递推关系的特征方程，很容易解方程得到两个解：X₁ =(1+5^1/2 )/2和X₂ =(1-5^1/2 )/2，根据所求特征方程的解，可以写出f_n 的通项公式，f_n =C₁*(X₁ )ⁿ +C₂*(X₂ )ⁿ ，将f₀ =0,f₁ =1带入可得C₁₌ (5^1/2 -2)/5^1/2 和C₂₌ (5^1/2 +2)/5^1/2 ，将C_{1 ,}C_2, X_1, X₂上式即可得到f_n的通项公式。你可能觉得奇怪，这个数列的每一项都是整数，怎么通项中有带根号的项，你可以自己验证一下，算出来的结果正好都是整数。

推广开来，我可以写出递推关系的一般表达式：h_n-a₁*h_n-1-a₂*h_n-2... ...-a_k*h_n-k=b_n(a_k不为0且n>=k),其中a_1,a_{2,... ...}a_n-k为常数，若b_n=0,程该递推方程是齐次的，否则就是非齐次的。

针对齐次方程（b_n=0), 设通解为h_n=c * Xⁿ ，带入递推方程可得x^k-a₁*x^k-1-a₂*x^k-2...... a_k=0,解方程得X_1,X_2,...... X_k,这些值成为原递推方程的特征解。若解得的各个跟互不相同，则h_n的通项公式为h_n=C₁*(X₁ )ⁿ +C₂*(X₂ )ⁿ+... ...C_k*(X_k )ⁿ，若存在重跟，则在通项公式中与它对应的项相应的乘n即可。例如，若X₁=X_2,则通项公式改为h_n=C₁*(X₁ )ⁿ +C₂*n*(X₂ )ⁿ+... ...C_k*(X_k )ⁿ；若X₁=X₂=X_3,则通项公式改为h_n=C₁*(X₁ )ⁿ +C₂*n*(X₂ )ⁿ+C₃*n*n*(X₃ )ⁿ+... ...C_k*(X_k )ⁿ；然后再根据初始条件，分别求出系数C_1，C_{2，... ...}C_k的值即可。

若是非齐次方程，可以先求出对应的齐次方程的通项公式，设为f_{n ，}然后找原递推方程的一个特解，设为s_{n ，}则非齐次方程的通解为h_n=f_n+s_n ，f_n可有上面的方法求出，那我们如何找到递推方程的特解呢？这是一个比较难的问题，我们只关注一下其中最简单的两种，即b_n为多项式或指数形式的情况，若b_n为多项式，递推方程的特解的形式是跟b_n齐次的多项式；若b_n为指数函数aⁿ，则特解的形式为c * aⁿ，这里要注意，如果a恰好与X_1,X_2,...... X_k的某个值相等，那么特解的形式要相应的乘以n。下面我们来看个简单的例子，加深一下了解：

设递推方程为：h_n=3*h_n-1+3ⁿ,初始条件h₀=2 ，求数列的通项公式。

根据上述解题步骤，其对应的齐次方程为h_n=3 * h_{n-1 ，}设通项f_n=c * Xⁿ ，带入齐次递推方程，得X=3，故齐次方程的通项形式为f_n=c * 3ⁿ；接下来求原递推方程的特解，由于b_n=3ⁿ，不妨设特解形式为s_n =b * 3ⁿ带入原递推方程化简可得b=b+1,这显然不可能，仔细一看，原来指数中的3恰好为齐次方程的特征解，故特解的形式应该为s_n =b * 3ⁿ*n,带入原递推方程，可得b=1，于是特解为s_n = 3ⁿ* n，所以原递推方程的解为h_n=f_n+s_n =c * 3ⁿ+3ⁿ* n ，由初始条件h₀=2可知，c=2，所以h_n=2 * 3ⁿ+3ⁿ* n=(2+n)*3ⁿ

下面来介绍生成函数的方法，这中方法适合求解计数问题。

令h₀,h_{1,... ...}h_n...... 为一无穷数列，则它的生成函数为一个无穷级数：

g(x)=h₀+h₁*X+h₂*X²+... ...+h_n*Xⁿ+... ...

而有限序列h₀,h_{1,... ...}h_n可以看成是形为h₀,h_{1,... ...}h_n,0,0.... 的无限序列，其生成函数为

g(x)=h₀+h₁*X+h₂*X²+... ...+h_n*Xⁿ

例如，每一项都为1的无穷序列1，1，1... ... 的生成函数为

g(x)=1+X+X²+X³+... ...，这是一个等比数列的和，在其收敛域内，其值为g(x)=1/(1-x)，我们可以理解为函数1/(1-x)中包含了一个每一项都为1的无穷序列的信息。

再举个例子，令m为一整数,关于二项式系数C_{m,0 ,}C_{m,1 ,}C_{m,2 ... ...}C_m,m的生成函数为

g(x)=C_m,0+C_m,1* X +C_m,2*X²_{... ...}C_m,m*X^m=(1+X)^m

首先，我们需要把以前学过的二项式公式进行一下扩展，公式如下所示：

(a+b)^m=C_m,0*b^m+C_m,1* a*b^m-1 +C_m,2*a²*b^m-2_{... ...}C_m,m*a^m,

以前我们限定m必须为正整数，现在我们把它扩展为整个实数范围内，扩展后的每一项计算方法不变

C_m,r=m*(m-1)*(m-2)*......*(m-r+1) / (r!) , r=0,1,2,3,... n ,......,

因此，对于任意一个实数m，二项式公式扩展为

(a+b)^m=C_m,0*b^m+C_m,1* a*b^m-1 +C_m,2*a²*b^m-2_{... ...}C_m,m*a^m+C_m,m+1*a^m+1*b^-1+... ...

即它有无穷多的项，以前所学的公式只是m为整数的条件下的一种特殊情况，在这种情况下，对所有的n>m,C_m,n=0.

根据扩展的二项式公式，可以求得当m为负整数时，C_m,r=m*(m-1)*(m-2)*......*(m-r+1) / (r!)

=(-1)^r*(-m)*(-m+1)*(-m+2)*......*(-m+r-1)/(r!)

=(-1)^r*C_-m+r-1,r

上述结果应尽量记住，因为后面会经常用到上述公式。

下面我们来看一个比较有启发性的例子，考虑方程X₁+X₂+... ...+X_k=n (其中，n=0,1,2... ...)，设其非负整数解的个数为h_n,显然h_0,h₁,h₂,... ...h_n... ...构成一个无穷数列，其对应的生成函数为：

g(x)=h₀+h₁*X+h₂*X²+... ...+h_n*Xⁿ+... ...

先把原方程的解答放一边，考虑一个多项式相乘的例子，如（1+X+X²+X³)*(1+X+X²)*(1+X+X²+X³+X⁴)，现在我想知道结果中X³的系数，就是我会列举所有可能的组合情况，分别是1*1*X³,1*X*X²,1*X²*X,X*1*X²,X*X*X,X*X²*1,X²*1*X,X²*X*1,X³*1*1,总共有九种情况，因此我可以判定X³的系数为9.再仔细想想这个过程，由于1可以表示成X⁰，上面九个子项中X的指数和加起来都为3，如果限定X₁只能取0，1，2，3，X₂只能取0，1，2，X₃只能取0，1，2，3，4，那么上述九个子项正好跟X₁+X₂+X₃=3 的所有非负整数解是一一对应的，也就是说我们可以把这个求方程非负整数接个数的问题转化为一个求多项式系数的问题。

现在，我们在回到原来的方程，现在要求满足方程X₁+X₂+... ...+X_k=n (其中，n=0,1,2... ...)的非负整数解的个数，我们首先将这个问题等价成一个求多项式系数的问题呢，很显然，由于这里没有对各项的上界进行限制，我们可以模仿上面的多项式，如下：

（1+X+X²+X³+...)*(1+X+X²+X³+...)*(1+X+X²+X³+...)...(一共k个），即（1+X+X²+X³+...)^k

而原来的问题就转化为求（1+X+X²+X³+...)^k中Xⁿ的系数，这个系数就是无穷数列h_0,h₁,h₂,... ...h_n... ...的通项公式。那我们怎样才能求出Xⁿ的系数呢，考虑我们已经掌握的技术，在X的收敛域内，1+X+X²+X³+...=1/(1-x),用一下这个公式，（1+X+X²+X³+...)^k=（1/(1-x))^k=(1-X)^-k,很简洁吧，刚才大量的信息都表示到(1-X)^-k里面了，现在我们只要求出(1-X)^-k中Xⁿ的系数即可，可是这怎么求呢，考虑一下我们现在掌握的知识，对了，刚才我已经把二项式公式扩展到整个实数范围了，不记得的话再仔细看看前面的公式，根据公式可知，包含Xⁿ的项为C_-k,n(-X)ⁿ=(-1)ⁿ*C_k+n-1,n*(-1)ⁿ*Xⁿ=C_k+n-1,n*Xⁿ,即Xⁿ的系数为C_k+n-1,n，所以原方程X₁+X₂+... ...+X_k=n (其中，n=0,1,2... ...)的非负整数解的个数为C_k+n-1,n