arxiv: 2604.23133 · v1 · submitted 2026-04-25 · 🧮 math.PR · cs.NA· math.CO· math.NA

Recognition: unknown

High-Precision Framework for Expected Hitting Times Analysis in the Dice-Sum Process

Tipaluck Krityakierne , Thotsaporn Aek Thanatipanonda

Authors on Pith no claims yet

Pith reviewed 2026-05-08 07:26 UTC · model grok-4.3

classification 🧮 math.PR cs.NAmath.COmath.NA

keywords expected hitting timedice sum processdynamic programmingtruncation error boundsperfect squaresstochastic processesovershoot correction

0 comments

The pith

A dynamic programming method with analytic truncation correction computes the expected dice-sum hitting time to perfect squares to 1017 decimal places.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a one-variable recurrence for the expected number of die rolls until the cumulative sum first enters a target set H. By truncating the state space at a cutoff N and adding a closed-form overshoot correction derived from the tail, the approach supplies rigorous two-sided bounds on the approximation error. This framework is applied to the set of perfect squares, producing an explicit numerical value for the expectation that is certified correct to more than a thousand decimal places. The method uses constant memory and linear time in the cutoff, making high-precision evaluation feasible for any fixed target set.

Core claim

The authors show that the expected hitting time for the cumulative sum of fair six-sided dice satisfies a single-variable dynamic-programming equation independent of the step index. Truncating the resulting infinite system at a large cutoff N and augmenting it with an analytically computed overshoot term yields strict two-sided error bounds. For the target set of perfect squares this procedure evaluates the expectation explicitly as 7.07976423755110510389555305690818489468..., correct to 1017 decimal places.

What carries the argument

The one-variable dynamic-programming formulation for expected hitting times together with an analytically derived overshoot correction term that supplies rigorous truncation-error bounds.

If this is right

The reported numerical value is the exact expectation for the perfect-square target set within the stated 1017-place guarantee.
The same truncation-plus-overshoot procedure applies to any other fixed target set and returns comparable two-sided error bounds.
All computations require only constant memory and run in time linear in the chosen cutoff.
The method removes explicit dependence on the number of rolls, replacing the usual two-variable formulation with a single-variable recurrence.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same truncation technique could be applied to compute higher moments or the full distribution of the hitting time rather than only its expectation.
The high-precision value for perfect squares offers a benchmark against which other approximation schemes for random-walk hitting times can be validated.
Extensions to dice with different numbers of faces or to non-uniform step distributions would follow by replacing the transition probabilities in the recurrence while retaining the overshoot correction.
The linear-time scaling suggests the framework remains practical even when the cutoff must be taken extremely large to reach machine precision.

Load-bearing premise

The tail behavior of the overshoot beyond any finite cutoff admits an exact analytic expression that produces strict two-sided bounds on the residual error.

What would settle it

A direct numerical solution of the full infinite system to a precision finer than 1017 decimals that differs from the reported value in any digit after the 1017th place.

Figures

Figures reproduced from arXiv: 2604.23133 by Thotsaporn Aek Thanatipanonda, Tipaluck Krityakierne.

**Figure 1.** Figure 1: Probability pn that the partial sums of fair die rolls ever hit the integer n, for 1 ≤ n ≤ 100. Proof. From the explicit formula, pn − 2 7 = 1 7 view at source ↗

read the original abstract

We study the expected number of rolls required for the cumulative sum of a fair six-sided die to first enter a prescribed target set $H\subset\mathbb{Z}_{\ge0}$. A one-variable dynamic-programming formulation is introduced that removes dependence on the roll count. Within this framework, the infinite process is truncated at a large cutoff $N$ and corrected by an analytically derived overshoot term that accounts for the rare event of exceeding $N$ before entering $H$. Explicit bounds on this residual yield a strict two-sided estimate of the truncation error. The method is numerically efficient, requiring constant memory and linear time in the cutoff. For the perfect-square target set $H=\{n^2:n\in\mathbb{N}\}$, all quantities are evaluated explicitly, yielding \[ \mathbb{E}[T]=7.07976423755110510389555305690818489468\ldots, \] provably correct to 1,017 decimal places. This constitutes the most precise result known to date and establishes a general framework for high-accuracy computation of discrete hitting times.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper gives a practical DP reduction plus bounded overshoot correction for high-precision expected hitting times on dice sums, with a concrete 1017-decimal result for perfect-square targets.

read the letter

This paper gives a practical DP reduction plus bounded overshoot correction for high-precision expected hitting times on dice sums, with a concrete 1017-decimal result for perfect-square targets. The core advance is a one-variable dynamic program that eliminates explicit dependence on the number of rolls, followed by truncation at a large cutoff N and an analytically derived overshoot term that supplies rigorous two-sided error bounds. For the perfect-square set they evaluate the expectation explicitly and report a value correct to 1017 places, which is a clear numerical improvement over prior approximations. The method runs in linear time with constant memory, which is attractive for repeated computations. The underlying Markov-chain and renewal arguments are standard, but the specific error-control construction appears carefully tailored to this setting and avoids the usual floating-point uncertainty. One soft spot is that the strength of the claim rests on the overshoot bounds being free of gaps; the abstract sketches the approach, so a referee would want to see the full derivation and any edge-case checks. The paper also focuses on a single target set, so the framework's performance on other H is asserted rather than demonstrated at length. Overall the central computation looks solid and non-circular. This is for readers working on exact or high-accuracy hitting-time calculations in discrete probability, such as random-walk or gambling problems. It deserves peer review because the error analysis is explicit and the numerical result is reproducible in principle.

Referee Report

0 major / 1 minor

Summary. The manuscript introduces a one-variable dynamic-programming formulation for the expected hitting time T of the cumulative sum of fair six-sided die rolls to a target set H. The infinite process is truncated at cutoff N and corrected by an analytically derived overshoot term that supplies strict two-sided bounds on the residual truncation error. For H the set of perfect squares, all quantities are evaluated explicitly to obtain the numerical value E[T] = 7.07976423755110510389555305690818489468… claimed to be provably correct to 1,017 decimal places.

Significance. If the overshoot bounds are rigorously established, the work supplies the highest-precision explicit value known for this hitting-time expectation and a general, parameter-free computational framework that is memory-efficient and linear-time in the cutoff. The direct evaluation (no fitting or self-referential definitions) and the provision of explicit two-sided error bounds constitute clear strengths for reproducibility and verification.

minor comments (1)

[Abstract] The abstract states that 'all quantities are evaluated explicitly'; the main text should clarify which components admit closed forms versus direct numerical evaluation under the DP.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary, recognition of the significance of the overshoot bounds and computational framework, and recommendation to accept the manuscript.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation relies on a standard one-variable dynamic programming recurrence for the expected hitting time of a Markov chain on the non-negative integers, combined with a truncation at finite cutoff N whose residual error is bounded using an analytically derived overshoot correction from renewal theory. Both the DP setup and the overshoot bounds are obtained directly from the problem definition and standard probabilistic arguments; the final high-precision numerical value for the perfect-square target set is produced by explicit evaluation of the truncated system rather than by fitting, renaming, or self-referential closure. No load-bearing step reduces to a prior result by the same authors, to a fitted parameter renamed as a prediction, or to an ansatz smuggled through citation. The method is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework relies on classical probability axioms for independent rolls and standard results on hitting times for random walks on the non-negative integers; no new free parameters are introduced and the overshoot correction is derived rather than postulated.

axioms (2)

standard math The six-sided die is fair and rolls are independent, following standard Kolmogorov probability axioms.
Used to define the transition probabilities of the sum process.
domain assumption An analytically derivable overshoot distribution exists for the truncation at N that yields strict two-sided error bounds.
This is the load-bearing step that converts the truncated computation into a rigorous approximation of the infinite-process expectation.

pith-pipeline@v0.9.0 · 5501 in / 1414 out tokens · 69171 ms · 2026-05-08T07:26:29.248250+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

When Does the Dice Sum Become Prime?
math.PR 2026-05 unverdicted novelty 6.0

The expected number of die rolls until the sum is prime is computed to more than 1000 decimal places via DP truncation with exponential-decay error bounds.

Reference graph

Works this paper leans on

15 extracted references · cited by 1 Pith paper

[1]

Feller,An Introduction to Probability Theory and Its Applications, Vol

W. Feller,An Introduction to Probability Theory and Its Applications, Vol. 1, John Wiley & Sons, New York, 1968

1968
[2]

Karlin,A First Course in Stochastic Processes, Academic Press, 1968

S. Karlin,A First Course in Stochastic Processes, Academic Press, 1968

1968
[3]

D. R. Cox,Renewal Theory, Methuen, London, 1962

1962
[4]

Feller,An Introduction to Probability Theory and Its Applications, Vol

W. Feller,An Introduction to Probability Theory and Its Applications, Vol. 2, John Wiley & Sons, New York, 1971

1971
[5]

Farahmand et al., Truncated approximate dynamic programming with task-dependent terminal value, inProceedings of the AAAI Conference on Artificial Intelligence30(1) (2016)

A. Farahmand et al., Truncated approximate dynamic programming with task-dependent terminal value, inProceedings of the AAAI Conference on Artificial Intelligence30(1) (2016)

2016
[6]

DasGupta, Solution to Puzzle 17,Institute of Mathematical Statistics Bulletin46(5) (2017), 9

A. DasGupta, Solution to Puzzle 17,Institute of Mathematical Statistics Bulletin46(5) (2017), 9. https://imstat.org/2017/07/14/ student-puzzle-corner-18-and-solution-to-puzzle-17/

2017
[7]

M. M. Conroy,A Collection of Dice Problems, Online note, 2018

2018
[8]

Alon and Y

N. Alon and Y. Malinovsky, Hitting a prime in 2.43 dice rolls (on average),The American Statistician77(3) (2023), 301–303

2023
[9]

Martinez and D

L. Martinez and D. Zeilberger, How many dice rolls would it take to reach your favorite kind of number?,Maple Transactions3(3) (2023), Autumn issue

2023
[10]

N. Alon, Y. Malinovsky, L. Martinez, and D. Zeilberger, Hitting k primes by dice rolls, The Electronic Journal of Combinatorics32(4) (2025), 4.16

2025
[11]

Billingsley,Probability and measure, John Wiley & Sons, 2017

P. Billingsley,Probability and measure, John Wiley & Sons, 2017

2017
[12]

Zeilberger, The C-finite Ansatz,The Ramanujan Journal31(1) (2013), 23–32

D. Zeilberger, The C-finite Ansatz,The Ramanujan Journal31(1) (2013), 23–32

2013
[13]

Kauers and P

M. Kauers and P. Paule,The Concrete Tetrahedron: Symbolic Sums, Recurrence Equations, Generating Functions, Asymptotic Estimates, Springer, Vienna, 2011. 15

2011
[14]

Krityakierne and T

T. Krityakierne and T. A. Thanatipanonda, C-finite, Holonomic, and C2-finite, inApplied Mathematical Analysis and Computations I: 1st SGMC, Statesboro, USA, April 2–3, 2021 (Virtual), Springer Proceedings in Mathematics & Statistics471(2024), 255–268

2021
[15]

Chern, Hitting a prime by rolling a die with infinitely many faces.The American Statistician78(3) (2024), 297–303

S. Chern, Hitting a prime by rolling a die with infinitely many faces.The American Statistician78(3) (2024), 297–303. 16

2024