[论文解读] Optimal Homologous Cycles, Total Unimodularity, and Linear Programming

Tamal K. Dey, Anil N. Hirani+1 more|arXiv (Cornell University)|Jan 2, 2010

Topological and Geometric Data Analysis参考文献 20被引用 16

一句话总结

该论文证明，在单形复形中，当边界矩阵为全单模时，使用整数系数计算最优同调循环可在多项式时间内求解，这恰好发生在所有相对同调群均为无挠群时。关键贡献在于通过同调条件对全单模性进行了拓扑表征，使得此类复形可使用线性规划而非整数规划求解。

ABSTRACT

Given a simplicial complex with weights on its simplices, and a nontrivial cycle on it, we are interested in finding the cycle with minimal weight which is homologous to the given one. Assuming that the homology is defined with integer coefficients, we show the following : For a finite simplicial complex $K$ of dimension greater than $p$, the boundary matrix $[\partial_{p+1}]$ is totally unimodular if and only if $H_p(L, L_0)$ is torsion-free, for all pure subcomplexes $L_0, L$ in $K$ of dimensions $p$ and $p+1$ respectively, where $L_0$ is a subset of $L$. Because of the total unimodularity of the boundary matrix, we can solve the optimization problem, which is inherently an integer programming problem, as a linear program and obtain integer solution. Thus the problem of finding optimal cycles in a given homology class can be solved in polynomial time. This result is surprising in the backdrop of a recent result which says that the problem is NP-hard under $\mathbb{Z}_2$ coefficients which, being a field, is in general easier to deal with. One consequence of our result, among others, is that one can compute in polynomial time an optimal 2-cycle in a given homology class for any finite simplicial complex embedded in $\mathbb{R}^3$. Our optimization approach can also be used for various related problems, such as finding an optimal chain homologous to a given one when these are not cycles.

研究动机与目标

• 确定在何种拓扑条件下，寻找最优同调循环（即在给定同调类中具有最小加权ℓ¹范数的循环）的问题可在多项式时间内求解。
• 解决一个看似矛盾的现象：该问题在Z₂上为NP难，但在Z上却为多项式时间可解，尽管Z的结构不如域丰富。
• 以相对同调群中无挠性的形式，提供边界矩阵全单模性的拓扑表征。
• 实现对高维复形中最优循环与链的高效计算，特别是针对流形和嵌入R^d中的复形。
• 通过识别解保持整数而无需整数规划的条件，将线性规划的应用扩展到计算拓扑中的整数约束优化问题。

提出的方法

• 通过放松整数约束，将最优同调循环问题（OHCP）表述为线性规划，依赖边界矩阵的全单模性以确保解为整数。
• 证明边界矩阵∂_{p+1}为全单模当且仅当对所有维度分别为p和p+1的纯子复形L₀ ⊂ L，均有H_p(L, L₀)无挠。
• 利用全单模性与线性规划解的整数性之间的等价性，将OHCP的整数规划形式化简为线性规划。
• 将线性规划形式化应用于计算三角剖分d-流形及嵌入R^d中的有限单形复形中的最优(d−1)-循环。
• 引入线性规划的变体，增加约束以强制解的系数属于{−1, 0, 1}，确保结果具有几何意义。
• 在2D环面和3D实心环面上对方法进行实验验证，表明最优循环的系数常为±1。

实验结果

研究问题

• RQ1 在何种拓扑条件下，整数系数下的最优同调循环问题可在多项式时间内求解？
• RQ2 尽管Z不是域，为何该问题在Z上可解而Z₂上为NP难？
• RQ3 边界矩阵的全单模性能否完全通过子复形同调的性质来表征？
• RQ4 当复形具有相对挠性时，线性规划形式化是否仍能产生整数解？若能，其条件为何？
• RQ5 该框架能否扩展至计算非循环的最优链，或最优同调基？

主要发现

• 边界矩阵∂_{p+1}为全单模当且仅当对所有维度分别为p和p+1的纯子复形L₀ ⊂ L，均有H_p(L, L₀)无挠。
• 当边界矩阵为全单模时，最优同调循环问题可作为线性规划求解，且无需整数规划即可获得整数解。
• 在任意三角剖分的紧致可定向d-流形或任意嵌入R^d的有限单形复形中，计算最优(d−1)-循环的问题可在多项式时间内求解。
• 实践中，该方法产生的解其系数属于{−1, 0, 1}，尽管这在理论上无法保证，但可通过线性规划中的附加约束强制实现。
• 即使复形具有相对挠性，线性规划仍可能产生整数解，但此情况并非保证，且取决于输入链。
• 该结果提供了全单模性的拓扑表征，将其与相对同调群中无挠性的存在性联系起来。