線形代数において、行列のパーマネントの計算は、定義が似ているように見えるにもかかわらず、行列式の計算よりも難しい問題であると考えられている。
パーマネントは、行列式と同様に、異なる行と列にある行列要素の集合の積の和として定義されます。ただし、行列式では集合の偶奇性に基づいて各積に±1の符号が付けられるのに対し、パーマネントではすべての積に+1の符号が付けられます。
行列式はガウス消去法によって多項式時間で計算できますが、パーマネントは多項式時間では計算できないと一般的に考えられています。計算複雑性理論において、ヴァリアントの定理は、パーマネントの計算は#P困難であり、すべての要素が0または1である行列の場合は#P完全であると述べています(ヴァリアント、1979) 。このことから、パーマネントの計算はNPよりもさらに計算が難しいと考えられている問題のクラスに分類されます。対数空間一様ACC 0回路ではパーマネントの計算は不可能であることが知られています(アレンダー&ゴア、1994 )。
行列のパーマネントを計算するための厳密なアルゴリズムと近似アルゴリズムの開発は、活発な研究分野である。
n × n行列A = ( a i,j )のパーマネントは次のように定義される。
ここでの総和は、対称群S nのすべての要素 σ 、すなわち1, 2, ..., nのすべての順列にわたっています。この式は、行列式の対応する式と、行列式では各積に順列σ の符号が掛けられるのに対し、この式では各積に符号がないという点だけが異なります。この式は、すべての順列にわたって総和を取り、その総和の中で各行列要素を掛け合わせるという、単純に式を展開するアルゴリズムに直接変換できます。これにはn! n回の算術演算が必要です。
最もよく知られている[ 1 ]一般的な厳密アルゴリズムは、 HJ Ryser (1963 )によるものである。Ryserの方法は、次のように与えられる包含排除公式に基づいている[ 2 ] 。 Aからk列を削除して得られるものとする。行和の積である、そしての値の合計全体的に可能な限り。 それから
行列の要素を用いて次のように書き直すことができる[ 3 ]
Ryser の式と同程度(あるいは2倍速い)と思われる別の式は、2 つの博士論文に見られます。( Balasubramanian 1980 )、( Bax 1998 )、 ( Bax & Franklin 1996 )を参照してください。この式を見つける方法はかなり異なり、それぞれ Muir 代数の組み合わせ論と有限差分理論に関連しています。不変理論に関連する別の方法は、対称テンソルの偏極恒等式を経由することです( Glynn 2010 )。この式は、これらの著者全員によって発見されたように、無限に多くの他の式に一般化されますが、それらが基本式よりも速いかどうかは明らかではありません。( Glynn 2013 )を参照してください。
この種の最も単純な既知の公式(場の特性が2でない場合)は
外側の合計は全体にわたるベクトル。
二部グラフにおける完全マッチングの数は、グラフの二部隣接行列のパーマネントによってカウントされ、任意の 0-1 行列のパーマネントは、このようにしてグラフにおける完全マッチングの数として解釈できます。平面グラフ(二部グラフであるかどうかに関わらず) の場合、FKT アルゴリズムは、グラフのTutte 行列のエントリの慎重に選択された部分集合の符号を変更することによって、多項式時間で完全マッチングの数を計算します。これにより、結果として得られる歪対称行列のPfaffian (その行列式の平方根) が完全マッチングの数になります。この手法は、完全な二部グラフK 3,3と同相な部分グラフを含まないグラフに一般化できます。[ 5 ]
ジョージ・ポリアは、 01行列Aの一部の要素の符号を変更することで、新しい行列の行列式がAのパーマネントとなるのはどのような場合かという疑問を提起した[ 6 ]。すべての01行列がこのように「変換可能」なわけではない。実際、線形写像は存在しないことが知られている(マーカス&ミンク(1961) )。そのためすべての人々のために行列「変換可能」な行列の特徴付けは、Little (1975)によって与えられ、そのような行列は、 Pfaffian 方向を持つ二部グラフの二部隣接行列であることを示しました。Pfaffian 方向とは、すべての偶数サイクルに対して、エッジの方向が となるような方向です。そのために完全マッチングがあり、C に沿って方向付けられたエッジの数が奇数である(したがって、反対方向のエッジの数も奇数である)。また、これらのグラフは、 と同相な部分グラフを含まないグラフと正確に一致することも示された。上記のとおり。
モジュロ2では、パーマネントは行列式と同じである。また、モジュロで計算することもできます。時間が経つにつれてのためにしかし、2のべき乗でない任意の数を法とするパーマネントを計算することは非常に困難である。ヴァリアント(1979)
グリン(2010)は、素数pを法とする計算のためのさまざまな公式を提示している。まず、偏微分を用いた記号計算を用いる公式がある。
第二に、p = 3 の場合、n×n 行列に対して次の式が成り立つ。行列の主小行列式を含む(Kogan(1996)):
どこは、行と列によって誘発される インデックス、 そしてはで一方、空の小行列の行列式は1と定義される。
上記の展開は、任意の特性pに関して、次の双対恒等式のペアとして一般化できる。 どちらの式においても、和はすべての ( p − 1) 組について取られる。集合の分割であるp − 1 個の部分集合に分割する。これらの部分集合の中には、空集合が含まれる場合もある。
前者の式は、対称的なハフニアンに対応する類似式を持つ。そして奇妙なp:
同じインデックスの集合について和をとる。さらに、特性がゼロの場合、パーマネントと行列式の両方を含む同様の畳み込み和表現は、ハミルトン閉路多項式(次のように定義される)を与える。どこは、サイクルを1つだけ持つn順列の集合です。
特性2では、後者の等式はしたがって、これは任意のユニタリのハミルトンサイクル多項式を多項式時間で計算する機会を提供する(つまり、どこ(n × n単位行列)は、そのような行列の各小行列式がその代数的補行列と一致するため、単位行列である。どここれは、 n × n単位行列のインデックス 1,1 の要素を 0 に置き換えたものです。さらに、これはn × n単位行列に一般化することもできます。としてどこは {1, ..., n }の部分集合です。は、すべてのkに属するインデックスk、kの要素を 0 に置き換えたn × n単位行列です。、そして定義しますどこは、各サイクルに少なくとも 1 つの要素が含まれる n 置換の集合である。。
この式は、標数3の体において以下の恒等式も意味する。
任意の可逆
任意のユニタリつまり、正方行列そのためどこは対応するサイズの単位行列です。
どこは、対応する要素の立方である要素を持つ行列です。。
また、( Kogan (1996) ) 正方行列を定義すると、k-半ユニタリとして1-半ユニタリ行列のパーマネントは、標数3の体上では多項式時間で計算可能ですが、k > 1 の場合は#3-P-完全問題となります。(標数2のハミルトン閉路多項式に関する並行理論があります。ユニタリ行列上での計算は多項式時間で実行可能ですが、k-半ユニタリ行列の場合は任意のk > 0に対して #2-P-完全問題となります。)後者の結果は2017年に本質的に拡張され(Knezevic & Cohen (2017))、標数3では正方行列のパーマネントとその部分逆行列(そして四角いので、可逆であること):
そして、k行またはk -1行のサブセットが別の(互いに素な)k行のサブセットの線形結合として表現できるn × n行列のパーマネントの計算を、それぞれ( n - k )×( n - k)行列または(n- k +1)× ( n - k +1)行列のパーマネントの計算に多項式時間で還元することを可能にする。これにより、特性3のパーマネントを「保存」する圧縮演算子(行列式の計算に適用されるガウス修正に類似)が導入される。(類似して、特性2のハミルトン閉路多項式も不変行列圧縮を持つことに注目すべきである。これは、3つの等しい行を持つ任意のn × n行列A、またはn >2の場合、i番目とj番目の行が同一で、 i番目とj番目の列が同一であるインデックスのペアi、jに対してham( A ) = 0となるという事実を考慮すると、である。 (同様に。)その演算子の閉包は、転置変換(演算子が行列をそのまま残すたびに使用される)とともに逐次適用された極限として定義され、行列のクラスに適用された場合、あるクラスから別のクラスへの演算子マッピングでもあります。圧縮演算子は 1-半ユニタリ行列のクラスをそれ自身とユニタリ行列および 2-半ユニタリ行列のクラスに写像しますが、1-半ユニタリ行列のクラス (およびユニタリ行列から任意の行ベクトルで 1 つの行を置き換えることによって得られる行列のクラス- このような行列のパーマネントは、ラプラス展開を介して 1-半ユニタリ行列のパーマネントの和であり、したがって多項式時間で計算可能です) の圧縮閉包はまだ不明であり、標数 3 におけるパーマネントの計算複雑性の一般的な問題とP 対 NPの主要な問題に密接に関連しています。 ( Knezevic & Cohen (2017)で示されたように)、このような圧縮閉包が標数 3 の体上のすべての正方行列の集合であるか、少なくともパーマネントの計算が#3-P 完全である行列クラス(2-半ユニタリ行列のクラスなど) を含む場合、パーマネントはこの特性に関して多項式時間で計算可能である。
さらに、特性3に存在する永続性を保つ圧縮の類似物を他の素特性について見つけて分類するという問題が定式化され(Knezevic & Cohen (2017) )、 n × n行列に対して次の恒等式が与えられた。そして、2つのn次元ベクトル(すべての要素が集合{0, ..., p − 1}から構成される)そしてそのため任意の素数特性pにおいて有効:
ここで、n × m行列の場合n次元ベクトルそしてmベクトル両方のベクトルは、すべての要素が集合 {0, ..., p − 1} から得られます。は、から受け取った行列を表します。繰り返しによってi = 1, ..., nのi番目の行を掛け、j = 1, ..., mのj番目の列の倍数(行または列の重複度がゼロの場合、その行または列が削除されたことを意味するため、この概念は部分行列の概念の一般化である)、は、すべての要素が 1 に等しい n 次元ベクトルを表します。この恒等式は、行列の小行列式をその逆行列の小行列式で表す古典的な公式と全く同じであり、したがって、行列式とパーマネントが相対的イマナントとして一種の双対性を持っていることを(再び)示しています。(実際には、対称行列のハフニアンに対するそれ自身の類似物です。)そして奇素数pは)
さらに、素数標数 p の部分逆の場合のより広い一般化として、、四角いので、反転可能でサイズもx、 そしてそこには同一性も存在する
ここで、共通の行/列多重度ベクトルそして行列の場合対応する行/列の多重度ベクトルを生成するそして、s,t = 1,2、そのブロックの場合(同じ懸念事項)(等式の右辺にある部分逆数)。
Aの要素が非負の場合、パーマネントは、パーマネントの値 M と任意の ε > 0 の誤差を除いて、確率的多項式時間で近似的に計算できます。言い換えれば、完全多項式時間ランダム化近似スキーム (FPRAS) が存在します( Jerrum 、 Sinclair 、Vigoda ( 2001 ) )。
計算において最も難しいステップは、与えられた二部グラフ内のすべての完全マッチングの集合からほぼ均一にサンプリングするアルゴリズム、すなわち完全多項式時間でほぼ均一にサンプリングするアルゴリズム(FPAUS)を構築することです。これは、分布がほぼ均一で混合時間が多項式時間となるマルコフ連鎖を定義および実行するためにメトロポリスルールを使用するマルコフ連鎖モンテカルロアルゴリズムを用いることで実現できます。
パーマネントの自己還元性を利用して、グラフ内の完全マッチングの数を近似的に数えることが可能です。これは、FPAUSと、Jerrum、Valiant 、 Vazirani (1986)によるサンプリングから計数へのよく知られた還元法を組み合わせることによって実現できます。完全一致の数を表すおおよそ、特定のエッジの場合で多数のマッチングをサンプリングすることによってそして、そのうちのいくつが一致しているかを数える比率の推定値を得ることができる数すると、 どこ同じ方法を再帰的に適用することで近似できる。
パーマネントが特に興味深いもう1つの行列クラスは、正半定値行列である。[ 7 ]ストックマイヤー計数法を用いることで、これらの行列をクラス内で計算することができる。しかし、これは一般的に実現不可能なクラスと考えられています。PSD行列のパーマネントを準指数因子で近似することはNP困難であり、-hard [ 8 ]スペクトルにさらに制約が課される場合、より効率的なアルゴリズムが知られています。 1 つのランダム化アルゴリズムはボソンサンプリングのモデルに基づいており、量子光学に固有のツールを使用して、正定値半正定値行列のパーマネントを特定の確率変数の期待値として表現します。 後者は、その標本平均によって近似されます。[ 9 ]このアルゴリズムは、特定の正定値半正定値行列のセットに対して、そのパーマネントを加法誤差を除いて多項式時間で近似します。これは、Gurvits による標準的な古典的な多項式時間アルゴリズムよりも信頼性が高いです。[ 10 ]
{{citation}}ISBN /日付の不一致(ヘルプ)