[論文レビュー] Circle separability queries in logarithmic time

Greg Aloupis, Luis Barba+1 more|arXiv (Cornell University)|Mar 28, 2012

Data Management and Algorithms参考文献 16被引用数 3

ひとこと要約

本稿では、n個の点の集合をO(n log n)時間で事前処理することで、与えられた凸m角形を除くすべての点を含む最小円をO(log n + log m)時間で照会するデータ構造を提示する。遠方点ボロノイ図と多角形の辺に対する二分探索を活用することで、最適な円の中心と接点を効率的に特定し、点集合および凸多角形障害物の両方を統一的な照会フレームワークで、最適な事前処理を伴って処理する。

ABSTRACT

Let $P$ be a set of $n$ points in the plane. In this paper we study a new variant of the circular separability problem in which a point set $P$ is preprocessed so that one can quickly answer queries of the following form: Given a geometric object $Q$, report the minimum circle containing $P$ and exluding $Q$. Our data structure can be constructed in $O(n\log n)$ time using O(n) space, and can be used to answer the query when $Q$ is either a circle or a convex $m$-gon in $O(\log n)$ or $O(\log n + \log m)$ time, respectively.

研究の動機と目的

• 事前処理後に、与えられた凸多角形Qを除く点集合Pの最小包含円を高速に照会することを可能にすること。
• 複数の照会が事前処理後に予想されるオンライン版円分離問題に対処すること。
• 幾何データ構造を用いて、従来のO(log n · log m)照会時間の上限を改善し、O(log n + log m)を達成すること。
• 点集合と凸多角形障害物の両方を統一的な照会フレームワークでサポートし、最適な事前処理を実現すること。

提案手法

• 点集合Pを遠方点ボロノイ図V(P)を用いて事前処理し、これは最小包含円の中心c_Pをルートとする根付き木として表現される。
• V(P)上に、辺を二分探索することで最小分離円の中心c'を特定するデータ構造を維持する。
• 各照会用の凸m角形Qに対して、PとQ間の内側共通接線を計算し、Qの境界上での探索区間を定義する。
• Qの多角形境界上で二分探索を実行し、最適な円C'とQの接点φ'を特定する。
• 各ステップで、V(P)に沿った候補の円の中心を評価し、線分と円の交差チェックを用いて、その結果得られる円がQを排除しているかをテストする。
• 二段階の二分探索を適用する：まずV(P)の辺を対象にc'を特定する探索、次に多角形の辺を対象に接点φ'を特定する探索。

実験結果

リサーチクエスチョン

• RQ1 事前処理済みの点集合Pと照会用の凸m角形Qの間の最小分離円を、照会あたりサブ線形時間で特定できるか？
• RQ2 多角形障害物を伴う円分離問題において、事前処理時間と照会時間の最適なトレードオフは何か？
• RQ3 遠方点ボロノイ図の幾何的性質と凸多角形の接線を組み合わせることで、どのように効率的な探索を達成できるか？
• RQ4 この問題において、照会時間をO(log n · log m)からO(log n + log m)に短縮することは可能か？

主な発見

• 提案されたデータ構造は、点集合Pと照会用の凸m角形Qの間の最小分離円をO(log n + log m)時間で照会可能である。
• 事前処理フェーズはO(n log n)時間で実行され、O(n)の記憶領域を要し、効率的かつスケーラブルである。
• アルゴリズムは2つの不変条件を正しく維持する：分離円の中心c'は遠方点ボロノイ図内のパス上にあり、接点φ'は多角形境界の部分チェーン上にある。
• ボロノイ構造と多角形の辺に対する探索を分離することで、従来のO(log n · log m)の上限を著しく改善した。
• アルゴリズムはO(log n + log m)回の反復で終了し、各ステップは事前に計算されたポインタと幾何プリミティブを用いてO(1)時間で実行される。
• 探索区間が1つの辺にまで短縮された場合でも、無限直線ではなく辺自体への直接距離計算により、解がロバストに保たれる。