拉格朗日反演

介绍

拉格朗日反演是求关于函数方程的幂级数展开系数非常重要的工具，可以用于组合计数函数的系数提取。

形式幂级数

对任意域$F$我们定义其上的形式幂级数为：

$$f(x)=a_0+a_1x+\cdots+a_nx^n+\cdots,\ a_i\in F.$$
记所有的形式幂级数为$F[[x]]$.（从交换代数的观点来说$F[[x]]$实际上是多项式环$F[x]$在理想$(x)$处的完备化。）

容易看到$F[[x]]$显然的构成整环，我们记其商域为$F((x))$，由形如：

$$f(x)=a_{-m}x^{-m}+\cdots+a_{-1}x^{-1}+a_0+\cdots +a_nx^n+\cdots$$ 的元素构成。

拉格朗日反演

定理：若$f(x),g(x)\in F[[x]]$且$f(g(x))=x$，那么

$$[x^n]g(x)=\frac{1}{n}[x^{-1}]\frac{1}{f(x)^n}.$$ 特别的，若 $f(x)=x/\phi(x)$那么 $$[x^n]g(x)=\frac{1}{n}[x^{n-1}]\phi(x)^n.$$ 这里 $[x^n]f(x)$表示取 $f(x)$中 $x^n$的系数。

证明：记$g(x)=\sum_{i\ge 1}b_ix^i$，由形式幂级数性质（无常数项的形式幂级数在复合运算下构成群）有

$$x=g(f(x))=\sum_{i\ge 1}b_if(x)^i.$$ 求导可得 $$1=\sum_{i\ge1}ib_if(x)^{i-1}f'(x).$$ 为了得到 $b_n$，我们两边除以 $f(x)^n$有 $$\frac{1}{f(x)^n}=\sum_{i=1}^{n-1}\frac{ib_i}{i-n}(f(x)^{i-n})'+nb_n\frac{f'(x)}{f(x)}+\sum_{i> n}\frac{ib_i}{i-n}(f(x)^{i-n})'.$$ 我们用算子 $[x^{-1}]$作用两边，容易知道上式右边前后两项系数均为 $0$（幂级数求导后不会出现$x^{-1}$的项），中间相计算如下 $$\begin{align}\frac{f'(x)}{f(x)}&=\frac{a_1+2a_2x+3a_3x^2+\cdots}{a_1x+a_2x^2+\cdots}\\ &=\frac{a_1+2a_2x+3a_3x^2+\cdots}{a_1x}\cdot \frac{1}{1+\left(\frac{a_2}{a_1}x+\frac{a_3}{a_1}x^2+\cdots\right)}\\ &=\left(x^{-1}+\frac{2a_2}{a_1}+\cdots\right)\left(1-z\left(\frac{a_2}{a_1}+\frac{a_3}{a_1}x+\cdots\right)\right) \end{align}.$$ 所以 $$[x^{-1}]\frac{f'(x)}{f(x)}=1.$$ 简单计算即得 $$b_n=\frac{1}{n}[x^{-1}]\frac{1}{f(x)^n}.$$ $\Box$

Reference

[1]. Lagrange Inversion Formula.