[논문 리뷰] Perfect sampling algorithm for Schur processes

Dan Betea, Cédric Boutillier외 4명|arXiv (Cornell University)|2014. 07. 14.

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와 아츠체 다이아몬드에 대한 도미노 셔플링—을 하나의 알고리즘 프레임워크 아래 통합하는 것.
• 대칭 스체르 과정과 무한한 매개변수를 포함하는 일반화에 알고리즘을 확장하여 대칭 타일링과 오버파트션의 샘플링을 가능하게 하는 것.
• 대상 분포에 필요한 최소한의 난수 비트만 읽는 다항식 복잡도를 갖는 계산적으로 효율적인 방법을 제공하는 것.
• 이 알고리즘에 의해 가능해지는 대규모 시뮬레이션을 통해 타일링 모델에서 새로운 극한 형태 현상들을 밝혀내는 것.

제안 방법

• 알고리즘은 폴름의 RSK 대응에 대한 기술을 기반으로 한 성장 다이어그램을 사용하며, 상호작용하는 분할 수열을 포함하는 스체르 과정 설정에 적응시켰다.
• 샘플링 과정의 정확성과 효율성을 보장하기 위해 스체르 함수에 대한 코시 및 리틀우드 항등식의 이항 증명에 의존한다.
• 입력은 기하 또는 베르누이 난수와 난수 비트로 구성되며, 그림 형태의 다이어그램적 방식으로 기초 확률 측도를 표현한다.
• 알고리즘이 엔트로피 최적화되어 있으며, 목표 분포에서 정확하게 샘플을 생성하기 위해 필요한 최소한의 난수 비트만 사용한다.
• 대칭 스체르 과정의 경우, 리틀우드 항등식에 기반한 적응이 이루어지며, 대칭 RSK 알고리즘을 일반화하고 자유 경계 타일링의 샘플링을 가능하게 한다.
• 프랑케르의 측도와 혼합 매개변수 모델과 같은 극한 경우로의 확장이 이루어지며, 이 경우 알고리즘은 지수 시계를 갖는 연속 시간 역학으로 수렴한다.

실험 결과

연구 질문

• RQ1 스체르 과정의 전체 클래스—평면 분할과 아츠체 다이아몬드 타일링을 포함—에 대해 단일의 효율적이고 엔트로피 최적의 샘플링 알고리즘을 구성할 수 있는가?
• RQ2 RSK 알고리즘과 도미노 셔플링 알고리즘은 어떻게 동일한 조합적 프레임워크 아래 통합될 수 있는가?
• RQ3 스체르 과정 구성의 극한 형태 행동은 무엇이며, 대규모 샘플링을 통해 경험적으로 밝혀질 수 있는가?
• RQ4 이 알고리즘은 대칭 스체르 과정과 무한한 매개변수를 포함하는 기타 변형으로 일반화될 수 있는가?
• RQ5 알고리즘의 역학과 알려진 연속 시간 과정, 예를 들어 보로딘-페라리 역학에서의 지수 시계에 의해 지배되는 과정 사이의 관계는 무엇인가?

주요 결과

• 제안된 SchurSample 알고리즘은 다항식 복잡도와 최소한의 난수를 갖는 완전 샘플링을 달성하여 정확하고 엔트로피 최적화되어 있다.
• 이 알고리즘은 평면 분할에 대한 RSK 기반 샘플링과 아츠체 다이아몬드에 대한 도미노 셔플링 알고리즘을 일반화하여 두 주요 조합적 샘플링 절차를 통합한다.
• 이 알고리즘을 사용한 시뮬레이션은 기존에 관측되지 않았던 기울인 도미노 타일링과 피라미드 분할에서의 새로운 극한 형태 현상들을 드러낸다.
• 플랑케르 극한에서 알고리즘은 고전적 RSK 대응의 포issonized 형태로 축소되며, 폴름의 성장 다이어그램 기술을 복원한다.
• 혼합 매개변수 모델의 경우 알고리즘은 고전적 RSK와 RSK 유형의 역학 사이를 보간하며, 올컨넬의 포issonized RS 대응과 보로딘-페라리의 연속 시간 역학과 일치한다.
• 이 방법은 자연스럽게 대칭 스체르 과정으로 확장되며, 리틀우드 항등식 기반의 적응을 통해 대칭 평면 분할과 오버파트션의 샘플링을 가능하게 한다.