[論文レビュー] An Approximation Algorithm for Shortest Descending Paths

Mustaq Ahmed, Anna Lubiw|arXiv (Cornell University)|May 9, 2007

Computational Geometry and Mesh Generation参考文献 20被引用数 7

ひとこと要約

本稿では、一般の多面体地形上の最短降下パス（SDP）を計算するための$(1+\epsilon)$-近似アルゴリズムを提示する。このアルゴリズムは、$O(n^2X/\epsilon)$個のステイナー点を用いた新しい離散化手法を採用している。アルゴリズムは、$O(n^3(X/\epsilon)^2)$の事前処理を経て、頂点またはステイナー点に対するクエリ時間は$O(n)$、それ以外の点に対しては$O(nX/\epsilon)$の時間で処理可能であり、地形の幾何的性質を捉える$X = \frac{L}{h} \sec\theta$を用いる。

ABSTRACT

A path from s to t on a polyhedral terrain is descending if the height of a point p never increases while we move p along the path from s to t. No efficient algorithm is known to find a shortest descending path (SDP) from s to t in a polyhedral terrain. We give a simple approximation algorithm that solves the SDP problem on general terrains. Our algorithm discretizes the terrain with O(n^2 X / e) Steiner points so that after an O(n^2 X / e * log(n X /e))-time preprocessing phase for a given vertex s, we can determine a (1+e)-approximate SDP from s to any point v in O(n) time if v is either a vertex of the terrain or a Steiner point, and in O(n X /e) time otherwise. Here n is the size of the terrain, and X is a parameter of the geometry of the terrain.

研究の動機と目的

• 一般の多面体地形上での最短降下パス（SDP）を計算するための効率的でないアルゴリズムの欠如に応えること。
• インフラ整備のための灌漑用水路ルーティングやスキー道計画などの応用分野におけるSDPの実用的近似解を提供すること。
• ステイナー点による地形の離散化を用いることで、SDP計算の計算量を低減しつつ、近似の保証を維持すること。
• 各ソース頂点に対して一度の事前処理フェーズを経て、地形上の任意の点に対して高速なクエリ時間（$O(n)$）を達成すること。

提案手法

• 辺に沿って高さ$z = j\delta$および面固有の高さ$h(x)$にステイナー点を追加することで地形を離散化し、$\delta = \frac{\epsilon h \cos\theta}{4n}$を用いる。
• 地形の頂点とステイナー点を頂点とする有向グラフ$G$を構築し、それらの間の降下セグメントを辺とする。
• ソース$s$から出発する最短降下パス木を、バッシュホックアルゴリズムを用いて$G$上で計算し、すべてのパスが降下的かつ$1+\epsilon$の近似で保証されるようにする。
• 得られた木$T$を用いて、頂点またはステイナー点へのパスを$O(n)$時間でクエリ可能とし、他の点については候補の直前点探索により$O(nX/\epsilon)$時間で処理する。
• 幾何的パラメータ$L$（最長辺）、$h$（頂点から辺までの最小距離）、$\theta$（最大鋭角）を用いて$X = \frac{L}{h} \sec\theta$を定義し、近似品質の上限を定める。
• 任意の真のSDPに対して、ステイナー点を通るパスが存在し、その長さが最適解の$1+\epsilon$倍以内であることを保証する。

実験結果

リサーチクエスチョン

• RQ1 真の多項式時間アルゴリズムが知られていない一般の多面体地形上でも、$(1+\epsilon)$-近似の最短降下パスを効率的に計算できるか？
• RQ2 ステイナー点による地形の離散化は、良い近似降下パスの存在を保つためにどのように行うべきか？
• RQ3 近似品質$\epsilon$、幾何的パラメータ$L, h, \theta$、および結果として得られる時間・空間計算量とのトレードオフは何か？
• RQ4 頂点およびステイナー点に対して$O(n)$のクエリ時間を維持しつつ、任意の点に対しては$O(nX/\epsilon)$時間で処理できるか？
• RQ5 複雑さの上限において、幾何的パラメータ$X = \frac{L}{h} \sec\theta$への依存を軽減することは可能か？

主な発見

• 本アルゴリズムは、一般の多面体地形上での最短降下パス問題に対して$(1+\epsilon)$-近似を達成する。
• 追加されるステイナー点の数は$O(n^2X/\epsilon)$であり、$X = \frac{L}{h} \sec\theta$で地形の幾何的性質を捉える。
• 事前処理時間は$O(n^3(X/\epsilon)^2)$であり、これはソース$s$からの最短パス木の構築コストを支配する。
• クエリ時間は頂点またはステイナー点に対して$O(n)$、それ以外の点に対しては$O(nX/\epsilon)$であり、効率的なパス取得を可能にする。
• 空間計算量は$O(n^2X/\epsilon)$であり、グラフ内の頂点とステイナー点の数に一致する。
• 本手法は単純かつ実用的であり、複雑さの上限における$X$への依存を低減することでさらなる改善が可能である。