{
"cells": [
{
"cell_type": "markdown",
"id": "17463a96",
"metadata": {},
"source": [
"---\n",
"title: \"Deutsch-Jozsaアルゴリズム\"\n",
"description: \"量子情報と量子計算に関する無料の IBM 講座\"\n",
"---\n",
"\n",
"\n",
"\n",
"# Deutsch-Jozsaアルゴリズム\n",
"\n",
"Deutschのアルゴリズムは、クエリ問題では古典的なアルゴリズムよりも優れているが、その優位性は1クエリ対2クエリという極めてささやかなものである。\n",
"Deutsch-Jozsaアルゴリズムは、この利点を拡張するもので、実際、2つの異なるクエリー問題を解くのに使うことができる。\n",
"\n",
"ドイチュ・ヨッサ・アルゴリズムの量子回路の説明である。\n",
"図には示されていないが、解決しようとする特定の問題によっては、さらに古典的な後処理ステップが必要になることもある。\n",
"\n",
"\n",
"\n",
"もちろん、このアルゴリズムがどのような問題を解決するのかについては、実際には議論していない。\n",
"\n",
"\n",
"\n",
"## ドイッチュ＝ヨージャ問題\n",
"\n",
"まず、Deutsch-Jozsaアルゴリズムが元々解こうとしていたクエリー問題（ *Deutsch-Jozsa問題として*知られている）から始めよう。\n",
"\n",
"この問題の入力関数は、任意の正の整数に対して $f:\\Sigma^n \\rightarrow \\Sigma$ の形をとる。 $n.$ Deutschの問題と同様に、 $f$ が定数なら $0$ を、 $f$ が釣り合いなら $1$ を出力することである。これは、関数が値 $0$ をとる入力文字列の数と、関数が値 $1$ をとる入力文字列の数が等しいことを意味する。\n",
"\n",
"$n$ が $1,$ より大きい場合、定数でも釣り合いでもない $f:\\Sigma^n \\rightarrow \\Sigma$ という形の関数が存在することに注意。\n",
"例えば、関数 $f:\\Sigma^2\\rightarrow\\Sigma$ は次のように定義される\n",
"\n",
"$$\n",
"\\begin{aligned}\n",
"f(00) & = 0 \\\\\n",
"f(01) & = 0 \\\\\n",
"f(10) & = 0 \\\\\n",
"f(11) & = 1\n",
"\\end{aligned}\n",
"$$\n",
"\n",
"この2つのカテゴリーには当てはまらない。\n",
"Deutsch-Jozsa問題では、このような関数は気にしない。\n",
"つまり、この問題では、 $f$ が一定か均衡のどちらかであるという*約束が*ある。\n",
"\n",
"\n",
" 入力：関数 $f:\\{0,1\\}^n\\rightarrow\\{0,1\\}$ \\ 約束： $f$ が一定か均衡のいずれかである。 \\ 出力： $f$ が定数なら $0$, $f$ が釣り合いなら $1$\n",
"\n",
"\n",
"Deutsch-Jozsaアルゴリズムは、単一のクエリを用いて、以下の意味でこの問題を解決する。\n",
"$n$ の測定結果のすべてが $0,$ である場合、関数 $f$ は定数である。\n",
"それ以外の場合、すなわち測定結果のうち少なくとも1つが $1,$ である場合、関数 $f$ は平衡である。\n",
"別の言い方をすれば、上述の回路の後に、測定結果の論理和（OR）を計算して、ドイッチュ・ジョザ問題の出力ビットを生成する、古典的な後処理ステップが続くということである。\n",
"\n",
"\n",
"\n",
"### アルゴリズム分析\n",
"\n",
"Deutsch-Jozsa問題に対するDeutsch-Jozsaアルゴリズムの性能を分析するためには、ハダマードゲートの単一層の作用について考えることから始めると役に立つ。\n",
"ハダマード演算は、通常の方法で行列として表すことができる、\n",
"\n",
"$$\n",
"H = \\begin{pmatrix}\n",
"\\frac{1}{\\sqrt{2}} & \\frac{1}{\\sqrt{2}} \\\\[2mm]\n",
"\\frac{1}{\\sqrt{2}} & -\\frac{1}{\\sqrt{2}}\n",
"\\end{pmatrix},\n",
"$$\n",
"\n",
"しかし、この操作を標準的な基底状態に対する作用という観点から表現することもできる：\n",
"\n",
"$$\n",
"\\begin{aligned}\n",
"H \\vert 0\\rangle & = \\frac{1}{\\sqrt{2}} \\vert 0 \\rangle + \\frac{1}{\\sqrt{2}} \\vert 1 \\rangle\\\\[3mm]\n",
"H \\vert 1\\rangle & = \\frac{1}{\\sqrt{2}} \\vert 0 \\rangle - \\frac{1}{\\sqrt{2}} \\vert 1 \\rangle.\n",
"\\end{aligned}\n",
"$$\n",
"\n",
"この2つの方程式は1つの式にまとめることができる、\n",
"\n",
"$$\n",
"H \\vert a \\rangle = \\frac{1}{\\sqrt{2}} \\vert 0 \\rangle + \\frac{1}{\\sqrt{2}} (-1)^a \\vert 1 \\rangle\n",
"= \\frac{1}{\\sqrt{2}} \\sum_{b\\in\\{0,1\\}} (-1)^{ab} \\vert b\\rangle,\n",
"$$\n",
"\n",
"のどちらを選んでも同じである。 $a\\in\\Sigma.$\n",
"\n",
"ここで、1つの量子ビットの代わりに $n$、それぞれの量子ビットに対してハダマード演算が行われたとする。\n",
"$n$ 量子ビットの複合演算は、テンソル積 $H\\otimes \\cdots \\otimes H$ （ $n$ 回）によって記述される。ここでは簡潔かつ明瞭にするために $H^{\\otimes n}$ と表記する。\n",
"上記の公式を使い、展開して単純化すると、 $n$ の標準基底状態に対するこの複合操作の作用を次のように表すことができる：\n",
"\n",
"$$\n",
"\\begin{aligned}\n",
" & H^{\\otimes n} \\vert x_{n-1} \\cdots x_1 x_0 \\rangle \\\\\n",
" & \\qquad = \\bigl(H \\vert x_{n-1} \\rangle \\bigr) \\otimes \\cdots \\otimes \\bigl(H \\vert x_{0} \\rangle \\bigr) \\\\\n",
" & \\qquad = \\Biggl( \\frac{1}{\\sqrt{2}} \\sum_{y_{n-1}\\in\\Sigma} (-1)^{x_{n-1} y_{n-1}} \\vert y_{n-1} \\rangle \\Biggr)\n",
" \\otimes \\cdots \\otimes\n",
" \\Biggl( \\frac{1}{\\sqrt{2}} \\sum_{y_{0}\\in\\Sigma} (-1)^{x_{0} y_{0}} \\vert y_{0} \\rangle \\Biggr) \\\\\n",
" & \\qquad = \\frac{1}{\\sqrt{2^n}} \\sum_{y_{n-1}\\cdots y_0 \\in \\Sigma^n}\n",
" (-1)^{x_{n-1}y_{n-1} + \\cdots + x_0 y_0} \\vert y_{n-1} \\cdots y_0 \\rangle.\n",
"\\end{aligned}\n",
"$$\n",
"\n",
"ちなみにここでは、 $n$ の長さのバイナリ文字列を、Qiskitのインデックス規則に従って、 $x_{n-1}\\cdots x_0$ と $y_{n-1}\\cdots y_0,$ のように書いています。\n",
"\n",
"この式は、上記の量子回路を解析するための便利なツールを提供してくれる。\n",
"ハダマードゲートの第1層が実行された後、 $n+1$ の量子ビット（左端/下端の量子ビットを含み、他とは別に扱われる）の状態は次のようになる\n",
"\n",
"$$\n",
"\\bigl( H \\vert 1 \\rangle \\bigr) \\bigl( H^{\\otimes n} \\vert 0 \\cdots 0 \\rangle \\bigr)\n",
"= \\vert - \\rangle \\otimes \\frac{1}{\\sqrt{2^n}} \\sum_{x_{n-1}\\cdots x_0 \\in \\Sigma^n} \\vert x_{n-1} \\cdots x_0 \\rangle.\n",
"$$\n",
"\n",
"$U_f$、この状態は次のように変換される\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{\\sqrt{2^n}}\n",
"\\sum_{x_{n-1}\\cdots x_0 \\in \\Sigma^n} (-1)^{f(x_{n-1}\\cdots x_0)} \\vert x_{n-1} \\cdots x_0 \\rangle\n",
"$$\n",
"\n",
"ドイチュのアルゴリズムの分析で見られたのとまったく同じ位相キックバック現象によって。\n",
"\n",
"次に、2層目のハダマードゲートが実行され、（上式によって）この状態は次のように変換される\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{2^n}\n",
"\\sum_{x_{n-1}\\cdots x_0 \\in \\Sigma^n}\n",
"\\sum_{y_{n-1}\\cdots y_0 \\in \\Sigma^n}\n",
"(-1)^{f(x_{n-1}\\cdots x_0) + x_{n-1}y_{n-1} + \\cdots + x_0 y_0}\n",
"\\vert y_{n-1} \\cdots y_0 \\rangle.\n",
"$$\n",
"\n",
"この式はやや複雑に見えるので、関数についてもっと詳しく知らなければ、異なる測定結果を得る確率についてあまり多くのことを結論づけることはできない。 $f.$\n",
"\n",
"幸いなことに、私たちが知る必要があるのは、すべての測定結果が $0$ である確率だけである。それは、 $f$ が一定であるとアルゴリズムが判断する確率だからである。\n",
"この確率には簡単な公式がある。\n",
"\n",
"$$\n",
"\\Biggl\\vert\n",
"\\frac{1}{2^n}\n",
"\\sum_{x_{n-1}\\cdots x_0 \\in \\Sigma^n}\n",
"(-1)^{f(x_{n-1}\\cdots x_0)}\n",
"\\Biggr\\vert^2\n",
"= \\begin{cases}\n",
"1 & \\text{if $f$ is constant}\\\\[1mm]\n",
"0 & \\text{if $f$ is balanced}\n",
"\\end{cases}\n",
"$$\n",
"\n",
"なお、これらの値は、\n",
"$\\vert 0^{\\otimes n} \\rangle$ という状態が観測される確率に対応するものであり、\n",
"ドイッチュ・ジョザ問題の最終的な古典的な出力ビットに直接対応するものではないことに注意してください。 このアルゴリズムは、すべての測定結果が\n",
"$0$（$f$が定数であることを示す）である場合に$0$を出力し、それ以外の場合は\n",
"$1$を出力する\n",
"（$f$が平衡状態にあることを示す）。\n",
"\n",
"より詳細には、 $f$ が定数である場合、すべての文字列について $f(x_{n-1}\\cdots x_0) = 0$ のどちらかとなる。 $x_{n-1}\\cdots x_0,$ $f(x_{n-1}\\cdots x_0) = 1$ となり、その場合の和の値は $2^n,$ となる。 $x_{n-1}\\cdots x_0,$ である。 $-2^n.$ $2^n$ で割り、絶対値の2乗をとると、次のようになる。 $1.$\n",
"\n",
"一方、 $f$ が釣り合っている場合、 $f$ は文字列の半分に値 $0$ をとり、 $x_{n-1}\\cdots x_0$ と残りの半分に値 $1$ をとる。したがって、合計の $+1$ の項と $-1$ の項は相殺され、次のような値が残る。 $0.$\n",
"\n",
"約束が果たされる限り、アルゴリズムは正しく動作すると結論づける。\n",
"\n",
"\n",
"\n",
"### 古典的難度\n",
"\n",
"Deutsch-Jozsaアルゴリズムは毎回機能し、約束が守られれば常に正しい答えを出し、クエリーは1回で済む。\n",
"Deutsch-Jozsa問題に対する古典的なクエリーアルゴリズムと比較してどうか？\n",
"\n",
"まず、Deutsch-Jozsa問題を正しく解く*決定論的*古典アルゴリズムは、指数関数的に多くの問い合わせを行わなければならない： $2^{n-1} + 1$ 最悪の場合、クエリーは指数関数的に多くなる。\n",
"その理由は、もし決定論的アルゴリズムが $f$、 $2^{n-1}$ 以下の異なる文字列を問い合わせ、毎回同じ関数値を得るなら、両方の答えが可能であるということだ。\n",
"関数は一定かもしれないし、バランスが取れているかもしれないが、運が悪いことにクエリーはすべて同じ関数値を返す。\n",
"\n",
"しかし、決定論的アルゴリズムにはランダム性や不確実性がないため、特定の関数でシステマティックに失敗することになる。\n",
"従って、この点では古典的アルゴリズムよりも量子の方が大きな利点がある。\n",
"\n",
"しかし、 *確率的な*古典的アルゴリズムは、わずか数回のクエリーで非常に高い確率でDeutsch-Jozsa問題を解くことができるというキャッチがある。\n",
"特に、長さ $n$ の数種類の文字列をランダムに選び、それらの文字列に対して $f$ 問い合わせた場合、 $f$ のバランスが取れたときに、すべての文字列に対して同じ関数値が得られる可能性は低い。\n",
"\n",
"具体的には、 $k$ 入力文字列 $x^1,\\ldots,x^k \\in \\Sigma^n$ を一様にランダムに選び、 $f(x^1),\\ldots,f(x^k),$ を評価し、関数値がすべて同じであれば $0$ に、そうでなければ $1$ に答えるとすると、 $f$ が一定の場合は常に正解となり、 $f$ がちょうど釣り合う場合は確率で不正解となる。 $2^{-k + 1}.$ $k = 11,$ を例にとると、このアルゴリズムは $99.9$ %以上の確率で正解となる。\n",
"\n",
"このような理由から、古典的アルゴリズムに対する量子の利点はまだ控えめであるが、それでもドイチュのアルゴリズムに対する改善を示す定量的な利点である。\n",
"\n"
]
},
{
"cell_type": "markdown",
"id": "83b0e8b5",
"metadata": {},
"source": [
"\n",
"\n",
"## Deutsch-Jozsa with Qiskit\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "c7839d8f",
"metadata": {},
"outputs": [],
"source": [
"from qiskit import QuantumCircuit\n",
"from qiskit_aer import AerSimulator\n",
"import numpy as np"
]
},
{
"cell_type": "markdown",
"id": "59738a5c",
"metadata": {},
"source": [
"QiskitでDeutsch-Jozsaアルゴリズムを実装するために、まずDeutsch-Jozsa問題の約束を満たす関数をランダムに選び、クエリーゲートを実装した量子回路を生成する関数 `dj_query` 。\n",
"50％の確率で機能は一定であり、50％の変化で機能は均衡する。\n",
"この2つの可能性のそれぞれについて、そのタイプの関数から一様に関数が選択される。\n",
"引数は関数の入力ビット数。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "b39e7af8",
"metadata": {},
"outputs": [],
"source": [
"def dj_query(num_qubits):\n",
" # Create a circuit implementing for a query gate for a random function\n",
" # satisfying the promise for the Deutsch-Jozsa problem.\n",
"\n",
" qc = QuantumCircuit(num_qubits + 1)\n",
"\n",
" if np.random.randint(0, 2):\n",
" # Flip output qubit with 50% chance\n",
" qc.x(num_qubits)\n",
" if np.random.randint(0, 2):\n",
" # return constant circuit with 50% chance\n",
" return qc\n",
"\n",
" # Choose half the possible input strings\n",
" on_states = np.random.choice(\n",
" range(2**num_qubits), # numbers to sample from\n",
" 2**num_qubits // 2, # number of samples\n",
" replace=False, # makes sure states are only sampled once\n",
" )\n",
"\n",
" def add_cx(qc, bit_string):\n",
" for qubit, bit in enumerate(reversed(bit_string)):\n",
" if bit == \"1\":\n",
" qc.x(qubit)\n",
" return qc\n",
"\n",
" for state in on_states:\n",
" qc.barrier() # Barriers are added to help visualize how the functions are created.\n",
" qc = add_cx(qc, f\"{state:0b}\")\n",
" qc.mcx(list(range(num_qubits)), num_qubits)\n",
" qc = add_cx(qc, f\"{state:0b}\")\n",
"\n",
" qc.barrier()\n",
"\n",
" return qc"
]
},
{
"cell_type": "markdown",
"id": "60b36f11",
"metadata": {},
"source": [
"クエリーゲートの量子回路実装は、通常通り `draw` 。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "b4a6df3d",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"display(dj_query(3).draw(output=\"mpl\"))"
]
},
{
"cell_type": "markdown",
"id": "6ed6dd43",
"metadata": {},
"source": [
"次に、クエリーゲートの量子回路実装を引数として、Deutsch-Jozsa回路を作成する関数を定義する。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "79ab8817",
"metadata": {},
"outputs": [],
"source": [
"def compile_circuit(function: QuantumCircuit):\n",
" # Compiles a circuit for use in the Deutsch-Jozsa algorithm.\n",
"\n",
" n = function.num_qubits - 1\n",
" qc = QuantumCircuit(n + 1, n)\n",
" qc.x(n)\n",
" qc.h(range(n + 1))\n",
" qc.compose(function, inplace=True)\n",
" qc.h(range(n))\n",
" qc.measure(range(n), range(n))\n",
" return qc"
]
},
{
"cell_type": "markdown",
"id": "71c07df9",
"metadata": {},
"source": [
"最後に、Deutsch-Jozsa回路を一度実行する関数が定義されている。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "8a6cf538",
"metadata": {},
"outputs": [],
"source": [
"def dj_algorithm(function: QuantumCircuit):\n",
" # Determine if a function is constant or balanced.\n",
"\n",
" qc = compile_circuit(function)\n",
"\n",
" result = AerSimulator().run(qc, shots=1, memory=True).result()\n",
" measurements = result.get_memory()\n",
" if \"1\" in measurements[0]:\n",
" return \"balanced\"\n",
" return \"constant\""
]
},
{
"cell_type": "markdown",
"id": "49113ebf",
"metadata": {},
"source": [
"ランダムに関数を選び、その関数に対するクエリゲートの量子回路実装を表示し、その関数に対してDeutsch-Jozsaアルゴリズムを実行することで、我々の実装をテストすることができる。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "aca4745e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"'balanced'"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"f = dj_query(3)\n",
"display(f.draw(\"mpl\"))\n",
"display(dj_algorithm(f))"
]
},
{
"cell_type": "markdown",
"id": "89e58116",
"metadata": {},
"source": [
"\n",
"\n",
"## バーンスタイン＝ヴァジラニ問題\n",
"\n",
"次に、 *バーンスタイン＝ヴァジラニ問題として*知られる問題について説明しよう。\n",
"これは*フーリエ・サンプリング問題とも*呼ばれるが、この問題にはもっと一般的な定式化もあり、その名前でも呼ばれている。\n",
"\n",
"まず、いくつかの表記法を紹介しよう。\n",
"長さ $n,$ の任意の2つのバイナリ文字列 $x = x_{n-1} \\cdots x_0$ と $y = y_{n-1}\\cdots y_0$ に対して、次のように定義する\n",
"\n",
"$$\n",
"x \\cdot y = x_{n-1} y_{n-1} \\oplus \\cdots \\oplus x_0 y_0.\n",
"$$\n",
"\n",
"この演算を *2進ドット積と*呼ぶことにする。\n",
"別の言い方をすれば、こうなる。\n",
"\n",
"$$\n",
"x \\cdot y =\n",
"\\begin{cases}\n",
"1 & x_{{n-1}} y_{n-1} + \\cdots + x_0 y_0 \\text{ is odd}\\\\[0.5mm]\n",
"0 & x_{{n-1}} y_{n-1} + \\cdots + x_0 y_0 \\text{ is even}\n",
"\\end{cases}\n",
"$$\n",
"\n",
"これは対称的な操作であり、 $x$ と $y,$ を入れ替えても結果は変わらない。\n",
"バイナリのドット積 $x \\cdot y$ を、文字列 $y$ が $1,$ を持つ位置における $x$ のビットのパリティ、あるいは同等に、文字列 $x$ が を持つ位置における $y$ のビットのパリティと考えると便利なことがある。 $1.$\n",
"\n",
"この表記法を手にして、我々は今、Bernstein-Vazirani問題を定義することができる。\n",
"\n",
"\n",
" 入力：関数 $f:\\{0,1\\}^n\\rightarrow\\{0,1\\}$ \\ 約束：すべての $x\\in\\Sigma^n$ に対して $f(x) = s\\cdot x$ が成り立つバイナリ文字列 $s = s_{n-1} \\cdots s_0$ が存在する。 \\ 出力：文字列 $s$\n",
"\n",
"\n",
"この問題には、実は新しい量子アルゴリズムは必要ない。ドイチュ・ヨッサ・アルゴリズムが解決してくれる。\n",
"わかりやすくするために、ORを計算するという古典的な後処理を含まない上記の量子回路を、 *ドイチュ・ヨッサ回路と*呼ぶことにしよう。\n",
"\n",
"\n",
"\n",
"### アルゴリズム分析\n",
"\n",
"Bernstein-Vazirani問題の約束を満たす関数に対してDeutsch-Jozsa回路がどのように働くかを分析するために、まず簡単な観察から始めよう。\n",
"バイナリー・ドット・プロダクトを用いると、 $n$ ハダマード・ゲートの標準基底状態（ $n$ クビット）に対する作用を次のように記述することができる。\n",
"\n",
"$$\n",
"H^{\\otimes n} \\vert x \\rangle = \\frac{1}{\\sqrt{2^n}} \\sum_{y\\in\\Sigma^n} (-1)^{x\\cdot y} \\vert y\\rangle\n",
"$$\n",
"\n",
"Deutschのアルゴリズムを分析したときに見たのと同様、これは、任意の整数 $k$ に対する値 $(-1)^k$ が、 $k$ が偶数か奇数かにのみ依存するからである。\n",
"\n",
"Deutsch-Jozsa回路に目を向けると、第1層のハダマードゲートが実行された後、 $n+1$ の量子ビットの状態は次のようになる\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{\\sqrt{2^n}} \\sum_{x \\in \\Sigma^n} \\vert x \\rangle.\n",
"$$\n",
"\n",
"その後、クエリーゲートが実行され、（位相キックバック現象によって）状態は次のように変換される\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{\\sqrt{2^n}} \\sum_{x \\in \\Sigma^n} (-1)^{f(x)} \\vert x \\rangle.\n",
"$$\n",
"\n",
"ハダマードゲートの層の作用に関する公式を使えば、ハダマードゲートの第2層はこの状態を次のように変換することがわかる\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{2^n}\n",
"\\sum_{x \\in \\Sigma^n} \\sum_{y \\in \\Sigma^n} (-1)^{f(x) + x \\cdot y} \\vert y \\rangle.\n",
"$$\n",
"\n",
"ここで、和の中の指数 $-1$、いくつかの簡略化を行うことができる。\n",
"$f(x) = s\\cdot x$、文字列 $s = s_{n-1} \\cdots s_0,$、次のように表現できる\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{2^n}\n",
"\\sum_{x \\in \\Sigma^n} \\sum_{y \\in \\Sigma^n} (-1)^{s\\cdot x + x \\cdot y} \\vert y \\rangle.\n",
"$$\n",
"\n",
"$s\\cdot x$ と $x\\cdot y$ は2進数値なので、加算を排他的論理和に置き換えることができる。 $-1$ の指数で整数にとって重要なのは、偶数か奇数かということだけだからである。\n",
"バイナリー・ドット・プロダクトの対称性を利用して、このような表現が得られる：\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{2^n}\n",
"\\sum_{x \\in \\Sigma^n} \\sum_{y \\in \\Sigma^n} (-1)^{(s\\cdot x) \\oplus (y \\cdot x)} \\vert y \\rangle.\n",
"$$\n",
"\n",
"(わかりやすくするために括弧がつけられているが、2進のドット積は排他的論理和よりも優先順位が高いものとして扱われるのが通例であるため、実際には必要ない)\n",
"\n",
"この時点では、以下の公式を使用する。\n",
"\n",
"$$\n",
"(s\\cdot x) \\oplus (y \\cdot x) = (s \\oplus y) \\cdot x\n",
"$$\n",
"\n",
"この公式は、ビットに関する同様の公式によって求めることができる、\n",
"\n",
"$$\n",
"(a c) \\oplus (b c) = (a \\oplus b) c,\n",
"$$\n",
"\n",
"バイナリのドット積とビット単位の排他的論理和の展開と一緒に：\n",
"\n",
"$$\n",
"\\begin{aligned}\n",
"(s\\cdot x) \\oplus (y \\cdot x)\n",
"& = (s_{n-1} x_{n-1}) \\oplus \\cdots \\oplus (s_{0} x_{0}) \\oplus\n",
"(y_{n-1} x_{n-1}) \\oplus \\cdots \\oplus (y_{0} x_{0}) \\\\\n",
"& = (s_{n-1} \\oplus y_{n-1}) x_{n-1} \\oplus \\cdots \\oplus (s_{0} \\oplus y_{0}) x_{0} \\\\\n",
"& = (s \\oplus y) \\cdot x\n",
"\\end{aligned}\n",
"$$\n",
"\n",
"これにより、測定直前の回路の状態をこのように表現することができる：\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{2^n}\n",
"\\sum_{x \\in \\Sigma^n} \\sum_{y \\in \\Sigma^n} (-1)^{(s\\oplus y)\\cdot x} \\vert y \\rangle.\n",
"$$\n",
"\n",
"最後のステップは、すべてのバイナリ文字列に対して機能する、さらに別の公式を使うことである。 $z = z_{n-1}\\cdots z_0.$\n",
"\n",
"$$\n",
"\\frac{1}{2^n}\n",
"\\sum_{x \\in \\Sigma^n} (-1)^{z \\cdot x}\n",
"= \\begin{cases}\n",
"1 & \\text{if $z = 0^n$}\\\\\n",
"0 & \\text{if $z\\neq 0^n$}\n",
"\\end{cases}\n",
"$$\n",
"\n",
"ここでは、このレッスンであと何度か使うことになる、文字列の簡単な表記法を使う。 $0^n$ は、長さがすべてゼロの文字列である。 $n.$\n",
"\n",
"この公式が機能することを論証する簡単な方法は、2つのケースを別々に考えることである。\n",
"もし $z = 0^n,$ ならば、すべての文字列に対して $z\\cdot x = 0$ $x\\in\\Sigma^n,$。したがって、和の各項の値は $1,$ となり、和と除算によって $1$ が得られる。 $2^n.$ $x$ 一方、 $z$ のいずれかのビットが $1,$ に等しい場合、バイナリドット積 $z\\cdot x$ は、 $x\\in\\Sigma^n$ の選択肢のちょうど半分について $0$ に等しく、残りの半分について $1$ に等しくなります。 $z\\cdot x$ の値が反転する（ $0$ から $1$ へ、または $1$ から $0$ へ）のは、 $z$ が $1.$\n",
"\n",
"ここで、この式を測定前の回路の状態を単純化するために適用すると、次のようになる\n",
"\n",
"$$\n",
"\\vert - \\rangle \\otimes \\frac{1}{2^n}\n",
"\\sum_{x \\in \\Sigma^n} \\sum_{y \\in \\Sigma^n} (-1)^{(s\\oplus y)\\cdot x} \\vert y \\rangle\n",
"= \\vert - \\rangle \\otimes \\vert s \\rangle,\n",
"$$\n",
"\n",
"という事実によって、 $s\\oplus y = 0^n$ if and only if $y = s.$ したがって、この測定は、我々が探しているストリング $s$ を正確に明らかにしている。\n",
"\n",
"\n",
"\n",
"### 古典的難度\n",
"\n",
"Deutsch-Jozsa回路は1回の問い合わせでBernstein-Vazirani問題を解くが、古典的な問い合わせアルゴリズムはこの問題を解くために少なくとも $n$。\n",
"\n",
"これは、いわゆる*情報理論的な*議論によって推論することができるが、この場合は非常に単純である。\n",
"古典的なクエリーは、それぞれ解に関する1ビットの情報を明らかにする。 $n$ ビットの情報を明らかにする必要があるので、少なくとも $n$ クエリーが必要になる。\n",
"\n",
"実際、バーンスタイン-ヴァジラニ問題を古典的に解くことは可能である。それは、 $n$ の文字列のそれぞれについて、 $1,$ を1つずつ、それ以外のすべてのビットについては $0$ を持つ関数を問い合わせることで、 $s$ のビットを1つずつ明らかにすることである。\n",
"したがって、この問題に対する古典的アルゴリズムに対する量子の利点は、 $1$ クエリ対 $n$ クエリである。\n",
"\n"
]
},
{
"cell_type": "markdown",
"id": "28583734",
"metadata": {},
"source": [
"\n",
"\n",
"## バーンスタイン＝ヴァジラニ問題とQiskit\n",
"\n",
"ドイチュ・ヨツァ回路はすでに実装したが、ここではそれを使ってバーンスタイン・ヴァジラニ問題を解く。\n",
"まず、任意の2進文字列が与えられた場合のBernstein-Vazirani問題のクエリーゲートを実装する関数を定義する。 $s.$\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 18,
"id": "11861a7e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"def bv_query(s):\n",
" # Create a quantum circuit implementing a query gate for the\n",
" # Bernstein-Vazirani problem.\n",
"\n",
" qc = QuantumCircuit(len(s) + 1)\n",
" for index, bit in enumerate(reversed(s)):\n",
" if bit == \"1\":\n",
" qc.cx(index, len(s))\n",
" return qc\n",
"\n",
"\n",
"display(bv_query(\"1011\").draw(output=\"mpl\"))"
]
},
{
"cell_type": "markdown",
"id": "c24e91e8",
"metadata": {},
"source": [
"ここで、先に定義した `compile_circuit` 関数を使って、関数上でDeutsch-Jozsa回路を実行する関数を作ることができる。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 19,
"id": "7db2ea99",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"'1011'"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"def bv_algorithm(function: QuantumCircuit):\n",
" qc = compile_circuit(function)\n",
" result = AerSimulator().run(qc, shots=1, memory=True).result()\n",
" return result.get_memory()[0]\n",
"\n",
"\n",
"display(bv_algorithm(bv_query(\"1011\")))"
]
},
{
"cell_type": "markdown",
"id": "86c52e37",
"metadata": {},
"source": [
"\n",
"\n",
"### 命名法に関する注記\n",
"\n",
"Bernstein-Vazirani問題の文脈では、Deutsch-Jozsaアルゴリズムが \"Bernstein-Vaziraniアルゴリズム \"と呼ばれるのが一般的である\n",
"というのも、BernsteinとVaziraniが彼らの研究で明確にしていたように\\*、\\* このアルゴリズムはDeutsch-Jozsaアルゴリズムだからだ。\n",
"\n",
"Deutsch-JozsaアルゴリズムがBernstein-Vazirani問題を解くことを示した後、BernsteinとVaziraniが行ったことは（上記の通り）、 *再帰的フーリエサンプリング*問題として知られる、より複雑な問題を定義することだった。\n",
"これは高度に仕組まれた問題で、問題のさまざまなインスタンスに対する解答が、ツリー状に配置された問題の新たなレベルを効果的に解き放つ。\n",
"Bernstein-Vazirani問題は本質的に、この複雑な問題の基本ケースに過ぎない。\n",
"\n",
"再帰的フーリエサンプリング問題は、量子アルゴリズムが確率的アルゴリズムに対していわゆる*超多項式的な*優位性を持ち、Deutsch-Jozsaアルゴリズムによる古典的アルゴリズムに対する量子の優位性を凌駕する問い合わせ問題の最初の例として知られている。\n",
"直観的に言えば、この問題の再帰的バージョンは、量子アルゴリズムの $1$ 対 $n$ の利点を、より大きなものに増幅する。\n",
"\n",
"この優位性を確立する数学的分析で最も困難な点は、古典的なクエリーアルゴリズムでは多くのクエリーを行わないと問題を解決できないことを示すことである。\n",
"これは極めて典型的な例である。多くの問題において、それらを効率的に解決する創造的な古典的アプローチを除外することは非常に難しい。\n",
"\n",
"サイモンの問題と、次のセクションで説明するそのアルゴリズムは、古典的アルゴリズムに対する量子の超多項式的な（実際には指数関数的な）利点を示す、より単純な例を提供している。\n",
"とはいえ、それ自体が興味深い計算問題ではある。\n",
"\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"id": "a1b8767d",
"source": "© IBM Corp., 2017-2026"
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3"
}
},
"nbformat": 4,
"nbformat_minor": 5
}