[論文レビュー] 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.

研究の動機と目的

• 与えられた同調類内での重み付き1ノルムが最小となるサイクル（最適な同調的サイクル）を求める問題が、どのような位相的条件下で多項式時間で解けるかを特定すること。
• Z₂上でNP困難であるのに対し、Z上で多項式時間で解けるという表面的矛盾を解明すること。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}に制限する追加制約を導入した線形計画法の変種を提案し、幾何学的に意味のある結果を得ることを保証する。
• 2次元トーラスおよび3次元固体アンヌラスを用いた実験的評価により、最適なサイクルが±1の係数を多く得られることを示した。

実験結果

リサーチクエスチョン

• RQ1 整数係数における最適な同調的サイクル問題が、どのような位相的条件下で多項式時間で解けるか。
• RQ2 Zは体ではないにもかかわらず、Z₂上でNP困難であるのに対し、Z上で多項式時間で解けるのはなぜか。
• RQ3 境界行列の完全単模性を、部分複体の同調群の性質のみで特徴づけられるか。
• RQ4 複体に相対的ねじれがある場合でも、線形計画法の定式化が整数解をもたらすことはあり得るが、そのような場合の条件は何か。
• RQ5 この枠組みを、サイクルでない最適なチェインを計算するのにも、あるいは最適な同調基底を計算するのにも拡張可能か。

主な発見

• 境界行列 ∂_{p+1} が完全単模であることは、すべての次元がpおよびp+1である純粋な部分複体 L₀ ⊂ L に対して、H_p(L, L₀) がねじれを持たないことに同値である。
• 境界行列が完全単模である場合、最適な同調的サイクル問題は線形計画法として解くことができ、整数計画法を用いなくても整数解が得られる。
• 任意の三角化された向き付け可能でコンパクトなd次元多様体、またはR^dに埋め込まれた任意の有限単体的複体において、最適な(d−1)-サイクルを計算する問題は、多項式時間で解ける。
• 実際の計算では、解の係数が{−1, 0, 1}に収まることが頻繁に得られるが、理論的には保証されない。この性質は線形計画法の追加制約により強制可能である。
• 複体に相対的ねじれがある場合でも、線形計画法が整数解をもたらすことはあり得るが、これは保証されず、入力チェインに依存する。
• 本結果により、完全単模性の位相的特徴づけが得られ、相対的同調群におけるねじれの不在と関連づけられた。