This article generalizes the study of ramified optimal transport with capacity constraint in transport multi-paths by generalizing the $\mathbf{M}_α$ cost to $\mathbf{M}_{α,c}$, which incorporates capacity constraints into the cost function. Equipped with $\mathbf{M}_{α,c}$ cost, we prove the existence of optimal transport path, $\mathbf{M}_{α,c}$ related inequalities, decomposition of any general transport paths, and occurrence of direct line segments in an optimal transport path.
论文ID : 2510.10557标题 : Optimal transport paths with capacity induced cost function作者 : Haotian Sun, Qinglan Xia分类 : math.OC (Optimization and Control)发表时间 : 2025年10月12日论文链接 : https://arxiv.org/abs/2510.10557v1 本文通过将M α M_α M α 成本推广为M α , c M_{α,c} M α , c 来推广对具有容量约束的分支最优传输的研究,该成本将容量约束纳入成本函数中。配备M α , c M_{α,c} M α , c 成本,我们证明了最优传输路径的存在性、M α , c M_{α,c} M α , c 相关不等式、任意一般传输路径的分解,以及最优传输路径中直线段的出现。
该研究要解决的核心问题是分支最优传输中的容量约束问题。传统的Monge-Kantorovich传输问题使用传输映射和传输计划来刻画测度之间的传输,其总成本独立于连接源点和目标点的实际"路径"。
现有方法的局限性 :在先前的工作11 中,作者使用"多路径"来处理容量约束的传输问题,但这种方法存在收敛性问题,如文中Example 1.2所示,满足θ ( x ) ≤ c θ(x) ≤ c θ ( x ) ≤ c 条件的传输路径无法保证收敛。多路径方法的缺陷 :尽管多路径方法解决了收敛问题,但仍存在缺陷。如Remark 1.3和Figure 3所示,存在可容许的传输路径,其每条边上的权重小于等于c,且其边界等于多路径边界的和,但其M α M_α M α 成本小于等于多路径的M α M_α M α 成本。方法创新的必要性 :需要更新容量约束分支传输问题的刻画方法,使得可容许传输路径的集合相比多路径"扩展",同时仍在某种意义上"包含"θ ( x ) ≤ c θ(x) ≤ c θ ( x ) ≤ c 条件。提出了新的成本函数M α , c M_{α,c} M α , c :将传统的M α M_α M α 成本推广为M α , c M_{α,c} M α , c 成本,将容量约束直接纳入成本函数中证明了最优传输路径的存在性 :在新的成本函数下建立了最优传输路径存在性定理建立了M α , c M_{α,c} M α , c 相关不等式 :包括次可加性、单调性等重要性质提供了传输路径的分解理论 :证明了任意传输路径可以分解为权重等于容量整数倍的路径和权重小于容量的路径之和分析了最优路径中直线段的出现 :证明了在最优传输路径中,权重等于容量整数倍的部分往往通过直线段进行传输给定两个等质量的Radon测度μ − μ^- μ − 和μ + μ^+ μ + ,支撑在R m \mathbb{R}^m R m 中的紧集上,参数α ∈ [ 0 , 1 ] α ∈ [0,1] α ∈ [ 0 , 1 ] 和容量约束c > 0 c > 0 c > 0 ,寻找最小化M α , c M_{α,c} M α , c 成本的传输路径T ∈ Path ( μ − , μ + ) T ∈ \text{Path}(μ^-, μ^+) T ∈ Path ( μ − , μ + ) 。
对于任意T = τ ( M , θ ( x ) , ξ ( x ) ) ∈ Path ( μ − , μ + ) T = τ(M, θ(x), ξ(x)) ∈ \text{Path}(μ^-, μ^+) T = τ ( M , θ ( x ) , ξ ( x )) ∈ Path ( μ − , μ + ) ,新的传输成本定义为:
M α , c ( T ) : = c α ⋅ ∫ M ⌊ θ ( x ) c ⌋ + ( θ ( x ) c − ⌊ θ ( x ) c ⌋ ) α d H 1 M_{α,c}(T) := c^α \cdot \int_M \left\lfloor \frac{θ(x)}{c} \right\rfloor + \left(\frac{θ(x)}{c} - \left\lfloor \frac{θ(x)}{c} \right\rfloor\right)^α dH^1 M α , c ( T ) := c α ⋅ ∫ M ⌊ c θ ( x ) ⌋ + ( c θ ( x ) − ⌊ c θ ( x ) ⌋ ) α d H 1
其中⌊ θ ( x ) / c ⌋ \lfloor θ(x)/c \rfloor ⌊ θ ( x ) / c ⌋ 表示不超过θ ( x ) / c θ(x)/c θ ( x ) / c 的最大整数。
当c → ∞ c → ∞ c → ∞ 时:lim c → ∞ M α , c ( T ) = M α ( T ) \lim_{c→∞} M_{α,c}(T) = M_α(T) lim c → ∞ M α , c ( T ) = M α ( T ) (回到传统M α M_α M α 成本) 当c → 0 c → 0 c → 0 时:成本主要由整数部分c α ∫ M ⌊ θ ( x ) / c ⌋ d H 1 c^α \int_M \lfloor θ(x)/c \rfloor dH^1 c α ∫ M ⌊ θ ( x ) / c ⌋ d H 1 决定 定义辅助函数H c , α ( x ) : = ⌊ x / c ⌋ + ( x / c − ⌊ x / c ⌋ ) α H_{c,α}(x) := \lfloor x/c \rfloor + (x/c - \lfloor x/c \rfloor)^α H c , α ( x ) := ⌊ x / c ⌋ + ( x / c − ⌊ x / c ⌋ ) α ,具有以下性质:
当α ∈ ( 0 , 1 ] α ∈ (0,1] α ∈ ( 0 , 1 ] 时,H c , α ( x ) H_{c,α}(x) H c , α ( x ) 严格递增、凹函数且连续 当α = 0 α = 0 α = 0 时,H c , α ( x ) H_{c,α}(x) H c , α ( x ) 递增、分段常数,在整数点有跳跃不连续性 次可加性:H c , α ( x 1 + x 2 ) ≤ H c , α ( x 1 ) + H c , α ( x 2 ) H_{c,α}(x_1 + x_2) ≤ H_{c,α}(x_1) + H_{c,α}(x_2) H c , α ( x 1 + x 2 ) ≤ H c , α ( x 1 ) + H c , α ( x 2 ) 与显式约束θ ( x ) ≤ c θ(x) ≤ c θ ( x ) ≤ c 不同,新方法通过成本函数隐式处理容量约束,避免了收敛性问题。
成本函数中的⌊ θ ( x ) / c ⌋ \lfloor θ(x)/c \rfloor ⌊ θ ( x ) / c ⌋ 项表示"总"权重在每点被细分为若干组件,每个组件权重不超过c c c 。
Proposition 3.6证明了对于多路径T ⃗ ∈ Path c ( μ − , μ + ) \vec{T} ∈ \text{Path}_c(μ^-, μ^+) T ∈ Path c ( μ − , μ + ) ,有M α , c ( T ) ≤ M α ( T ⃗ ) M_{α,c}(T) ≤ M_α(\vec{T}) M α , c ( T ) ≤ M α ( T ) ,其中T = ∑ k = 1 ∞ T k T = \sum_{k=1}^∞ T_k T = ∑ k = 1 ∞ T k 。
定理3.4 :给定支撑在紧集上的等质量Radon测度μ − , μ + μ^-, μ^+ μ − , μ + ,对于任意α ∈ [ 0 , 1 ] α ∈ [0,1] α ∈ [ 0 , 1 ] 和c > 0 c > 0 c > 0 ,存在T ∗ ∈ Path ( μ − , μ + ) T^* ∈ \text{Path}(μ^-, μ^+) T ∗ ∈ Path ( μ − , μ + ) 使得M α , c ( T ∗ ) ≤ M α , c ( T ) M_{α,c}(T^*) ≤ M_{α,c}(T) M α , c ( T ∗ ) ≤ M α , c ( T ) 对所有T ∈ Path ( μ − , μ + ) T ∈ \text{Path}(μ^-, μ^+) T ∈ Path ( μ − , μ + ) 成立。
引理3.3 :对于任意可求长1-电流T T T ,有
M ( T ) ≤ c 1 − α M α , c ( T ) ≤ M ( T ) + c ⋅ size ( T ) M(T) ≤ c^{1-α}M_{α,c}(T) ≤ M(T) + c \cdot \text{size}(T) M ( T ) ≤ c 1 − α M α , c ( T ) ≤ M ( T ) + c ⋅ size ( T )
命题3.5 :对于h ∈ R h ∈ \mathbb{R} h ∈ R ,有
当0 ≤ h ≤ 1 0 ≤ h ≤ 1 0 ≤ h ≤ 1 时:M α , c ( h T ) ≤ h α M α , c ( T ) M_{α,c}(hT) ≤ h^α M_{α,c}(T) M α , c ( h T ) ≤ h α M α , c ( T ) 当h ≥ 1 h ≥ 1 h ≥ 1 时:h α M α , c ( T ) ≤ M α , c ( h T ) h^α M_{α,c}(T) ≤ M_{α,c}(hT) h α M α , c ( T ) ≤ M α , c ( h T ) 定理4.3 :给定传输路径T T T 和其上的任意环R R R ,若满足条件
max x ∈ N ( θ ( x ) c − ⌊ θ ( x ) c ⌋ ) + min x ∈ N ( θ ( x ) c − ⌊ θ ( x ) c ⌋ ) ≤ 1 \max_{x∈N} \left(\frac{θ(x)}{c} - \left\lfloor \frac{θ(x)}{c} \right\rfloor\right) + \min_{x∈N} \left(\frac{θ(x)}{c} - \left\lfloor \frac{θ(x)}{c} \right\rfloor\right) ≤ 1 max x ∈ N ( c θ ( x ) − ⌊ c θ ( x ) ⌋ ) + min x ∈ N ( c θ ( x ) − ⌊ c θ ( x ) ⌋ ) ≤ 1
则可找到无环传输路径T 1 , T 2 T_1, T_2 T 1 , T 2 使得:
∂ T = ∂ ( T 1 + T 2 ) ∂T = ∂(T_1 + T_2) ∂ T = ∂ ( T 1 + T 2 ) θ 1 ( x ) = c ⋅ n ( x ) θ_1(x) = c \cdot n(x) θ 1 ( x ) = c ⋅ n ( x ) ,其中n ( x ) ∈ Z + n(x) ∈ \mathbb{Z}^+ n ( x ) ∈ Z + θ 2 ( x ) ≤ c θ_2(x) ≤ c θ 2 ( x ) ≤ c M α , c ( T 1 + T 2 ) = M α , c ( T 1 ) + M α , c ( T 2 ) ≤ M α , c ( T ) M_{α,c}(T_1 + T_2) = M_{α,c}(T_1) + M_{α,c}(T_2) ≤ M_{α,c}(T) M α , c ( T 1 + T 2 ) = M α , c ( T 1 ) + M α , c ( T 2 ) ≤ M α , c ( T ) 命题4.5 和推论4.6 :在最优传输路径中,对于权重等于容量整数倍的部分,如果存在连接两点的路径,则该路径必须是直线段。
在简单的"2点到1点"传输情况下,文章详细计算了交汇点处的角度。设交汇点为t ′ t' t ′ ,三个方向的单位向量为n ⃗ 1 , n ⃗ 2 , n ⃗ 3 \vec{n}_1, \vec{n}_2, \vec{n}_3 n 1 , n 2 , n 3 ,则最优条件为:
k 1 n ⃗ 1 + k 2 n ⃗ 2 + k 3 n ⃗ 3 = 0 k_1\vec{n}_1 + k_2\vec{n}_2 + k_3\vec{n}_3 = 0 k 1 n 1 + k 2 n 2 + k 3 n 3 = 0
其中k i k_i k i 是对应的修正成本。角度公式为:
cos ( θ 1 ) = k 1 2 + k 3 2 − k 2 2 2 k 1 k 3 \cos(θ_1) = \frac{k_1^2 + k_3^2 - k_2^2}{2k_1k_3} cos ( θ 1 ) = 2 k 1 k 3 k 1 2 + k 3 2 − k 2 2
当c → ∞ c → ∞ c → ∞ 时:角度行为与传统M α M_α M α 成本相同 当c → 0 c → 0 c → 0 时:lim c → 0 k 1 / k 3 = m 1 / ( m 1 + m 2 ) \lim_{c→0} k_1/k_3 = m_1/(m_1 + m_2) lim c → 0 k 1 / k 3 = m 1 / ( m 1 + m 2 ) ,导致所有点共线 本文建立在以下重要工作基础上:
Monge-Kantorovich传输理论 1,7 :经典最优传输理论分支最优传输 8,9 :使用传输路径而非传输映射几何测度论 4,6 :可求长电流和平坦范数理论容量约束传输 11 :作者之前关于多路径方法的工作下半连续性理论 2 :用于证明存在性的关键工具新的成本函数M α , c M_{α,c} M α , c 成功地将容量约束纳入传输成本,避免了显式约束带来的收敛问题 在新成本下最优传输路径存在,且具有良好的数学性质 任意传输路径可以分解为"整数容量"部分和"分数容量"部分 最优路径中的整数容量部分倾向于使用直线段传输 计算复杂性 :新成本函数涉及取整操作,可能增加数值计算的复杂性参数敏感性 :成本函数对容量参数c c c 较为敏感,特别是在c c c 很小时高维推广 :角度计算等几何分析主要在二维情况下进行,高维情况更复杂数值算法开发 :设计高效的数值方法求解M α , c M_{α,c} M α , c 最优化问题应用拓展 :将理论应用于实际的网络设计、物流优化等问题随机情况 :考虑需求和供给具有随机性的情况理论创新性强 :巧妙地通过修改成本函数解决容量约束问题,避免了显式约束的技术困难数学严谨性高 :证明完整,涵盖存在性、唯一性、分解性质等多个方面几何直观清晰 :通过具体的角度计算和例子,提供了良好的几何直观理论体系完整 :从基本定义到高级性质,构建了完整的理论框架实际应用验证不足 :论文主要关注理论发展,缺乏实际应用场景的验证计算方法缺失 :没有提供具体的数值算法来求解优化问题与现有方法的定量比较 :缺乏与多路径方法等现有方法的定量性能比较理论贡献 :为容量约束最优传输理论提供了新的数学框架方法论价值 :通过成本函数设计处理约束的思路具有一般性价值应用潜力 :在网络优化、供应链管理等领域具有应用前景网络设计 :考虑边容量限制的最优网络设计物流优化 :具有运输容量约束的供应链优化城市规划 :考虑道路容量的交通网络规划生物网络 :血管网络、神经网络等生物系统建模论文引用了23篇重要参考文献,主要包括:
最优传输经典著作:Ambrosio等1 ,Villani7 分支最优传输:Xia8,9 几何测度论:Lin & Yang4 ,Simon6 作者前期工作:Xia & Sun10,11 本论文在理论层面取得了重要进展,为容量约束最优传输问题提供了新的数学框架。虽然在实际应用验证方面还有待加强,但其理论创新性和数学严谨性使其成为该领域的重要贡献。