[论文解读] Random Permutations, Random Sudoku Matrices and Randomized Algorithms

Krasimir Yordzhev|arXiv (Cornell University)|Dec 1, 2013

graph theory and CDMA systems被引用 8

一句话总结

本文提出了一种新颖的随机化算法，通过利用Πₙ矩阵的结构（即(2n)×n矩阵，其行是ℤₙ的排列），来生成随机的n²×n²数独矩阵。该算法引入了一个双射映射φ: Πₙ → Σₙ²，将数独约束表示为二值矩阵，从而实现高效的不相交性检查。其主要贡献是一项O(n⁴)时间复杂度的算法，通过逐步组合不相交的Πₙ矩阵来构建有效的数独矩阵，显著优于朴素的蒙特卡洛采样方法。

ABSTRACT

Some randomized algorithms, used to obtain a random $n^2 imes n^2$ Sudoku matrix, where $n$ is a natural number, is reviewed in this study. Below is described the set $Π_n$ of all $(2n) imes n$ matrices, consisting of elements of the set $\mathbb{Z}_n =\{ 1,2,\ldots ,n\}$, such that every row is a permutation. It is proved that such matrices would be particularly useful in developing efficient algorithms in generating Sudoku matrices. An algorithm to obtain random $Π_n$ matrices is presented in this paper. The algorithms are evaluated according to two criteria - probability evaluation, and time evaluation. This type of criteria is interesting from both theoretical and practical point of view because they are particularly useful in the analysis of computer programs.

研究动机与目标

• 开发一种高效的随机化算法，用于生成随机的n²×n²数独矩阵。
• 形式化集合Πₙ，即(2n)×n矩阵的集合，其行是ℤₙ的排列，作为数独矩阵的构建模块。
• 建立双射映射φ: Πₙ → Σₙ²，将Πₙ矩阵转换为编码数独区块约束的二值矩阵。
• 通过概率和时间复杂度标准评估该算法，证明其在效率上优于朴素的蒙特卡洛采样方法。

提出的方法

• 将Πₙ定义为(2n)×n矩阵的集合，其中每一行是ℤₙ的排列，作为数独构造的基础。
• 引入算法3.1，通过为每一行置换ℤₙ的元素来生成Πₙ中的随机矩阵。
• 开发算法3.3，将每个Pₖ ∈ Πₙ映射到一个二值矩阵A = φ(Pₖ) ∈ Σₙ²，其中A通过单位矩阵编码区块位置。
• 利用双射映射φ验证不相交性：当且仅当它们对应的Σₙ²矩阵无重叠的1时，两个Πₙ矩阵是不相交的。
• 通过仅在总和B = ∑Aₖ保持为二值矩阵（即元素为0或1）时，将k·Aₖ累加到最终的数独矩阵S中，逐步构建S。
• 应用算法4.1，该算法使用while循环不断重采样Pₖ，直到找到一个不相交的Πₙ矩阵，从而确保满足有效的数独约束。

实验结果

研究问题

• RQ1 Πₙ矩阵能否作为生成随机数独矩阵的结构化基础？
• RQ2 是否存在Πₙ矩阵与二值矩阵Σₙ²之间的双射映射，且能保持数独约束？
• RQ3 能否通过其对应的Σₙ²表示，高效验证Πₙ矩阵的不相交性？
• RQ4 所提出的算法在生成随机数独矩阵时，时间复杂度是否优于朴素的蒙特卡洛采样？
• RQ5 使用该方法构建数独矩阵时，每次迭代的成功概率是多少？

主要发现

• 集合Πₙ包含n!²ⁿ个元素，且每个P ∈ Πₙ唯一对应一个二值矩阵A = φ(P) ∈ Σₙ²，从而建立了双射关系。
• 当且仅当它们对应的Σₙ²矩阵无重叠的1时，两个Πₙ矩阵是不相交的，这使得通过矩阵加法实现高效的验证成为可能。
• 所提出的算法4.1的时间复杂度为O(n⁴)，由于优化了约束检查，显著优于朴素的O(n⁶)蒙特卡洛方法。
• 每次迭代的成功概率为p(n) = σₙ / n²ⁿ²，其中σₙ是有效n²×n²数独矩阵的数量，σ₂ = 288，σ₃ ≈ 6.671 × 10²¹。
• 该算法通过构造保证正确性：仅满足行、列和区块排列约束的矩阵才会被添加到最终的数独矩阵中。
• 该方法通过在二值矩阵中编码约束，避免了对完整数独约束的暴力验证，将每次验证的时间从O(n³)降低至O(n)。