递推关系与母函数

本文最后更新于:2022年8月3日 晚上

『组合数学』2 递推关系与母函数

2.1 递推与递归关系

  • 递归
  • 汉诺塔
  • Fib

2.2母函数

  • 母函数

    G(x)=h(1)x+h(2)x++h(n)xnG(x) = h(1)x+h(2)x+\cdots+h(n)x^n

G(x)是序列h的母函数

Hanoi 应用:

G(x)=2xG(x)+x1x:(x<1)G(x)=x(12x)(1x)=112x11x=i=0n2ixii=0nxih(n)=2n1G(x) = 2xG(x) + \frac{x}{1-x}:(x<1)\\ \\G(x) = \frac{x}{(1-2x)(1-x)} = \frac 1{1-2x} -\frac 1{1-x} = \sum^n_{i=0}2^{i}x^i -\sum^n_{i=0}x^i\\h(n) = 2^n - 1

  • 组合数的证明

    C(n,0)C(m,0)++C(n,m)C(m,m)=C(n+m,m)(1+x)n×(1+(1x)m)=1xm(1+x)m+nC(n,0)C(m,0)+……+C(n,m)C(m,m) = C(n+m,m)\to (1+x)^n\times (1+(\frac 1x) ^m) =\frac {1}{x^m}(1+x)^{m+n}

  • 2n2^n

  • 求导函数

2.3 Fibonacci

对Fibonacci等式两边相加易得

G(x)=x+x2+xG(x)x2+x2G(x)G(x)=x1xx2α=1+52,  β=152hn=αnβn5limn+=15(1+52)nG(x) = x+x^2 + xG(x)-x^2+x^2G(x)\\ G(x) =\frac{x}{1-x-x^2}\\ \alpha = \frac{1+\sqrt 5}{2} , \ \ \beta = \frac{1-\sqrt 5} 2\\ h_n = \frac {\alpha^n-\beta^n}{\sqrt 5} \\lim_{n\to +\infty}= \frac 1{\sqrt 5 }(\frac{1+\sqrt 5}{2})^n

常见等式:

  • F1+F2++Fn=Fn+21F_1+F_2 +……+F_n = F_{n+2}-1

  • F1+F3++F2n1=F2nF_1+F_3 +……+F_{2n-1} = F_{2n}

  • F12++Fn2=FnFn+1F^2_1+……+F^2_n = F_nF_{n+1}

    前两个直接列式即可,第三个变换:F2(n)=F(n)(F(n+1)F(n1))F^2(n) = F(n)(F(n+1)-F(n-1))累加即可

2.4 母函数的性质

  • 现推型:

    bk={0,k<laklklb_k = \begin{cases} 0,\quad \quad \qquad k<l \\ a_{k-l} \qquad k\ge l \end{cases}B(x)=xlA(x)B(x)= x^lA(x)

    bk=ak+lb_k = a_{k+l}B(x)=[A(x)k=0l1akxk]/xlB(x) = [A(x)-\sum^{l-1}_{k=0}a_kx^k]/x_l

    bk=i=0kaib_k = \sum^k_{i=0}a_iB(x)=A(x)1xB(x) = \frac{A(x)}{1-x}

    bk=i=kaib_k = \sum^\infty_{i=k} a_iB(x)=A(1)xA(x)1xB(x) = \frac{A(1)-xA(x)}{1-x}

    bk=kakb_k = ka_kB(x)=xA(x)B(x)=xA^\prime(x)

    bk=ak1+kb_k = \frac{a_k}{1+k}B(x)=1x0xA(x)dxB(x) = \frac{1}{x} \int^x_0A(x)dx

    ck=i=0kaibkic_k =\sum^k_{i=0}a_ib_{k-i}C(x)=A(x)B(x)C(x) = A(x)B(x)

基本列出式子累加即可,(好像最后一个还是记忆一下。。)

2.5 线性常系数齐次递推关系

  • 一般性的线性常系数齐次

    特征多项式:C(x)=xk+c1xk1+....+ck=0C(x) = x^k+ c_1 x^{k-1}+....+c_k = 0

    求得特征方程C(x)=0C(x) = 0

    由此得到的确定序列的G(x)=a0+a1x +a2x2+...G(x) = a_0 + a_1x  +a_2 x^2 + ...的母函数为

    相加有G(x)=P(x)R(x)G(x) = \frac{P(x)}{R(x)}

    且有R(x)=xkC(1x)R(x) =x^kC(\frac 1x)

C(x) 复数域中分解为n个根    由此将

R(x)=(1α1x)k1(1α2x)k2...(1αnx)knR(x) = (1-\alpha_1 x)^{k_1}(1-\alpha_2 x)^{k_2}...(1-\alpha_n x)^{k_n}

定理2-1

有理分式P(x)R(x)\frac{P(x)}{R(x)}有P(x)<R(x),则可化为唯一的部分分式表示

  • 有k重根

    • 数学归纳法,n次时变为更小次项得证
  • 一般化:

    • 特征方程没有重根

      • 可直接求得G( x )
      • 证明:
        • 范德蒙行列式不为零,\because无重根\to唯一
      • 特征方程的根有复数
        • 对于共轭复根,特殊表示形式:
      • 特征方程有重根:
      • 对每一次的n次项各项求和得:

      an=j=1kAjC(n+j1,j1)ana_n =\sum_{j=1}^kA_jC(n+j-1,j-1)a^n\\

  • TODO: 回去证明下

非齐次递推关系

基础做法:

  1. 可解的情况非齐次为bnsnb_ns^n, 依次代入原方程中得到一特解

  2. 非齐次为rnb(n)r^n b(n)则有特解的形式为:

    rn(k0nm+k1nm+1+...+kpnm+p)r^n(k_0n^m+k_1n^{m+1}+...+ k_pn^{m+p})

    1. 其中r为方程的m重根,
    2. b(n)是n的p次多项式,
    3. 若r不是根,令m = 0
    4. 练习
      • 1+2++n=Sn1+2+……+n = S_n

则有Sn3Sn1+3Sn21=0S_n -3S_{n-1}+3S_{n-2}-1= 0

  • 不断降阶有(x1)4=0特征方程为(x-1)^4 = 0的方程

  • 由此方法求得

  • 例题:

    G(x)=1(1x)(1x2)(1x3)G(x) = \frac 1{(1-x)(1-x^2)(1-x^3)}xn中x^n的系数ana_n

    • G(x)=1(1x)3×(1+x)×(1+x+x2)G(x) = \frac{1}{(1-x)^3\times(1+x)\times(1+x+x^2)}

2.7 整数的拆分

  • 整数分解为若干整数的和

    以一题为例更好说明:一克的砝码33枚,22克的砝码44枚,44克的砝码22枚,试问能称出多少种重量,各有多少种称法

    G(x)=(1+x+x2+x3)(1+x2+x4+x6+x8)(1+x4+x8)==1+x+2x3+.....+5x11+...+x19G(x) = (1+x+x^2+x^3)(1+x^2+x^4+x^6+x^8)(1+x^4+x^8) = \\ \qquad \qquad \qquad = 1+x+2x^3+.....+5x^{11}+...+x^{19}

    其中5x115x^{11}表示组合成11共有五种方法

2.8 指数型母函数

仍是以例题为例更好说明:88个元素中a1a_1重复了33次,a2a_2重复了22次,a3a_3重复了33次,从中取rr个组合,其组合数为crc_r.

若求母函数时将xx分为x1,x2,x3x_1,x_2,x_3则可知道对应的四次项(取四个组合)中的来源,从而根据去重原理计算。但是目前有更好的方法:

Ge(x)=(1+1x+x22!+x33!)(1+1x+x22!)(1+1x+x22!+x33!)=1+...+704!x4...G_e(x) = (1+\frac 1x + \frac{x^2}{2!}+\frac{x^3}{3!})(1+\frac 1x + \frac{x^2}{2!})(1+\frac 1x + \frac{x^2}{2!}+\frac{x^3}{3!}) = 1+...+\frac{70}{4!}x^4...

而取44个球的组合数正是系数。巧妙的解决了问题

2.9.1 Ferrers 图像

共轭Ferrers图:

  • def:
    • 将行列互换,关于 y = -x 对称的一对图形成为一对共轭Ferrers图形

推出结论:

  • 整数nn拆分为kk个数的和的拆分 ==nn差分成最大数为k的拆分数相等。

  • 整数nn拆分成最多不超过mm个数的和的拆分数,和nn拆分成最大不超过mm的拆分数相等

    并且有正好拆分成mm个数的母函数为:

    xm(1x)(1x2)(1xm)(1)\frac{x^m}{(1-x)(1-x^2)…(1-x^m)} \qquad \qquad \qquad (1)

  • 整数nn拆分成互不相同的若干奇数的和的拆分书,和nn拆分成自共轭的Ferrers图像的拆分数相等

2.9.2 错排

  • 递推关系:Dn=(n1)(Dn1+Dn2)D_n = (n-1) (D_{n-1}+D_{n-2})

    基本解释:(n1)Dn1=(n-1) D_{n-1}= 先将n1n-1元素错排,选一交换

    (n1)Dn2=(n-1) D_{n-2}= 先选一交换,再对剩余n2n-2元素错排

    更进一步有:

    DnnDn1=(Dn1(n1)Dn2)DnnDn1=(1)nD_n-nD_{n-1} = -(D_{n-1}-(n-1)D_{n-2}) \\D_n-nD_{n-1} = (-1)^n

    其幂次母函数:Ge(x)=D0+D1x1!++Dnxn!G_e(x) = D_0+\frac{D_1x}{1!}+……+\frac{D_nx}{n!}

    并且由DnnDn1=(1)nD_n-nD_{n-1} = (-1)^n

    Ge(x)xGe(x)=exGe(x)=ex1xG_e(x)-xG_e(x) = e^{-x}\\\quad G_e(x) = \frac{e^{-x}}{1-x}

2.9.3 广义二项式定理

(1y)α=1αy+α(α1)2!y2α(α1)(α2)3!y3++(1)rα(α1)(α2)(αr+1)r!yr+(1-y)^\alpha = 1-\alpha y+\frac{\alpha(\alpha-1)}{2!}y^2-\frac{\alpha(\alpha -1)(\alpha-2)}{3!}y^3+……+(-1)^r\frac{\alpha(\alpha -1)(\alpha-2)……(\alpha -r+1)}{r!}y^r+……

例题:序列{1,1×3,1×3×5,}\left\{ 1,1\times3,1\times3\times5,…… \right\}的母函数

2.10 非线性递推关系

  • Stirling 数
  • 双参数

引入:对nn个有区别的球放到两个有区分的盒子里,若要求第一个盒子方kk个球,第二个盒子放nkn-k个球k=0,1,2,,n(k = 0,1,2,……,n)则方案书应为x1+x2n(x_1+x_2)^nx1kx2nkx_1^kx_2^{n-k}的系数C(n,k)C(n,k)

若推广至nn个不同的球放置在mm个不同的盒子中,则不同的方案数为:

(nn1n2nm)=n!n1!n2!nm!\dbinom{n}{n_1n_2\cdots n_m} = \frac{n!}{n_1!n_2!\cdots n_m!}

  • 定理2-10-1(x1+x2++xm)n(x_1+x_2+\cdots+x_m)^n展开式项数=(n+m1n)\dbinom {n+m-1}{n}

  • 证明:

    取m个允许重复的组合的全体,再去重

【定义1:】:[x]n=x(x1)(xn+1)=s(n,0)+s(n,1)x+s(n,2)x2++s(n,n)xn[x]_n = x(x-1)\cdots(x-n+1) = s(n,0)+s(n,1)x+s(n,2)x^2+\cdots+s(n,n)x^n

称其中为第一类SritlingSritling

易知:s(n+1,k)=s(n,k1)ns(n,k)s(n+1,k) = s(n,k-1)-ns(n,k)

【定义2:】n个有区别的球放到m个相同的盒子中,并要求无一空盒,其不同的方案用S(n,m)S(n,m)表示,成为第二类StirlingStirling

  • 性质如下:
    • S(n,0)=0,S(n,1)=1S(n,2)=2n11,S(n,n1)=C(n,2),S(n,n)=1S(n,0) = 0,S(n,1) = 1\\S(n,2) = 2^{n-1}-1,S(n,n-1) = C(n,2),S(n,n) = 1
    • S(n,2)证明:排除掉全部在一个盒子的情况
  • 【定理2-10-3】第二类StirlingStirling数,满足如下递推公式
    • S(n,m)=mS(n1,m)+S(n1,m1),(n>1,m1)S(n,m) = mS(n-1,m)+S(n-1,m-1),\qquad (n>1,m\ge 1)

最后给出StirlingStirling公式的计算公式

S(n,m)=1m!k=0m1(1)k(mk)(mk)nS(n,m) = \frac 1{m!}\sum^{m-1}_{k=0} (-1)^k\dbinom mk(m-k)^n

  • CatalanCatalan

    • def
      • 凸n变形,不相交于n变形内部的对角线,拆分为三角形
      • 其不同拆分数CnC_n为卡特兰数
    1. CatalanCatalan递推关系

      Cn+1=C2Cn+C3Cn1++CnC2C_{n+1} = C_2C_n+C_3C_{n-1}+……+C_nC_2

      • 证明

        v1vn+1v_1 v_{n+1}为固定边,vkv_k为另一顶点,分割成三角形+k边形与nk+2n-k+2边形 ,可知彼此间不重复

      (n3)Cn=n2(C3Cn1++Cn1C3)(n-3)C_n = \frac{n}{2} (C_3C_{n-1}+……+C_{n-1}C_3)

      • 证明

        V1V2VnV_1 V_2 V_n为四边形三个顶点,

        • 等式右侧:以某点的两边为邻边,选一个顶点作对角线,有n3n-3个,可重复执行nn次,且对于无向直线划分的区域相同,÷ C22\div\ C^2_2

        • 等式左侧:计算中,每次nn边形进行剖分后都使用了n3n-3条对角线,也即是重复了n3n-3 。(其实单看了好多遍这一句话不太理解。。)

          但是可以想到的是对于每次剖分需要n3n-3条边,每一次使用对角线划分成的两个区域中,剩余的n4n-4条连线早晚也会被用作“对角线”,因此算上原本的一次便是重复了n3n-3

          自认为终于是理解了QAQ

    2. CatalanCatalan数计算公式

      Cn+1=1n(2n2n1)C_{n+1} = \frac 1n \dbinom{2n-2}{n-1}

      • 证明:

        (C3Cn1++Cn1C3)=Cn+12Cn(n3)Cn=n2(Cn+12Cn)nCn+1=(4n6)Cn(C_3C_{n-1}+……+C_{n-1}C_3) = C_{n+1}-2C_n\\ \Rightarrow (n-3)C_n = \frac n2(C_{n+1}-2C_n) \\ \Rightarrow nC_{n+1} = (4n-6)C_n

        nCn+1=Hn+1nC_{n+1} = H_{n+1}


递推关系与母函数
http://example.com/2021/05/12/Combinatorial Mathematics/CM-递推关系与母函数/
作者
BFlame
发布于
2021年5月12日
许可协议