[論文レビュー] Perfect sampling algorithm for Schur processes

Dan Betea, Cédric Boutillier+4 more|arXiv (Cornell University)|Jul 14, 2014

Random Matrices and Applications参考文献 23被引用数 15

ひとこと要約

本稿では、シュール関数のカウチおよびリトルウッド恒等式の統一的双対証明とフォーミンの成長図枠組みを用いて、平面分割、アズテック・ダイアモンドのタイル張り、パラミッド分割を含む組合せ的対象であるシュール過程の多項式時間の完全サンプリングアルゴリズムを提示する。このアルゴリズムはエントロピー最適であり、最小限のランダムビットを必要とし、RSKアルゴリズムとドミノシャッフルの両方を一般化し、効率的なサンプリングを可能にするとともに、タイリングモデルにおける新たな極限形状現象を明らかにする。

ABSTRACT

We describe random generation algorithms for a large class of random combinatorial objects called Schur processes, which are sequences of random (integer) partitions subject to certain interlacing conditions. This class contains several fundamental combinatorial objects as special cases, such as plane partitions, tilings of Aztec diamonds, pyramid partitions and more generally steep domino tilings of the plane. Our algorithm, which is of polynomial complexity, is both exact (i.e. the output follows exactly the target probability law, which is either Boltzmann or uniform in our case), and entropy optimal (i.e. it reads a minimal number of random bits as an input). The algorithm encompasses previous growth procedures for special Schur processes related to the primal and dual RSK algorithm, as well as the famous domino shuffling algorithm for domino tilings of the Aztec diamond. It can be easily adapted to deal with symmetric Schur processes and general Schur processes involving infinitely many parameters. It is more concrete and easier to implement than Borodin's algorithm, and it is entropy optimal. At a technical level, it relies on unified bijective proofs of the different types of Cauchy and Littlewood identities for Schur functions, and on an adaptation of Fomin's growth diagram description of the RSK algorithm to that setting. Simulations performed with this algorithm suggest interesting limit shape phenomena for the corresponding tiling models, some of which are new.

研究の動機と目的

• 平面分割やタイリングを含む、タイリングや整数分割といった広範なランダムな組合せ的対象のクラスであるシュール過程の一般的で正確かつエントロピー最適なサンプリングアルゴリズムの開発。
• 平面分割のRSKやアズテック・ダイアモンドのドミノシャッフルといった既存のサンプリング手順を、単一のアルゴリズム的枠組みに統合すること。
• 無限個のパラメータを含む一般化を伴う対称シュール過程へのアルゴリズムの拡張により、対称タイリングやオーバーパartitionsのサンプリングを可能にすること。
• ターゲット分布に必要な最小限のランダムビットのみを読み取る多項式時間の計算効率の良い方法を提供すること。
• アルゴリズムによって可能になる大規模シミュレーションを通じて、タイリングモデルにおける新たな極限形状現象の同定。

提案手法

• フォーミンによるRSK対応の記述に基づく成長図を用い、入れ子の分割列を伴うシュール過程の設定に適応したアルゴリズム。
• シュール関数のカウチおよびリトルウッド恒等式の双対証明に依拠し、サンプリングプロセスの正しさと効率性を保証。
• 入力は幾何的またはベルヌーイ的確率変数とランダムビットから構成され、図的・図解的形で確率測度を表現。
• アルゴリズムはエントロピー最適であり、ターゲット分布から正確にサンプルを生成するために必要な最小限のランダムビットのみを用いる。
• 対称シュール過程の場合、リトルウッド恒等式に適応され、対称RSKアルゴリズムを一般化し、自由境界タイリングのサンプリングを可能にする。
• プラナチェル測度や混合パrameterモデルなどの極限ケースに拡張され、指数分布の時計を伴う連続時間ダイナミクスに収束する。

実験結果

リサーチクエスチョン

• RQ1 平面分割やアズテック・ダイアモンドのタイル張りを含む、シュール過程の全クラスに対して、単一で効率的かつエントロピー最適なサンプリングアルゴリズムを構築可能か？
• RQ2 RSKアルゴリズムとドミノシャッフルアルゴリズムは、どのように同一の組合せ的枠組みに統合できるか？
• RQ3 シュール過程の配置における極限形状の挙動は何か？ そして、大規模なサンプリングによってそれらは実証的に発見可能か？
• RQ4 アルゴリズムは、無限個のパラメータを含む対称シュール過程やその他の変種に一般化可能か？
• RQ5 アルゴリズムのダイナミクスと、ボロディン＝フェラリのダイナミクスに従う指数時計を伴う既知の連続時間プロセスとの関係は何か？

主な発見

• 提案されたSchurSampleアルゴリズムは、多項式時間の複雑性と最小限のランダム性を達成し、正確かつエントロピー最適である。
• このアルゴリズムは、平面分割のRSKに基づくサンプリングとアズテック・ダイアモンドのドミノシャッフルの両方を一般化し、2つの主要な組合せ的サンプリング手順を統合する。
• アルゴリズムを用いたシミュレーションにより、急勾配ドミノタイリングやパラミッド分割における、以前に観察されていなかった挙動を含む、タイリングモデルにおける新たな極限形状現象が明らかになった。
• プラナチェル極限において、アルゴリズムは古典的RSK対応のポisson化版に還元され、フォーミンの成長図記述を回復する。
• 混合パrameterモデルでは、古典的RSKとRSK型ダイナミクスの間を補間し、オー・コンネルのポisson化RS双対写像およびボロディン＝フェラリの連続時間ダイナミクスと整合する。
• この手法は自然に対称シュール過程へ拡張され、リトルウッド恒等式に基づく適応により、対称平面分割やオーバーパーティションのサンプリングが可能になる。