形式陈述
奇点分析公理库奇点分析Singularity analysis · Singularity analysis of generating functions在带缺口的复邻域内把主导奇点的局部幂对数展开逐项转移为生成函数系数渐近。中最基本的转移尺度是
广义二项式定理先给出精确系数
再由 Gamma 比值的渐近展开得到
转移定理把这条标准计算推广到一般函数。设 能解析延拓到以 为尖点的 -域,并在该域内一致地满足
其中 为实数,,且暂不取非正整数。则逐项转移给出
若把 Gamma 比值只保留首项,会额外产生 ;当它大于给定余项时,须把比值继续展开到足够阶。
幂乘对数的尺度可由对 求导获得;例如在 避开 Gamma 极点时, 的首项多出 。
直觉
“转移”不是把 随手换成 ,而是把一族已经算清的局部模型当作字典。由Cauchy 系数公式公理库Cauchy 积分定理与公式Cauchy integral theorem · Cauchy integral formula全纯函数在可缩闭路上的积分为零,并可由边界积分重建内部值与所有导数。得到的围道被压到距离 约 的区域后,尺度变换 把局部幂变为 Hankel 围道上的 ;该积分正好等于 。因此 决定指数因子 ,奇性指数决定 ,局部振幅 决定常数。
-解析和一致余项是定理能够逐项工作的理由。只知道实轴上 ,无法排除函数在复方向出现尖峰;只写出前几项而不控制余项,也无法保证被省略部分的系数更小。
例子与边界
取
广义二项式展开给出 。这里 、,且 ,所以转移公式直接得到
这不仅给出增长阶,还正确恢复首个修正项;它与对阶乘使用 Stirling 展开的结果一致。
单一主导点的假设不可暗中省略。对
而言, 与 都是主导奇点。两点的局部贡献相加后表现为奇偶分离:
只转移 附近的一项,会得到错误的非零奇数系数。周期结构应列出全部同模奇点,再按相位 求和;主项抵消时须分别陈述各剩余类。
当 是非正整数时,纯幂可能只是多项式,Gamma 分母的零点提示首项已经消失;若同时出现对数,必须使用相应的对数转移式。自然边界、本性奇点以及随 移动的贡献区也不属于这套固定代数奇点字典。
推论与应用
若 在 附近解析,先把 作 Taylor 展开再乘到奇异尺度上,便能系统地产生任意阶系数渐近。代数生成函数常在隐式方程的临界点形成平方根展开,于是大量树类共享 的形状;有理函数的极点则给出指数乘多项式。
多元生成函数公理库多元生成函数Multivariate generating function · Multivariable generating function用多个形式变量同时记录组合对象的大小与参数,并按方向提取联合系数的生成函数。若可压成一元对角线,有时也能使用本定理。当贡献点随 移动或没有可隔离的主导奇点时,应改用鞍点法公理库鞍点法Saddle-point method · Method of steepest descent选择穿过相位驻点的系数积分轮廓,以局部高斯积分提取主项并控制远离驻点部分的渐近方法。;转移定理并非无条件代换法。
参考资料
- Philippe Flajolet and Andrew Odlyzko, “Singularity Analysis of Generating Functions,” SIAM Journal on Discrete Mathematics 3(2), 1990, pp. 216–240。
- Philippe Flajolet and Robert Sedgewick, Analytic Combinatorics, Cambridge University Press, 2009, Chapter VI, Theorems VI.1–VI.4。
- Peter Henrici, Applied and Computational Complex Analysis, Vol. 2, Wiley, 1977, Chapter 11。