[论文解读] Full asymptotic expansion for Polya structures

Antoine Genitrini|arXiv (Cornell University)|May 3, 2016

Stochastic processes and statistical mechanics参考文献 12被引用 4

一句话总结

本文提出一种通用方法，通过将生成函数视为Cayley树函数的扰动，推导出无标号非平面树（特别是Pólya树、根身份树和层级结构）的完整渐近展开。利用Weierstrass准备定理与Puiseux展开，计算了至18阶的显式奇点与渐近展开，展示了在树规模较大时相对误差低于10⁻¹⁵的高精度近似。

ABSTRACT

In order to obtain the full asymptotic expansion for Polya trees, i.e. rooted unlabelled and non-plane trees, Flajolet and Sedgewick observed that their specification could be seen as a slight disturbance of the functional equation satisfied by the Cayley tree function. Such an approach highlights the complicated formal expressions with some combinatorial explanation. They initiated this process in their book but they spared the technical part by only exhibiting the first- order approximation. In this paper we exhibit the university of the method and obtain the full asymptotic expansions for several varieties of trees. We then focus on three different varieties of rooted, unlabelled and non-plane trees, Polya trees, rooted identity trees and hierarchies, in order to calculate explicitly their full singular expansions and asymptotic expansions.

研究动机与目标

• 开发一种通用、系统化的方法，用于计算无标号非平面树类别的完整渐近展开。
• 通过解析组合学，将Flajolet与Sedgewick的一阶近似扩展至完整的奇点与渐近展开。
• 为三种关键树类（Pólya树、根身份树与层级结构）计算显式系数。
• 通过展示在树规模≥500时，使用8阶近似可使相对误差低于10⁻¹⁵，验证该方法的精度。
• 建立一个适用于其他组合结构（如序列-并联网络与系统发育树）的框架。

提出的方法

• 将树的生成函数建模为T(z) = ζ(z) exp(T(z))，其中ζ(z)为扰动函数，将其视为对Cayley树函数C(z) = z exp(C(z))的扰动。
• 应用Weierstrass准备定理分析函数方程，并在z = 1/e处主导奇点附近推导出T(z)的完整Puiseux（奇点）展开。
• 使用Faa di Bruno公式与Bell多项式，计算Cayley函数C(z)的Puiseux展开系数。
• 通过引入ζ(z)在主导奇点处的导数，将C(z)的奇点展开推广至T(z)。
• 通过转移定理将Puiseux展开转化为系数渐近性，推导出T_n的完整渐近展开。
• 通过计算并比较每类树的τ_i至18阶系数，对方法进行数值验证。

实验结果

研究问题

• RQ1 能否开发一种通用方法，用于计算无标号非平面树在首阶近似之外的完整渐近展开？
• RQ2 Pólya树、根身份树与层级结构的渐近展开系数行为如何？是否存在可识别的模式？
• RQ3 ζ(z)对Cayley函数的扰动在多大程度上影响奇点与渐近展开的结构？
• RQ4 对于大尺寸树，高阶近似（如8阶）相比首阶近似的精度增益如何？
• RQ5 该框架能否扩展至无根树及其他组合结构（如序列-并联网络或系统发育树）？

主要发现

• 显式推导出Cayley树函数C(z)的完整Puiseux展开，其系数涉及Bell多项式与(1 - ez)^(n/2)的幂次。
• 对于Pólya树，8阶渐近近似可使树规模≥500时的相对误差低于10⁻¹⁵。
• Pólya树与层级结构的τ_i系数在截断至n阶后均被下界控制，而根身份树的τ_i系数则被上界控制。
• 该方法实现高达20位有效数字的数值精度，由于计算限制，大尺寸树的有效近似阶数被限制为8。
• 通过函数方程P(z) = T(z) - ½T²(z) + ½T(z²)，可将该框架扩展至无根树，前提是已知T(z)的完整展开。
• 该方法可泛化应用于其他结构，如序列-并联网络（OEIS A058385–A058387）与系统发育树（OEIS A000311）。