Abstract:Sudoku is a widely popular $\mathcal{NP}$-Complete combinatorial puzzle whose prospects for studying human computation have recently received attention, but the algorithmic hardness of Sudoku solving is yet largely unexplored. In this paper, we study the statistical mechanical properties of random Sudoku grids, showing that puzzles of varying sizes attain a hardness peak associated with a critical behavior in the constrainedness of random instances. In doing so, we provide the first description of a Sudoku \emph{freezing} transition, showing that the fraction of backbone variables undergoes a phase transition as the density of pre-filled cells is calibrated. We also uncover a variety of critical phenomena in the applicability of Sudoku elimination strategies, providing explanations as to why puzzles become boring outside the typical range of clue densities adopted by Sudoku publishers. We further show that the constrainedness of Sudoku puzzles can be understood in terms of the informational (Shannon) entropy of their solutions, which only increases up to the critical point where variables become frozen. Our findings shed light on the nature of the $k$-coloring transition when the graph topology is fixed, and are an invitation to the study of phase transition phenomena in problems defined over \emph{alldifferent} constraints. They also suggest advantages to studying the statistical mechanics of popular $\mathcal{NP}$-Hard puzzles, which can both aid the design of hard instances and help understand the difficulty of human problem solving.

Random Latin squares and Sudoku designs generation

Latin Puzzles

Fun with Latin Squares

Random sampling of Latin squares via binary contingency tables and probabilistic divide-and-conquer

Problem Solving at the Edge of Chaos: Entropy, Puzzles and the Sudoku Freezing Transition

On some randomized algorithms and their evaluation

A Simple Approach to Constructing Quasi-Sudoku-based Sliced Space-Filling Designs

Subsquares in random Latin squares and rectangles

Mutually orthogonal Sudoku Latin squares and their graphs

On an algorithm for receiving Sudoku matrices

A Novel Evolutionary Algorithm with Column and Sub-Block Local Search for Sudoku Puzzles

Extremal and monotone behaviour of the Sudoku number and related critical set parameters

Rating and generating Sudoku puzzles

Sudoku Number of Graphs

Investigating the Existence of Costas Latin Squares via Satisfiability Testing

Genarate and Solve Sudoku Puzzle

Generating Random Graphic Sequences

Random generation of optimal saturated designs

Sliced Latin hypercube designs with arbitrary run sizes

Designs based on the cycle structure of a Latin square autotopism

Towards a Latin-Square Search Engine.