法国数学家拉普拉斯(对!就是那个拉普拉斯变换的发明者)在他的著作中提出了生成函数的概念
我认为生成函数解决组合问题就是把问题转换成二项式或递推关系问题的方法
问题
用1元、5元、10元三种纸币能拼出多少种100元的方案?
这其实是一个组合问题,用三种纸币组合成100元
这种题用生成函数来做就很方便
下面我们来分析这道题
分析1元硬币的使用情况(五元十元可类比)
一元硬币在最终的方案中可能用到0次、1次、2次.......
好,接下来,我们尝试将这个问题转换成二项式问题
我们将0次写为 x 0 x^0 x0,1次写为 x 1 x^1 x1,2次写成 x 2 x^2 x2
则,针对1元硬币的所有情况就可以表示为 A = x 0 + x 1 + x 2 . . . A=x^0+x^1+x^2... A=x0+x1+x2...
类比到其他两种硬币:
五元硬币: B = x 0 + x 5 + x 10 . . . B=x^0+x^5+x^{10}... B=x0+x5+x10...
十元硬币: C = x 0 + x 10 + x 20 . . . C=x^0+x^{10}+x^{20}... C=x0+x10+x20...
所以硬币组合的所有情况就是将以上三个多项式相乘
A l l C o u n t = ( x 0 + x 1 + x 2 . . . ) ∗ ( x 0 + x 5 + x 10 . . . ) ∗ ( x 0 + x 10 + x 20 . . . ) AllCount=(x^0+x^1+x^2...)*(x^0+x^5+x^{10}...)*(x^0+x^{10}+x^{20}...) AllCount=(x0+x1+x2...)∗(x0+x5+x10...)∗(x0+x10+x20...)
稍微观察可知,给出一个 和n,则多项式展开后的 x n x^n xn项的系数就是组合成n的组合的数量,我们现在想把它化成一个 ( m + x ) n (m+x)^n (m+x)n的形式,这样我们就可以用二项式的系数公式解决问题了。
但是,我们不可能全部打开计算,那就是穷举了,所以要利用数学知识来解决这个问题
观察一元硬币的表达式,明显是等比数列,则可利用等比数列求和公式得到
A = 1 − x n 1 − x A=\displaystyle\frac{1-x^n}{1-x} A=1−x1−xn
我们再来化简这个式子,我们知道x的取值实际上是没有意义的,所以为了化简,我们规定 0 < x < 1 0<x<1 0<x<1
所以可得:
A = 1 1 − x A=\displaystyle\frac{1}{1-x} A=1−x1
类比得: B = 1 1 − x 5 C = 1 1 − x 10 B=\displaystyle\frac{1}{1-x^5}\ \ \ \ C=\displaystyle\frac{1}{1-x^{10}} B=1−x51 C=1−x101
这个时候我们发现ABC不能化简成一个二项式的形式,也就是说不能用组合数来表示第N项前面的系数
所以我们转而求其递推式 D n D_n Dn
A式不动,仍旧表示只用1元硬币的情况,我们来分析一下A中的系数表示
由于 A ∗ ( 1 − x ) = 1 A*(1-x) = 1 A∗(1−x)=1得到 A − A x = 1 A-Ax = 1 A−Ax=1说明 A n = A n − 1 ( A n 表 示 A 式 指 数 为 n 项 的 系 数 ) A_n=A_{n-1} (A_n表示A式指数为n项的系数) An=An−1(An表示A式指数为n项的系数)
这是为什么呢?由于 A x Ax Ax可以看做将A的每一项指数加1,相减后剩下第一项 x 0 x^0 x0,后面的项都减没了,说明n-1项的系数与n项的系数是一样的即: A n = A n − 1 A_n=A_{n-1} An=An−1
接下来用同样的方法处理,我们用N式表示用1元和5元的情况,D式表示全用的情况即:
N = ( x 0 + x 5 + x 10 . . . ) A N = (x^0+x^5+x^{10}...)A N=(x0+x5+x10...)A
D = ( x 0 + x 10 + x 20 . . . ) N D = (x^0+x^{10}+x^{20}...)N D=(x0+x10+x20...)N
我们处理一下上式:
KaTeX parse error: Can't use function '$' in math mode at position 14: N(1-x^5)-A=0 $̲ 打开得 $N-Nx^5-A…
得: N n − N n − 5 − A n = 0 N_n -N_{n-5}-A_n=0 Nn−Nn−5−An=0
同理可得: D n − D n − 10 − N n = 0 D_n -D_{n-10}-N_n=0 Dn−Dn−10−Nn=0
这样我们就得到了一个递推关系式:
D n = D N − 10 + N n D_n = D_{N-10}+N_n Dn=DN−10+Nn
N n = N n − 5 + A n N_n =N_{n-5}+A_n Nn=Nn−5+An
A n = A n − 1 = 1 A_n = A_{n-1} = 1 An=An−1=1
按照此递归关系式可以得到线性时间的 D n D_n Dn算法
代码:
#include "stdafx.h"
#include
int gene5(int n) {
if (n < 5) return 1;
else return gene5(n - 5) + 1;
}
int generation(int n) {
if (n < 5) { return 1; }
else if(n<10)return gene5(n);
else return generation(n - 10) + gene5(n);
}
int main()
{
std::cout<
我们搞定了组合问题,但是组合问题的解决方法不能解决排列问题,比如说:
拿6个球,有红球绿球,有多少种方法?
我们来看一下区别,之前钱的问题没有顺序一块和五块哪个先挑出来是无所谓的的因素,这次红绿之间有顺序的问题
这样我们就没办法用之前的方法来操作了,这时候,我们用到了拉普拉斯提出的另一种生成函数——指数生成函数
指数生成函数的一般形式是 x + 1 2 ! x 2 + 1 3 ! x 3 . . . x+\displaystyle\frac{1}{2!}x^2+\displaystyle\frac{1}{3!}x^3... x+2!1x2+3!1x3...
处理最后得出的生成函数表达式时要化成这样的基本形式,每项前写上它的系数(不理解看后面)
我们来看一下与常规生成函数的区别,式子的每一项的指数 n 还是代表选择 n 个绿(红)球
如果没有条件约束,两个式子乘起来就是绿球和红球能抽出的所有情况的生成函数
但是跟常规生成函数的区别在于,每一项前面多出了一个常数项 1 r ! \displaystyle\frac{1}{r!} r!1
这个数有什么用呢?
我猜想它跟组合数有关系,你想要把排列的方案数显示在生成函数中,必须做点什么,我们来看看它到底做了什么才能把组合数巧妙地融合进生成函数中
先来看看我们想要什么,对于某个确定的n(球的总数),当绿球取r个,红球取n-r个时,我们想要的结果是 ∑ r = 0 n C n r \sum\limits_{r=0}^n{C^r_n} r=0∑nCnr
明确答案后我们来看一个例子
如果是三个球的组合排列,两个红球和三个绿球组成的排列中,对应的项 1 2 ! x 2 \displaystyle\frac{1}{2!}x^2 2!1x2与 1 3 ! x 3 \displaystyle\frac{1}{3!}x^3 3!1x3乘积再化成 x 5 x^5 x5对应的形式
10 5 ! x 5 \displaystyle\frac{10}{5!}x^5 5!10x5
那么10就是两个红球和三个绿球的方案数,初见很神奇,看一下其中的原理,我们知道:
C n r = n ! r ! ( n − r ) ! C^r_n = \displaystyle\frac{n!}{r!(n-r)!} Cnr=r!(n−r)!n!
而 1 2 ! x 2 \displaystyle\frac{1}{2!}x^2 2!1x2和 1 3 ! x 3 \displaystyle\frac{1}{3!}x^3 3!1x3分别是 1 r ! x r \displaystyle\frac{1}{r!}x^r r!1xr和 1 ( n − r ) ! x n − r \displaystyle\frac{1}{(n-r)!}x^{n-r} (n−r)!1xn−r(n=5,r=2)
他们的乘积可以化成:
n ! r ! ( n − r ) ! n ! x n \displaystyle\frac{\displaystyle\frac{n!}{r!(n-r)!}}{n!}x^n n!r!(n−r)!n!xn
分子就是我们要求的方案数。
这时候我们验证了从数学上来看这个公式的合理性,我们再想一想,它的物理意义究竟是什么?
最终方案数的多少,可以用组合的和来表示,n!代表了针对一个组合所有排列的合计,而 1 r ! x r \displaystyle\frac{1}{r!}x^r r!1xr中的 1 r ! \displaystyle\frac{1}{r!} r!1用来消除重复的方案数,也就是绿色或红色自己本身的排序。
##三、生成函数求解任意常系数k阶齐次递推关系
##四、一些例子和总结
…未完