【算法讲16：生成函数（理论部分）】

学习资源：

• 南京大学的学习视频 [算法竞赛入门] 生成函数：函数与数列之间的桥梁 (蒋炎岩)
  视频里的内容有些错误，做了初步的修改。
  本文可能有一些内容会错误，恳请斧正。

前置

• 多项式，如 $f(x)=a{x}^2+bx+c$
• 递推/数列，如 $f(n)=f(n-1)+f(n-2)$
• 无穷求和，如 $0.\dot{9}=1$

三个问题的关联

• “无穷多项”的多项式 $=$ 多项式+数列+无穷求和
  $$\begin{array}{rl}f(x)& =a_0+a_{1}x+a_2{x}^2+a_3{x}^3+\cdots \\ & =\sum _{i\ge 0}{a}_i{x}^i\end{array}$$
  这里面有一个无穷数列 $[a_0,a_1,a_2,\cdots ]$
  我们把这个函数 $f(x)$ 通过多项式的无穷求和和数列建立起了联系。

神奇的多项式

• 我们令 $a=[1,1,1,\cdots ]$
  此时 $f(x)=1+x^2+x^3+\cdots$
  我们带入一些值进去，看看等于什么
  $f(1)=1+1+1+\cdots =\infty$
  $f(\frac{1}{2})=1+\frac{1}{2}+\frac{1}{4}+\frac{1}{8}+\cdots$
  当公比 $q\ne 1$ 时这就是一个等比数列求和公式：
  $$f(x)=\frac{1-q^{\infty }}{1-q}=\{\begin{array}{ll}\frac{1}{q-1}& |q|<1\\ \infty & |q|>1\end{array}$$
  换句话说，我们得到了：
  $$\frac{1}{1-x}\leftrightarrow [1,1,1,\cdots ]$$
• 我们之前也学过其他一些函数的多项式展开：
  $$\begin{gathered} sin(x)\leftrightarrow [0,1,0,-\frac{1}{3!},0,\frac{1}{5!},0,-\frac{1}{7!}\cdots ] \\ {e}^x\leftrightarrow [1,\frac{1}{1!},\frac{1}{2!},\frac{1}{3!},\cdots ] \end{gathered}$$
• 挑战：
  $$F(x)=x+x^2+2x^3+3x^4+5x^5+8x^6+\cdots =?$$
  什么样的 $F(x)$ 生成斐波那契数列？
  $$\begin{array}{rl}F(x)& =f(0)x^0+f(1)x^1+f(2)x^2+\cdots \\ xF(x)& =f(0)x^1+f(1)x^2+f(2)x^3+\cdots \\ x^{2}F(x)& =f(0)x^2+f(1)x^3+f(2)x^4+\cdots \end{array}$$
  可以看到，第三行加第二行可以把第一行后面的大部分内容给消掉，得到：
  $$(1-x-x^2)F(x)=x$$
  故，最终答案为：
  $$F(x)=\frac{x}{1-x-x^2}$$
  因此，得到了：
  $$\frac{x}{1-x-x^2}\leftrightarrow [0,1,1,2,3,5,8,\cdots ]$$

数列的$\lceil$生成函数$\rfloor$

• 生成函数：函数与数列之间的桥梁。
  我们可以用函数表示各种数列：
  $$\begin{array}{rl}c\cdot F(x)& \to 数列倍增\\ x\cdot F(x)& \to 数列右移\\ [F(x)-f(0)]/x& \to 数列左移\\ F_1(x)\pm F_2(x)& \to 数列求和/差\\ F^{\prime }(x)& \to 数列乘幂次+左移\\ F_1(x)\cdot F_2(x)& \to 数列卷积\end{array}$$
  举些例子：
  $$\begin{array}{rl}\frac{1}{1-x}& \leftrightarrow [1,1,1,\cdots ]\\ \frac{x}{1-x-x^2}& \leftrightarrow [0,1,1,2,3,5,8,\cdots ]\\ \frac{1}{1-x}+\frac{x}{1-x-x^2}& \leftrightarrow [1,2,2,3,4,6,9,\cdots ,f(n)+1,\cdots ]\\ x(\frac{1}{1-x}+\frac{x}{1-x-x^2})& \leftrightarrow [0,1,2,2,3,\cdots ,f(n+1)+1,\cdots ]\\ \frac{2}{1-x}& \leftrightarrow [2,2,2,\cdots ]\\ \frac{1}{1-ax}& \leftrightarrow [1,a,a^2,a^3,a^4,\cdots ]\end{array}$$
• 我们这里主要关注 数列卷积
  假设我们有两个无穷多项式：
  $$\begin{array}{rl}F(x)& =a_0+a_{1}x+a_2{x}^2+a_3{x}^3+\cdots \\ G(x)& =b_0+b_{1}x+b_2{x}^2+b_3{x}^3+\cdots \end{array}$$
  我们想知道 $F\cdot G$ 是什么？
  $$F(x)\cdot G(x)=a_0{b}_0+(a_1{b}_0+a_0{b}_1)x+\cdots +\left(\sum _{i=0}^n{a}_i{b}_{n-i}\right)x^n+\cdots$$
• 试一试以下的例子
  $$F(x)\cdot \frac{1}{1-x}$$
  我们已经知道了
  $$\begin{array}{rl}F(x)& =a_0+a_{1}x+a_2{x}^2+a_3{x}^3+\cdots \\ \frac{1}{1-x}& =1+x+x^2+x^3+\cdots \end{array}$$
  我们套卷积的式子，得到：
  $$F(x)\cdot \frac{1}{1-x}=a_0+(a_0+a_1)x+\cdots +\left(\sum _{i=0}^n{a}_i\right)x^n+\cdots$$
• 例子：
  求
  $$S_n=\sum _{i=1}^n{i}^2$$
  我们用生成函数做：
  $$\begin{array}{rl}\frac{1}{1-x}& \leftrightarrow [1,1,1,1,1,1\cdots ]\\ \frac{1}{(1-x{)}^2}& \leftrightarrow [1,2,3,4,5,6,7,\cdots ]\\ \frac{1}{(1-x{)}^3}& \leftrightarrow [1,3,6,10,15,\cdots ]三角形数列\\ \frac{x}{(1-x{)}^3}& \leftrightarrow [0,1,3,6,10,15,\cdots ]向右移项\\ \frac{1+x}{(1-x{)}^3}& \leftrightarrow [1,4,9,16,25,\cdots ]相加得到正方形数列\\ \frac{x(1+x)}{(1-x{)}^4}& \leftrightarrow [0,1,5,14,28,53,\cdots ]所求\end{array}$$
  我们得到了一个生成函数：
  $$\frac{x}{(1-x{)}^4}+\frac{x^2}{(1-x{)}^4}$$
  利用超纲知识 ??? 我们知道了 $(1-x{)}^{-4}$ 是 $[C_{n+3}^3]$
  所以最终
  $$S_n=C_{n+2}^3+C_{n+1}^3=\frac{n(n+1)(2n+1)}{6}$$

更多的生成函数

• 指数生成函数：
  $$F(x)=a_0+a_{1}x+a_2\frac{x^2}{2!}+a_3\frac{x^3}{3!}+\cdots$$
• 狄利克雷级数生成函数
  $$F(s)=a_1+\frac{a_2}{2^s}+\frac{a_3}{3^s}+\frac{a_4}{4^s}+\cdots$$
• $F\cdot G$ 的推导
  指数生成函数卷积推导：
  $$\begin{array}{rl}F(x)& =a_0+a_{1}x+a_2\frac{x^2}{2!}+a_3\frac{x^3}{3!}+\cdots \\ G(x)& =b_0+b_{1}x+b_2\frac{x^2}{2!}+b_3\frac{x^3}{3!}+\cdots \\ F(x)\cdot G(x)& =a_0{b}_0+(a_0{b}_1+a_1{b}_0)x+第三项+(\frac{a_0{b}_4}{4!0!}+\frac{a_1{b}_3}{3!1!}+\frac{a_2{b}_2}{2!2!}+\frac{a_3{b}_1}{3!1!}+\frac{a_4{b}_0}{4!0!})x^4+\cdots \\ F(x)\cdot G(x)& =a_0{b}_0+(a_0{b}_1+a_1{b}_0)x+\cdots +\left(\sum _{i=0}^n\frac{a_i{b}_{n-i}\cdot n!}{i!(n-i)!}\right)\frac{x^n}{n!}+\cdots \\ F(x)\cdot G(x)& =a_0{b}_0+(a_0{b}_1+a_1{b}_0)x+\cdots +\left(\sum _{i=0}^n{C}_n^i{a}_i{b}_{n-i}\right)\frac{x^n}{n!}+\cdots \end{array}$$
  狄利克雷级数卷积推导：
  $$\begin{array}{rl}F(s)& =a_1+\frac{a_2}{2^s}+\frac{a_3}{3^s}+\frac{a_4}{4^s}+\cdots \\ G(s)& =b_1+\frac{b_2}{2^s}+\frac{b_3}{3^s}+\frac{b_4}{4^s}+\cdots \\ F(s)\cdot G(s)& =a_1{b}_1+(a_1{b}_2+a_2{b}_1)/2^s+\cdots +(a_1{b}_{18}+a_2{b}_9+a_3{b}_6+a_6{b}_3+a_9{b}_2+a_{18}{b}_1)/18^s+\cdots \\ F(s)\cdot G(s)& =a_1{b}_1+(a_1{b}_2+a_2{b}_1)/2^s+\cdots +\left(\sum _{i\times j=n}{a}_i{b}_j\right)/n^s+\cdots \\ F(s)\cdot G(s)& =\sum _{n\ge 1}\left(\sum _{d|n}{a}_d{b}_{n/d}\right)/n^s\end{array}$$
• $Rieman{n}^{\prime }sZetaFunction$ 黎曼函数
  $$\zeta (s)=1+\frac{1}{2^s}+\frac{1}{3^s}+\cdots$$
  其实就是狄利克雷级数，数列为 $[1,1,1,\cdots ]$
  还有一些性质：
  $$\begin{gathered} {\zeta }^2(s)=d(1)+\frac{d(2)}{2^s}+\frac{d(3)}{3^s}+\cdots \\ \frac{1}{\zeta (s)}=\mu (1)+\frac{\mu (2)}{2^s}+\frac{\mu (3)}{3^s}+\cdots \end{gathered}$$

唯一的例题

• 求长度为 $n$ 的 $01$ 字符串中不可分解字符串的个数
  不可分解指不能写成两个或多个相同字符串的拼接。
• 设 $f(n)$ 是不能分解的个数。
  因为如果一个串有循环节 $2$ 和循环节 $3$ ，那么其最小循环节为 $1$。
  $2^n={\sum }_{d|n}f(d)$
  然后根据莫比乌斯反演去得到 (感觉视频里的式子写错了)
  $$f(n)=\sum _{d|n}\mu (d)2^{\frac{n}{d}}$$