Skip to content

解析生成函数

Analytic generating function · Generating function as an analytic function

将具有正收敛半径的生成级数视为复解析函数,以围道积分和奇点结构研究系数。

条目类型
定义

形式陈述 ​

给定数列 (an)n≥0,它的普通生成级数是

A(z)=∑n=0∞anzn.

若这条幂级数有正收敛半径 R>0,便在复圆盘 |z|<R 内定义一个解析函数,也称全纯函数。所谓“解析生成函数”,强调的是使用这个函数的解析性质研究系数;生成级数可以来自递推、直接计数或组合构造,不必先经过某一种特定的符号方法。

Cauchy–Hadamard 公式给出

1R=lim supn→∞|an|1/n,

并采用 1/∞=0、1/0=∞ 的约定。它先确定系数的指数增长尺度,并不自动给出每一项的渐近等价式。

对 0<r<R,沿正向圆周应用 Cauchy 系数公式,得到

an=[zn]A(z)=12πi∫|z|=rA(z)zn+1dz,|an|≤M(r)rn,

其中 M(r)=max|z|=r|A(z)|。右边的估计来自“圆周长度 2πr 乘以被积函数的最大模”,所以无需先计算积分,也能得到系数上界。[1,2]

直觉

形式生成函数把一串系数写进一个可操作的表达式;解析生成函数再让这个表达式成为复平面上的函数。前者可以逐项比较系数,后者还能考察函数在哪些区域有界、在哪里出现极点,以及能否继续延拓。

Cauchy 公式把两层连接起来:一个系数可以从围住原点的整条曲线上恢复。选小圆还是大圆,算出的系数不变,但估计的松紧会改变。增大半径会使 r−n 迅速减小;与此同时,靠近奇点时 M(r) 可能增大。这一拉扯解释了为什么“选择围道”是分析的一部分。

当函数可以延拓到更大的区域,且只有可控的孤立奇点时,可以把积分围道向外移动,分离出奇点的贡献。奇点距原点多远,常决定指数因子;奇点的阶数或局部形状,进一步决定多项式因子和常数。若同一距离上有多个奇点,必须把它们一起计算。分母的零点还可能被分子约去,应先确认它是真正的奇点。

例子与边界

从递推完整走到系数 ​

设 F0=0,F1=1,对 n≥2 有 Fn=Fn−1+Fn−2。逐项乘以 zn 并相加,得到

F(z)−z=zF(z)+z2F(z),F(z)=z1−z−z2.

这一恒等式先在形式幂级数中成立。令

φ=1+52,ψ=1−52,

则 1−z−z2=(1−φz)(1−ψz),并且

F(z)=15(11−φz−11−ψz).

展开两个几何级数便有 Fn=(φn−ψn)/5。解析上,两个极点分别是 1/φ 与 1/ψ=−φ;前者离原点更近,而 |ψ|<1<φ,所以

Fn∼φn5.

这里不仅知道增长“像指数”,还知道了底数、常数,以及被舍去的第二个指数项。最近极点给出的贡献能够成为主项,是因为另一个极点严格更远。

相同半径不决定多项式因子 ​

几何级数 1/(1−z) 的系数恒为 1,而二阶极点给出

1(1−z)2=∑n≥0(n+1)zn.

两者的收敛半径都是 1,但后一组系数多了一个线性因子。收敛半径只控制指数尺度,不能替代对奇点阶数的分析。

两个同模奇点会留下奇偶结构 ​

考察

A(z)=11−z2=12(11−z+11+z).

1 和 −1 都是模为 1 的极点,它们分别贡献 1/2 和 (−1)n/2,因此

an=1+(−1)n2.

偶数项为 1,奇数项为 0。若只看正实轴上的极点,就会错误地预测每项约为 1/2。收敛半径仍然是 1,与 lim sup|an|1/n=1 完全一致;“指数尺度正确”和“逐项主项正确”是不同要求。

形式存在,不代表解析存在 ​

∑n!zn 是合法的形式生成级数,但对任意 z≠0,相邻项绝对值之比为 (n+1)|z|→∞,故收敛半径为零。此时仍可使用形式系数恒等式,却不能围绕原点套用上述解析积分。

同样的对象数 an=n!,换成指数生成函数后得到 ∑anzn/n!=1/(1−z),收敛半径变成 1。OGF 与 EGF 的解析行为可以完全不同;使用 EGF 提取系数后,必须乘回 n! 才是对象数。

反过来,同一个解析函数也可对应不同计数规范:ez=∑zn/n! 作为 OGF 编码 an=1/n!,作为 EGF 却编码 an=1。必须先说明所用编码,再解释系数。

另一端是整函数 ez:它没有有限奇点,而系数是 1/n!。因此“找到最近奇点”不是适用于所有生成函数的完整方案。某些函数还会以整个收敛圆周为自然边界,不能任意向外延拓。

若表达式带平方根或对数,原点附近的幂级数会选定一个局部分支。后续延拓和围道变形必须与这个分支相容;跨越支割线后重新使用原公式,会改变正在分析的函数。

推论与应用

解析生成函数把计数中的代数关系变成增长率问题。主导奇点分析适合有可控边界奇点的函数;孤立极点可用部分分式或留数处理,平方根、对数等奇点则需要相应的奇点分析定理。把一个局部展开转成系数渐近式,还要核查延拓区域、误差项和其他同模奇点,不能仅凭某种熟悉的外形猜答案。[1]

从一个标准奇点读出多项式因子 ​

对实数 α>0、ρ>0,广义二项式展开精确给出

[zn](1−z/ρ)−α=ρ−nΓ(n+α)Γ(α)Γ(n+1)∼ρ−nnα−1Γ(α).

这里 Γ 是把阶乘推广到非整数参数的 Gamma 函数,Γ(k)=(k−1)!。当 α=1,2 时,公式分别回到系数 ρ−n 与 (n+1)ρ−n。一般情形的渐近式由 Gamma 比值或 Stirling 展开得到。

把这条标准公式迁移到更一般的 A(z),需要控制奇点附近的解析延拓及误差项。典型的奇点转移定理要求在穿过收敛圆、但避开奇点向外射线的一类裂口邻域内解析,并控制其他同模长奇点;只有沿实轴的一个近似式,通常不足以逐系数转移。上面的双极点例子说明了为什么全局位置条件不可省略。

从系数上界到算法分析 ​

即使不计算精确主项,也可以优化 Cauchy 上界。若 an≥0 且系数不全为零,则 |A(reit)|≤A(r),从而 an≤A(r)r−n。在可微的内部最优半径处,对其对数求导得到

rA′(r)A(r)=n.

这正是常见的鞍点方程。它来自让上界尽可能小,而不是凭空选择一个特殊点;求得半径仍不等于已经证明精确的鞍点渐近式。

符号方法等组合构造负责给出正确的生成函数,解析工具负责解释系数在规模增大时的行为。两部分应分别核对:再精确的围道计算,也不会纠正前面把 OGF、EGF 或递推初值写错的问题。

参考资料
  • [1] Philippe Flajolet and Robert Sedgewick, Analytic Combinatorics, Cambridge University Press, 2009,Chapters IV–VIII;作者配套站点,其中第 IV 章说明与勘误定位涵盖解析对象、有理函数和主导极点。
  • [2] Herbert S. Wilf, generatingfunctionology, 2nd ed., Academic Press, 1994,Chapter 5;作者提供的合法开放版本入口。在线版与 2006 年第三版是不同版本。
  • [3] Robin Pemantle and Mark C. Wilson, Analytic Combinatorics in Several Variables, Cambridge University Press, 2013,Chapter 1;从单变量系数提取到多变量方法的进一步阅读。
  • [4] NIST Digital Library of Mathematical Functions,§1.10 的广义二项式展开及 §5.11(iii) 的 Gamma 比值渐近,在线参考手册,访问于 2026-09-21。
关系图谱13 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系