跳转至

Lagrange 反演

形式 Laurent 级数

我们已经知道形式幂级数环 C[[x]] 了,定义形式 Laurent 级数环:

C((x)):={kNakxk:NZ,akC}

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

若对于 f:=kNfkxkfN0 存在 g=kNgkxk 满足 fg=1 那么

gk:={fN1, if k=N,fN1i>Nfigki, otherwise

与形式幂级数类似的,我们也对非零的 f(x)=kNfkxk 定义:

ordf:=min{k:fk0}

显然对于 g0

ord(fg)=ord(f)+ord(g)

形式留数

形式留数是形式 Laurent 级数中 x1 项的系数。记 resf:=[x1]f

引理:对于任何形式 Laurent 级数 fresf=0

证明:考虑形式导数的定义 (xk)=kxk1

引理:对于任何形式 Laurent 级数 f,gres(fg)=res(fg)

证明:考虑乘法法则 (fg)=fg+fg 所以 0=res((fg))=res(fg)+res(fg)

引理:对于形式 Laurent 级数 f(x)0res(f/f)=ordf

证明:设 ordf=k 那么

res(ff)=res(kfkxk1+fkxk+fk+1xk+1+)=res(kfkx1+fk+fk+1x+)=k

引理:对于形式 Laurent 级数 f 和形式幂级数 g0res(f)ord(g)=res(f(g)g)

证明:考虑线性性,我们只需证明 f=xk 其中 kZ 的情况即可,若 k1 那么

resxk=0res(gkg)=res(1k+1(gk+1))=1k+1res((gk+1))=0

k=1 那么

resf=res(x1)=1res(f(g)g)=res(g/g)=ord(g)=res(f)ord(g)

复合逆

A(x)B(x):=A(B(x))

命题f(x):=k1fkxk 存在复合逆 f1(x) 当且仅当 f(0)=0f(0),此时 f1(x) 是唯一的。进一步说:若 g(x)=k1gkxk 满足 f(g(x))=xg(f(x))=x 那么 g(x)=f1(x)

证明:考虑

g(f(x))=g1(f1x+f2x2+f3x3+)+g2(f1x+f2x2+)2+g3(f1x+)3+=g1f1x+(g1f2+g2f12)x2+(g1f3+2g2f1f2+g3f13)x3+

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

{g1f1=1g1f2+g2f12=0g1f3+2g2f1f2+g3f13=0

我们只能在 f10 时才能解出第一个等式,然后依次可以解出 g2,

特别的,考虑 f(h(x))=x 那么 g(f(h(x)))=g(x),进而 g(x)=gfh(x)=xh(x)=h(x)

Lagrange 反演公式

f(x),g(x)C[[x]] 满足 f(g(x))=g(f(x))=x。取 Φ(x)C[[x]](或 Φ(x)C((x))),那么

[xn]Φ(f(x))=[xn1]Φ(x)g(x)g(x)(xg(x))n=[x1]Φ(x)g(x)g(x)n+1

证明

[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)

一些读者可能会更加熟悉下面的版本:对于 kZ0,nZ>0

[xn]f(x)k=kn[xnk](xg(x))n

或者

[xn]Φ(f(x))=1n[xn1]Φ(x)(xg(x))n=1n[x1]Φ(x)g(x)n

发现

res(Φ(x)g(x)nnΦ(x)g(x)g(x)n+1)=res((Φ(x)g(x)n))=0

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

参考文献

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