[论文解读] On the Balance of Unrooted Trees

Mareike Fischer, Volkmar Liebscher|arXiv (Cornell University)|Oct 27, 2015

Graph theory and applications参考文献 7被引用 12

一句话总结

本文提出了一套统一框架，用于分析无根系统发育树上的优化问题，通过研究树 τ 上所有非平凡分裂 σ 的泛函 Φ_f(τ) = Σ f(‖σ‖)，其中 f 是分裂大小上的非负函数。研究证明，当 f 为递增函数时，羽状树（caterpillars）使 Φ_f 最小化；当 f 为凹函数或凸函数时，半正则树（semi-regular trees）使平衡性最大化且 Φ_f 最小化，揭示了平衡树形并非唯一，而是依赖于 f 的性质。

ABSTRACT

We solve a class of optimization problems for (phylogenetic) $X$-trees or their shapes. These problems have recently appeared in different contexts, e.g. in the context of the impact of tree shapes on the size of TBR neighborhoods, but so far these problems have not been characterized and solved in a systematic way. In this work we generalize the concept and also present several applications. Moreover, our results give rise to a nice notion of balance for trees. Unsurprisingly, so-called caterpillars are the most unbalanced tree shapes, but it turns out that balanced tree shapes cannot be described so easily as they need not even be unique.

研究动机与目标

• 系统表征在无根系统发育树的所有非平凡分裂上，泛函 Φ_f(τ) = Σ f(‖σ‖) 的极值树形。
• 通过证明平衡树形并非唯一且依赖于 f 的选择，解决长期存在的树平衡定义模糊问题。
• 将多种应用（如 TBR 邻域大小、Gromov 型距离和拓扑指数）统一于单一分析框架之下。
• 基于 f 的单调性、凹性或凸性，建立 Φ_f(τ) 最小值与最大值的必要且充分条件。
• 通过将拓扑指数与树形上的加权度量关联，将该理论推广至更广泛的拓扑指数类别。

提出的方法

• 将 Φ_f(τ) 定义为所有非平凡分裂 σ ∈ Σ*(τ) 上 f(‖σ‖) 的和，其中 ‖σ‖ = min(|A|,|B|) 对于分裂 A|B。
• 使用分裂大小序列 (s_2, ..., s_{⌊n/2⌋}) 表示树形，并分析其帕累托最小值。
• 证明当 f 为递增函数时，Φ_f(τ) 的最小值唯一由羽状树实现；当 f 为凹函数时，最小值由半正则树实现。
• 将该框架应用于已知泛函：f(k) = k(n−k) 对应 TBR 邻域大小，f(k) = k 对应 ℓ¹-Gromov 距离，f(k) = √(k(n−k)) 对应 ℓ²-Gromov 距离。
• 利用 Φ_f(τ) 推导树空间在 Gromov 型距离下的直径上界，表明 diam_D₁ ≤ n²−2n−4 且 diam_D₂ ∼ πn^{3/2}（当 n 为偶数时）。
• 通过定义 w(σ,τ) = g(‖σ‖) 并构造加权度量 d_τ^w，将方法推广至其他指数，得到 Γ^w(τ) = Φ_f(τ)，其中 f(k) = g(k)k(n−k)。

实验结果

研究问题

• RQ1 对于给定函数 f，哪些树形使泛函 Φ_f(τ) = Σ f(‖σ‖) 在所有非平凡分裂上取得最小值或最大值？
• RQ2 为何羽状树与半正则树在不同系统发育泛函中反复作为极值树形出现？
• RQ3 是否可以形式化树平衡的概念，使其依赖于泛函 f，而非假设存在单一通用定义？
• RQ4 对于给定 f，是否存在多个非同构的树形使 Φ_f(τ) 取得相同的最小值，还是最小化者唯一？
• RQ5 现有拓扑指数（如 Wiener 指数或 Γ 指数）能否嵌入此分裂大小泛函的统一框架中？

主要发现

• 对于任意递增函数 f，所有具有 n 个叶的无根系统发育树中，Φ_f(τ) 的最小值唯一由羽状树实现。
• 对于凹函数或凸函数 f，Φ_f(τ) 的最小值由半正则树实现，这些树形不一定唯一，且可能与其他帕累托最小分裂大小序列共存。
• 当 f(k) = k(n−k) 时，泛函 Φ_f(τ) 对应于树 τ 的 TBR 邻域大小，其最大值由羽状树实现。
• 在 ℓ¹-Gromov 距离下，树空间的直径上界为 n²−2n−4（当 n 为偶数时）；在 ℓ²-Gromov 距离下，直径约为 πn^{3/2}。
• 具有 ‖σ‖ = 2 的分裂数（即 'cherries' 数量）由半正则树最大化，由羽状树最小化，确认了它们作为极值树形的角色。
• 该框架可推广至加权拓扑指数：Γ^w(τ) = Φ_f(τ)，其中 f(k) = g(k)k(n−k)，表明许多现有指数均为此统一形式的特例。