Dead ends in square-free digit walks
Evan Chen, Chris Cummins, Ben Eltschig, Dejan Grubisic, Leopold Haller, Letong Hong, Andranik Kurghinyan, Kenny Lau, Hugh Leather, Seewoo Lee, Aram Markosyan, Ken Ono, Manooshree Patel, Gaurang Pendharkar, Vedant Rathi, Alex Schneidman, Volker Seeker, Shubho Sengupta, Ishan Sinha, Jimmy Xin, Jujian Zhang
TL;DR
This work studies dead ends in base-$b$ square-free digit walks, defining a dead end as a square-free $N$ with every one-digit extension $bN+d$ non-square-free. It shows the stochastic prediction of dead-end density is drastically overestimated in base-10, proving $D(X)=c_{ ext{dead}}X+O(X/\sqrt{\log X})$ with $c_{ ext{dead}}\approx 1.317\times 10^{-9}$, far smaller than the predicted $\approx 5.218\times 10^{-5}$. The main method is a rigorous inclusion-exclusion sieve that expresses $c_{ ext{dead}}(b)$ as a finite alternating sum of Euler products, valid for all bases $b\ge 2$, and yields explicit numerical values for small prime bases. The results are fully formalized in Lean/Mathlib via AxiomProver, and the paper provides explicit computational details and a numerical strategy for $c_{ ext{dead}}$—highlighting arithmetic dependencies that suppress dead ends relative to naive stochastic models.
Abstract
We study "dead ends" in square-free digit walks: square-free integers $N$ such that, in base $b$, every one-digit extension $bN+d$ is non-square-free. In base $10$, the stochastic independence model of Miller et al. suggests that infinite square-free walks occur with probability near $1$, corresponding to an asymptotic dead-end density of $\approx 5.218\times 10^{-5}$. We prove that the true asymptotic dead-end density satisfies \[ c_{\mathrm{dead}} \approx 1.317\times 10^{-9}, \] roughly a factor of $\sim 4\times 10^4$ smaller than the prediction. For every base $b\geq 2$, we prove that dead-end densities exist and are given by a closed-form expression (as a finite alternating sum of Euler products). The argument is fully formalized in Lean/Mathlib, and was produced automatically by AxiomProver from a natural-language statement of the problem.
