要旨：タイトルに対する回答：大きくYESとは言えないが、ある程度可能、である。教育的観点から、量子コンピューティングを身近な問題に適用してみたいと考えていた。その一例として、1777年に発表された「ビュフォンの針（πの近似計算）」に着目した。たくさん平行線が引かれた紙の上に、多数の針をばらまき、それらが平行線と交わる確率pを求める。すると（平行線の間隔と針の長さが共に１の場合）、π=2/Pとなる！とても不思議に思えるが、それは高校数学で証明できる。そのシミュレーションは、通常モンテカルロ法による。だが、ここでは量子的方法で実行を試みた。結論としては、完全な量子計算にすることは困難だが、量子コンピューティングの重要な側面に触れるという教育的意義はあった。
🟢「ビュフォンの針」の古典版と量子版の概要
円周率πの近似値を「ビュフォンの針」によって求める。従来のモンテカルロ法と、量子振幅推定（QAE: Quantum Amplitude Estimation）を用いる方法の両方で試した。どちらの方法でも、多数の針に対して、針が平行線と交差するかどうかを、針の中心位置と角度から判定することが必要である。
モンテカルロ法では、各針について交差の有無を調べ、交差した針の本数を数えて合計し、それを全体の本数で割ることで、交差確率の近似値が得られる。
一方QAEでは、状態準備回路によって、針の中心位置と角度に対応する量子状態の全てを、量子重ね合わせ状態として準備する。その後、オラクルと呼ばれる量子回路によって交差状態のみに位相反転を与え、量子干渉を利用して交差状態に対応する振幅を増幅させる。さらに、その振幅に対応する回転角（固有位相）を推定することで、交差確率pを直接推定することができる。
どちらの方法でも交差確率pが求まれば、平行線の間隔と針の長さがともに1の場合、（ビュフォンの針の肝である）p=2/πが成立するので、πを推定することができる。
🟢古典版「ビュフォンの針」の実行
まず、モンテカルロ法によるシミュレーションのアニメをご覧頂きたい。
乱数系列に強く依存するので、何とも言えないが、図１の通り、8,900本の針を投入した場合、πの近似値として５桁正解の3.1415まで得られたことに驚いた。
🟢量子版「ビュフォンの針」の実行
上で述べた本来のオラクルは、sinθの計算と比較演算による交差判定を行い、交差状態にある針を意味する量子状態の位相を反転させる。だが、それを完全な量子回路として実現するには多くの困難を伴う。そのため、ここでは、交差判定は、あらかじめ古典計算で済ませることにした。それに基づく位相反転機能だけをオラクル量子回路に持たせた。
そのようなオラクルを図２に示した。針の中心位置yと角度θをそれぞれ、分解性能を３量子ビットで表現しているので、それらで決まる8×8=64個の針の量子状態を全て重ねた状態を一度に持つことができる。このオラクルは、そのうちの交差状態に対して最終段の標識量子ビットq6を反転させている。
例えば、最初のMCX（多制御付きNOT）では、針の量子状態|000000〉が交差と判定された後、標識ビットq6が反転するようになっているので、|000000〉⊗|1〉という状態となる。その後、標識ビットが|1〉の場合は、|000000〉の位相が反転される。
次に、このオラクルを、量子振幅推定器QAEに渡す。QAEは、重要な量子回路ライブラリとして使えるようになっている。QAEは、量子干渉を利用して交差状態に対応する振幅をGrover演算で増幅させ、その振幅に対応する回転角（固有位相）を推定することで、交差確率pを直接推定することができる。
ともかく、このようにして、量子的ビュフォンをシミュレーションで実行させた。その結果を以下に示す。yとθの解像度と、QAEで使う評価用量子ビット数に強く依存するが、今回はπ=3.12というある程度の精度を得た。上記の6量子ビットと標識用１量子ビットに加えて、AQEで必要な評価用９量子ビットの合計16量子ビットを使っている。
==簡易OracleとQAEによる、量子的「ビュフォンの針」の実行結果==
使用した全量子ビット数 : 16
QAE推定確率 P (交差確率) : 0.639260
計算された円周率 π (推定値) : 3.128618
実際の円周率 (理論値) : 3.141593
誤差 : 0.012974
================================================
🟢補 足
精度の面では、モンテカルロ法では試行回数Nを増やすと推定誤差はO(1/√N)で減少する。したがって精度を10倍向上させるためには、理論上100倍の試行が必要となる。一方、QAEではオラクルの呼び出し回数Nに対して推定誤差がO(1/N)で減少するため、精度を10倍向上させるためには理論上10倍のオラクル呼び出しで済む。この結果、QAEは古典モンテカルロ法に対して二次加速（quadratic speedup）を実現できることが理論的に示されている。
ただし、ここで用いた回数N意味は両者で異なることに注意が必要である。モンテカルロ法におけるNは試行した針の本数を表すのに対し、QAEにおける回数Nはオラクル（交差判定回路）の呼び出し回数を表す。したがって、この二次加速は実行時間そのものではなく、交差判定へのアクセス回数（クエリ複雑度）の観点での理論的な優位性を意味している。
🟢感 想
「ビュフォンの針」による円周率近似において、簡易的、部分的にではあるが、オラクルという量子回路を設定し、それをQAE（量子振幅推定）に繋げて結果を出すことができた。古典的にはモンテカルロ法で求める交差確率を、個々の試行での加算を行わずに、振幅増幅によって直接読み出すという、量子コンピューティングの重要な側面に触れることはできた。ただし、本来の（本物の）オラクル回路にはなっていない点は明確にしておきたい。今後、教育的により適切な量子計算の例題を探求していきたい。
0 件のコメント:
コメントを投稿