[論文レビュー] Approximating Multiplicatively Weighted Voronoi Diagrams: Efficient Construction with Linear Size

Joachim Gudmundsson, Martin P. Seybold+1 more|arXiv (Cornell University)|Dec 23, 2021

Computational Geometry and Mesh Generation被引用数 1

ひとこと要約

本稿では、最適なサイズとほぼ線形の構築時間を持つε-近似乗法的重み付きボロノイ図（ε-AMWVD）を構築するための新規アルゴリズムを提示する。コア、適応的リファインメント、バイセクタコアセットを導入することで、出力サイズがOd(n log(1/ε)/εd−1) となり、この境界のOD(log(n)/ε(d+5)/2) 以内の構築時間を達成する。これは、Θ(1)d要因を除いて、MWVDに対して初めて最適サイズを達成する構築法であることを証明する。

ABSTRACT

Given a set of $n$ sites from $\mathbb{R}^d$, each having some positive weight factor, the Multiplicatively Weighted Voronoi Diagram is a subdivision of space that associates each cell to the site whose weighted Euclidean distance is minimal for all points in the cell. We give novel approximation algorithms that output a cube-based subdivision such that the weighted distance of a point with respect to the associated site is at most $(1+\varepsilon)$ times the minimum weighted distance, for any fixed parameter $\varepsilon \in (0,1)$. The diagram size is $O_d(n \log(1/\varepsilon)/\varepsilon^{d-1})$ and the construction time is within an $O_D(\log(n)/\varepsilon^{(d+5)/2})$-factor of the size bound. We also prove a matching lower bound for the size, showing that the proposed method is the first to achieve \emph{optimal size}, up to $Θ(1)^d$-factors. In particular, the obscure $\log(1/\varepsilon)$ factor is unavoidable. As a by-product, we obtain a factor $d^{O(d)}$ improvement in size for the unweighted case and $O(d \log(n) + d^2 \log(1/\varepsilon))$ point-location time in the subdivision, improving the known query bound by one $d$-factor. The key ingredients of our approximation algorithms are the study of convex regions that we call cores, an adaptive refinement algorithm to obtain optimal size, and a novel notion of \emph{bisector coresets}, which may be of independent interest. In particular, we show that coresets with $O_d(1/\varepsilon^{(d+3)/2})$ worst-case size can be computed in near-linear time.

研究の動機と目的

• 乗法的重み付きボロノイ図（MWVD）のための効率的で最適サイズの構築法が不足していること、特に2次元ではサイズが二次的になり、非凸領域や三角不等式の違反により困難であることを解決すること。
• 各セル内のすべての点に対して(1+ε)-近似最近接サイトクエリを保証する立方体ベースのε-近似MWVD（ε-AMWVD）を設計すること。
• 出力サイズを最適化し、一致する下界によって証明すること。また、Har-PeledとKumar（2004）のような従来の一般枠組みよりも優れた近似的最適な構築時間を達成すること。
• 高次元における点位置クエリ時間の低減をO(d log n)からO(d log n + d² log(1/ε))に改善し、未重み付きボロノイ図のサイズをdO(d)の要因で縮小すること。
• コア、バイセクタコアセット、適応的リファインメントといった新たな幾何的プリミティブを導入・活用し、効率的かつ証明可能に最適な近似を実現すること。

提案手法

• コアは、最大n−1個のアポロニウスのバイセクタの交差によって定義される凸領域であり、ボロノイセルの本質的構造を捉える。
• 適応的リファインメントアルゴリズムは、コアをd次元の立方体に再帰的に分割し、出力サイズを最小限に抑えつつε-近似を保証する。
• バイセクタコアセットは、各コーンごとにO(1/ε(d+3)/2)個のサイトを選択することで構築され、角度的および半径的平行移動の境界といった幾何的性質を用いて近似誤差を制御する。
• アルゴリズムはアフィン変換とコーンベースのサイト選択を用い、直径と重みに基づく4つのサイトペair（LL、LH、HL、HH）のすべてのケースにおいて(1+ε)近似要因を保証する。
• O(1/ε²C)個のサイトを1コーンあたりに持つ新規コアセット構築法を適用し、(1+εC)近似を維持する。また、重みの単調性を用いて冗長なサイトのチェックを削除する。
• 最終的なε-AMWVDは、これらのコアセットを組み合わせ、効率的な点位置検索のための圧縮されたQuadTreeベースの検索構造を用いて構築する。

実験結果

リサーチクエスチョン

• RQ1 任意のε ∈ (0,1)に対して、出力サイズがOd(n log(1/ε)/εd−1) であるε-近似MWVDを構築することは可能か？
• RQ2 出力サイズのOD(log(n)/ε(d+5)/2) 以内の構築時間を達成することは可能か？ これは理論的下界に対してΘ(1)d要因を除いて一致する。
• RQ3 コアとコアセットの使用により、この方法の副産物として未重み付きボロノイ図のサイズをdO(d)の要因で縮小できるか？
• RQ4 高次元（d = O(log n / ε²)）において、(1+ε)近似を保ちながら効率的な点位置クエリ時間（O(d log n + d² log(1/ε))）を達成することは可能か？
• RQ5 新規のバイセクタコアセットの概念を用いて、証明可能かつ(1+ε)近似保証を持つほぼ線形時間のコアセット計算を達成できるか？

主な発見

• 提案されたアルゴリズムは、出力サイズがOd(n log(1/ε)/εd−1) であるε-AMWVDを構築し、これはΘ(1)d要因を除いて下界と一致するため、サイズの最適性が確立される。
• 構築時間は出力サイズのOD(log(n)/ε(d+5)/2) 倍であり、MWVDに対して初めてこの近似的最適な構築時間を達成する方法である。
• 高次元における点位置クエリ時間はO(d log n + d² log(1/ε)) に改善され、従来の境界よりもd倍良い。
• このアプローチの副産物として、未重み付きボロノイ図のサイズがdO(d)の要因で改善され、クエリ時間も1つのd要因で短縮される。
• バイセクタコアセットは、出力サイズに対してほぼ線形時間OD(n log n / ε3(d+1)/2) でサイズOd(1/ε(d+3)/2) で計算可能である。
• 本稿では、サイズの境界におけるlog(1/ε)要因が避けられないことを証明し、MWVD近似における重要な未解決問題を解決した。