相対エントロピーから情報と物理を組み立てる 第6部 第4章 §6.4.4 ★★★ 1回の問い合わせで dkv を増やせる量は 2∥P^v∣ψk⟩∥ 以下で、足し上げると Δk≤4k2。成功には Δk≳N が要るので k≳N/2
1回の問い合わせで引き離せる量
6.4.1 節の一般のアルゴリズム (6.4.5) に戻る。候補のレジスタが ∣v⟩ にある部分への射影を、作業用のレジスタも含めてP^v=∣v⟩⟨v∣⊗1 と書く。∑vP^v=1 で、O^v=1−2P^v である。
定理 6.4.3(問い合わせの回数と状態の離れ方 ★★★)
どの k 回の量子アルゴリズムについても
Δk=v∑∥∣ψkv⟩−∣ψk⟩∥2≤4k2.(6.4.20)
証明
k+1 回目の問い合わせの前後を比べる。∣ψk+1v⟩=U^k+1O^v∣ψkv⟩、∣ψk+1⟩=U^k+1∣ψk⟩ で、ユニタリはノルムを変えないので
dk+1v=∥O^v∣ψkv⟩−∣ψk⟩∥≤∥O^v(∣ψkv⟩−∣ψk⟩)∥+∥(O^v−1)∣ψk⟩∥=dkv+2∥P^v∣ψk⟩∥.
1つ目の項はもう一度ユニタリ O^v を外し、2つ目は O^v−1=−2P^v を使った。d0v=0 から始めて足し上げるとdkv≤2∑j=0k−1ajv(ajv=∥P^v∣ψj⟩∥)である。Cauchy–Schwarz の不等式 (∑j=0k−1ajv)2≤k∑j=0k−1(ajv)2 を使い、v について和を取ると、∑v(ajv)2=∑v∥P^v∣ψj⟩∥2=1 なので
Δk≤4v∑(j=0∑k−1ajv)2≤4kj=0∑k−1v∑(ajv)2=4k2.□
証明の要は、問い合わせ1回で ∣ψv⟩ を ∣ψ⟩ から引き離せる量が、空の状態 ∣ψj⟩ が ∣v⟩ に置いている振幅 ajv の2倍以下だという点にある。空の状態は v を知らないので、その振幅を N 個の候補に配るしかなく、∑v(ajv)2=1 から平均の振幅は 1/N 程度にとどまる。k 回の問い合わせで dkv は高々 2k/N 程度しか増えず、その2乗の和 Δk は 4k2 を超えない。1つ1つの dkv は振幅(確率の平方根)の大きさで k に比例して増えるので、確率の大きさで見ると k2 に比例して増える。この形は前の節で見た Grover の成功確率 (2k+1)2/N と同じである。
成功には大きく引き離すことが要る
定理 6.4.4(成功の条件 ★★★)
どの v についても、k 回の後に候補のレジスタを測って v が出る確率が 1−ε 以上(0≤ε<1)なら、
Δk≥N(1−ε)−2(1−ε)N.(6.4.21)
証明
成功の確率は ∥P^v∣ψkv⟩∥2≥1−ε である。射影はノルムを増やさないので
dkv≥∥P^v(∣ψkv⟩−∣ψk⟩)∥≥∥P^v∣ψkv⟩∥−∥P^v∣ψk⟩∥≥1−ε−bv,bv=∥P^v∣ψk⟩∥.
c=1−ε とおく。c−bv≥0 なら2乗して (dkv)2≥(c−bv)2≥c2−2cbv である。c−bv<0 なら c2−2cbv<c2−2c2≤0≤(dkv)2 である。どちらの場合も (dkv)2≥c2−2cbv が成り立つ。v について和を取り、∑vbv≤N∑v(bv)2=N(Cauchy–Schwarz)を使えばΔk≥Nc2−2cN を得る。
系 6.4.5(探索の下界(BBBV) ★★★)
すべての v について成功確率 1−ε 以上で当たりを見つける量子アルゴリズムは、
k≥21(1−ε)N−2(1−ε)N=21−εN(1+O(N−1/2))(6.4.22)
回の問い合わせを必要とする。
定理 6.4.3 と定理 6.4.4 を合わせればよい。Grover のアルゴリズムは (π/4)N≈0.785N 回で成功確率 1−O(1/N) を達成し、下界は 0.5N なので、両者は N の比例係数の違い(π/2 倍)を除いて一致する。C. Zalka は、より細かい議論で、成功確率が 1 に近い場合の最小の回数が N→∞ で (π/4)N に一致すること、つまり Grover のアルゴリズムが比例係数まで含めて最適であることを示した1C. Zalka (1999), Phys. Rev. A 60, 2746。。
Grover のアルゴリズム自身で Δk を計算しておく。空のオラクルでは G^=R^s で、R^s∣s⟩=∣s⟩ なので ∣ψk⟩=∣s⟩ のままである。式 (6.4.13) から ⟨s∣ψkv⟩=sinθsin((2k+1)θ)+cosθcos((2k+1)θ)=cos(2kθ) で、
ΔkGrover=N(2−2cos2kθ)=4Nsin2(kθ)≤4Nk2sin2θ=4k2(6.4.23)
となる(∣sinkθ∣≤k∣sinθ∣ を使った)。kθ が小さいうちは sin(kθ)≈kθ≈k/N で、Grover のアルゴリズムは不等式 (6.4.20) をほぼ等号で満たしている。1回の問い合わせで引き離せる量の上限を、毎回ほぼ使い切っているということである。
下界が言っていること
下界 (6.4.22) は、オラクルの中身を問わず、問い合わせの回数だけで述べた主張である。したがって、問題に構造があってオラクルの中身を覗ける場合には当てはまらない。Shor のアルゴリズム(第6部第3章)が指数的に速いのは、位数発見という問題に周期という構造があり、量子 Fourier 変換がそれを読むからである。一方、充足可能性問題(SAT)のような NP 完全問題に対して、候補の割り当てを1つずつ試すことしかできないなら、2n 個の候補に対して 2n/2 回ほどの問い合わせが要り、指数的な手間は残る。BBBV の論文は、この意味で「オラクルを相手にする限り、NP の問題が量子の多項式時間で解けるとは限らない」ことを示した。これは NP と BQP(量子の多項式時間で解ける問題のクラス)の関係を決めた結果ではなく、構造を使わない方法では足りないことを示した結果である。
例題 6.4.4 100万個の候補での上界と下界 ★★
オリジナルN=106 とする。
- Grover のアルゴリズムの最適な回数 k∗ と、そのときの失敗の確率を求めよ。
- 誤りなく(ε=0)当てるのに要る問い合わせの回数の下界 (6.4.22) を求めよ。
- k=k∗ での ΔkGrover と上限 4k2 を比べよ。
解答
(1) θ=arcsin(10−3)=1.0000002×10−3、π/(4θ)=785.40 なので k∗=785。失敗の確率は 1−sin2(1571θ)=4.2×10−8 で、下限 sin2θ=10−6 よりずっと小さい(角 (2k∗+1)θ がたまたま π/2 に近い)。
(2) 21106−2×103=21998000=499.5 なので、少なくとも 500 回。古典の手順では誤りなく当てるのに N−1=999999 回かかる。
(3) Δ785Grover=4×106sin2(785θ)=2.00×106、上限は 4×7852=2.46×106 である。kθ=0.785≈π/4 まで来ると sin(kθ) は kθ より目に見えて小さくなるので、等号から離れる。成功のためには定理 6.4.4 から Δk≥106−2×103=9.98×105 が要り、Grover はその約2倍まで引き離している。
1 C. Zalka (1999), Phys. Rev. A 60, 2746。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com