この記事は、量子サポートベクターマシーンを利用して、MNISTから配布されているデータベースの分類を行う方法をある程度略して説明するものです。変分量子計算の方法論は、機械学習における分類問題にも応用されています。その1つがQuantum Support Vector Machine(QSVM)です[1]。これはクラス分類問題において、評価関数に現れるカーネルを量子状態の内積で代用して、それを解くというものです。この方法において、過学習が抑えられるということが報告されています。本来ならば、実装してデータを採りたいところですが、手書き文字の認識分類をさせるだけでも画像を取り込んだ上に数字とアルファベットの全36種類の分類を行わなければならないため、私のパソコン環境では実装出来ないと判断しました。そのため、理論の話にとどめます。まずはタイトル通り、機械学習における2クラス分類の基礎から説明します。

機械学習の基礎:

機械学習の基本は、与えられたデータxに対するクラス分類を行うことです。もっとも簡単な例はn次元平面上における一次関数における評価関数の値の範囲の中に有るか無いかを判定する問題です。入力するデータを $\bm{x}$ として式で表せば、

$0=\bm{w}^\dagger \bm{x} + b$ -(1)

となります。分類のうちどちらかを表す空間 $C\in\{+1,-1\}$ 、i番目の入力データに対応するその分類を表す符号変数 $y_i\in\{+1,-1\}$ を導入し、 $\bm{w}$ のノルムで割りマージンを取ることで、

$y_i(\bm{w}^\dagger \bm{x}_i + b)\geq 1, ~ i\in\{1,2,\cdots,t\}$ -(2)

の関係が導かれます。ここで、tは入力されるトレーニング用データの数です。こうして、データセット $\bm{x}(\bm{w},b)$ に対する評価関数の結果はかっこの中身の符号を返す関数sign()を使って、

$m(\bm{x}(\bm{w}, b))=sign(\bm{w}^\dagger \bm{x} + b)$ -(3)

となります。これは異なる分類に属する距離（マージン)の最も大きい2点間をつなぐサポートベクトルのノルムwを最小化する問題です。そのため、ラグランジュの未定乗数法と同じ形に表すことができ、

$L_p=\frac{1}{2}\mid\mid \bm{w} \mid\mid - \sum_{i=1}^t \alpha_i y_i(\bm{w}^\dagger \bm{x}_i + b) + \sum_{i=1}^t \alpha_i$ -(4)

となります。 $\partial L_p /\partial b = \partial L_p / \partial w = 0$ となるため計算機での処理が可能なように、

$L_D=\sum_{i=1}^t\alpha_i-\frac{1}{2}\sum_{i,j=1}^t \alpha_i \alpha_j y_i y_j \bm{x}_i^\dagger \bm{x}_j$ -(5)

と変換できます。この関数は分類数空間を $C\in \{1,2,\cdots,c\}$ として束縛条件 $0\leq \alpha \leq c$ , $\sum_{i=1}\alpha_i y_i = 0$ に従います。この関数はコスト関数と呼ばれる、多クラス分類において最も重要な式です。この関数を最大化する $\alpha$ の組を見つけることでデータセットを正しく分類することがクラス分類の目的です。この中の $\bm{x}_i^\dagger \bm{x}_j$ は二値カーネルと呼ばれるもので、 $K(\bm{x}_i,\bm{x}_j)$ と表されます。これを問題の特徴に応じて適切なものを選ぶことで、計算の精度が格段に上がることが知られています[2]。QSVMにおいては、これは異なるデータセットから作った状態間の内積 $\mid \langle \Phi(\bm{x}_i)\mid \Phi(\bm{x}_j)\rangle \mid^2$ となります。

量子多クラス分類:

この節では、いよいよ量子多クラス分類の説明をしていきます。多少複雑ですが、それぞれのプロセスで行うことはそれほど複雑ではありません。このアルゴリズムの工程は主に3つです。

1, 量子多クラス分類機にトレーニング用データセットt子を学習させます。

2, 式(5)のコスト関数を計算します。

3, 量子多クラス分類機における1,で最適化された変数θと経験的コスト関数 $R_{emp}$ を基に、テストデータセットを分類させます。

1,と3,の大まかな工程を図1,2にそれぞれ示します。

図1 トレーニング用データセットを学習させるプログラムの大まかな流れ。 $\theta$ は最適化変数、 $W(\theta)$ はハミルトニアンに相当する量子回路です。

図2 学習した最適化変数と経験的コスト関数を利用して $s\in S$ の分類を行う流れ。

2クラス分類の場合はそれらの図において、

$y\in \{+1,-1\},~f:\{0,1\}\rightarrow \{+1,-1\}, M_y=\frac{1}{2}(I+y\sum_{z\in\{0,1\}^n}f(z)\mid z \rangle \langle z \mid )$ -(6)

となります。 $\mid \Phi(\bm{x})\rangle$ を実現する回路は $R_z(\theta)$ ゲート、 $exp(i\theta Z_i Z_j)$ の回路で実装します。 $W(\theta)$ も同様です。

$p_y(\bm{x})=\langle \Phi(\bm{x}) \mid W^\dagger(\theta)M_yW(\theta)\mid \Phi(\bm{x}) \rangle$ -(7)

から、

$R_{emp}(\theta)=\frac{1}{T}\sum_{x\in T}Pr(\tilde{m}(\bm{x})\neq m(\bm{x}))$ -(8)

を計算し、経験的コスト関数 $R_{emp}(\theta)$ と変数 $\theta$ を最適化します。2, は2クラス分類の場合は少々異なり、量子ビット数をnとしてパウリゲート $P_{\alpha},\alpha \in \{1,2,\cdots,4n\}$ に対して、

$w_{\alpha}(\theta)=Tr[W^\dagger(\theta)\bm{f}W(\theta)P_{\alpha}]$ -(9)

と、

$\Phi_{\alpha}(\theta)=\langle \Phi(\bm{x})\mid P_\alpha \mid \Phi(\bm{x})\rangle$ -(10)

から、

$m(\bm{x})=sign(\frac{1}{2^n}\sum_{\alpha}w_{\alpha}(\theta)\Phi_{\alpha}(\bm{x})+b)$ -(11)

を量子ビット数nに対して計算します。また、 $m(\bm{x}_i)=y_i$ とします。3,の量子カーネルはSWAP-test[3]で計算し、式(11)を $\alpha_j^{*}$ とbの連立方程式として解くことで、計算出来ます。多クラス分類の場合は式(6)を、 $M_y$ を $y$ を2進数で表した場合のi番目のビット $y^i$ を用いて $[g_i,g_j]=0$ の関係を満たすパウリゲートによって表し、

$y\in \{1,2,\cdots,c \},~f:\{0,1\}\rightarrow \{+1,-1\}, M_y=\frac{1}{2}(I - \Pi_{i=1}^{log_2 c}g_i^{y^i})$ -(12)

としてあとは1対他あるいは1対1戦略で逐次分類することで図1,2の通りの方法で学習が完了します。 $R_{emp}=P_{err}$ となります。MNISTデータセットも、計算時間と仕様ゲート数は2値分類に比べて大幅に増加しますが、この方法で分類が可能です。

編集後記:

今回は私自身機械学習はタイトル通りゼロから実装を学んだため、非常にざっくりとした説明になってしまいました。もしかしたら、内容に誤りがあるかもしれません。しかしながら、変分量子計算の応用の中で一番難しい部類に入るアルゴリズムの論文を解読し、その内容を分かりやすく説明することは、私にとっても有意義でした。量子機械学習はFault-Tolerant-Quantum-Computer(FTQC)における最有力応用候補として期待されており、更なる発展が期待されています。そしてすでに、量子化学計算においてもクラスターを機械学習で最適化する手法が出ています。なので、私もこれを機に、量子機械学習の研究を始めようと思います。

[1] Vojtech Havlicek, Antonio D. Córcoles, Kristan Temme, Aram W. Harrow, Abhinav Kandala, Jerry M. Chow, Jay M. Gambetta, arxiv, 1804.11326[quant-ph](2018)

[2] Qiskitで量子SVMを実装して性能評価してみた - Qiita

[3] Juan Carlos Garcia-Escartin and Pedro Chamorro Posada, arxiv, 1304.4923[quant-ph](2013)

Re:ゼロから始める量子多クラス分類

Hikaru Wakaura