精读笔记
Problem Setting
这篇论文实际解决的是 hard-constraint average distance problem 的拓扑结构定理，而不是平均距离泛函本身的定义或 minimizer 存在性。问题可压缩为一句：在长度预算固定、集合必须紧致连通的约束下，最优 free Dirichlet region 是否必然是有限端点、无环、仅三叉分支的树状网络。
真正困难点在于 hard constraint。soft-penalty 下，长度惩罚项直接给出“缩短网络会降低 penalty”的一阶收益；hard-constraint 下，任何删减长度都必须再把预算用回来，否则不能和最优解公平比较。已有条件性拓扑理论已经说明：如果 minimizer 有 atom，则无环、有限端点、三叉分支基本可推出。因此瓶颈不是拓扑论证本身，而是如何在任意维 hard-constraint 情形证明 atom 必然存在。
以前方法卡在两个地方：一是二维 phi(t)=t 的 atom 证明利用了较特殊的平面结构；二是高维 barycentre/self-consistency 路线在 Obrie25 中依赖较强 cost 条件和全局扰动，难以覆盖 phi(t)=t。关键矛盾是：没有 atom 时，局部投影质量看似可以无限分散；但 minimizer 又必须抵抗任何“删一点、加一点”的长度预算交换。论文正是利用这个矛盾逼出 atom。
Motivation
已有路线不够的根本原因是它们把 hard-constraint 的长度预算处理得太全局。全局 rescaling 虽然自然，但在 phi(t)=t 这类原型 cost 下，高阶项难控，尤其靠近集合本身的点会破坏光滑展开。soft-penalty 的证明又不能直接迁移，因为 hard-constraint 缺少 penalty multiplier 的显式一阶项。
作者的核心观察是：barycentre field 不只是 principal curves 语境中的 self-consistency 工具，它实际上记录了 average distance minimizer 是否存在可用一阶下降方向。若 field 非零，额外长度能线性改善目标；若 noncut point 附近可以删除长度保持连通，那么这种线性改善会强迫该点拥有正投影质量，即 atom。
关键缺口因此变成：如何从“没有 atom”推出矛盾。本文的答案是局部化：不要全局缩放整个网络，而是在一个 noncut-neighbourhood 删除长度，在另一个点附近加 cross。这样一阶项由 barycentre field 控制；无 atom 又等价于 field 平凡，于是主要剩下可估计的高阶项。
Core Idea
核心思想是建立一个闭环：barycentre field 非平凡 ⇒ noncut points 是 atoms；存在 atom ⇒ scaling constant 正 ⇒ barycentre field 非平凡；然后证明 minimizer 不可能处在“无 atom / field 平凡”的状态。这个闭环把拓扑问题转成一阶变分场的问题。
和 prior 的本质区别在于，本文没有把 atom 当成一个需要直接从 Voronoi cell 几何中找出的对象，而是把它看作 hard-constraint optimality 的变分后果。noncut point 提供可删除的长度，barycentre field 衡量新增长度的边际收益；如果所有投影质量都是 diffuse 的，那么局部删除造成的损失太小，而 cross insertion 的收益仍由非 1 维测度条件保证，最终违反最优性。
这种建模方式改变了信息组织：过去是“先有 atom，再做拓扑”；这里是“通过一阶场和局部预算交换推出 atom”。这使证明更 generalizable，因为它不依赖二维结构，也不依赖 soft-penalty 的外生 multiplier。
Method
1. 最近点投影与 barycentre field：作者先把集合问题提升到 projection map 上。这样平均距离泛函的扰动可以在 fiber 层面分解，barycentre field 成为每个投影点的平均外力。它解决的是集合变分不可微和投影多值的问题。
2. 一阶近似与高阶项控制：Proposition 2.6 给出一般一阶展开，Proposition 2.18 专门控制局部改动时的 remainder。后者很关键，因为 hard-constraint 证明最终靠微小局部删改，不能只知道一阶项。
3. ambiguous locus 为零：对 minimizer，最近点投影在 mu-a.e. 意义下唯一。这一步不是技术洁癖，而是确保 atom、nu=(pi_Sigma)#mu、barycentre field 都不是选择依赖的伪对象。
4. noncut-neighbourhood 删除：noncut point 的作用是提供可以删掉的小片段，同时保持集合连通。删除长度后得到预算余量，这是 hard-constraint 里模拟 soft-penalty 一阶收益的唯一自然方式。
5. cross insertion：加入 d-dimensional cross 是局部使用回收预算的机制。它对周围质量产生可量化改进，且收益可和 rho-mass 联系起来。论文的关键是把这个收益和删除造成的高阶损失放到同一尺度比较。
Key Insight / Why It Works
最核心贡献不是定义 barycentre field，而是证明它和 atom 在 hard-constraint minimizer 中几乎是同一个结构信号。非平凡 barycentre field 意味着存在一阶下降方向；Theorem 3.3 表明，如果某个 noncut point 没有 atom，那么从它附近回收长度再沿 barycentre 方向花掉长度会改进目标，矛盾。因此 field 非零会强迫所有 noncut points 原子化。
反方向更微妙：atom 会让 scaling constant 为正。直觉上，一个 atom 对应一团质量全部投影到同一点；围绕该点加 cross 能以线性尺度降低其距离成本，因此任何保持最优性的网络必须在全局 rescaling 方向上有正的边际阻力。这把 atom 转成 barycentre 非平凡性。
Theorem 4.9 的有效性来自局部预算交换，而不是 scaling。假设没有 atom，则由等价性得到 barycentre field 平凡，一阶项消失。此时删除 noncut-neighbourhood 的代价主要是二阶/小量；但另一个 noncut point 附近的 rho-mass 由 Proposition 3.2 至少线性存在，cross insertion 产生足够收益。选择 delta 与 epsilon 成比例后，正项压过误差，矛盾成立。
我认为最实质的 insight 是：hard-constraint 的 atom 不是 measure-theoretic accident，而是 optimality under local length exchange 的必然后果。辅助性贡献是高阶项估计和 cross 常数处理；这些是证明闭合所需，但不是概念核心。这里没有 engineering/scaling 成分，所谓 scaling constant 只是数学变分量，不是经验 scaling law。
Relation To Prior Work
这篇论文处在 Buttazzo-Oudet-Stepanov average distance problem、Paolini-Stepanov 拓扑结构理论、以及 principal curves/barycentre field 三条线的交汇处。
与 Buttazzo-Stepanov 2003 的区别：后者在 d=2、phi(t)=t、mu<<H^2 下证明完整拓扑描述，核心依赖二维情形的 atom 机制；本文把 atom 证明推广到任意维，并且处理原型 phi(t)=t。
与 Paolini-Stepanov / Stepanov 2006 的关系：本文不是重新证明完整拓扑分类，而是补上 Stepanov 条件定理所需的 atom 假设。真正新增的信息是 atom existence，而非三叉分支结论本身。
与 Lu-Slepcev soft-penalty 路线的区别：soft-penalty 的长度项直接提供可用下降机制，hard-constraint 没有这个外源项。本文的局部删除-插入机制可以看成在 hard-constraint 内部构造一个“隐式 multiplier”，但它不是简单重写 soft-penalty。
与 Obrie25 / Kobayashi-Hayase-Kim 的区别：barycentre field 已有，但先前 hard-constraint atom 证明对 phi=t^p 的指数范围有限，且依赖更全局的 perturbation。本文的实质创新是把 proof localize，从而覆盖 phi(t)=t 和 p>=2 的任意维情形。
Dataset / Evaluation
这是一篇纯数学论文，没有 dataset、benchmark 或实验 evaluation。它的“评估”应理解为定理覆盖范围与假设强度是否支持主张。
主 claim 是 hard-constraint minimizer 在任意维下满足完整拓扑描述。论文确实在 phi∈C^{1,1}_{loc}、strictly increasing、满足 alpha_2、且 mu(B_epsilon(x))=o(epsilon) 的条件下证明了这一点，并特别覆盖 phi(t)=t 和 phi(t)=t^p, p>=2。对原始 Buttazzo-Oudet-Stepanov 的 phi=t 原型情形，这是强证据。
但 evaluation 也有边界：结果不覆盖 mu 支撑在 1 维集合上的情形，而这不是技术小问题，文中 Remark 3.6 的例子说明 hard-constraint 在 1 维支撑下确实可能出现高阶分支反例。因此 claim 的真实范围是“非 1 维集中测度下的 hard-constraint average distance minimizers”，不是任意 measure。
Limitation
1. measure 假设是核心前提，不是可轻易去掉的 regularity convenience。mu(B_epsilon(x))=o(epsilon) 本质上保证局部 cross 附近的坏质量项足够小，并排除 1 维支撑导致的反例。
2. cost 假设仍有明显上限。C^{1,1}_{loc}、strict monotonicity、alpha_2 支撑了高阶项控制和 Paolini-Stepanov 拓扑定理调用；p<1 被明确排除，1<p<2 的某些边界正则性问题没有完全解决。
3. 拓扑结论强，但几何信息有限。无环、有限端点、三叉分支并不告诉我们边的正则性、夹角、长度分布、唯一性、稳定性，也不描述 ell 变化时网络如何分叉。
4. endpoint 数量的 bound 依赖 Lambda_Sigma 这种 barycentre-field 量，理论上有意义，但文中未充分说明它是否可估、是否 sharp、是否能预测实际 endpoint rate。
5. 方法把问题从“找 atom”转移为“控制局部预算交换的误差”。这次在给定假设下成功闭合，但对 Wasserstein-H^1 或更弱 cost/measure 条件是否仍成立，文中未充分说明。
Takeaway
- 1. hard-constraint average distance problem 的拓扑结构，关键不是直接研究 branching，而是先证明 atom；atom 是整个拓扑定理的门锁。
- 2. barycentre field 的价值被重新定位：它不是 principal curves 的 self-consistency 附属工具，而是 hard-constraint minimizer 的一阶 optimality certificate。
- 3. 局部预算交换比全局 rescaling 更稳健。
- 这个 insight 很可能可迁移到其他长度约束网络优化问题，尤其是存在连通性约束和距离型目标的场景。
一句话总结
这篇论文把 hard-constraint average distance minimizer 的高维拓扑猜想归约并解决为一个 barycentre-field/atom 等价与局部长度预算交换问题，是从条件性拓扑理论走向完整结构定理的关键一步。