[論文レビュー] Random Latin squares and Sudoku designs generation

Roberto Fontana|arXiv (Cornell University)|May 16, 2013

graph theory and CDMA systems参考文献 6被引用数 4

ひとこと要約

本稿では、特別に構築されたグラフにおける最大クリークとしてモデル化することで、ラテン方陣およびスウドクの均一なランダム生成のための新規アルゴリズムを提示する。Cliquer などのグラフ解析ツールを用いて、クリーク列挙と置換に基づく変換を活用することで、すべての可能な構成から均一にサンプリングする。n=7 までに成功裏に実装され、スケーラブルなグラフ解析を用いることで、より大きな n に対しても実行可能である。

ABSTRACT

Uniform random generation of Latin squares is a classical problem. In this paper we prove that both Latin squares and Sudoku designs are maximum cliques of properly defined graphs. We have developed a simple algorithm for uniform random sampling of Latin squares and Sudoku designs. It makes use of recent tools for graph analysis. The corresponding SAS code is annexed.

研究の動機と目的

• 既存のマルコフ連鎖ベースのアプローチに起因する制限を克服するため、ラテン方陣およびスウドクの設計に対する均一なランダムサンプリング手法の開発。
• 定義されたグラフ構造におけるラテン方陣と最大クリークの理論的同等性を確立し、新たなアルゴリズム的手法を可能にする。
• SAS を用いた実装を含め、正確な均一サンプリングを n=7 まで効率的かつ実用的に実現するため、グラフ解析ツール（例：Cliquer）を用いた実装。
• スウドク設計へのフレームワークの拡張を図る。これは、ブロック制約を満たすようにグラフ構造を変更することで実現する。
• サブグラフサンプリングと現代のグラフ解析ソフトウェアを用いて、より高い順序（例：n=9）におけるスケーラビリティを検討する。

提案手法

• ノードが固定点を持たない置換（derangements）を表すグラフにおいて、ラテン方陣を最大クリークとしてモデル化する。エッジは、すべての位置で互いに素である置換同士を結ぶ。
• {1,…,n} の置換としての derangements をノードとするグラフ 𝒢ₙ を構築し、すべての位置でペアワイズに互いに素である置換の間にエッジを設ける。
• スウドク設計のため、3×3 のサブグリッド構造を保証するように変更されたグラフ 𝔾̃ₙ を定義する。9×9 の場合に有効である。
• グラフ解析ソフトウェア（例：Cliquer）を用いて、𝒢ₙ や 𝔾̃ₙ のサイズ n のすべての最大クリークを列挙し、それらが有効なラテン方陣またはスウドク設計に対応することを確認する。
• ランダムに選択されたクリークを用いて、均一なランダムラテン方陣を生成する。これには、記号 2 から n に対するランダムな置換を適用し、列を再配置することで完全な均一性を達成する。
• スウドクの場合は、ブロック構造を保つために、列および行の置換をバンド／スタックの置換に制限する。

実験結果

リサーチクエスチョン

• RQ1 適切に定義されたグラフにおいて、ラテン方陣およびスウドク設計が最大クリークとして正確に同定可能であるか？
• RQ2 n=7 までの順序について、このようなグラフにおけるすべての最大クリークの列挙が計算的に実行可能であるか？
• RQ3 Jacobson と Matthews（1996）のマルコフ連鎖モンテカルロ手法と比較して、提案手法の効率性および均一性はどの程度優れているか？
• RQ4 n=9 のようなより高い順序において、サブグラフサンプリングを用いて均一サンプリングを近似できる範囲はどの程度か？
• RQ5 現代のグラフ解析ツールは、大規模なラテン方陣の均一サンプリングに必要なクリーク列挙の規模を処理できるか？

主な発見

• 本稿では、ラテン方陣およびスウドク設計が特別に構築されたグラフにおける最大クリークに正確に対応することを証明し、それらの生成に新たな理論的基盤を提供した。
• n=7 の場合、導出された置換グラフにおける最大クリーク 16,942,080 個を列挙することで、61,479,419,904,000 個の有効なラテン方陣を生成した。
• SAS で実装された本手法は、全 17,972 個のスウドク置換のうち 809 個をランダムに選択したサンプルから、有効な 9×9 スウドク設計を成功裏に生成した。
• 9×9 スウドク部分グラフにおける最大クリークのサイズは 8 であり、そのような最大クリークが 73 個同定された。これにより構造的妥当性が確認された。
• 本手法は原則としてスケーラブルである。最近のグラフ解析ソフトウェアの進歩により、より大きな n に対しても均一サンプリングが可能であるが、計算コストは n に伴い増加する。
• クリーク列挙とその後の記号および列のランダム置換を活用することで、本アルゴリズムは均一なサンプリングを達成し、最終出力にバイアスを生じさせない。