[論文レビュー] Answers to Some Questions about Explicit Sinkhorn Limits posed by Mel Nathanson

Shalosh B. Ekhad, Doron Zeilberger|arXiv (Cornell University)|Feb 27, 2019

Computability, Logic, AI Algorithms被引用数 4

ひとこと要約

この論文は、正の要素を持つ3×3対称行列のSinkhorn極限に対して、グレブナー基底を用いて対角スケーリング行列を解くことにより、明示的な代数的公式を提供する。双随机行列の極限行列およびそのスケーリング因子は、元の行列の要素を変数とする4次方程式の根として閉形式で得られる。これは、メル・ナタンソンが提起した、明示的なSinkhorn極限とアルゴリズムの有限回収束に関する未解決問題を解消する。

ABSTRACT

At the Jan. 2018 Joint Mathematics Meetings, Avi Wigderson gave a series of three fascinating lectures whose starting point was the Sinkhorn algorithm. One of the people in the audience was Mel Nathanson, and this lead him to ask some intriguing questions, some of which we answer in this note.

研究の動機と目的

• メル・ナタンソンが要請した、正の要素を持つ3×3対称行列のSinkhorn極限に対する明示的な代数的表現を提供すること。
• ナタンソンの問い、すなわちSinkhornアルゴリズムが有限回のステップで終了できるかどうかを、そのような例を明示的に構成することで解明すること。
• グレブナー基底技術を用いて、Sinkhorn極限およびスケーリング行列の正確な代数的表現を導出すること。
• パラメータを含む特別な3×3対称行列のケースに一般化し、そのSinkhorn極限に対して閉形式の解を得ること。

提案手法

• 著者らは、ブッヒバーガーのアルゴリズムを実装したMapleパッケージSINKHORN.txtを用い、Sinkhorn極限を定義する方程式系を記号的に解く。
• 一般の3×3対称行列に対して、変数$ z $に関する4次方程式が導出され、その正の根が対角行列$ X $のスケーリング係数を決定する。ここで$ s_{11} = a_{11}z $である。
• Sinkhornアルゴリズムの交互な行と列の正規化操作を記号的に適用し、双随机性の条件を課えることで多項式制約を導出する。
• 特定の例を構築することで、Sinkhornアルゴリズムが2ステップで収束することを確認し、有限収束を検証する。
• パラメータ付き要素（例：$ K, L, M $）を含む特別なケースに対しては、根がSinkhorn極限およびスケーリング行列を導く低次元化された4次方程式が得られる。
• この手法は、より大きな$ n \times n $対称行列へ一般化可能であるが、定義方程式の次数は$ 2^{n-1} $に比例するため、$ n \geq 4 $では記号計算が非現実的になる。

実験結果

リサーチクエスチョン

• RQ1 正の要素を持つ3×3対称行列のSinkhorn極限は、行列要素の代数的式として閉形式で表現可能か？
• RQ2 Sinkhornアルゴリズムが有限回のステップで終了する行列は存在するか？ そのような行列は明示的に構成可能か？
• RQ3 グレブナー基底技術を用いて、Sinkhorn極限およびスケーリング行列の正確な代数的表現を導出可能か？
• RQ4 スケーリング係数の定義多項式の構造は何か？ その次数は行列サイズにどのように依存するか？
• RQ5 要素が$ \frac{r(r+1)}{2}, 1, 1 $である行列$ A(r) $に対して、一般の$ r $に対してSinkhornアルゴリズムは有限回で終了するか？

主な発見

• 一般の3×3対称行列$ A $のSinkhorn極限$ S $は、$ s_{11} = a_{11}z $として明示的に与えられ、ここで$ z $は$ A $の要素を係数とする4次方程式の正の根であり、その係数は次数8の同次多項式である。
• 対角スケーリング行列$ X $は$ x_{11} = \sqrt{z} $を満たし、$ X $および$ S $の残りの要素は公開URLで提供可能であるが、論文に記載するには長すぎる。
• 行列$ A = \begin{pmatrix} K & 1 & 1 \\ 1 & L & 1 \\ 1 & 1 & 1 \end{pmatrix} $に対して、Sinkhorn極限は$ z $に関する4次方程式で決定され、係数は$ K $および$ L $に依存する。特定の$ K, L $の選択では、2ステップでアルゴリズムが終了する。
• 行列$ A = \begin{pmatrix} K & 1 & 1 \\ 1 & L & 1 \\ 1 & 1 & M \end{pmatrix} $に対しても、Sinkhorn極限は再び$ z $に関する4次方程式で与えられ、係数は$ K, L, M $を含み、解は明示的に計算可能である。
• 著者らは、2回の反復でSinkhorn極限に到達する行列$ M = \begin{pmatrix} 2 & 2 & 6 \\ 2 & 1 & 2 \\ 9 & 3 & 3 \end{pmatrix} $を構築し、非双随机行列に対しても有限収束が成立することを確認した。
• 行列$ A(r) $で$ a_{11} = \frac{r(r+1)}{2} $のとき、Sinkhorn極限は$ S(r) $として明示的に計算可能であり、$ r = 1 $または$ r = -2 $でない限り、記号的反復において非ゼロの行和差が持続するため、有限回収束は成立しない。