What happens when Sudoku scales?

What happens when Sudoku scales?

3 min read
Part of the papers collection. Introductions to papers I have worked on, with some context for the research question and links to the paper and supporting material. See all 9 posts in the collection

There is now a public repository of 434,201 generated Sudoku instances, ranging from 6x6 to 36x36. That collection is the main result of my paper Scaling Sudoku as a Constraint Problem, to be presented at ModRef 2026. The paper describes how the collection was generated, and its experiments provide a map for navigating it.

The change in scale is easier to see than describe. The widget below shows one published base instance at each size; the board stays the same width while the number of cells grows.

Generated instances at five sizes

Grid size

9×9 grid · 3×3 boxes

24 clues · 29.6% filled

View this instance in the repository

Building the collection#

In his 2005 paper Sudoku as a Constraint Problem, Helmut Simonis used published 9x9 Sudoku instances as a compact laboratory for constraint modelling. He compared increasingly strong propagation schemes by asking which ones could solve a puzzle without search. My paper keeps that experimental setup, generates the instances instead, and changes the grid size.

The repository starts with 32,000 unique-solution puzzles at sizes 6x6, 9x9, 16x16, 25x25, and 36x36. Another 402,201 instances are saved along 170,000 hardness walks.

A hardness walk starts from one of the generated puzzles and adds clues from its known solution. Whenever an easier propagation scheme becomes sufficient, that version of the puzzle is saved. Alongside the difficult starting points, the repository therefore contains paths from those puzzles towards easier propagation categories.

Each puzzle comes as a text file with JSON metadata, and the repository includes tools for validating the corpus and selecting smaller subsets.

What changes with size#

Dominant hardness buckets change from Val at 6x6 to ValBS at 9x9, DomBS at 16x16, and Search at 25x25 and 36x36; propagation-only coverage is 100%, 100%, 84.9%, 0%, and 0%.

The figure adapts the cross-size summary from the paper. Each bucket names the weakest successful propagation scheme. Val uses value propagation, Dom uses the stronger domain-consistent model, and BS adds bounds shaving. Search means that none of the tested schemes solved the puzzle by propagation alone.

The shift with size is quite sharp. All generated 6x6 and 9x9 starting puzzles are covered by the tested family. At 16x16, propagation alone solves 84.9%. None of the 1,000 generated 25x25 or 1,000 generated 36x36 starting puzzles are solved without search, even with domain-consistent propagation and domain shaving.

That does not make the larger puzzles isolated. Adding clues moves them towards easier propagation categories, but the walk gets longer: the average number of clues needed for a verified category drop grows from 1.20 at 6x6 to 24.12 at 36x36.

Within this generated corpus, size changes which propagation schemes suffice, rather than merely making the same kind of puzzle slower. There is an important qualification: clue density is a feature of the generated puzzles, not an experimental parameter held fixed across sizes. Median fill rises from 27.8% at 6x6 to 52.2% at 36x36 as a result of the uniqueness-preserving clue-removal process. The comparison describes this collection; it is not a claim about every generalized Sudoku puzzle.

The hardness labels and walks give several ways to choose a useful subset from the public repository. The instances can then be reclassified or tested with other models and solvers. The repository is the result to take away from this work; the generation method and the map of its contents are in the paper.