[論文レビュー] A Polynomial Time Algorithm to Compute Geodesics in CAT(0) Cubical Complexes

Koyo Hayashi|arXiv (Cornell University)|Oct 26, 2017

Topological and Geometric Data Analysis参考文献 21被引用数 1

ひとこと要約

この論文は、一般次元のCAT(0)立方複体における測地線を計算する最初の多項式時間アルゴリズムを提示する。この手法は、ミラー、オービアン、およびプロヴァンのアルゴリズムをサブルーチンとして用い、反復的ブレークポイント更新法に基づく。計算時間は |P| および log(1/ϵ) の多項式時間で、P は複体の不一致対を用いたポセット表現である。(1+ϵ)-近似測地線が得られる。

ABSTRACT

This paper presents the first polynomial time algorithm to compute geodesics in a CAT(0) cubical complex in general dimension. The algorithm is a simple iterative method to update breakpoints of a path joining two points using Miller, Owen and Provan's algorithm (2015) as a subroutine. Our algorithm is applicable to any CAT(0) space in which geodesics between two close points can be computed, not limited to CAT(0) cubical complexes.

研究の動機と目的

• 一般のCAT(0)立方複体における測地線が多項式時間で計算可能かどうかという未解決問題を解消すること。
• 特定の部分クラス（木空間や2次元複体など）に限定されない、実用的で効率的なアルゴリズムの開発。
• 任意のCAT(0)立方複体における測地線の多項式時間近似スキーム（EPTAS）の提供。
• 測地線が近接する点の間で計算可能な任意のCAT(0)空間において有効な一般枠組みの確立。

提案手法

• CAT(0)立方複体 K を、不一致対（PIP）を用いたポセット P で表現し、[0,1]^|P| への幾何的埋め込みを可能にする。
• 出発点と到達点を結ぶ n 個のブレークポイントを持つ多角形パスを初期化する。
• 現在のパス・セグメントを含む最小のセルへの直交射影を用いて、反復的にブレークポイントを更新する。
• ミラー、オービアン、プロヴァンのアルゴリズムをサブルーチンとして用い、切断されたCAT(0)直角錐空間（頂点のスターバイアス）における測地線を計算する。
• スターバイアスの凸性と直交射影の性質を活用し、各ステップでより短いパスに収束することを保証する。
• 反復回数と精度を、座標のビット複雑度を用いて制限し、|P| および log(1/ϵ) の多項式時間で実行可能であることを保証する。

実験結果

リサーチクエスチョン

• RQ1 任意のCAT(0)立方複体における測地線が多項式時間で計算可能か？
• RQ2 半定値計画法や2次錐計画法に依存しない一般アルゴリズム枠組みは存在するか？
• RQ3 閉形式解が存在しない場合、CAT(0)立方複体における測地線の近似の計算複雑度はいかほどか？
• RQ4 PIPの構造をどのように活用して、効率的な反復的測地線計算アルゴリズムを設計できるか？
• RQ5 このアルゴリズムは、類似した局所的凸性と射影性を持つ他のCAT(0)空間へ拡張可能か？

主な発見

• 提案されたアルゴリズムは、|P| および log(1/ϵ) の多項式時間で、長さが d(p, q) + ϵ 未満のパスを計算可能であり、一般のCAT(0)立方複体における測地線計算の多項式時間解法という未解決問題を解決する。
• このアルゴリズムは単純であり、半定値計画法や2次錐計画法といった複雑な最適化技術を避けており、代わりに反復的更新と直交射影に依存する。
• 正当性は、頂点のスターバイアスの凸性と、点のセルへの直交射影が1次スケルトングラフにおけるゲートと一致することに依存する。
• 2点が距離1未満の距離にある場合、それらの最小セルが交差することを保証しており、収束性と正当性に不可欠である。
• δ-ミッドポイント計算における座標のビット複雑度は O(log(m/δ)) で有界であり、効率的な数値近似を可能にする。
• この手法は、近接する点間の測地線が計算可能な任意のCAT(0)空間に一般化可能であり、立方複体を超えて広く適用可能である。