Lagrange 反演

形式 Laurent 级数

我们已经知道形式幂级数环 C[[x]]\mathbb{C}\lbrack\lbrack x\rbrack\rbrack 了,定义形式 Laurent 级数环:

C((x)):={∑k≥Nakxk:N∈Z,ak∈C}\mathbb{C}\left(\left(x\right)\right):=\left\lbrace \sum_{k\geq N}a_kx^k : N\in\mathbb{Z},a_k\in \mathbb{C}\right\rbrace

我们可以仿照形式幂级数的乘法逆元定义来定义 C((x))\mathbb{C}\left(\left(x\right)\right) 上元素的乘法逆元:

若对于 f:=∑k≥Nfkxkf:=\sum_{k\geq N}f_kx^k 且 fN≠0f_N\neq 0 存在 g=∑k≥−Ngkxkg=\sum_{k\geq -N}g_kx^k 满足 fg=1fg=1 那么

gk:={fN−1, if k=−N,−fN−1∑i>Nfigk−i, otherwiseg_k:= \begin{cases} f_N^{-1}, &\text{ if }k=-N\text{,} \\ -f_N^{-1}\sum_{i> N}f_ig_{k-i}, &\text{ otherwise} \end{cases}

与形式幂级数类似的,我们也对非零的 f(x)=∑k≥Nfkxkf(x)=\sum_{k\geq N}f_kx^k 定义:

ord⁡f:=min⁡{k:fk≠0}\operatorname{ord} f:=\min\lbrace k:f_k\neq 0\rbrace

显然对于 g≠0g\neq 0 有

ord⁡(fg)=ord⁡(f)+ord⁡(g)\operatorname{ord} (fg)=\operatorname{ord}(f)+\operatorname{ord}(g)

形式留数

形式留数是形式 Laurent 级数中 x−1x^{-1} 项的系数.记 res⁡f:=[x−1]f\operatorname{res} f:=\lbrack x^{-1}\rbrack f.

引理:对于任何形式 Laurent 级数 ff 有 res⁡f′=0\operatorname{res} f'=0.

证明:考虑形式导数的定义 (xk)′=kxk−1\left(x^k\right)'=kx^{k-1}.

引理:对于任何形式 Laurent 级数 f,gf,g 有 res⁡(f′g)=−res⁡(fg′)\operatorname{res}(f'g)=-\operatorname{res}(fg').

证明:考虑乘法法则 (fg)′=f′g+fg′(fg)'=f'g+fg' 所以 0=res⁡((fg)′)=res⁡(f′g)+res⁡(fg′)0=\operatorname{res}((fg)')=\operatorname{res}(f'g)+\operatorname{res}(fg').

引理:对于形式 Laurent 级数 f(x)≠0f(x)\neq 0 有 res⁡(f′/f)=ord⁡f\operatorname{res}(f'/f)=\operatorname{ord}f.

证明:设 ord⁡f=k\operatorname{ord}f=k 那么

res⁡(f′f)=res⁡(kfkxk−1+⋯fkxk+fk+1xk+1+⋯)=res⁡(kfkx−1+⋯fk+fk+1x+⋯)=k\begin{aligned} \operatorname{res}\left(\frac{f'}{f}\right)&=\operatorname{res}\left(\frac{kf_kx^{k-1}+\cdots}{f_kx^k+f_{k+1}x^{k+1}+\cdots}\right) \\ &=\operatorname{res}\left(\frac{kf_kx^{-1}+\cdots}{f_k+f_{k+1}x+\cdots}\right) \\ &=k \end{aligned}

引理:对于形式 Laurent 级数 ff 和形式幂级数 g≠0g\neq 0 有 res⁡(f)ord⁡(g)=res⁡(f(g)g′)\operatorname{res}(f)\operatorname{ord}(g)=\operatorname{res}(f(g)g').

证明:考虑线性性,我们只需证明 f=xkf=x^k 其中 k∈Zk\in\mathbb{Z} 的情况即可,若 k≠−1k\neq -1 那么

res⁡xk=0res⁡(gkg′)=res⁡(1k+1(gk+1)′)=1k+1res⁡((gk+1)′)=0\begin{aligned} \operatorname{res}x^k&=0 \\ \operatorname{res}(g^kg')&=\operatorname{res}\left(\frac{1}{k+1}\left(g^{k+1}\right)'\right) \\ &=\frac{1}{k+1}\operatorname{res}\left(\left(g^{k+1}\right)'\right) \\ &=0 \end{aligned}

若 k=−1k=-1 那么

res⁡f=res⁡(x−1)=1res⁡(f(g)g′)=res⁡(g′/g)=ord⁡(g)=res⁡(f)ord⁡(g)\begin{aligned} \operatorname{res}f&=\operatorname{res}\left(x^{-1}\right)=1 \\ \operatorname{res}(f(g)g')&=\operatorname{res}(g'/g) \\ &=\operatorname{ord}(g) \\ &=\operatorname{res}(f)\operatorname{ord}(g) \end{aligned}

复合逆

记 A(x)∘B(x):=A(B(x))A(x)\circ B(x):=A(B(x)).

命题:f(x):=∑k≥1fkxkf(x):=\sum_{k\geq 1}f_kx^k 存在复合逆 f⟨−1⟩(x)f^{\langle -1\rangle}(x) 当且仅当 f(0)=0≠f′(0)f(0)=0\neq f'(0),此时 f⟨−1⟩(x)f^{\langle -1\rangle}(x) 是唯一的.进一步说:若 g(x)=∑k≥1gkxkg(x)=\sum_{k\geq 1}g_kx^k 满足 f(g(x))=xf(g(x))=x 或 g(f(x))=xg(f(x))=x 那么 g(x)=f⟨−1⟩(x)g(x)=f^{\langle -1\rangle}(x).

证明:考虑

g(f(x))=g1(f1x+f2x2+f3x3+⋯ )+g2(f1x+f2x2+⋯ )2+g3(f1x+⋯ )3+⋯=g1f1x+(g1f2+g2f12)x2+(g1f3+2g2f1f2+g3f13)x3+⋯\begin{aligned} g(f(x))&=g_1(f_1x+f_2x^2+f_3x^3+\cdots ) \\ &+g_2(f_1x+f_2x^2+\cdots )^2 \\ &+g_3(f_1x+\cdots )^3 \\ &+\cdots \\ &=g_1f_1x+(g_1f_2+g_2f_1^2)x^2+(g_1f_3+2g_2f_1f_2+g_3f_1^3)x^3+\cdots \end{aligned}

因为 g(f(x))=xg(f(x))=x 所以有下面的方程组

{g1f1=1g1f2+g2f12=0g1f3+2g2f1f2+g3f13=0⋮\begin{cases} g_1f_1&=1 \\ g_1f_2+g_2f_1^2&=0 \\ g_1f_3+2g_2f_1f_2+g_3f_1^3&=0 \\ \vdots \end{cases}

我们只能在 f1≠0f_1\neq 0 时才能解出第一个等式,然后依次可以解出 g2,…g_2,\dots.

特别的,考虑 f(h(x))=xf(h(x))=x 那么 g(f(h(x)))=g(x)g(f(h(x)))=g(x),进而 g(x)=g∘f∘h(x)=x∘h(x)=h(x)g(x)=g\circ f\circ h(x)=x\circ h(x)=h(x).

Lagrange 反演公式

令 f(x),g(x)∈C[[x]]f(x),g(x)\in\mathbb{C}\lbrack\lbrack x\rbrack\rbrack 满足 f(g(x))=g(f(x))=xf(g(x))=g(f(x))=x.取 Φ(x)∈C[[x]]\Phi(x)\in\mathbb{C}\lbrack\lbrack x\rbrack\rbrack(或 Φ(x)∈C((x))\Phi(x)\in\mathbb{C}\left(\left(x\right)\right)),那么

[xn]Φ(f(x))=[xn−1]Φ(x)g′(x)g(x)(xg(x))n=[x−1]Φ(x)g′(x)g(x)n+1\begin{aligned} \lbrack x^n\rbrack\Phi(f(x))&=\lbrack x^{n-1}\rbrack\Phi(x)\frac{g'(x)}{g(x)}\left(\frac{x}{g(x)}\right)^n \\ &=\lbrack x^{-1}\rbrack\frac{\Phi(x)g'(x)}{g(x)^{n+1}} \end{aligned}

证明:

[xn]Φ(f(x))=res⁡(Φ(f(x))xn+1)=res⁡(Φ(f(g(x)))g′(x)g(x)n+1)⋅(ord⁡(g(x)))−1=res⁡(Φ(x)g′(x)g(x)n+1)\begin{aligned} \lbrack x^n\rbrack\Phi(f(x))&=\operatorname{res}\left(\frac{\Phi(f(x))}{x^{n+1}}\right) \\ &=\operatorname{res}\left(\frac{\Phi(f(g(x)))g'(x)}{g(x)^{n+1}}\right)\cdot \left(\operatorname{ord}(g(x))\right)^{-1} \\ &=\operatorname{res}\left(\frac{\Phi(x)g'(x)}{g(x)^{n+1}}\right) \end{aligned}

一些读者可能会更加熟悉下面的版本:对于 k∈Z≥0,n∈Z>0k\in\mathbb{Z}_{\geq 0},n\in\mathbb{Z}_{>0} 有

[xn]f(x)k=kn[xn−k](xg(x))n\lbrack x^n\rbrack f(x)^k=\frac{k}{n}\lbrack x^{n-k}\rbrack\left(\frac{x}{g(x)}\right)^n

或者

[xn]Φ(f(x))=1n[xn−1]Φ′(x)(xg(x))n=1n[x−1]Φ′(x)g(x)n\begin{aligned} \lbrack x^n\rbrack \Phi(f(x))&=\frac{1}{n}\lbrack x^{n-1}\rbrack \Phi'(x)\left(\frac{x}{g(x)}\right)^n \\ &=\frac{1}{n}\lbrack x^{-1}\rbrack\frac{\Phi'(x)}{g(x)^n} \end{aligned}

发现

res⁡(Φ′(x)g(x)n−nΦ(x)g′(x)g(x)n+1)=res⁡((Φ(x)g(x)n)′)=0\begin{aligned} \operatorname{res}\left(\frac{\Phi'(x)}{g(x)^n}-n\frac{\Phi(x)g'(x)}{g(x)^{n+1}}\right)&=\operatorname{res}\left(\left(\frac{\Phi(x)}{g(x)^n}\right)'\right) \\ &=0 \end{aligned}

可以通过我们已经证明的部分导出.

参考文献

  1. Richard P. Stanley and Sergey P. Fomin. Enumerative Combinatorics Volume 2 (Edition 1).
  2. Ira M. Gessel. Lagrange Inversion.

本文内容来自 OI Wiki《Lagrange 反演》;作者与修订历史见原文。AIOJ 对排版、站内链接和资源路径做了格式适配。

内容依 CC BY-SA 4.0 与 SATA 许可共享;附加条款以原文声明为准。感谢 OI Wiki 作者与贡献者;请按 SATA 约定前往 OI Wiki 仓库 Star 支持。