[論文レビュー] A fast quantum mechanical algorithm for estimating the median

Lov K. Grover|ArXiv.org|Jul 29, 1996

Quantum Computing Algorithms and Architecture参考文献 5被引用数 31

ひとこと要約

この論文は、精度 ε で N 個の要素の中央値を O(1/ε) ステップで推定する量子アルゴリズムを提示している。これは、古典的手法が Ω(1/ε²) 個のサンプルを必要とするのに対し、2乗の高速化を達成する。この手法は、量子アモニチュード拡大と位相推定を活用して、しきい値以下の要素を効率的に数え上げ、量子サブルーチンを用いてしきい値を繰り返し精錬することで、高速な中央値推定を実現する。

ABSTRACT

Consider the problem of estimating the median of N items to a precision epsilon, i.e., the estimate should be such that, with a high probability, the number of items, with values both smaller than and larger than this estimate, is less than N*(1+epsilon)/2. Any classical algorithm to do this will need at least O(1/epsilon^2) samples. Quantum mechanical systems can simultaneously carry out multiple computations due to their wave like properties. This paper describes an O(1/epsilon) step algorithm for the above estimation.

研究の動機と目的

• 古典的手法よりも効率的に、N 個の要素の中央値を精度 ε で推定する量子アルゴリズムを開発すること。
• 量子重ね合わせと干渉を活用して、古典的境界を超える少ないクエリ回数で必要な情報を得ること。
• 量子アルゴリズムが統計的推定問題（特に中央値推定）を、古典的手法よりも速く解けることを示すこと。
• 与えられたしきい値以下の要素の割合を推定する実用的な量子サブルーチンを提供し、これを二分探索に応用して中央値を特定すること。

提案手法

• 与えられたしきい値 µ より小さい N 個の要素の割合を推定するために、量子アモニチュード拡大を用いる。
• 高い信頼性で、µ より上と下の要素数の不均衡 ε を推定するための量子サブルーチンを採用する。
• 位相推定技術を用いて、量子アモニチュード測定から ε の符号と大きさを抽出する。
• サブルーチンが ε を推定する反復的呼び出しにより中央値に収束する二分探索フレームワークを適用する。
• アルゴリズムの初期段階で、すべての N ステートに均一な重ね合わせを作成するために、フーリエ型変換（アダマール変換）を用いる。
• しきい値以下の状態をマークするための選択的位相回転を用い、干渉に基づくアモニチュードの増幅を実現する。

実験結果

リサーチクエスチョン

• RQ1 古典的手法と比較して、量子アルゴリズムは、未ソートのデータセットの中央値推定において顕著な高速化を達成できるか？
• RQ2 与えられた精度 ε で中央値を推定するために、最小でどの程度の量子クエリが必要か？
• RQ3 量子アモニチュード拡大と位相推定をどのように組み合わせれば、しきい値以下の要素の割合を高精度に推定できるか？
• RQ4 しきい値の割合推定用の量子サブルーチンを、中央値の特定に向けた二分探索において再帰的に利用できるか？

主な発見

• アルゴリズムは O(1/ε) の量子ステップで中央値を推定し、古典的手法の Ω(1/ε²) 複雑度と比較して2乗の高速化を達成する。
• アモニチュード推定を用いて、中央値しきい値以下の要素の割合を高い信頼性で推定し、誤差は O(θε) で抑えられる。ここで θ は制御パラメータである。
• 不均衡 ε を推定するサブルーチンは、最大期待誤差 O(θε) の値を返し、アルゴリズムは真の ε がこの範囲内にある確率が非常に高いことを保証する。
• 不均衡の符号は、しきい値をわずかに変更し、推定された割合の変化を観測することで特定され、中央値の正しく特定が可能になる。
• 中央値推定は、量子サブルーチンを O(log N) 回呼び出す二分探索により実現され、全体の効率が保たれる。
• この手法は、中央値しきい値以下の要素に対応する状態のアモニチュードを増幅するために、量子干渉とユニタリー発展に依存している。