第一次接触生成函数在
OI 中的组合计数题,然后在网上找了许多讲解生成函数的 blog 学习,从胡小兔大佬的
【趣谈生成函数】入门,再到从 cmd 大佬
【多项式计数杂谈】里一览生成函数的全貌,非常感谢这些大佬们让我这个无名之辈窥见这一美丽的数学结构,也了解到许多有趣的组合计数。
网上介绍生成函数的 blog 里,一般直接给出关于生成函数的运算法则以及基本公式,将它们默认为正确,在此基础上探讨生成函数在组合计数、求和式计算中的各种技巧与应用,却并未给出这些运算法则和基本公式的正确性证明。
自然的,在学习过程中遇到的各种天花乱坠技巧使我好奇它们的正确性证明是什么。
于是我也写了一篇更详细 blog 作为学习笔记,但是这篇 blog 纯粹是从数学的角度上写的,对学习生成函数各种应用没有什么帮助,纯粹只是理性愉悦。
况且讲生成函数应用的好 blog 太多太多了,我再去写应用介绍不过是重复造轮子,写的也不如其它 blog 好。
网上介绍生成函数正确性证明的 blog 几乎见不到(因为大家都不关注它吧),也许这篇 blog 能有点意义。
生成函数的水深不见底,这篇 blog 不可能一次性写完,只能慢慢完善它吧。
目前谈及的内容有
- 生成函数的定义、为什么生成函数是无穷求和
- 生成函数的加法运算、乘法运算、除法运算
- 生成函数的取模运算(多项式取模)
- 生成函数的极限
- 无穷求和交换次序
- 生成函数的复合运算、复合逆
- 生成函数求导
- 生成函数泰勒展开
上面这些内容都给了完整的证明,但还有很多生成函数里的内容没有涉及到,将来可能会填下面的坑。
- 生成函数积分
- exp 函数、ln 函数
- 生成函数开根
- 生成函数的分式域、拉格朗日反演
- 多元生成函数
事实上生成函数在组合计数之外的领域也有很多应用,它是一种化简计算的利器,在数学分析里也有各种应用,了解它的数学基础也是一件挺有趣的事情。
因为本人只学过数学分析,没学过抽象代数,写的东西都很初级,遇到啰嗦的地方也可以跳过不看。
语文太差了,这篇 blog 肯定写得很枯燥,请见谅。
如果有什么公式错误或者证明错误之类的可以私信或评论区,我尽量改正。
术语约定
- 自然数集 N={0,1,2,3,...},N+={1,2,3,...}。
- 设 P 是某一数域,如整数域、有理数域、实数域,P∗=P−{0} 代表去掉 0 元后的数域。
- 用 {ai:i∈N} 来表示如下序列
a0,a1,a2,a3,...,ak,...
其中
ai 可以是
P 里的数字,也可以是生成函数,或者其它的东西,总之用序列来表示一组有序对象。
-
本文中的序列均指无穷序列。
-
对于一个数字序列
{ai∈P},若
∃k∈N,使得
∀n>k 有
an=0,则称该序列为有限序列,意思这个无穷序列在某一项后均为
0。
正文
在这篇文章中,首先会抛弃掉生成函数作为 “函数” 的意义,认为它只序列的另一种表示方式,然后再逐步证明这个 “ 序列” 具有的种种 “函数” 性质。
生成函数,也称形式幂级数,多项式(评论区有人认为这有问题,我不太确定,但这应该不影响阅读)。
设有序列
{ai:ai∈P}。
根据序列
{ai},生成函数
f(x) 被定义为
f(x)=a0+a1x+a2x2+a3x3+...+anxn+...
也可以简记为
f(x)=i=0∑∞aixi
在生成函数
f(x) 中,
x 不是未知数,它只是一个占位符,
xi 用来标记它是第几个占位符。
对于
f(x) 第
i 位上的
ai,我们记作
f[i],
f[i] 也被称作
f(x) 的第
i 次项系数,简称
i 次项。
特别的有
f(x) 的
0 次项称作常数项。
于是
f(x)=i=0∑∞f[i]xi
f(x)=f[0]+f[1]x+f[2]x2+f[3]x3+...+f[n]xn+...
∑∞ 在这里的意义并不是无穷项求和,上述定义看似是无穷项的求和,实质上只是一个记号,只是对序列的另一种表达方式,因为目前已经抛弃掉生成函数作为函数的意义,这一点尤为值得注意。
在下文,会逐步证明
∑∞ 代表序列的意义与生成函数无穷项求和的意义是等价的。
为了提取序列中的某一项系数方便,也记
[xi]f(x)=f[i]
这并不是多此一举,比如提取
1+x−2x2+3x3 的二次项系数可以简记为
[x2](1+x−2x2+3x3) 。
对于两个生成函数
f(x) 与
g(x),
f(x)=g(x) 当且仅当
∀i∈N 均满足
f[i]=g[i],这阐明了生成函数相等的概念。
若产生生成函数
f(x) 的序列
{f[i]} 是有限序列,则称
f(x) 是有限多项式、有限幂级数,否则称
f(x) 是无穷多项式、无穷幂级数。
若
f(x) 是有限多项式,它的最高次项为
n,也可以将
f(x) 记为
f(x)=i=0∑nf[i]xi
或者
f(x)=f[0]+f[1]x+f[2]x2+...+f[n]xn
注意,这仍然不意味着上述两种记法中,符号
i=0∑n 或
+ 意味着有限项求和,它们仍然只是对有限序列的另一种表达方式,因为我们还未定义关于生成函数的加法运算。
例如:
1. 除第
0 项为
a∈P 外的项皆为
0 的序列
{a,0,0,0,...},它产生的生成函数为
f(x)=i=0∑0f[i]xi=a
2. 序列
{1,−2,1,−3,0,0,0,...} 产生的生成函数为
f(x)=i=0∑3f[i]xi=1−2x+x2−3x3
生成函数的运算
加法运算
设
f(x) 与
g(x) 是两个生成函数,它们之间进行加法运算的结果也是一个生成函数,记作
f(x)+g(x),满足
[xi](f(x)+g(x))=f[i]+g[i](i∈N)
于是
f(x)+g(x)=i=0∑∞(f[i]+g[i])xi
P 内的数字对加法满足交换律和结合律,因此可以直接推出生成函数对加法运算满足交换律和结合律,也就是说
f(x)+g(x)=g(x)+f(x)
(f(x)+g(x))+h(x)=f(x)+(g(x)+h(x))
根据生成函数的加法运算,我们可以定义对生成函数的求和式。
设
fi(x) (i=1,2,...,n) 是一组生成函数,定义它们的有限项求和
i=1∑nfi(x)=f1(x)+f2(x)+f3(x)+...+fn(x)=j=0∑∞(i=1∑nfi[j])xj
这里仍然强调一下,这并不是说
∑∞ 与
∑n 这两个求和可交换次序,前者是对序列的一种表达方式,后者是对生成函数或普通数字的有限项求和,也就是说对生成函数求和后体现的是下面的一个新序列
{(i=1∑nfi[j]):j∈N}
乘法运算\卷积运算
对于两个生成函数
f(x) 与
g(x),对它们做乘法运算的结果也是一个生成函数,记作
f(x)⋅g(x) 或简写成
f(x)g(x),满足
[xi]f(x)g(x)=j=0∑if[j]g[i−j]
相当于
f(x)g(x)=i=0∑∞(j=0∑if[j]g[i−j])xi
这也就是卷积运算。
与生成函数的加法运算类似,乘法运算也满足交换律和结合律,同时乘法运算对加法运算还满足分配律,也就是说
f(x)g(x)=g(x)f(x)
(f(x)g(x))⋅h(x)=f(x)⋅(g(x)h(x))
(f(x)+g(x))⋅h(x)=f(x)h(x)+g(x)h(x)
(证明待补)
ps : 太懒了不想写啊qwq
定义了乘法运算后顺便可以定义幂次运算。
∀n∈N+,
f(x) 的
n 次幂记作
fn(x) 或
f(x)n,定义为
fn(x)=i=1∏nf(x)
也即
f(x) 自乘
n 次的结果。
特别的,定义
f(x) 的
0 次幂为
1,
f0(x)=1,这说明无论
f(x) 是什么生成函数,
f0(x) 都为
1。
若对
f(x) 的
n 次幂提取
m 次项系数,可以记为
[xm]fn(x)=fn[m]
注意要把
fn[m] 与
f[m]n 加以区分,前者是
fn(x) 的第
m 次项系数,后这是
f[m] 的
n 次幂。
∀a∈P,我们证明
af(x)=i=0∑∞af[i]xi
解释一下,等式左边的
a 是一个生成函数
g(x),满足
g[0]=a, g[i]=0 (i≥1),等式右边的
a 是一个数字。
根据生成函数乘法的定义,
∀i∈N,有
[xi]g(x)f(x)=j=0∑ig[j]f[i−j]=g[0]f[i−0]=af[i]
于是结论成立。
对于两个有限生成函数
f(x),g(x),设它们的最高次分别是
n,m,可以证明
f(x)g(x)=i=0∑nj=0∑mf[i]g[j]xij
即等式左边的生成函数等于右边一堆生成函数的求和。
这分配律的直接推论。
证明
设有一组生成函数
fi(x)=f[i]xi(i=0,...,n) 和一组生成函数
gi(x)=g[i]xi(i=0,...,m)
易知
f(x)=i=0∑nfi(x)g(x)=i=0∑mgi(x)
且
∀i∈{0,...,n},j∈{0,...,m} 均有
fi(x)gj(x)=f[i]g[j]xij
于是
f(x)g(x)=(i=0∑nfi(x))(j=0∑mgj(x))=(i=0∑nfi(x)(j=0∑mgj(x)))=i=0∑nj=0∑mfi(x)gj(x)=i=0∑nj=0∑mf[i]g[j]xij
除法运算\乘法逆元
设有一个生成函数
f(x),定义
f(x) 的乘法逆元是另一个生成函数
g(x),满足
f(x)g(x)=1
解释一下,
f(x) 与
g(x) 的乘积是一个新的生成函数,记这个生成函数为
h(x),则
h(x) 满足
h[0]=1 且
h[i]=0 (i≥1)。
f(x) 的乘法逆元简称乘逆,记作
f(x)1 或
f(x)−1、
f−1(x)。
光定义了乘法逆元还不行,还要探讨一下它是否存在,注意到
1=h[0]=f[0]g[0]
于是
g(x) 存在的一个必要条件是
f[0]=0,反之当
f[0]=0 时必有
g[0]=f[0]1。
我们还可以证明,
f[0]=0 还是
g(x) 存在的充分条件,方法是利用递推公式构造出一个
g(x) 来。
注意到
∀i∈N+ 有
h[i]=0j=0∑ig[j]f[i−j]=0g[i]=−f[0]1j=0∑i−1g[j]f[i−j]
于是我们找到了一个关于
g[i] 的递推公式,也就构造出了
g(x),且这个构造方法产生的
g(x) 唯一。
所以对于某生成函数
f(x),
f(x)1 存在的充要条件是
f[0]=0,若
f(x) 的乘逆存在,则
f(x)1 唯一。
例1
设有生成函数
f(x)=1−x,请求出它的乘法逆元。
解
f[0]=1=0,于是
f(x) 有乘法逆元,记作
g(x),于是
g[0]=f[0]1=1。
∀i∈N+,有
g[i]=−f[0]1j=0∑i−1g[j]f[i−j]=−j=i−1∑i−1g[j]f[i−j]=−g[i−1]f[1]=g[i−1]
其中利用了
f(x) 的
2 次以上的项系数均为
0,且
f[0]=1, f[1]=−1 的性质。
由上面的递推式可以看出
g(x) 的每一项系数都相同,等于
g[0]=1,于是
f(x)1=g(x)=i=0∑∞xi
也就是
1−x1=i=0∑∞xi
这也是 oiers 喜闻乐见的式子。
设
f(x) 与
g(x) 两个存在乘法逆元的多项式,求证
f(x)g(x)1=f(x)1⋅g(x)1
证明
根据生成函数乘法的交换律与结合律有
f(x)1⋅g(x)1⋅f(x)g(x)=(f(x)1⋅f(x))⋅(g(x)1⋅g(x))=1⋅1=1
由生成函数乘法逆元的定义知
f(x)1⋅g(x)1=f(x)g(x)1
这个性质说明了乘逆运算与乘法运算可交换次序。
例2
∀n∈N,求证
(1−x)n+11=k=0∑∞(nn+k)xk
证明
我们使用数学归纳法,当
n=0 的情况已经证明过了。
现在通过
n 证明
n+1 时的结论也成立。
(1−x)n+21=1−x1⋅(1−x)n+11=(i=0∑∞xi)⋅(j=0∑∞(nn+j)xj)=k=0∑∞(i=0∑k(nn+i))xk=k=0∑∞(n+1n+1+k)xk
其中用到了关于组合数的一个公式
i=0∑k(nn+i)=(n+1n+1+k)
证明仍然用数学归纳法。
现在假设对于
k−1 (k≥1) 结论是成立的,于是有
(n+1n+1+k)=(nn+k)+(n+1n+k)=(nn+k)+(n+1n+1+(k−1))=(nn+k)+i=0∑k−1(nn+i)=i=0∑k(nn+i)
证毕。
例3
设有生成函数
f(x)=1+x,
∀n∈N 有
fn(x)=k=0∑n(kn)xk
这是二项式定理的直接推论。
若
g(x) 有乘法逆元,则
f(x) 除
g(x) 的结果记作
g(x)f(x),定义为
g(x)f(x)=f(x)⋅g(x)1
令
h(x)=g(x)f(x),则除法的含义是
h(x) 与
g(x) 的乘积结果为
f(x)。
生成函数取模运算
对于生成函数
f(x),设它的前
n (n≥1) 项(次数
≤n−1 的项)构成的有限生成函数是
f∗(x),定义
f(x) 对
xn 取模的结果为
f∗(x),记作
f∗(x)=f(x)%xn
简而言之,生成函数对
xn 取模就是只保留前
n 项,忽略掉
n 次以上的项。
比如当
f(x)=1−x1 时,
f(x)%x5=1+x+x2+x3+x4。
一般而言,若关于生成函数的等式是在取
xn 模的意义下,即等式两边的前
n 项对应相等,则记作
leftside=rightside(modxn)
容易验证生成函数取模运算满足以下几种性质。
(f(x)+g(x))%xn=f(x)%xn+g(x)%xn
(f(x)⋅g(x))%xn=((f(x)%xn)⋅(g(x)%xn))%xn
f(x)1%xn=f(x)%xn1%xn
证明
记
f(x) 的乘逆为
g(x),注意到构造乘法逆元的递推公式中,得到
g(x) 的前
n 项只需要
f(x) 的前
n 项。
于是
f(x)%xn1 的前
n 项就是
g(x) 的前
n 项。
若记
f∗(x),g∗(x) 分别为
f(x) 和
g(x) 的前
n 项,则一定有
f∗(x)⋅g∗(x)=1(modxn)
极限运算
设
{fn(x)} (n∈N) 是一个生成函数序列,定义生成函数序列的极限是一个生成函数
f(x),记作
f(x)=n→∞limfn(x)
满足
∀m∈N, ∃N∈N,使得任取
k∈[0,m]∧N,
∀n>N 都有
fn(x) 的第
k 次项系数相等。
也即当
n 充分大的时候
f(x) 的前
m+1 项系数都固定了下来。
生成函数序列极限存在还有另一种等价的定义。
f(x) 是
{fn(x)} 的极限等价于
∀m∈N, ∃N∈N+,使得
∀n>N 都有
fn(x)%xm=f(x)%xm
这个对生成函数极限的定义与数学分析中对极限的定义不同,不沿用数学分析中对极限的定义是因为这个极限定义在生成函数这个代数结构中更优良。
生成函数序列的极限不一定存在,若它的极限存在,也称该生成函数序列极限收敛或称它是收敛的生成函数序列。
首先需要阐明一些关于生成函数序列极限的性质。
(1) 极限运算与加法运算可交换次序。
设
{fn(x)}, {gn(x)} 是两个收敛的生成函数序列,它们的极限分别为
f(x),g(x),则有
{fn(x)+gn(x)} 的极限收敛,并且满足
f(x)+g(x)=n→∞lim(fn(x)+gn(x))
证明
由极限的定义知
∀m∈N。
∃N1∈N,使得
∀n>N1 都有
fn(x)%xm=f(x)%xm。
同时
∃N2∈N,使得
∀n>N2 都有
gn(x)%xm=g(x)%xm。
取
N=max{N1,N2},则
∀n>N 有
fn(x)%xm=f(x)%xm, gn(x)%xm=g(x)%xm。
由生成函数取模运算与加法运算交换可次序得
(f(x)+g(x))%xm=f(x)%xm+g(x)%xm=fn(x)%xm+gn%xm=(fn(x)+gn(x))%xm
根据生成函数序列极限的定义可知
n→∞lim(fn(x)+gn(x))=f(x)+g(x)
(2) 极限运算与乘法运算可交换次序。
设
{fn(x)}, {gn(x)} 是两个收敛的生成函数序列,它们的极限分别为
f(x),g(x),则有
{fn(x)⋅gn(x)} 的极限收敛,并且满足
f(x)⋅g(x)=n→∞limfn(x)⋅gn(x)
证明
由极限的定义知
∀m∈N, ∃N∈N,使得
∀n>N 都有
fn(x)%xm=f(x)%xm, gn(x)%xm=g(x)%xm。
由生成函数取模运算可与乘法运算交换次序得
(f(x)⋅g(x))%xm=(f(x)%xm⋅g(x)%xm)%xm=(fn(x)%xm⋅gn(x)%xm)%xm=(fn(x)⋅gn(x))%xm
于是
n→∞limfn(x)⋅gn(x)=f(x)⋅g(x)
(3) 极限运算与除法运算可交换次序。
设
{fn(x)},{gn(x)} 是两个收敛的生成函数序列,它们的极限分别为
f(x),g(x),其中
gn(x) 有乘逆,则有
g(x)f(x)=n→∞limgn(x)fn(x)
证明
首先
∀n∈N 均有
gn[0]=0,
gn(x) 的极限
g(x) 也一定满足
g[0]=0,所以
g(x) 有乘逆,等式的左边有意义。
∀m∈N, ∃N∈N,使得
∀n>N 都有
fn(x)%xm=f(x)%xm, gn(x)%xm=g(x)%xm。
由生成函数取模运算可与除法运算交换次序得
g(x)f(x)%xm=g(x)%xmf(x)%xm%xm=gn(x)%xmfn(x)%xm%xm=gn(x)fn(x)%xm
于是
n→∞limgn(x)fn(x)=g(x)f(x)
(4) 极限运算可与取模运算交换次序。
(证明待补,较容易)
生成函数的本质是序列,
∑∞ 只是对序列的一种表达方式,但在定义好了关于极限运算后,终于可以解释为什么要用无穷求和的符号
∑∞ 来表示生成函数。
在此之前先定义什么是对生成函数的无穷项求和。
与在数学分析中分析数字序列的无穷项求和一样,我们同样用部分和序列的极限来定义无穷项求和。
对于一个生成函数序列
F={fn(x)},定义它的部分和序列也是个生成函数序列
S={sn(x)},满足
sn(x)=i=0∑nfi(x)
即
sn(x) 是序列
F 的前
n+1 项和(从
0 开始计数)。
若
n→∞limsn(x) 存在,记为
s(x),则称
s(x) 为生成函数序列
F 的无穷项求和(或称无穷级数),记为
s(x)=i=0∑∞fi(x)
注意这里的
∑∞ 记号不再指对序列的表达方式,而是货真价实的无穷项求和。
根据生成函数极限的定义,
n→∞limsn(x) 存在的充要条件是
∀m∈N+, ∃N∈N 使得当
n>N 时就有
sn(x)%xm=s(x)%xm
就是说
s(x) 的前
m 项就固定下来了。
所以
∀n>N 有
fn(x)%xm=0,也就是说在
fN(x) 之后的
fn(x),不再对
s(x) 的前
m 项有影响,所以
fN(x) 之后的
fn(x) 的前
m 项为
0。
所以
i=0∑∞fi(x) 存在等价于,
∀m∈N+, ∃N∈N,使得
∀n>N 均有
fn(x)%xm=0。
回到一开始的问题:为什么用
∑∞ 来表示生成函数?
设有某序列
{ai},以及生成函数序列
F={fi(x)} 满足
fi(x)=aixi。
令
F 的部分和序列为
S={si(x)},于是
sn(x)=i=0∑nfi(x)=i=0∑naixi
作生成函数
f(x),满足
f[i]=ai。
∀m∈N+,取
N=m−1,则
∀n>N 均有
fn(x)%xm=0。
根据极限定义可知
S 的极限存在,且
∀i∈N 均有
[xi]s(x)=n→∞lim[xi]sn(x)=ai=f[i]
所以
S 的极限是
f(x)。
i=0∑∞fi(x)=n→∞limsn(x)=f(x)
于是有
i=0∑∞fi(x)=i=0∑∞f[i]xi
至此,我们可以发现等式左边的无穷项求和与等式右边的对序列的表达方式其实是等价的。
f(x)=i=0∑∞f[i]xi
从现在开始,上面这个式子既可以解读为序列
{f[i]} 的一种表达方式,又可以解读为对生成函数序列
{fi(x)=f[i]xi} 的无穷项求和。
生成函数的无穷求和求和也满足线性性质(称为无穷级数的线性性质)。
设有一组生成函数的序列
Fj={fi,j(x)},满足
Fj 的无穷级数存在,即
∀j 有
i∑∞fi,j(x)
存在。
再设
aj(x) 是一组有限生成函数,则有
j∑naj(x)i∑∞fi,j(x)=i∑∞j∑naj(x)fi,j(x)
证明
j∑naj(x)i∑∞fi,j(x)=m→∞limj∑naj(x)i∑mfi,j(x)=m→∞limi∑mj∑naj(x)fi,j(x)
然后
∀m∈N+
∀j,因为
Fj 有极限,所以
∃Nj 使得
∀i>Nj,有
fi,j(x)%xm=0。
取
N=jmax{Nj},则
∀i>N,
∀j 均有
fi,j(x)%xm=0。
于是
(j∑naj(x)fi,j(x))%xm=(j∑naj(x)(fi,j(x)%xm))%xm=(j∑naj(x)⋅0)%xm=0
于是
m→∞limi∑mj∑naj(x)fi,j(x) 存在,为
i∑∞j∑naj(x)fi,j(x)。
同样的我们可以定义关于生成函数的无穷乘积。
(待补充)
复合运算
生成函数的名字里带“函数”两个字,说明它本身带有函数运算的性质,也就是生成函数的复合运算,只有建立了生成函数复合运算的相关理论,才能让生成函数发挥更大的威力。
设
f(x),g(x) 是两个生成函数,定义它们复合运算为一个对生成函数的无穷项求和,记其结果为
h(x),则
h(x)=n→∞limi=0∑nf[i]gi(x)
或者
h(x)=i=0∑∞f[i]gi(x)
h(x) 也记作
f(g(x)) 或
h=f∘g。
来看看这个关于生成函数的无穷级数在什么情况下收敛吧。
若
f(g(x)) 存在,则一定至少满足下面两种情况其中之一。
- f(x) 是有限多项式。
- g(x) 满足 g[0]=0。
第一条就不需要解释了吧,只解释第二条。
若
f(x) 是无穷多项式,且
g[0]=0,则
∀n∈N 均
∃i>n 使得
[x0]f[i]=0,注意到
[x0]gi(x)=0,则
[x0]f[i]gi(x)=0。
于是
∀m∈N+ 均有
f[i]gi(x)%xm=0,所以对
f[i]gi(x) 的无穷求和的极限一定不存在。
反之当两者情况的其中之一满足时,
f(g(x)) 一定存在且唯一。
与普通的函数的复合运算类似,生成函数的复合运算也可以与加法运算、乘法运算、除法运算交换,以下我们一一证明。
(1) 复合运算可与加法运算交换次序,即
(f+g)(h(x))=f(h(x))+g(h(x))
证明
f(h(x))+g(h(x))=i=0∑∞f[i]hi(x)+i=0∑∞g[i]hi(x)=i=0∑∞(f[i]+g[i])hi(x)=(f+g)(h(x))
其中利用了级数运算的线性性质。
(2) 复合运算与乘法运算可交换次序。
(f⋅g)(h(x))=f(h(x))⋅g(h(x))
证明
若
f(x) 与
g(x) 为有限多项式,根据无穷级数的线性性质,则结论当然成立。
否则,要使等式两边有意义则必须满足
h[0]=0。
∀n∈N 我们证明
[xn](f⋅g)(h(x))=[xn]f(h(x))⋅g(h(x))。
因为
h[0]=0 所以
f[i]hi(x)%xn+1=0, g[i]hi(x)%xn+1=0(i>n)。
于是
[xn]f(h(x))⋅g(h(x))=[xn](i=0∑nf[i]hi(x))⋅(j=0∑ng[j]hj(x))=[xn]i=0∑n(j=0∑if[i]g[i−j])hi(x)=[xn](f⋅g)(h(x))
(3) 复合运算与乘逆运算可交换次序。
设
f(x) 存在乘逆,
g(x)=f(x)1,于是
g(h(x))=f(h(x))1
证明
首先有
g⋅f=1,于是
(g⋅f)(h(x))=1。
然后利用复合运算与乘法运算可交换次序
1=(g⋅f)(h(x))=g(h(x))⋅f(h(x))
所以
g(h(x))=f(h(x))1
(4) 在
g[0]=0 的情况下,计算
f(g(x)) 时复合运算可与取模运算交换次序。
记
f∗(x)=f(x)%xn,
g∗(x)=g(x)%xn,则有
f(g(x))=f∗(g∗(x))(modxn)
证明
因为
g[0]=0,当
i≥n 时有
gi(x)%xn=0,于是
f(g(x))%xn=(i=0∑n−1f[i]gi(x))%xn=(i=0∑n−1f[i]g∗i(x))%xn=f∗(g∗(x))%xn
(5) 对于常项为
0 一次项非
0 的生成函数,复合运算满足结合律。
设
f(x),g(x),h(x) 是三个不同的生成函数,满足
f[0]=g[0]=h[0]=0, f[1]=0, g[1]=0, h[1]=0。
证明
设
p(x)=f(g(x)), q(x)=g(h(x))。
任取
n∈N+,在
(modxn) 的情况下,
f(x),g(x),h(x) 均是有限多项式,有限多项式对复合运算满足结合律,于是有
p(h(x))=f(q(x))(modxn)
且
∀m≥n,均满足
∀i∈[0,n−1]∧N 有
[xi]p(h(x))=[xi]f(q(x))(modxm)
而
n 是任意选取的,故
∀k∈N,都能找到
n>k,于是
p(h(x)) 的前
n 项等于
f(q(x)) 的前
n 项,于是
[xk]p(h(x))=[xk]f(q(x))。
也就是说随着
n 的增大,
p(h(x)) 与
f(q(x)) 前面会有越来越多的项被固定,且这些被固定的项逐一对应相等。
故对于无穷幂级数
f(x),g(x),h(x) 而言,
f(g(h(x))) 也满足结合律,即
p(h(x))=f(q(x))
证毕。
在讨论复合运算时,不得不讨论逆元,就像定义完乘法运算后就讨论的乘法逆元一样。
设
f(x) 是某生成函数,如果存在某生成函数
g(x) 使得
g(f(x))=x,则称
g(x) 为
f(x) 的复合逆元,记作
g(x)=f−(x)。
注意要将复合逆元
f−(x) 与乘法逆元
f−1(x) 的记号区分开来。
容易验证
f(x) 存在复合逆元的一个必要条件是
f[0]=0,f[1]=0,也即
f(x) 必须是一个常项为
0 且
1 次项非
0 的生成函数,它的复合逆元
f−(x) 若存在也一定满足
f−[0]=0,f−[1]=0。
以下我们证明
f[0]=0,f[1]=0 也是复合逆元存在的充分条件,并且可以通过递推的方式得到
f(x) 的复合逆元
f−(x)。
为了方便记
g(x)=f−(x)。
根据
g(f(x))=x 可知
[x1]g(f(x))=1[xn]g(f(x))=0 (n>1)
首先有
1=[x1]g(f(x))=g[1]f[1]
故
g[1]=f[1]1
∀n>1,因为
f[0]=0,所以
fi(x)%xn+1=0 (i≥n+1),于是关于
g[i]fi(x) 的求和只需要求和
i≤n 的项,也就是说
0=[xn]g(f(x))=[xn]i=1∑ng[i]fi(x)=g[n]f[1]n+[xn]i=1∑n−1g[i]fi(x)
于是
g[n]=−f[1]n1i=1∑n−1g[i]fi[n]
注意到等式的右边只涉及到
g[1],g[2],...,g[n−1],于是
g(x) 的递推公式就求出来了。
根据此递推公式构造出的
g(x) 是唯一的,故
f(x) 的复合逆元也是存在且唯一的。
尽管
g(f(x))=x,但我们并不知道是否也有
f(g(x))=x,但现在我们可以证明
f(g(x))=x 也是成立的。
证明
记
ε(x)=x 于是
g∘f=ε。
因为
g[0]=0,g[1]=0,可知
g(x) 存在复合逆元
g−(x),
g−∘g=ε。
因为生成函数的复合运算满足结合律,于是
gε∘g(g∘f)∘gg∘(f∘g)g−∘(g∘(f∘g))(g−∘g)∘(f∘g)ε∘(f∘g)f∘g=g=g=g=g=g−∘g=g−∘g=ε=ε
于是
f(g(x))=x。
这说明了
f(x) 与
f−(x) 互为复合逆元,
f−∘f=f∘f−=ε。
例
设
f(x)=1−x1,则
f(x−1) 不存在。
如果直接将
x−1 带入到
1−x1 里会得到
2−x1 的结果,然而这是错的。
f(x) 是无穷多项式,但是
[x0](x−1)=0,所以下面的无穷求和是不存在的。
i=0∑∞(x−1)i
如果硬是要展开则得到
j=0∑∞xji=0∑∞(ji)(−1)i−j
且不说能不能交换求和次序,即使能交换,交换后的结果也不存在。
生成函数求导运算
设生成函数
f(x)=i=0∑∞f[i]xi,则它的导数记作
f′(x),定义为
f′(x)=i=0∑∞(i+1)f[i+1]xi
也可以是下面这种形式
f′(x)=i=1∑∞if[i]xi
注意求导后常数项
f[0] 被丢失了。
生成函数的导数也被记作
(f(x))′。
在利用生成函数求递推式时,求导是一种有力的工具,在介绍生成函数导数的应用之前需要阐明一些性质。
生成函数的求导公式与数学分析中的普通函数的求导公式类似,接下来会一一证明。
(1) 设
f(x),g(x) 是某两个生成函数,
a,b∈P,则有
(af(x)+bg(x))′=af′(x)+bg′(x)
证明
这是显然的。
(2) 设
f(x),g(x) 是某两个生成函数,则有
(f(x)⋅g(x))′=f′(x)g(x)+f(x)g′(x)
证明
等式的左边。
(f(x)⋅g(x))′=i=0∑∞(i+1)(j=0∑i+1f[j]g[i+1−j])xi
等式右边的第一项。
f′(x)g(x)=(i=0∑∞(i+1)f[i+1]xi)⋅(j=0∑∞g[j]xj)=i=0∑∞(j=0∑i(j+1)f[j+1]g[i−j])xi=i=0∑∞(j=1∑i+1jf[j]g[i−(j−1)])xi=i=0∑∞(j=0∑j+1jf[j]g[i+1−j])xi
等式右边的第二项。
f(x)g′(x)=(i=0∑∞f[i]xi)⋅(j=0∑∞(j+1)g[j+1]xj)=i=0∑∞(j=0∑if[j+1](i+1−j)g[i+1−j])xi=i=0∑∞(j=0∑i+1f[j+1](i+1−j)g[i+1−j])xi
将利用无穷级数的线性性质,等式右边的两项合并后就得到等式的左边。
(3) 设
f(x) 是某生成函数,
n∈N+,则有
(fn(x))′=nfn−1(x)f′(x)
证明
使用数学归纳法。
(fn−1(x))′=(n−1)fn−2(x)f′(x)
于是
(fn(x))′=(fn−1(x)⋅f(x))′=(fn−1(x))′f(x)+fn−1(x)f′(x)=(n−1)fn−2(x)f′(x)f(x)+fn−1(x)f′(x)=nfn−1(x)f′(x)
(4) 求导运算与生成函数取模运算的次序交换满足
f′(x)%xn=(f(x)%xn+1)′
证明是显然的。
(5) 极限运算可与求导运算交换次序。
设有收敛的生成函数序列
{fn(x)},其极限为
f(x),则有
f′(x)=n→∞limfn′(x)
证明
由生成函数序列极限的定义可知
∀m∈N+, ∃N∈N,使得
∀n>N 有
fn(x)=f(x)%xm+1。
于是
f′(x)%xm=(f(x)%xm+1)′=(fn(x))′%xm=fn′(x)%xm
于是生成函数序列
{fn′(x)} 收敛于
f′(x)。
推论
生成函数的无穷求和(若存在的话)可与求导运算交换次序,即
(i=0∑∞fi(x))′=i=0∑∞fi′(x)
(6) 设
f(x),g(x) 是某两个生成函数且
f(g(x)) 存在,则有
(f(g(x)))′=f′(g(x))g′(x)
证明
根据
(4) 和
(5) 的结论有
(f(g(x)))′=(i=0∑∞f[i]gi(x))′=i=0∑∞(f[i]gi(x))′=i=1∑∞if[i]gi−1(x)g′(x)=g′(x)i=1∑∞if[i]gi−1(x)=f′(g(x))g′(x)
(7) 设
f(x) 是一个有乘逆的生成函数,则有
(f(x)1)′=−f2(x)f′(x)
证明
首先可以证明,
∀a∈P∗,有
(a+x1)′=−(a+x)21
将等式的左右两边分别展开为无穷求和即可证明。
于是
(f(x)1)′=(f[0]+(f(x)−f[0])1)′=−(f[0]+(f(x)−f[0]))2(f[0]+(f(x)−f[0]))′=−f2(x)f′(x)
推论
设
f(x),g(x) 是两个生成函数,
g[0]=0,于是有
(g(x)f(x))′=g2(x)f′(x)g(x)−f(x)g′(x)
综上可以整理出以下四条求导法则。
(af(x)+bg(x))′(f(x)g(x))′(g(x)f(x))′f(g(x))′=af′(x)+bg′(x)=f′(x)g(x)+f(x)g′(x)=f2(x)f′(x)g(x)−f(x)g′(x)=f′(g(x))g′(x)
泰勒展开公式
设
g(x) 是某生成函数,
a(x),b(x) 是某生成函数,且满足
g(a(x)),g(b(x)) 均存在。
则有
g(b(x))=i=0∑∞g(i)(a(x))i!(b(x)−a(x))i
证明
为了表示方便将
a(x),b(x) 简写成
a,b。
i=0∑∞g(i)(a)i!(b−a)i=i=0∑∞j=0∑i(ji)bjai−j(−1)i−j=j=0∑∞bji=j∑∞(ji)ai−j(−1)i−ji!g(i)(a)(*1)
注意到
i=j∑∞(ji)ai−j(−1)i−ji!g(i)(a)=i=j∑∞(ji)ai−j(−1)i−ji!1k=i∑∞kig[k]xk−i=i=j∑∞k=i∑∞(ji)ai−j(−1)i−j(ik)g[k]xk−i=k=j∑∞g[k]ak−ji=j∑k(ik)(ji)(−1)i−j(*2)
因为
(ik)(ji)=(jk)(i−jk−j)
于是
i=j∑k(ik)(ji)(−1)i−j=(1−1)k−j=[k=j]
代回式
(*2) 得到
i=j∑∞(ji)ai−j(−1)i−ji!g(i)(a)=k=j∑∞g[k]ak−j[k=j]=g[j]
再代回式
(*1) 得到
i=0∑∞g(i)(a)i!(b−a)i=j=0∑∞g[j]bj=g(b)
证毕。