特許 7451606 学習システム及び学習方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-03-08

(45)【発行日】2024-03-18

(54)【発明の名称】学習システム及び学習方法

(51)【国際特許分類】

   G06N 10/60 20220101AFI20240311BHJP

   G06F 18/2433 20230101ALI20240311BHJP

   G06N 3/0475 20230101ALI20240311BHJP

   G06N 20/00 20190101ALI20240311BHJP

【ＦＩ】

G06N10/60

G06F18/2433

G06N3/0475

G06N20/00 130

【請求項の数】 5

(21)【出願番号】P 2022101100

(22)【出願日】2022-06-23

(65)【公開番号】P2024002105

(43)【公開日】2024-01-11

【審査請求日】2022-06-23

(73)【特許権者】

【識別番号】592131906

【氏名又は名称】みずほリサーチ＆テクノロジーズ株式会社

(73)【特許権者】

【識別番号】303028745

【氏名又は名称】株式会社みずほフィナンシャルグループ

(74)【代理人】

【識別番号】100105957

【弁理士】

【氏名又は名称】恩田 誠

(74)【代理人】

【識別番号】100068755

【弁理士】

【氏名又は名称】恩田 博宣

(72)【発明者】

【氏名】宇野 隼平

(72)【発明者】

【氏名】手塚 宙之

【審査官】渡辺 一帆

(56)【参考文献】

【文献】特開２０２２－００２０８０（ＪＰ，Ａ）

【文献】特許第７００５８０６（ＪＰ，Ｂ１）

【文献】米国特許出願公開第２０２２／０１２１５４０（ＵＳ，Ａ１）

【文献】DALLAIRE-DEMERS, P et al.，"Quantum generative adversarial networks"，arXiv.org [online]，2018年，pp. 1-10，[retrieved on 2023.08.23], Retrieved from the Internet: <URL: https://arxiv.org/abs/1804.08641v2>，<DOI: 10.48550/arXiv.1804.08641>

【文献】御手洗 光祐，"変分量子アルゴリズムと量子機械学習"，電子情報通信学会誌，一般社団法人電子情報通信学会，2021年，第104巻, 第11号，pp. 1166-1173，ISSN 0913-5693

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｎ １０／００－１０／８０

Ｇ０６Ｆ １８／２４３３

Ｇ０６Ｎ ３／０２－３／１０

Ｇ０６Ｎ ２０／００－２０／２０

(57)【特許請求の範囲】

【請求項1】

量子状態の操作部と、前記量子状態を保持する状態保持部と、前記量子状態を観測する計測部とを含む量子計算部と、古典計算部と、前記量子計算部及び前記古典計算部を制御する管理部と、を備え、量子データセットの分布を学習する学習システムであって、
前記古典計算部が、潜在空間上の複数のサンプルを、超平面上に写像するジェネレータで用いる量子回路パラメータの値を保持し、
前記量子計算部が、ユニタリ演算子により構成された前記ジェネレータを用いて、前記潜在空間上の潜在変数の複数のサンプルと前記量子データセットとから、前記量子回路パラメータにおけるground cost及びその勾配を算出し、
前記古典計算部が、前記ground cost及び前記勾配を用いて、最適輸送損失及びその勾配を算出し、前記最適輸送損失の勾配を用いて、前記最適輸送損失が小さくなるように、前記量子回路パラメータを更新する更新処理を実行し、
前記管理部が、前記更新処理を前記最適輸送損失が収束するまで繰り返すことを特徴とする学習システム。

【請求項2】

前記管理部が、前記分布に応じて、新たに取得した量子データの異常判定を行なうことを特徴とする請求項１に記載の学習システム。

【請求項3】

前記管理部が、教師データとなる量子データセットにおける分布と、新たに取得した量子データとを比較することにより、前記異常判定を行なうことを特徴とする請求項２に記載の学習システム。

【請求項4】

前記ユニタリ演算子は、前記最適輸送損失の勾配を算出するために、前記量子回路パラメータで微分可能であることを特徴とする請求項１～３の何れか一項に記載の学習システム。

【請求項5】

量子状態の操作部と、前記量子状態を保持する状態保持部と、前記量子状態を観測する計測部とを含む量子計算部と、古典計算部と、前記量子計算部及び前記古典計算部を制御する管理部と、を備えた学習システムを用いて、量子データセットの分布を学習する方法であって、
前記古典計算部が、潜在空間上の複数のサンプルを、超平面上に写像するジェネレータで用いる量子回路パラメータの値を保持し、
前記量子計算部が、ユニタリ演算子により構成された前記ジェネレータを用いて、前記潜在空間上の潜在変数の複数のサンプルと前記量子データセットとから、前記量子回路パラメータにおけるground cost及びその勾配を算出し、
前記古典計算部が、前記ground cost及び前記勾配を用いて、最適輸送損失及びその勾配を算出し、前記最適輸送損失の勾配を用いて、前記最適輸送損失が小さくなるように、前記量子回路パラメータを更新する更新処理を実行し、
前記管理部が、前記更新処理を前記最適輸送損失が収束するまで繰り返すことを特徴とする学習方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、量子データセットの分布を学習するための学習システム及び学習方法に関する。

【背景技術】

【0002】

量子コンピュータを用いた機械学習(量子機械学習)はＮＩＳＱ（Noisy Intermediate-Scale Quantum Computer）において有力視されている適用領域の一つである。量子機械学習には、量子データを扱うものと、古典データを扱うものがあり、いずれも積極的に研究が行なわれている。量子機械学習は、特に、量子データを用いた場合、既存のコンピュータを用いた機械学習(古典機械学習)と比較して大きな優位性を有することが期待されている。量子機械学習は、古典機械学習と同様に、教師あり学習や教師なし学習や強化学習等に分類される。

【0003】

生成モデルは、教師なし学習の主要なタスクの一つである。与えられたデータセットをもとに、そのデータの背後にある分布を近似することを目的とする。より具体的には、与えられたデータセットの背後にある確率分布をα(x)と仮定する。そして、確率モデルをθでパラメトライズされた確率分布β_θ(x)としたときに、適切な損失関数Ｌ(α,β_θ)をなるべく小さくするようにθを学習する。古典機械学習において、生成モデルは異常検知、データの生成、データ構造の学習、次元削減などの目的で広く研究されているが、量子コンピュータへの適用事例はほとんどない。

【0004】

生成モデルは大きく、明示的モデル（prescribed model）と暗黙的モデル（implicit model）に分類される。明示的モデルでは、パラメータθを用いて、データセットの背後にある確率密度を関数として明示的に記述する。このため、明示的モデルでは、対数尤度を計算することが可能である。そして、最尤推定などにより、パラメータθを求めることが可能となる。一方で、暗黙的モデルでは、確率密度を明示的に記述せずに、生成モデルに従うサンプルを生成する。暗黙的モデルの利点としては、広大な全標本空間のうち、比較的低次元の超平面上にデータセットが分布するような現象を容易に表すことができること、データを生成する物理的なプロセスとして解釈できること等がある。

【0005】

暗黙的モデルでは、通常、全標本空間に比べて、十分に次元が低い潜在空間上のランダム変数(潜在変数)を仮定する。そして、潜在変数を標本空間に移す写像(ジェネレータ)を用いて生成モデルを表す。ここで、潜在変数はガウス分布や一様分布等の既知の分布に従うと仮定される。暗黙的モデルの学習では、適当な損失関数を用いて、データセットとジェネレータの生成する分布とを適合させるように、写像のパラメータθの学習を行なう。

【0006】

古典機械学習における生成モデルでは、損失関数として最適輸送損失／ワッサーシュタイン距離（Optimal Transportation loss/Wasserstein distance）が注目されており、画像解析、自然言語処理、金融等様々な分野で活用されている。最適輸送損失は、特に、尤度が計算できない場合や、確率分布のサポートが一致しない場合にも適用可能なこと、標本空間上の距離を自然に取り込めること等の利点がある。最適輸送損失は、確率分布α(x)を別の確率分布β(x)に移す最小のコストとして、以下で定義される。

【0007】

【数1】

ここで、c(x,y)≧0は、ｘからｙへの輸送コストを表しており、ground costと呼ばれる。また、ground costを最小化するπ（経路）の集合を最適輸送計画と呼ぶ。

【0008】

量子生成モデルを構築するにあたり、（１）ジェネレータ及び（２）損失関数に関して検討する必要がある。
量子データのジェネレータとして、パラメータ付きの量子回路(Parameterized Quantum Circuit、PQC)を用いることが考えられる。量子コンピュータを用いた既存の生成モデル(例えば量子ＧＡＮ等（非特許文献１を参照）)では、PQC U(θ)を用いて、以下のようにして、ジェネレータを定義することが提案されている。

【0009】

【数2】

ここで、｜Ψ_θ(z)>は、量子状態の分布である。
この回路では、わずかに潜在変数が異なるような量子状態同士が直交する。このため、量子空間上における連続的な分布を表すことが困難である。

【0010】

この問題を解決するため、ジェネレータとしては以下のように潜在変数も含めたPQC G(θ,z)を用いることとする。

【0011】

【数3】

例えば、図５に示す量子回路Ｃ１をジェネレータG(θ,z)として設定することが可能である。図５において、「N_L」は、ジェネレータG(θ,z)を構成する階層の数である。このジェネレータG(θ,z)では、同じ構造を持つ階層をN_L回繰り返し作用する。i番目の階層内ではまず、j番目のビットに角度θ_i,j× zη_i,j、方向k_i,jの単一量子ビットのパウリ回転を作用し、その後、全ビットに対するladder controlled Zゲートを作用する。ここで、回転方向k_i,jおよびzの方向η_i,jはランダムに選ばれる。

【0012】

次に、損失関数について述べる。最適輸送損失に関して、ジェネレータにより生成される量子状態のサンプルと与えられた量子データの間のground cost c(x,y)について検討する必要がある。最も単純には、ground costとして、量子状態間の重なりをもとにしたトレース距離を採用することが考えられるが、この距離を用いるとビット数が大きい大規模な系では損失関数の最小化が進まない問題(barren plateau問題)がある。

【0013】

この問題を解決する方法として、量子支援量子コンパイル（ＱＡＱＣ）の中で提案されているコスト関数がある（非特許文献２を参照）。ＱＡＱＣでは、個別ビットの出力結果を用いてコスト関数を構成することで、barren plateau問題を回避する方法が提案されている。

【0014】

学習においては、量子回路学習（非特許文献３を参照）と同様に、PQCのパラメータを調整することで、損失関数の最小化を行なう。
損失関数を最小化するためには、PQCのパラメータに対する損失関数の勾配を得ることが重要である。損失関数の勾配を、元の回路とほぼ同じ構造の回路を用いて推定する方法(パラメータシフトルール)が提案されている（非特許文献３、４を参照）。

【先行技術文献】

【非特許文献】

【0015】

【文献】Pierre-Luc Dallaire-Demers and Nathan Killoran，「Quantum generative adversarial networks」、［online］、American Physical Society、２０１８年７月２３日、［令和４年５月１９日検索］、インターネット＜https://journals.aps.org/pra/abstract/10.1103/PhysRevA.98.012324＞

【文献】Sumeet Khatri，Ryan LaRose，Alexander Poremba，Lukasz Cincio，Andrew T Sornborger，Patrick J Coles，「Quantum-assisted quantum compiling」、［online］、arXiv.org、Cornell University、２０１９年５月７日、［令和４年５月１日検索］、インターネット＜https://arxiv.org/abs/1807.00800＞

【文献】Kosuke Mitarai，Makoto Negoro，Masahiro Kitagawa，Keisuke Fujii、「Quantum circuit learning」、［online］、arXiv.org、Cornell University、２０１９年４月２４日、［令和４年５月１日検索］、インターネット＜https://arxiv.org/abs/1803.00745＞

【文献】Maria Schuld，Ville Bergholm，Christian Gogolin，Josh Izaac，Nathan Killoran、「Evaluating analytic gradients on quantum hardware」、［online］、arXiv.org、Cornell University、２０１８年１１月２７日、［令和４年５月１日検索］、インターネット＜https://arxiv.org/abs/1811.11184＞

【発明の概要】

【発明が解決しようとする課題】

【0016】

量子データを生成する量子生成モデルを構築するにあたり、適切なジェネレータ及び損失関数が検討されておらず、量子データセットの分布の学習が困難であった。

【課題を解決するための手段】

【0017】

本開示の学習システムは、量子状態の操作部と、前記量子状態を保持する状態保持部と、前記量子状態を観測する計測部とを含む量子計算部と、古典計算部と、前記量子計算部及び前記古典計算部を制御する管理部と、を備え、量子データセットの分布を学習する。前記古典計算部が、潜在変数の複数のサンプルと量子回路パラメータの値からジェネレータの構造を決定する処理を実行し、前記量子計算部が、前記ジェネレータと前記量子データセットからground cost及びその勾配を算出し、前記古典計算部が、前記ground cost及び前記勾配を用いて、最適輸送損失及びその勾配を算出し、前記最適輸送損失の勾配を用いて、前記最適輸送損失が小さくなるように、前記量子回路パラメータを更新する更新処理を実行し、前記管理部が、前記更新処理を前記最適輸送損失が収束するまで繰り返す。

【発明の効果】

【0018】

本発明によれば、量子データセットの分布を学習することができる。

【図面の簡単な説明】

【0019】

【図1】実施形態のシステム概略図である。

【図2】実施形態のハードウェア構成の説明図である。

【図3】実施形態のデータ空間における超平面の説明図である。

【図4】実施形態の潜在空間からデータ空間への写像の説明図である。

【図5】実施形態の量子回路の説明図である。

【図6】実施形態の量子データセットの分布の学習処理手順の説明図である。

【図7】実施形態の量子状態の異常検知処理手順の説明図である。

【発明を実施するための形態】

【0020】

以下、図１～図７を用いて、学習システム及び学習方法の一実施形態を説明する。本実施形態では、外部から入力された量子データセットを基に、量子データセットの背後にある分布を学習する。
このため、図１に示すように、ネットワークで接続されたユーザ端末１０、演算装置２０を用いる。

【0021】

（ハードウェア構成の説明）
図２を用いて、ユーザ端末１０、演算装置２０（管理部２１１、古典計算部２１２）を構成する情報処理装置Ｈ１０のハードウェア構成を説明する。情報処理装置Ｈ１０は、通信装置Ｈ１１、入力装置Ｈ１２、表示装置Ｈ１３、記憶装置Ｈ１４、プロセッサＨ１５を備える。なお、このハードウェア構成は一例であり、他のハードウェアにより実現することも可能である。

【0022】

通信装置Ｈ１１は、他の装置との間で通信経路を確立して、データの送受信を実行するインタフェースであり、例えばネットワークインタフェースや無線インタフェース等である。

【0023】

入力装置Ｈ１２は、利用者等からの入力を受け付ける装置であり、例えばマウスやキーボード等である。表示装置Ｈ１３は、各種情報を表示するディスプレイ等である。
記憶装置Ｈ１４は、ユーザ端末１０、演算装置２０の各種機能を実行するためのデータや各種プログラムを格納する装置である。記憶装置Ｈ１４の一例としては、ＲＯＭ、ＲＡＭ、ハードディスク等がある。

【0024】

プロセッサＨ１５は、記憶装置Ｈ１４に記憶されるプログラムやデータを用いて、ユーザ端末１０、演算装置２０における各処理を制御する。プロセッサＨ１５の一例としては、例えばＣＰＵやＭＰＵ等がある。このプロセッサＨ１５は、ＲＯＭ等に記憶されるプログラムをＲＡＭに展開して、各サービスのための各種プロセスを実行する。

【0025】

プロセッサＨ１５は、自身が実行するすべての処理についてソフトウェア処理を行なうものに限られない。例えば、プロセッサＨ１５は、自身が実行する処理の少なくとも一部についてハードウェア処理を行なう専用のハードウェア回路（例えば、特定用途向け集積回路）を備えてもよい。すなわち、プロセッサＨ１５は、以下を含む回路として構成し得る。

【0026】

（１）コンピュータプログラムに従って動作する１つ以上のプロセッサ
（２）各種処理のうち少なくとも一部の処理を実行する１つ以上の専用のハードウェア回路
（３）それらの組み合わせ
プロセッサは、ＣＰＵ並びに、ＲＡＭ及びＲＯＭ等のメモリを含み、メモリは、処理をＣＰＵに実行させるように構成されたプログラムコード又は指令を格納している。メモリすなわちコンピュータ可読媒体は、汎用又は専用のコンピュータでアクセスできるあらゆる利用可能な媒体を含む。

【0027】

（システム構成）
次に、図１を用いて、ユーザ端末１０、演算装置２０のシステム構成を説明する。
ユーザ端末１０は、利用者が用いるコンピュータ端末である。

【0028】

演算装置２０は、量子データセットの分布を学習するためのコンピュータである。
この演算装置２０は、後述する処理（管理段階、量子計算段階、古典計算段階等の各処理等）を行なう。そのためのプログラムを実行することにより、演算装置２０は、管理部２１１、古典計算部２１２、量子計算部２２として機能する。

【0029】

管理部２１１は、量子計算及び古典計算の管理を行なう。
古典計算部２１２は、潜在変数の複数のサンプルと量子回路パラメータの値からジェネレータの構造を決定する。このジェネレータは、微分可能なユニタリ演算子により構成する。更に、古典計算部２１２は、量子回路パラメータを更新する更新処理を実行する。

【0030】

量子計算部２２は、量子計算を行なう。本実施形態では、量子計算部２２は、潜在空間の潜在変数に依存する写像（ジェネレータ）を実行する。更に、量子計算部２２は、外部から入力された量子データセット及びジェネレータからground cost及びその勾配を算出する。

【0031】

この量子計算部２２は、操作部２２１、状態保持部２２２、計測部２２３を含む。
操作部２２１は、状態保持部２２２の量子ビットに対して量子演算に応じた量子操作を行なう。この場合、操作部２２１は、量子ゲート等により構成される量子回路を用いることによって、状態保持部２２２が保持する状態を操作（作成）する。

【0032】

状態保持部２２２は、複数の量子ビットを備え、任意の量子状態を保持する。各量子ビットは、電子準位、電子スピン、イオン準位、各スピン、光子等、任意の物理状態で複数の値の重ね合わせ状態を保持する。重ね合わせ状態を保持できれば、量子ビットは上記に限定されるものではない。

【0033】

計測部２２３は、状態保持部２２２の量子ビットの重ね合わせ状態を計算基底で観測する。計測部２２３は、状態保持部２２２の量子ビットの状態に応じて、ヒット回数を記録する。

【0034】

（量子データセットの分布の学習処理）
図６を用いて、量子データセットの分布の学習処理を説明する。
ここでは、N_r個の量子データからなる量子データセット{|ψ₁>,|ψ₂>,…,|ψ_Nr>}が与えられた時に、このデータセットの分布を近似するような生成モデルを学習する。

【0035】

図３に示すように、データ空間ＳＰ１上の量子データセットに、できる限り近くなるように、低次元多様体である超平面ＭＦ１を学習する。ここで、黒点で示したデータセットＤ１はデータ空間ＳＰ１上の量子データセットである。また、×点で示したデータセットＤ２は超平面ＭＦ１上の量子データセットである。

【0036】

図４に示すように、データセットＤ２は、潜在空間ＬＳ１上のサンプルを、ジェネレータＧ１を用いて、超平面ＭＦ１上へ写像したものである。
例えば、図５に示す量子回路Ｃ１をジェネレータG(θ,z)として設定する。このジェネレータG(θ,z)では、同じ構造を持つ階層をN_L回繰り返し作用する。i番目の階層内ではまず、j=1～n番目のすべてのビットに角度θ_i,j×zη_i,j、方向k_i,jの単一量子ビットのパウリ回転を作用し、その後、全ビットに対するladder controlled Zゲートを作用する。ここで、回転方向k_i,jおよびzの方向η_i,jはランダムに選ばれる。

【0037】

上述したように、このようなジェネレータG(θ,z)は、潜在空間ＬＳ１の分布を、超平面ＭＦ１上の分布に写像する役割を果たす。そして、ジェネレータG(θ,z)の出力が、与えられた量子データセットに、可能な範囲で近づくように、量子回路パラメータ｛θ｝を学習する。

【0038】

まず、演算装置２０は、潜在空間において既知の分布から潜在変数を生成する（ステップＳ１１）。具体的には、管理部２１１は、古典計算部２１２により、潜在空間ＬＳ１における既知の分布から、潜在変数としてサンプル数Ng個のサンプルz_jを生成する。

【0039】

【数4】

次に、演算装置２０は、ground costの推定を行なう（ステップＳ１２）。ここでは、管理部２１１は、量子計算部２２により、量子回路パラメータθと潜在変数のサンプルzjとから、量子状態のサンプル{G(θ,z₁)|0>, G(θ,z₂)|0>,…, G(θ,z_Ng)|0>}を生成する。次に、量子計算部２２により、サンプル{G(θ,z₁)|0>, G(θ,z₂)|0>,…, G(θ,z_Ng)|0>}と、量子データセット{|ψ₁>,|ψ₂>,…,|ψ_Nr>}とから、ground cost c_localを計算する。この計算には、非特許文献１に記載された下記式を用いる。

【0040】

【数5】

ここで、「ｎ」は量子ビットの個数、「p^(k) _i,j」はｋ番目のビットが「0」になる確率である。「Tr」は対角成分の和（トレース）、「I_i」はｉ番目のビットに作用する恒等演算子である。

【0041】

そして、量子計算部２２は、このground cost c_localをN_sショットにより近似する。ここでは、ground cost c_localの近似値を、以下のように表記する。

【0042】

【数6】

ここで、「N_r」は量子データセットのデータ数、「N_g」は潜在変数のサンプル数、「N_s」はショット数（同一の状態を用意する個数）である。
このground costの近似値は、確率（1-p^(k) _i,j）のベルヌーイ分布Be（1-p^(k) _i,j）に従う確率変数X^(s) _i,j,kを用いて以下のように定義される。

【0043】

【数7】

【0044】

そして、管理部２１１は、古典計算部２１２により、ground cost c_localの近似値を用いて、量子回路モデル「G(θ,z)」のパラメータ「θ」を固定したときの、最適輸送計画及び最適輸送損失の算出を行なう（ステップＳ１３）。

【0045】

【数8】

「Ｌ」は、最適輸送損失（Optimal Transport Loss）、「π_ij」は、データセットのi番目のデータ|ψ_i>とサンプルG(θ,z_j)|0>の間の輸送計画である。

【0046】

次に、管理部２１１は、古典計算部２１２により、上記式の線形計画法を解いて、最適輸送計画{π_i,j}及び最適輸送損失Ｌを求める（第１の最小化）。
次に、演算装置２０は、量子回路のパラメータの勾配を計算するためのパラメータシフトルールを用いてground costの勾配の算出を行なう（ステップＳ１４）。ここでは、例えば、非特許文献２、３に記載されたパラメータシフトルールを用いる。そして、ground cost c_localの勾配をN_sショットで近似する。ここで、ground cost c_localの勾配の近似値を以下のように表記する。

【0047】

【数9】

ここで、「N_p」は、量子回路パラメータ「θ」の個数である。

【0048】

次に、演算装置２０は、学習率及びground cost c_localの勾配を利用して、量子回路パラメータ｛θ｝の更新を行なう（ステップＳ１５）。ここでは、管理部２１１は、古典計算部２１２により、ground cost c_localの勾配及び最適輸送計画{π_i,j}を使って、最適輸送損失の勾配を計算する。

【0049】

【数10】

そして、管理部２１１は、計算した勾配及び学習率を用いて、量子回路パラメータ｛θ}をアップデートする（第２の最小化）。

【0050】

次に、演算装置２０は、収束するまで、ステップＳ１１からステップＳ１５の処理を繰り返す（ステップＳ１６）。ここでは、管理部２１１は、最適輸送損失の値に応じて収束を判定する。例えば、最適輸送損失が所定範囲内に含まれていない場合、ステップＳ１１以降の処理を繰り返す。
一方、管理部２１１は、最適輸送損失が所定範囲内に含まれており、収束と判定した場合、学習時処理を終了する。

【0051】

（量子状態の異常検知処理）
次に、図７を用いて、量子データ分布の学習処理の応用例として、量子状態の異常検知処理について説明する。ここで、新たに取得した量子データと、教師データセットから学習した分布とを比較することにより、異常かどうかの判断を行なう。

【0052】

まず、演算装置２０は、新たな量子データの取得を行なう（ステップＳ２１）。具体的には、管理部２１１は、ユーザ端末１０の指示に応じて、学習に用いたデータセットとは別の量子データ（新たな量子データ）を取得する。この新たな量子データとしては、観測により外部から入力されるデータや、量子計算により出力される量子状態等が考えられる。

【0053】

次に、演算装置２０は、ジェネレータを取得する（ステップＳ２２）。ここでは、管理部２１１は、通常時（正常時）の量子データセット（教師データセット）を用いて、予め学習されたジェネレータを取得する。

【0054】

次に、演算装置２０は、異常判定を行なう（ステップＳ２３）。ここでは、管理部２１１は、量子回路パラメータ｛θ}を用いて、異常判定を行なう。例えば、上記のジェネレータを用いて出力される量子状態の内、新たな量子データに最も近いデータを探索する。そして、管理部２１１は、新たな量子データと、最も近いデータとの距離(異常度)を算出する。そして、管理部２１１は、算出した異常度と、予め定めた異常判定値とを比較する。そして、管理部２１１は、比較結果に基づいて、異常の有無を判定し、判定結果をユーザ端末１０に出力する。

【0055】

本実施形態により、直感的には、以下のような効果を得ることができることが予想される。
（１）本実施形態では、潜在空間のデータを、データ空間ＳＰ１内の低次元の超平面ＭＦ１に写像する。これにより、写像を用いて、データセットと同様の分布に従う新たな量子データセットを生成することができる。

【0056】

（２）本実施形態では、局所的な測定から構成されるground costを用いることで、勾配消失を抑制できる。このため、最適輸送損失を用いた学習が可能となり、これを損失関数として量子データセットの分布を学習できる可能性がある。すなわち、最適輸送損失を損失関数として採用することで、暗黙的モデルを用いた量子生成モデルを構成できる。

【0057】

（３）本実施形態では、演算装置２０は、ground costの推定（ステップＳ１２）、最適輸送損失の算出（ステップＳ１３）、学習率及び勾配を利用した量子回路パラメータ｛θ｝の更新（ステップＳ１５）を行なう。これにより、最適輸送損失を最小化するような量子回路パラメータ｛θ｝を算出することができる。

【0058】

本実施形態は、以下のように変更して実施することができる。本実施形態及び以下の変更例は、技術的に矛盾しない範囲で互いに組み合わせて実施することができる。
・本実施形態では、量子状態の異常検知を行なう。すなわち、通常でない状態の検知に適用できる。例えば、量子コンピュータの異常検知、発生頻度が少ない状態の検知（例えば、重力波検知）、量子センシングに適用できる。

【0059】

また、量子コンピュータの計算結果として得られる量子データを用いた学習にも適用可能である。例えば、変分量子固有値ソルバー（ＶＱＥ：Variational Quantum Eigensolver）により得られた量子状態の学習、Born machineの波動関数の転移学習等に用いることができる。また、量子ダイナミクスの量子コンピュータでの計算として得られる量子状態を用いて、計算を簡易化するようなサロゲートモデルとしての利用も可能である。

【0060】

その他、例えば、古典データセットを埋め込んだ量子データセットにおいて、異常検知を行なうようにしてもよい。
・上記実施形態では、ground costを〔数７〕を用いて計算する。ground costを計算する計算式は、〔数７〕の数式に限定されるものではない。

【0061】

・上記実施形態では、ジェネレータG(θ,z)として、図５に示す公知の量子回路Ｃ１を用いる。ジェネレータG(θ,z)は、公知の量子回路Ｃ１に限定されるものではない。

【符号の説明】

【0062】

１０…ユーザ端末、２０…演算装置、２１１…管理部、２１２…古典計算部、２２…量子計算部、２２１…操作部、２２２…状態保持部、２２３…計測部。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】