REVIEW 4 major objections 6 minor 21 references
Solution of the 15 puzzle problem
T0 review · 4 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper establishes that the generalized 15 puzzle on an $n\times n$ torus mixes in order $n^4$ steps, with a single tile's position converging to Brownian motion on the torus at time $c_{\mathrm{puz}} n^4 t$ for…
desk verdict A serious, novel attack on Diaconis's 15-puzzle problem that likely gets the n^4 order right, but the sharp Brownian constant rests on a two-sentence characteristic-function bound (Lemma 22) that needs a real proof. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the renewal process of a marked tile: stopping times $t_i$ alternate between the empty square swapping with the tile from left or right and from above or below, with $H_i$ and $V_i$ the tile's horizontal and vertical displacements in those intervals and $r_i$, $s_i$ the elapsed times. The key identity is Lemma 21, which computes the limiting variance $s^2$ and mean time $\mu$ from the four directional return probabilities of Lemma 14, giving $c_{\mathrm{puz}}=2\mu/s^2=\frac52(\pi-1)$. The characteristic function of a single burst is expressed exactly as a resolvent formula (Lemma 23), then expanded at low frequencies; the resulting local limit theorem (Theorem 27) turns the renewal process into the Brownian $\theta$ kernel $\theta_t$. For several pieces, the mechanism shifts to a comparison of Dirichlet forms between the original puzzle, a symmetrized chain with overlapping pieces, and an abelian random walk with independent coordinates (Theorem 9), which yields the $d^2$ mixing bound $O(n^4)$.
What would settle it
Compute the limiting return probabilities in Lemma 14 by evaluating the $\mathbb{Z}^2$ lattice Green's functions $G_{1,0}$, $G_{1,1}$, and $G_{2,0}$ symbolically or to high precision; if the limits are not $\frac12$, $\frac12-\frac1\pi$, and $\frac2\pi-\frac12$, then $c_{\mathrm{puz}}$ changes, although the $n^4$ order would remain.
Extended reading notes
Core claim
The central discovery is that the mixing of a single tile in the $n^2-1$ puzzle reduces to a renewal process with a Brownian scaling limit. The tile moves only when the empty square is adjacent to it, and between such moves the empty square performs lazy simple random walk on the torus; the tile's displacement accumulates in alternating horizontal and vertical bursts. A local limit theorem for the triple (horizontal displacement, vertical displacement, elapsed time) of these bursts shows that the elapsed-time fluctuations are absorbed into the Brownian clock. Consequently the total variation distance of the tile at time $c_{\mathrm{puz}} n^4 t$ tends to $d_{\mathrm{Br}}(t)$, the total variation distance of Brownian motion on $(\mathbb{R}/\mathbb{Z})^2$, uniformly for $t$ in compact sets. For the whole board, the same renewal structure yields Poisson(1) fixed points after $n^4 f(n)$ steps and, via comparison with a 3-cycle walk, an upper bound of $O(n^4\log n)$.
Load-bearing premise
The sharp constant $\frac52(\pi-1)$ rests on four directional return probabilities in Lemma 14—stated as $\frac12$, $\frac12-\frac1\pi$, $\frac2\pi-\frac12$ with only a computer-algebra calculation as support—being exactly correct.
Editorial extensions
If this is right
- A single tile in an $n\times n$ puzzle needs order $n^4$ moves to mix: at time $c_{\mathrm{puz}} n^4 t$ its distance to uniformity approaches that of torus Brownian motion at time $t$.
- After $n^4 f(n)$ random moves with $f(n)\to\infty$, the number of tiles in their original positions converges to $\mathrm{Poisson}(1)$; if $f(n)$ stays bounded, this convergence fails.
- The fixed-point statistic shows no cut-off phenomenon: mixing happens gradually over order $n^4$ steps.
- The full board's total-variation and $\epsilon/|G|$-$\ell^\infty$ mixing time is $O(n^4\log n)$.
- For each fixed $d$, the empty square and any $d$ labeled pieces can be coupled with their stationary copies in expected time $O(n^4\log n)$.
Reading between the lines
- The renewal-and-local-limit structure should transfer to other 'moving-hole' walks, such as higher-dimensional grids or rectangular boards, where only the return-probability calculation changes.
- The Poisson(1) fixed-point result suggests the induced permutation is locally uniform on low-complexity representations already at order $n^4$; a further test would be whether the full board's mixing time can be sharpened from $O(n^4\log n)$ to $O(n^4)$.
- One could probe the sharp constant empirically by simulating small boards and comparing the single-tile distance to $d_{\mathrm{Br}}(c_{\mathrm{puz}}^{-1} t)$ at times $c_{\mathrm{puz}} n^4 t$; a mismatch would pinpoint the return-probability calculation as the fragile step.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to solve Diaconis' '15 puzzle problem' for the n x n puzzle on the torus with lazy random moves. Theorem 1 asserts that a single labeled piece converges in total variation, at time c_puz n^4 t with c_puz = (5/2)(π - 1), to the corresponding Brownian motion limit, uniformly on compact time intervals. Theorem 2 states that after n^4 f(n) steps with f(n) -> infinity, the number of fixed points converges to Poisson(1), while this convergence fails if f is bounded; Corollary 3 concludes there is no total variation cutoff for this statistic. Theorem 4 gives an O(n^4 log n) total variation mixing time upper bound for the full puzzle, and Theorem 5 states an O(n^4 log n) expected coupling time for the empty square and any fixed number of labeled pieces. The main technical ingredients are a renewal description of a tracked piece's motion, a local limit theorem for renewal sums (Theorem 27), spectral comparison arguments, and a comparison with the 3-cycle walk for the upper bound.
Significance. If the proofs can be completed, the paper resolves a long-standing problem of Diaconis and establishes a sharp, parameter-free Brownian constant c_puz for the single-piece mixing. The Poisson fixed-point result and the absence of cutoff for that statistic are natural and interesting consequences, and the coupling theorem and spectral comparison framework provide tools for further study. The derivation of the constant is not circular: it is obtained from limiting return probabilities of the empty square rather than fitted to the target Brownian distance. The paper also presents a new renewal/local-limit method for this type of constrained random walk on the symmetric group. These strengths make the paper a potentially significant contribution to the probability literature, provided the under-derived technical estimates are supplied.
major comments (4)
- [§5, Lemma 14 (Eq. (53))] The uniform bound |chi(xi1, xi2)| <= 1 - c max(xi1^2/n^2, xi2^2) is the load-bearing estimate for the truncation step in Theorem 27, and hence for the exact Brownian constant in Theorem 1. The proof is only two sentences: it invokes paths of bounded length or fixed displacement with positive probability and asserts that their phase variation has the given magnitude. That argument establishes, at best, that a single path has phase close to 1 in some regimes; it does not prove a uniform quadratic lower bound for 1 - |chi|, which requires controlling cancellation among many paths. The claim is especially delicate for xi2 near 1/2 and for xi1 near 0, where the phase scale is small relative to the stated quadratic bound. The authors should provide a rigorous derivation of Lemma 22, or replace it with a proved estimate that supports the frequency truncation in Theorem 27.
- [§9, Lemma 42] The limiting return probabilities p_(1,0) = 1/2, p_(0,±1) = 1/2 - 1/pi, and p_(-1,0) = 2/pi - 1/2 are stated after 'The exact values were calculated in Mathematica,' with no derivation shown. These values enter Lemma 21 and determine the constant c_puz = (5/2)(pi - 1) in Theorem 1. The values are checkable from the potential kernel and harmonic relations, but the reader should not be expected to trust a black-box computation at a load-bearing point of the proof. Please include an analytic derivation, or at least a complete and reproducible calculation, of the Green-function limits G_(1,0), G_(1,1), and G_(2,0) used in equations (60).
- [§6, Theorem 27] The proof of Lemma 42, which claims the d2 mixing time of the measure mu_S is O(n^2 log n), is only a sketch. In the one-dimensional representation case, the assertion that the spectrum is 'mixed in O(n^2 log n) steps' is stated without a quantitative argument covering the O(n^2) such representations. For higher-dimensional representations, the claim that an arbitrary factor n^2 in multiplicity can be saved by increasing the constant in the n^2 log n mixing time requires a justification, since the Plancherel sum includes multiplicities and the spectral gaps are only known up to constants. This lemma supports Theorem 4, so a complete proof should be supplied.
- [§6, Theorem 27] In the contour-shift step between equations (117) and (118), the horizontal integral is dismissed as bounded by O_A(n^{-A}) without a demonstration. The shift moves the integration contour by i(t - 2N mu_n)/(2 N^{1/2} v_n), and the resulting horizontal contribution is controlled only if the Gaussian factor decays uniformly in the summation range and in the error region; this decay is not shown. Since Theorem 27 is the main local limit theorem used in the proof of Theorem 1, this estimate should be stated and proved explicitly.
minor comments (6)
- [§8, Lemma 39] The notation in the lemma statement is garbled: 'let py,n,y P tp 1, 0q, p´1, 0q, p0, 1q, p0, ´1qu' should read 'let p_{y,n}, y in {...},' and the subscript typesetting throughout Section 5 should be cleaned up.
- [§8, Lemma 40] The proof of Lemma 39 is very brief and uses an 'auxilliary step' that is not precisely defined. Since the lemma is used to justify the Gaussian approximation in Lemmas 40 and 41, please spell out the construction and the induction more carefully.
- [§1, Abstract] In the proof of Lemma 40, the notation 'D d_m' in display (172) is unclear; it should be D^d_m or a similar indexed family. The same issue appears in the induction argument in Lemma 41. This is a presentation issue, but it makes already dense arguments harder to check.
- [§1, Introduction] The abstract contains the typo 'annˆ n numbered grid'; please correct to 'an n x n numbered grid'.
- [§6, Proof of Theorem 1] The related-work paragraph cites the authors' earlier paper [15], but that paper does not appear to be used in the proofs. If it is not needed, this is fine; if it is intended to motivate the renewal technique, a sentence clarifying the connection would help.
- [§8, Deduction of Theorem 2] The upper-bound part of the proof states that 'convolution with the remaining distributions can only decrease the total variation distance, as can removing the conditioning.' This is not immediate for the conditional distribution described, and a brief justification or reference would improve the exposition.
Circularity Check
No circularity: the Brownian constant is derived from independent return probabilities, not fitted to the target mixing time.
full rationale
The derivation is self-contained against external benchmarks. Theorem 1's Brownian constant is computed from the return probabilities of simple random walk on Z^2 (Lemma 14), via the renewal-theoretic variance and mean in Lemma 21; no term in c_puz = 2mu/s^2 is fitted to the Brownian distance d_Br(t). The local limit theorem (Theorem 27) supplies the asymptotic independence of the renewal sum and is proved by Fourier inversion with explicit error bounds rather than by importing the target result. Theorem 2's Poisson(1) limit follows from moment estimates and the O(n^4) d2 mixing bound of Theorem 30, which is established by comparison with symmetrized chains whose spectral estimates are proven directly in Lemmas 40-41. The only author self-citation is [15], which appears in a related-work survey and is not used as an input to any theorem. Lemma 14's Mathematica computation is an unshown calculation, not a circular step, and Lemma 22's brevity is a completeness concern rather than a circularity. No derivation step reduces by construction to its own inputs.
Assumptions & free parameters
assumptions (4)
- standard math Comparison theorem for Dirichlet forms of reversible Markov chains (Diaconis-Saloff-Coste, Theorem 8 here).
- domain assumption Spectral gap 3/(n−1) and mixing time of order n log n for the random walk generated by 3-cycles on Alt(n) (Lemma 11, from Helfgott-Seress-Zuk [14]).
- standard math Standard recurrence and Green's function facts for simple random walk on Z^2 (e.g., return probabilities and hitting time asymptotics).
- standard math Plancherel formula and Fourier inversion on finite groups.
Cite this review
Pith. "Pith review of Solution of the 15 puzzle problem." pith.science (2026). https://pith.science/paper/LQGV2LNF
@misc{pith2026190807106,
author = {Pith},
title = {Pith review of: Solution of the 15 puzzle problem},
year = {2026},
howpublished = {\url{https://pith.science/paper/LQGV2LNF}},
note = {Machine review of arXiv:1908.07106}
}
abstract
A generalized `$15$ puzzle' consists of an $n \times n$ numbered grid, with one missing number. A move in the game switches the position of the empty square with the position of one of its neighbors. We solve Diaconis' `15 puzzle problem' by proving that the asymptotic total variation mixing time of the board is at least order $ n^4 $ when the board is given periodic boundary conditions and when random moves are made. We demonstrate that for any $f(n) \to \infty$ with $n$, the number of fixed points after $n^4 f(n)$ moves converges to a Poisson distribution of parameter 1. The order of total variation mixing time for this convergence is $n^4$ without cut-off. We also prove an upper bound of order $n^{4 }\log n$ for the total variation mixing time.
Figures
Reference graph
Works this paper leans on
-
[1]
Cut-off for conjugacy-invariant random walks on the permutation group
Berestycki, Nathana¨ el, and B. S ¸eng¨ ul. “Cut-off for conjugacy-invariant random walks on the permutation group.” Probability Theory and Related Fields 173.3-4 (2019): 1197-1241
work page 2019
-
[2]
Cut-off for random to random c ard shuffle
Bernstein, M. and E. Nestoridi. “Cut-off for random to random c ard shuffle.” Annals of Probability , to appear
-
[3]
On the spectral gap for fin itely-generated subgroups of SU p2q
Bourgain, Jean, and Alex Gamburd. “On the spectral gap for fin itely-generated subgroups of SU p2q.” Inventiones mathematicae 171.1 (2008): 83-121
work page 2008
-
[4]
Uniform expansion bounds f or Cayley graphs of SL 2pFpq
Bourgain, Jean, and Alex Gamburd. “Uniform expansion bounds f or Cayley graphs of SL 2pFpq.” Annals of Mathematics (2008): 625-642
work page 2008
-
[5]
Group representations in probability and statistics
Diaconis, Persi. Group representations in probability and statistics. Lecture notes-monograph series 11 (1988): i-192
work page 1988
-
[6]
Some things we’ve learned (about Markov chain M onte Carlo)
Diaconis, Persi. “Some things we’ve learned (about Markov chain M onte Carlo).” Bernoulli 19.4 (2013): 1294-1305
work page 2013
-
[7]
Comparison techniq ues for random walk on finite groups
Diaconis, Persi, and Laurent Saloff-Coste. “Comparison techniq ues for random walk on finite groups.” The Annals of Probability (1993): 2131-2156. 56 YANG CHU AND ROBERT HOUGH
work page 1993
-
[8]
Comparison theore ms for reversible Markov chains
Diaconis, Persi, and Laurent Saloff-Coste. “Comparison theore ms for reversible Markov chains.” The Annals of Applied Probability (1993): 696-730
work page 1993
Show all 21 references
-
[9]
Generating a rando m permutation with random transpositions
Diaconis, Persi, and Mehrdad Shahshahani. “Generating a rando m permutation with random transpositions.” Probability Theory and Related Fields 57.2 (1981): 159-179
1981
-
[10]
Probability: theory and examples
Durrett, Rick. Probability: theory and examples. Vol. 49. Cambridge University Press, 2019
2019
-
[11]
Finding the shortest move-sequence in the g raph-generalized 15-puzzle is NP-hard
Goldreich, Oded. “Finding the shortest move-sequence in the g raph-generalized 15-puzzle is NP-hard.” Studies in Complexity and Cryptography. Miscellanea on the Interplay between Randomness and Computation . Springer, Berlin, Heidel- berg, 2011. 1-5
2011
-
[12]
Growth and generation in SL 2pZ{pZq
Helfgott, H. A. “Growth and generation in SL 2pZ{pZq.” Annals of Mathemat- ics, vol. 167, no. 2, 2008, pp. 601-623
2008
-
[13]
On the diameter of permutation groups
Helfgott, Harald A., and ´Akos Seress. “On the diameter of permutation groups.” Annals of Mathematics (2014): 611-658
2014
-
[14]
Random generators of the symmetric group: diameter, mixing time and spectral gap
Helfgott, Harald A., ´Akos Seress, and Andrzej Zuk. “Random generators of the symmetric group: diameter, mixing time and spectral gap.” Journal of Algebra 421 (2015): 349-368
2015
-
[15]
The random k cycle walk on the symmetric group
Hough, Bob. “The random k cycle walk on the symmetric group.” Probability Theory and Related Fields 165.1-2 (2016): 447-482
2016
-
[16]
Sandp iles on the square lattice
Hough, Robert D., Daniel C. Jerison, and Lionel Levine. “Sandp iles on the square lattice.” Communications in Mathematical Physics 367.1 (2019): 33-87
2019
-
[17]
Notes on the ‘15 ’ puzzle
Johnson, Wm Woolsey, and William Edward Story. “Notes on the ‘15 ’ puzzle.” American Journal of Mathematics 2.4 (1879): 397-404
-
[18]
Peres, and E.L
Levin, D.A., Y. Peres, and E.L. Wilmer. Markov chains and mixing times. American Mathematical Soc., 2009
2009
-
[19]
A real-time algorithm for the pn2 ´ 1q-puzzle
Parberry, Ian. “A real-time algorithm for the pn2 ´ 1q-puzzle.” Information Processing Letters 56.1 (1995): 23-28
1995
-
[20]
Finding a shortest s olution for the N ˆ N extension of the 15-puzzle is intractable
Ratner, Daniel, and Manfred K. Warmuth. “Finding a shortest s olution for the N ˆ N extension of the 15-puzzle is intractable.” AAAI. 1986
1986
-
[21]
Tao, Terence, and Van H. Vu. Additive combinatorics. Vol. 105. Cambridge University Press, 2006. Department of Mathematics, Stony Brook University, Stony B rook, NY, 11794 E-mail address, Yang Chu: yang.chu@stonybrook.edu E-mail address, Robert Hough: robert.hough@stonybrook.edu
2006
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.