中文

蝴蝶树的 Horton-Strahler 数

概率论 2026-04-21 v2 数据结构与算法 组合数学

摘要

Horton-Strahler(HS)数是源于水文学与寄存器分配的经典分支复杂度度量,本文将其研究对象限定为蝴蝶树——一种通过块合并操作生成的递归二叉树族。这些树作为蝴蝶排列的二叉搜索树出现,蝴蝶排列构成 N=2nN = 2^n 个元素上的对称群的 22-Sylow 子群,并在并行计算模型和结构化高斯消元中出现。对于作用于两棵独立卡特兰树(各含 mm 个节点)的单步合并操作,我们证明 HS(T1T2)/log2(2m)(\mathcal T_1 \oplus \mathcal T_2)/\log_2(2m) 依概率收敛于 1/21/2,因此经典的卡特兰标度在此受限构造下得以保持。在简单蝴蝶模型中,每一层由相同副本构成并由 nn 位字符串 x\mathbf{x} 编码,HS 数可表示为驱动于独立同分布比特 xjBern(p)x_j \sim \mathrm{Bern}(p) 的显式 88-状态马尔可夫链的加性泛函,并可在 O(n)\mathcal O(n) 时间内由 x\mathbf{x} 计算得出。这导出了一个完整的极限理论,包括强大数定律 HS(TnB)/nμp=pq/(1pq)(\mathcal T_n^B)/n \to \mu_p = pq/(1-pq) 几乎处处成立,以及方差为 σp2=pq(13pq2p2q2)/(1pq)3\sigma_p^2 = pq(1 - 3pq - 2p^2q^2)/(1-pq)^3 的泛函中心极限定理。对于通过递归合并独立子树得到的广义蝴蝶树,增量依赖于不断扩展的边轮廓,且该过程不存在有限状态约化。我们给出了一个 O(N)\mathcal O(N) 算法直接由 (N1)(N-1) 位编码计算 HS 数,刻画了零 HS 类,并结合小 nn 的精确枚举与 n=25n=25 前的蒙特卡洛模拟,支持均匀蝴蝶树中 HS(TnB)/nα0.4450(\mathcal T_n^B)/n \to \alpha \approx 0.4450 依概率收敛,将广义模型严格置于简单蝴蝶极限 1/31/3 与卡特兰极限 1/21/2 之间。

关键词

引用

@article{arxiv.2509.11384,
  title  = {The Horton-Strahler number of butterfly trees},
  author = {John Peca-Medlin},
  journal= {arXiv preprint arXiv:2509.11384},
  year   = {2026}
}