[論文レビュー] Towards Optimal Range Medians

Beat Gfeller, Peter Sanders|ArXiv.org|Jan 13, 2009

Video Surveillance and Tracking Methods参考文献 9被引用数 4

ひとこと要約

この論文は、クエリ区間で定義される部分配列の中央値を計算するためのRange Median Problem (RMP) の新しいアルゴリズムを提示する。値に基づく分割と分数連鎖を用いて、$O(n\log k + k\log n)$ 時間、$O(n\log k)$ 空間を達成し、従来の結果より対数的要因で改善され、$k = O(n)$ の場合に既知の下界と一致する。また、RAMモデルにおいて$O(\log n)$のクエリ時間で線形空間ソリューションを実現する。

ABSTRACT

We consider the following problem: given an unsorted array of $n$ elements, and a sequence of intervals in the array, compute the median in each of the subarrays defined by the intervals. We describe a simple algorithm which uses O(n) space and needs $O(n\log k + k\log n)$ time to answer the first $k$ queries. This improves previous algorithms by a logarithmic factor and matches a lower bound for $k=O(n)$. Since the algorithm decomposes the range of element values rather than the array, it has natural generalizations to higher dimensional problems -- it reduces a range median query to a logarithmic number of range counting queries.

研究の動機と目的

• 部分配列が区間によって定義される中央値を効率的に計算するためのRange Median Problem (RMP) を取り扱う。
• 従来の $O(n\log k + k\log n\log k)$ 時間を要するアルゴリズムを改善し、不要な $\log k$ 要因を排除する。
• $k = O(n)$ の場合に最適またはほぼ最適なクエリ時間を達成し、比較ベースのアルゴリズムの既知の下界と一致させる。
• 簡潔なデータ構造とビット配列上のランク操作を用いて、RAMモデルにおける線形空間への応用を可能にする。
• 高次元問題への一般化を図り、均し時間 $O(\log^2 n)$ 時間/操作で動的更新をサポートする。

提案手法

• アルゴリズムは、クイックソートに類似したアプローチを用いて、配列要素の値の範囲を再帰的に分割し、値に基づく木構造を構築する。
• 各クエリは、分数連鎖を活用することで、値のパーティション上で対数的数の範囲カウントクエリに還元される。
• RAMモデルでは、圧縮されたビット配列上のランク・セレクト操作を用いて範囲カウントを実装し、配列の縮小技術によって密度を維持する。
• データ構造は $O(n\log n)$ 時間で構築され、任意のランク中央値（中央値に限らない）に対して $O\left(\log n\right)$ クエリ時間を達成する。
• 動的バージョンでは、値で順序付けられた主構造としてBB(α)木を用い、位置で順序付けられた補助構造を備えて効率的なカウントと更新を実現する。
• 動的操作（挿入/削除）は、均し時間 $O(\log^2 n)$ でサポートされ、補助構造の回転と再構築を用い、均し時間コスト解析を適用する。

実験結果

リサーチクエスチョン

• RQ1 $k = O(n)$ の場合に、$O(n\log k + k\log n\log k)$ よりも良い時間計算量の範囲中央値クエリが達成可能か？
• RQ2 $O(n\log k + k\log n)$ 時間と $O(n\log k)$ 空間で実現可能か？ これは比較ベースのアルゴリズムの下界と一致するか？
• RQ3 RAMモデルにおいて、$O(n)$ 空間で $O(\log n)$ クエリ時間のまま維持できるか？
• RQ4 このアプローチを高次元の範囲中央値クエリに一般化可能か？
• RQ5 挿入・削除の動的操作を多項式対数時間/操作でサポート可能か？

主な発見

• 提案されたアルゴリズムは、$O(n\log k + k\log n)$ 時間と $O(n\log k)$ 空間を達成し、従来の $O(n\log k + k\log n\log k)$ 時間から対数的要因で改善された。
• このアルゴリズムは $k = O(n)$ の場合に漸近的に最適であり、比較ベースのアルゴリズムの $Ω(n\log k)$ 下界と一致する。
• $k = ω(n)$ の場合、最適時間から $\log n$ 要因以内にあり、制限付きポインタマシンモデルでは最適である。
• RAMモデルでは、$O(n)$ 空間と $O(\log n)$ クエリ時間を達成し、$O(n^\epsilon)$ クエリ時間を要する従来の線形空間ソリューションを改善する。
• この方法は自然に高次元に一般化可能であり、2次元では $O(\log n)$ クエリ時間を達成するのは依然困難であることが示されており、配列の縮小を2次元に拡張するのが困難である。
• 動的バージョンは、$O(\log^2 n)$ 均し時間で挿入・削除をサポートし、範囲中央値クエリを $O\left(\log^2 n\right)$ 最悪ケース時間で実行可能であり、画像中央値フィルタリングの既知の境界と一致する。