機会があり、チュートリアルをやってみました。<a href="https://qiskit.org/textbook/ch-algorithms/grover.html" rel="noopener noreferrer" target="_blank">https:</a><a href="https://qiskit.org/textbook/ch-algorithms/grover.html" rel="noopener noreferrer" target="_blank" class="hljs-regexp">//</a><a href="https://qiskit.org/textbook/ch-algorithms/grover.html" rel="noopener noreferrer" target="_blank">qiskit.org</a><a href="https://qiskit.org/textbook/ch-algorithms/grover.html" rel="noopener noreferrer" target="_blank" class="hljs-regexp">/textbook/</a><a href="https://qiskit.org/textbook/ch-algorithms/grover.html" rel="noopener noreferrer" target="_blank">ch-algorithms/grover.html</a> Groverのアルゴリズムは汎用量子アルゴリズムでハイブリッドではないタイプのものです。欲しいデータを探すために整理されていないデータでの検索での利用が想定されています。 １、初期の準備２、マーキング３、振幅増幅 の３ステップでできます。１と３は決まった形の回路になっているので、私たちが設計するのは２のマーキング回路で、オラクルとも呼ばれます。 <h2>インストール</h2>Qiskitをpythonで利用するには、pip経由でインストールします。量子回路の表記を豊かにするためにおまけも一緒にインストールします。 <pre class="ql-syntax" spellcheck="false">pip install qiskit pylatexenc
</pre> <h2>ツール読み込み</h2>ツールを読み込みます。qiskitとnumpyがメインです。 <pre class="ql-syntax" spellcheck="false">import numpy as np
from qiskit import *
%matplotlib inline
</pre> <h2>０、回路の初期化</h2>始める前に使う量子ビットを準備します。今回は2量子ビット利用します。 <pre class="ql-syntax" spellcheck="false">#2量子ビットの回路を作成する
grover = QuantumCircuit(2)
</pre> <h2>１、初期の準備</h2>Groverでは初期状態は重ね合わせで準備します。これは理論的な理由でそうなっています。別に他の状態を準備してもいいと思いますが、これが一番簡単です。回路は全ての量子ビットにアダマールゲートと呼ばれる重ね合わせの演算を適用します。 <pre class="ql-syntax" spellcheck="false">#初期状態を作る
grover.h([0,1])
</pre> hゲートを適用しました。一旦回路を確認してみます。下記で描画できます。 <pre class="ql-syntax" spellcheck="false">#量子回路を描画
grover.draw(&#39;mpl&#39;)
</pre> <img src="https://pro.blueqat.com/public/uploads/us-east-2:032aab15-1641-4f97-b43d-6d490c139451/2020/11/08/download.png"/> 初期状態は00,01,10,11の重ね合わせからスタートします。 <h2>２、マーキング</h2>マーキング回路は今回は11のデータを探すものとして作ってみます。状態ベクトルの具体的な説明は今回は省きますが、11のデータは状態ベクトルでは、00,01,10,11の最後の要素になります。最後の要素を見つけるには、CZゲートを使うとマーキングできます。CZゲートは[[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,-1]]のように、状態ベクトルの最後の要素にマイナスをつける行列になっています。 <pre class="ql-syntax" spellcheck="false">#マーキング回路を作成
grover.cz(0,1)
</pre> 早速作った回路を確認してみます。 <pre class="ql-syntax" spellcheck="false">#量子回路を描画
grover.draw(&#39;mpl&#39;)
</pre> <img src="https://pro.blueqat.com/public/uploads/us-east-2:032aab15-1641-4f97-b43d-6d490c139451/2020/11/08/download-1.png"/> <h2>３、振幅増幅</h2>こちらは定型作業です。マーキングを増幅させます。対角成分がマイナスの行列が作れればOKです。今回は下記を使いました。こちらのHH-ZZ-CZ-HHの仕組みについては、記事の最後におまけで計算を書いてみます。 <pre class="ql-syntax" spellcheck="false">#振幅増幅回路を作成
grover.h([0,1])
grover.z([0,1])
grover.cz(0,1)
grover.h([0,1])
grover.draw(&#39;mpl&#39;)
</pre> <img src="https://pro.blueqat.com/public/uploads/us-east-2:032aab15-1641-4f97-b43d-6d490c139451/2020/11/08/download-2.png"/> これで完成です！あとは測定と言って、計算結果を取り出す回路を最後につけます！ <h2>４、測定</h2>測定回路をつけます。 <pre class="ql-syntax" spellcheck="false"># 量子回路に測定を入れる
meas = QuantumCircuit(2, 2)
meas.measure(range(2), range(2))

# 測定込みの回路にする
qc = grover + meas

#量子回路を描画
qc.draw(&#39;mpl&#39;)
</pre> <img src="https://pro.blueqat.com/public/uploads/us-east-2:032aab15-1641-4f97-b43d-6d490c139451/2020/11/08/download-3.png"/> <h2>５、計算</h2>早速計算してみましょう。 今回はシミュレータを準備します。実機でgroverを計算した記事もありますので、別途参照してください。shotsと言って1024回同じ計算をします。今回はほぼ無意味ですが、、、そしてジョブを実行して、計算結果を取り出し、カウントします。 <pre class="ql-syntax" spellcheck="false">#シミュレータを準備
backend_sim = Aer.get_backend(&#39;qasm_simulator&#39;)

#計算を1024回行う
job_sim = execute(qc, backend_sim, shots=1024)

#計算結果を取り出す
result_sim = job_sim.result()

#計算結果をカウントする
counts = result_sim.get_counts(qc)
print(counts)
</pre> 予定通り11が1024回出ました。 <pre class="ql-syntax" spellcheck="false">{&#39;11&#39;: 1024}
</pre> ヒストグラムの作り方は、 <pre class="ql-syntax" spellcheck="false">#ヒストグラムへ落とし込む
from qiskit.visualization import plot_histogram
plot_histogram(counts)
</pre> <img src="https://pro.blueqat.com/public/uploads/us-east-2:032aab15-1641-4f97-b43d-6d490c139451/2020/11/08/download-4.png"/> こちらの方が見やすいです。オラクルは難しいものを作ればいろんな問題が解けますのでぜひ色々やってみてください。 <h2>おまけ</h2>途中の振幅増幅回路のHH-ZZ-CZ-HHがなぜこれなのか、こちらで確認します。マーキングは特定の状態ベクトルだけの符号を反転させます。その後に対角成分が全てマイナスの行列をかけるとマーキングした状態ベクトルの要素だけが増幅します。 <pre class="ql-syntax" spellcheck="false">import numpy as np

#アダマールゲートを作る
h = np.array([[1,1],[1,-1]])/np.sqrt(2)
h

#array([[ 0.70710678, 0.70710678],
# [ 0.70710678, -0.70710678]])

#アダマールゲートのテンソル積をとります
hh = np.kron(h,h)
hh

#array([[ 0.5, 0.5, 0.5, 0.5],
# [ 0.5, -0.5, 0.5, -0.5],
# [ 0.5, 0.5, -0.5, -0.5],
# [ 0.5, -0.5, -0.5, 0.5]])

#Zゲートを作る
z = np.array([[1,0],[0,-1]])
z

#array([[ 1, 0],
# [ 0, -1]])

#Zゲートのテンソル積
zz = np.kron(z,z)
zz

#array([[ 1, 0, 0, 0],
# [ 0, -1, 0, 0],
# [ 0, 0, -1, 0],
# [ 0, 0, 0, 1]])

#CZゲート
cz = np.array([[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,-1]])
cz

#array([[ 1, 0, 0, 0],
# [ 0, 1, 0, 0],
# [ 0, 0, 1, 0],
# [ 0, 0, 0, -1]])

#振幅増幅回路の行列を確認
hh@zz@cz@hh

#array([[-0.5, 0.5, 0.5, 0.5],
# [ 0.5, -0.5, 0.5, 0.5],
# [ 0.5, 0.5, -0.5, 0.5],
# [ 0.5, 0.5, 0.5, -0.5]])
</pre> これで、対角成分が全てマイナスになっているのが確認できました。以上です！では、量子と共に！

IBMQのQiskitでgroverのアルゴリズム

Yuichiro Minato