算法模板（5）：数学（3）：多项式2

生成函数

普通型母函数

普通型母函数的乘法和“多重集的组合问题”有关。

生成函数即母函数，就是离散数学课讲的，给一个无穷序列 { a 0 , a 1 , . . . a n , . . . } \{a_0,a_1,...a_n,...\} {a0​,a1​,...an​,...}，称 $f(x)=a_0+a_{1}x+a_2{x}^2+...+a_n{x}^n+...=\sum _{i=0}^{\infty }{a}_i{x}^i.$

举个例子， { C n 0 , C n 1 , . . . C n k , . . . C n n . . . } \{C_{n}^{0}, C_n^1,...C_n^k,...C_n^n...\} {Cn0​,Cn1​,...Cnk​,...Cnn​...} 的普通母函数的为 $(1+x{)}^n$​.
$$\{\begin{array}{l}\frac{1}{1-x}=1+x+x^2+x^3+\dots \\ \frac{1}{(1-x{)}^{N+1}}={\sum }_{i=0}^{\infty }{C}_{N+i}^N{x}^i\\ \frac{1-x^{N+1}}{1-x}=1+x+x^2+\cdots +x^N\\ \frac{1}{1-x^2}=1+x^2+x^4+\dots \\ \frac{x}{1-x^2}=x+x^3+x^5+\dots \end{array}$$

解释一下，前两个式子可以看作 $\frac{1}{G(x)}=F(x)$，即 $F(x)G(x)=1$​，已知 $G(x)$，然后凑出来 $F(x)$.

再举个例子。假设有 $k$ 个物品，每个物品无限多个，那么选 $n$ 个物品的组合数用隔板法或不定方程解的个数 ( $x_1+x_2+...+x_k=n$ ) 可以求出是 $C_{n+k-1}^{k-1}$. 而对于每个物品的多项式为 $f(x)=1+x+x^2+...+x^n+...$ 那么，这 $k$ 个多项式的卷积为 $\frac{1}{(1-x{)}^k}$，把这个展开为多项式，$x^n$ 的系数，其实就是上述问题的答案，即 $C_{n+k-1}^{k-1}.$ 这个也可以看作是 $\frac{1}{(1-x{)}^{N+1}}={\sum }_{i=0}^{\infty }{C}_{N+i}^N{x}^i$​ 的解释。

生成函数方法一般都是把原问题分成若干子问题。就比如这个选法码，子问题就是每个法码怎么选。而定义的多项式 $g(x)=a_0+a_{1}x+a_2{x}^2+...+a_n{x}^n$，x 的取值通常是 $(-1,1)$，这样可以保证多项式是收敛的。

例子1：

3132. 食物

问组成 $n$ 个物品有多少种方案.

• 承德汉堡：偶数个。
• 可乐：0 个或 1 个。
• 鸡腿：0 个，1 个或 2 个。
• 蜜桃多：奇数个。
• 鸡块：4 的倍数个。
• 包子：0 个，1 个，2 个或 3 个。
• 土豆片炒肉：不超过一个。
• 面包：3 的倍数个。

$$\{\begin{array}{l}1+x^2+x^4+...=\frac{1}{1-x^2}\\ 1+x\\ 1+x+x^2=\frac{1-x^3}{1-x}\\ x+x^3+x^5+...=\frac{x}{1-x^2}\\ 1+x^4+x^8+...=\frac{1}{1-x^4}\\ 1+x+x^2+x^3=\frac{1-x^4}{1-x}\\ 1+x\\ 1+x^3+x^6+...=\frac{1}{1-x^3}\end{array}$$

八个多项式相乘，最后得到的多项式是 $f(x)=\frac{x}{(1-x{)}^4}=\sum _{i=0}^{\infty }{C}_{i+3}^i{x}^i$

答案就是 $C_{n+2}^3$.

例子2：

求在十进制下，$N$ 位的（不允许有前导 $0$），写出来有偶数个 $5$ 的数字的个数。

设答案为 $a_n$​，把“偶数”改成“奇数”的话答案为 $b_n$​。显然地 $a_1=8,b_1=1$​。

对于一个 $n$ ​​位的，有偶数个 $5$ 的数，它只有两种可能：

• 一个 $n-1$ ​位的，有偶数个 $5$ 的数的末尾添了一个非 $5$数。
• 一个 $n-1$​ 位的，有奇数个 $5$ 的数的末尾添了一个 $5$。

$$\{\begin{array}{l}{a}_n=9a_{n-1}+b_{n-1}\\ b_n=9b_{n-1}+a_{n-1}\end{array}$$

$$\{\begin{array}{l}A(x)-9xA(x)-xB(x)=a_1=8\\ B(x)-9xB(x)-xA(x)=b_1=1\end{array}$$

解得
$$\{\begin{array}{l}A(x)=\frac{-71x+8}{(1-8x)(1-10x)}\\ A(x)=\frac{1}{2}(\frac{7}{1-8x}+\frac{9}{1-10x})\end{array}$$

$$A(x)=\frac{1}{2}\sum _{i=0}^{\infty }7\ast (8x{)}^i+9\ast (10x{)}^i$$

即 $a_n=\frac{1}{2}(7\ast 8^n+9\ast 10^n)$

例子3：给定一个长为 $n$​ 的序列 $a$​，求出其 $k$​ 阶差分或前缀和。

序列对应的生成函数是 $a(x)=\sum _{i=0}^{\infty }{a}_i{x}^i$​. 那么求一次前缀和相当于乘 $\sum _{i=0}^{\infty }{x}^i$​​，求一次差分相当于乘 $(1-x)$​.

那么 $k$​​​ 阶前缀和相当于 $a(x)\ast \frac{1}{(1-x{)}^k}=\sum _{i=0}^{\infty }\sum _{j=0}^i{C}_{k+j}^j{a}_{i-j}{x}^i$​​​.

求 $k$​ 阶差分相当于 $a(x)\ast (1-x{)}^k=\sum _{i=0}^{\infty }\sum _{j=0}^i(-1{)}^j{C}_i^j{a}_{i-j}{x}^i$​

生成函数入门

指数型生成函数

指数型母函数的乘法和“多重集的排列问题”有关。

$$\{\begin{array}{l}{e}^x=1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+\dots \\ \frac{e^x+e^{-x}}{2}=1+\frac{x^2}{2!}+\frac{x^4}{4!}+\frac{x^6}{6!}+\dots \\ \frac{e^x-e^{-x}}{2}=x+\frac{x^3}{3!}+\frac{x^5}{5!}+\frac{x^7}{7!}+\dots \end{array}$$​

式子 $1$ 对应的可以选任意多个，式子 $2$​ 对应只能选偶数个，式子 $3$ 对应只能选奇数个.

例题1

用两个1，两个3，一个5，能排成多少个不同的四位数？

展开 $(1+x+\frac{x^2}{2!})(1+x+\frac{x^2}{2!})$​。$\frac{x^4}{4!}$ 的系数就是答案。

可以从多重集的排列数考虑：设 $S$ 有 $k$ 个不同类型的元素，第 $i$ 个类型的元素的重数为 $n_i$ 的重集 $(i=1,2,...,k)$，即 S = { n 1 ⋅ a 1 , n 2 ⋅ a 2 , . . . , n k ⋅ a k } S = \{n_1 \cdot a_1, n_2 \cdot a_2,...,n_k \cdot a_k \} S={n1​⋅a1​,n2​⋅a2​,...,nk​⋅ak​}，则 $S$ 的排列数为
$$\frac{(n_1+n_2+...+n_k)!}{n_1!n_2!...n_k!}$$
例题2

求由 ACTG 组成的 $n$​ 位字符串的个数，其中 A 和 C 必须出现偶数次。

A 和 C 可以用之前的$\frac{e^x+e^{-x}}{2}$​表示，于是得出生成函数为

$$\frac{(e^x+e^{-x}{)}^2{e}^{2x}}{4}$$

括号拆成

$$\frac{e^{4x}+2e^{2x}+1}{4}$$

还是还原回sigma的形式：

$$\frac{1}{4}+\sum _{i=0}^{\infty }\frac{4^i+2\ast 2^i}{4}\ast \frac{x^i}{i!}$$​

$$a_n=\{\begin{array}{l}\frac{4^i+2\ast 2^i}{4}(i\ne 0)\\ 1(i=0)\end{array}$$