[論文レビュー] Algebraic Unimodular Counting

Jesús A. De Loera, Bernd Sturmfels|ArXiv.org|Apr 30, 2001

Commutative Algebra and Its Applications参考文献 11被引用数 8

ひとこと要約

本稿では、ユニモジュラー線形方程式系の非負整数解の個数を効率的に計算する2つの代数的アルゴリズム——グレブナー基底とBBKLP手法——を提示する。これにより、従来の計算限界を超えて連続表やコスタントの分割関数の正確な個数計算が可能になる。主な貢献は、$5\times5$の連続表や$A_6$根系に対応するユニモジュラー行列のためのチャネル多項式の計算が可能になったことであり、明示的な公式とオンライン計算機が提供されている。

ABSTRACT

We study algebraic algorithms for expressing the number of non-negative integer solutions to a unimodular system of linear equations as a function of the right hand side. Our methods include Todd classes of toric varieties via Gröbner bases, and rational generating functions as in Barvinok's algorithm. We report polyhedral and computational results for two special cases: counting contingency tables and Kostant's partition function.

研究の動機と目的

• ユニモジュラー行列に対する、$Ax = b$ の非負整数解の個数であるベクトル分割関数 $\phi_A(b)$ を計算する実用的な代数的アルゴリズムの開発。
• マウントの手法のような従来の補間ベース手法が、ハードウェア制限や行列固有のヒューリスティクスに制約を受けるのを克服すること。
• 特に連続表およびコスタントの分割関数に対して、ユニモジュラー系における $\phi_A(b)$ の完全なチャネル分解と多項式表現の計算。
• キリロフの $A_{m-1}$ 根系に対するチャネル数に関する未解決問題（$m=7$ まで）を解き、コスタント関数のオンライン計算機を提供すること。

提案手法

• 最初の手法は、多面体 $\{x \geq 0 : Ax = b\}$ に関連するトーリック多様体のTodd類をグレブナー基底を用いて計算し、代数幾何を活用して分割関数を抽出する。
• 2番目の手法は、バーヴィノクの有理生成関数とブリオンの定理に基づくBBKLPアルゴリズムを適用し、多面体内の格子点の生成関数を計算することで、$\phi_A(b)$ の効率的評価を可能にする。
• アルゴリズムはMacaulay2（グレブナー基底）およびMaple（BBKLP）で実装されており、コードは公開利用可能である。
• チャネル複体は、$\text{cone}(A)$ を、$\phi_A(b)$ が多項式になる凸多面体錐に分解することで計算され、チャネルと三角形分割の双対性を利用している。
• $K_n$（完全有向グラフの接続行列）に対しては、バランスの取れたカットと根ベクトルの同一直線上性を用いてチャネル複体を分析し、ファセット数が指数関数的に増加することが判明した。
• ユニモジュラー性を活用することで整数頂点と多項式的挙動を保証し、$b$ のサイズに依存しない正確な計算が可能になる。

実験結果

リサーチクエスチョン

• RQ1 グレブナー基底や有理生成関数に基づく代数的手法が、マウントの補間ベース手法よりもユニモジュラー数え上げにおいて優れているかどうか。
• RQ2 コスタントの分割関数の $A_{m-1}$ 根系に対するチャネル数とチャネル複体の構造は何か。
• RQ3 $4\times5$ および $5\times5$ の連続表に対して、任意のマージンを持つ非負整数解はいくつ存在し、これらを正確に計算できるか。
• RQ4 $K_n$ のチャネル複体におけるチャネル数の漸近的成長は何か。指数関数的挙動を示すか。
• RQ5 $K_n$ のチャネル多項式を明示的に計算できるか。また、それらを用いてコスタントの分割関数のオンライン計算機を構築できるか。

主な発見

• 本稿では、$A_6$ 根系に対するコスタントの分割関数のすべてのチャネル多項式を計算し、$m \leq 6$ におけるキリロフの未解決問題を解決した。$m \leq 7$ まででチャネル数を特定し、$m=7$ の場合に44,288個のチャネルがあることが判明した。
• チャネル複体 $\Gamma(K_n)$ は少なくとも $2^{\lfloor n/2 \rfloor}$ 個のファセットを持つことが示され、$n$ の増加に伴い指数関数的に複雑性が増加することが、バランスの取れたカットからの再帰的二重化によって確認された。
• $4\times4$ 連続表のケースでは、既知の7つのチャネルという結果を再現し、マウントの計算と整合的であるが、$5\times5$ 表へと拡張可能である。
• すべてのチャネルにおいて $\phi_A(b)$ の明示的な区分多項式公式が得られ、たとえば $\min\{b_1,b_2,b_3\} \geq 0$ のチャネルでは $\phi_{K_4}(b) = (b_1+1)(b_1+2)(b_1+3b_2+3)/6$ となり、次数3の多項式である。
• www.math.ucdavis.edu/~deloera/kostant.html に設置されたオンライン計算機は、計算済みの多項式を用いてコスタントの分割関数をリアルタイムで評価可能である。
• バーヴィノクのアルゴリズムに基づくBBKLP手法は、ユニモジュラー行列に対して非常に効率的であり、補間ベース手法を上回り、マウントの手法の実現可能性に関する疑問に明確に答えている。