【施工ing】生成函数与多项式——学习笔记

生成函数大概是一个无穷幂级数形式的函数，我们只关心它的形式，而不会去带入 $x$ 求值。可以看做是多项式，只是带入没有意义。它的一些运算可以对应组合意义，所以能通过它解决一些组合问题。

一般生成函数(OGF)：

$f(x)=a_0+a_1x^1+a_2x^2+a_3x^3+a_4x^4...$

指数型生成函数（EGF），之后会看到为什么要定义这样的函数：

$f(x)=a_0+\frac {a_1} {1!} x^1+\frac {a_2} {2!} x^2+\frac {a_3} {3!} x^3...$

重要的是要把生成函数的代数形式和组合意义对应起来：

考虑一类组合对象组成的集合 $A$ ,给每个 $A$ 中的元素定义一个 “大小”,记为 $|a|$。一个生成函数的系数 $a_n$ 等于大小为 $n$ 的元素的个数，这个生成函数是一类组合对象对应的生成函数。

加法就是组合对象的并。

乘法 $C=A*B$ ，$C$ 中的 $a_n=| \{(x,y)|x \in A, y\in B,|x|+|y|=n\}|$ ，好像叫做笛卡尔积，即是满足体积和为 $n$ 的有序二元组数。

定义 $SEQ(A)$ 是一个由生成函数 $A$ 得到的生成函数。

$$SEQ(A)=\sum_{i=0}^{\infty}A^i$$

对应的组合意义是，$SEQ(A)$ 是 $A$ 中的所有元素的有序组合。

例：$A=\{\text{'}0\text{'},\text{'}1\text{'} \}$, $|0|=|1|=1$，则 $SEQ(A)=\{每种长度的01串\}$

$A=\{1,2,3,... \}$ , $|i|=i$ ,则 $SEQ(A)=\{正整数有序拆分\}$

可以等比数列求和：

$$SEQ(A)=\sum_{i=0}^{\infty}A^i= \frac{1}{1-A}$$

关于指数生成函数（EGF）：

除一个阶乘，就可以对应带标号的组合对象，其他和普通生成函数类型。

定义 $SET(A)=\sum_{i=0}^{\infty} \frac {A^i} {i!}$ , 和 $SEQ(A)$ 不同的是算得的是无序的方案数。

由 $Taylor$ 展开，有 $e^x=1+\frac{x^1}{1!}+\frac{x^2}{2!}+\frac{x^3}{3!}+\frac{x^4}{4!}...$ ，所以 $SET(A)=\sum_{i=0}^{\infty} \frac {A^i} {i!}=e^A$

差不多就这些…

重要的是如何实现多项式的某些运算：

• 多项式求逆：

  ​

• $Newton$ 迭代：

  已知 $g$ , $g(f)=0$ ,可以迭代出 $f$。

  设已知 $f$ 前 $n$ 项，即 $f\equiv f0 \pmod {x^n}$ ，$f-f0 \equiv 0 \pmod{x^n}$

  $Taylor$ 展开得 $0=g(f)=g(f0)+\frac {g'(f0)(f-f0)}{1!}+\frac {g''(f0)(f-f0)^2}{2!}+... \equiv g(f0)+g'(f0)(f-f0) \pmod{x^{2n}}$

  所以 $g(f0)+g'(f0)(f-f0) \equiv 0 \pmod{x^{2n}} \to f=f0- \frac {g(f0)} {g'(f0)} \pmod {x^{2n}}$

  也就是说由前 $n$ 项能得到前 $2n$ 项，每次要多项式求逆，复杂度 $T(n)=T(n/2)+O(n\log n) =O(n\log n)$

• 给定 $A$ 求 $\ln A$ ：

  $f=ln A$ ，$f'=\frac {A'} {A}$ ，就好了。

• 给定 $A$ 求 $e^A$ ：

  $f=e^A$ ,设 $g(x)=\ln x - A$ ，则 $g(f)=0$ ，牛顿迭代即可。

【施工ing……..】