arXiv:2605.21098v2

A strange continued fraction associated with the Romik map

Yufei Chen, Karma Dajani, Yanyan Hu, Cor Kraaikamp

math.DSmath.NTmath.PR11J7028D05

Abstract

In 2008, Dan Romik studied in this journal Primitive Pythagorean Triples, or PPTs. In order to do so, he introduced a modified slow (subtractive) Euclidean algorithm, and showed that the underlying dynamical system of this Euclidean algorithm (the ``Romik system''), is ergodic and has a σσ-finite, infinite measure, of which is explicitly given. In this paper, the Romik system is further studied. Various basic properties are determined, such as the expansion of rational numbers and quadratic irrationals. Also (a version of) the planar natural extension of the Romik system is obtained, and the σσ-finite, invariant measure is explicitly given, and it is shown that it is ergodic. Furthermore, for Lebesgue almost every xx asymptotically half of the regular continued fraction (RCF) convergents of xx are among the Romik convergents. We also show that related to the Romik map a ``strange'' continued fraction can be given. ``Strange,'' as the set of possible partial quotients (i.e., digits) for any x[0,1]x\in [0,1] in this expansion is {0,±2}\{ 0, \pm 2\}. Various properties of this ``Romik expansion'' are given.

AI-generated audit

Audit summary

Audited against arXiv v2

Not a correctness certificate. A “Correct” result may include yellow typos or minor formal corrections that do not affect substantive soundness. It means this audit found no unresolved substantive error under the stated criteria; it does not replace expert scrutiny or formal verification.

Current report

Detailed mathematical audit

Generated August 15, 2026
01Statements3 reported findingsCorrect

The central statements checked are correct. The rational and quadratic-orbit results, the natural extension and invariant measure, the formal Romik expansion and convergence results, the conversion algorithm, and the almost-everywhere one-half deletion law follow from the supplied arguments. The literal notation defects listed in Part 2 have mechanically determined repairs and do not change these conclusions.

Propositions 2.1–2.5Correct

Arithmetic and dynamical foundations of the Romik map

Pages 5–11 · Propositions 2.1–2.5 · arXiv:2605.21098v2

For rational inputs, the branch formulas strictly reduce an appropriate denominator until the orbit reaches an endpoint. The regular-continued-fraction branch identities then imply eventual periodicity for quadratic irrationals. For the planar extension, the inverse branches, Jacobians, and density calculation give the invariant density (x+y2xy)2(x+y-2xy)^{-2} and its marginal 1/(x(1x))1/(x(1-x)); the branch coding supplies the asserted natural-extension and ergodicity conclusions.

Full paper, version 2
Theorem 3.1 and Proposition 4.3Correct

Formal expansion, convergents, and approximation identities

Pages 17–24 · Theorem 3.1, Lemma 4.1, Proposition 4.3, and Lemma 4.4 · arXiv:2605.21098v2

The branch-by-branch induction produces the stated digits and signs, and the matrix recurrence gives the convergents and determinant identity. Substitution of the tail Rk(x)R^k(x) yields the two exact approximation formulas after the absolute-value typo below is corrected. The analysis of repeated-convergent blocks proves that, for irrational xx, the denominators eventually exceed every prescribed bound.

Full paper, version 2
Algorithm 1 and Theorem 5.2Correct

Conversion from regular continued fractions and the metric deletion law

Pages 25–28 · Lemma 5.1, Algorithm 1, and Theorem 5.2 · arXiv:2605.21098v2

The insertion and singularization identities convert each regular partial quotient into the allowed Romik blocks and preserve the represented number. After the mechanically determined index corrections listed below, the algorithm covers the three possible first non-22 cases. Hopf's ratio theorem applies to the two indicator functions; the invariant-density integrals are log2\log 2 and 2log22\log 2, giving the asserted ratio 1/21/2.

Full paper, version 2
02Proofs5 reported findingsCorrect

The proofs of the central results are correct and complete after immediate-consequence closure. The findings below are genuine literal notation defects whose intended corrections are fixed by the adjacent matrix calculations or by the algorithm's own case definitions. They require textual correction but do not lower the proof status.

Central proof chainCorrect and complete

Branch identities, matrix recurrences, and ergodic ratio argument

Pages 5–28 · Sections 2–5 · arXiv:2605.21098v2

The denominator descent, finite-state periodicity argument, inverse-branch computation, continued-fraction matrix products, convergence analysis, and ratio-ergodic calculation supply every substantive step used by the main statements. The short natural-extension argument is recoverable from the displayed inverse branches and the generating branch coding; no additional estimate, compactness argument, or change of quantifiers is needed.

Full paper, version 2
Proposition 4.3(ii)Typo

The approximation formula is missing absolute-value bars

Pages 22–23 · Equation (38) and its proof · arXiv:2605.21098v2

Equation (38) prints xpn/qnx-p_n/q_n equal to a positive expression. Replace the left-hand side by xpn/qn\lvert x-p_n/q_n\rvert. The determinant identity in Lemma 4.1 fixes the omitted sign, and the proof immediately below computes the magnitude. For example, the displayed sample has x<p1/q1x< p_1/q_1, so the unsigned printed equality cannot hold. The corrected identity is exactly the estimate used later.

Full paper, version 2
Lemma 5.1Typo

The inserted-convergent indices are duplicated

Page 25 · Statement and matrix calculation in Lemma 5.1 · arXiv:2605.21098v2

The first inserted convergent is printed with denominator 2sn1+sn12s_{n-1}+s_{n-1}; replace the second sn1s_{n-1} by sn2s_{n-2}. In the subsequent product, replace the duplicated terms such as 2rn2+rn22r_{n-2}+r_{n-2} by 2rn1+rn22r_{n-1}+r_{n-2}, and make the analogous denominator correction. Direct multiplication of the displayed matrices uniquely gives these indices, and the final matrix already uses the correct recurrence.

Full paper, version 2
Algorithm 1Typo

Five symbols and indices contradict the surrounding case definitions

Pages 26–27 · Steps I–II and Equations (42)–(44) · arXiv:2605.21098v2

Step I defines nn as the first index with an2a_n\ne2, so its following infimum and empty-set condition must also use an2a_n\ne2, not an=2a_n=2. In the odd case, the exponent in Equation (42) must be (an1)/2(a_n-1)/2, matching the stated number of insertions. In case (iii), the number of insertions is an+11a_{n+1}-1, not an1a_n-1. In Step II, the terminal case is dn=2d_n=2 for every nn, not dn2d_n\ne2 for every nn. Finally, the continued fraction produced in case (iii) is Equation (44), not Equation (43). Each correction is forced by the immediately adjacent transformation and leaves the represented continued fraction unchanged.

Full paper, version 2
Lemma 4.1 and its setupTypo

The matrix group and greatest-common-divisor notation have the wrong signs

Pages 19–20 · Paragraph before Remark 4.1 and Lemma 4.1(ii) · arXiv:2605.21098v2

Matrices of determinant ±1\pm1 belong to GL2(Z)\mathrm{GL}_2(\mathbb Z), not necessarily SL2(Z)\mathrm{SL}_2(\mathbb Z), so replace the latter group name in this setup by the former. Also replace gcd(pn,qn)=±1\gcd(p_n,q_n)=\pm1 by gcd(pn,qn)=1\gcd(p_n,q_n)=1. The displayed determinant formula uniquely determines both corrections, and every subsequent Möbius and coprimality argument uses only determinant ±1\pm1.

Full paper, version 2
03Novelty0 reported findingsNo non-novelty findings

No non-novelty findings.

Detailed audit reportFull reasoning, manuscript locations, and references.
Open report PDF ↗

Author response

Challenge an audit finding

Local workflow preview

A listed author may submit formal evidence that an audit is inaccurate. The response would be considered in a fresh AI re-evaluation; it would not edit the audit automatically.

Paper
arXiv:2605.21098v2
Authors listed
Yufei Chen, Karma Dajani, Yanyan Hu, Cor Kraaikamp
Audit date
August 15, 2026
  1. 01Establish identityMatch an authenticated scholarly identity to this paper.
  2. 02Submit evidenceIdentify the finding and give a formal mathematical response.
  3. 03Re-evaluateA separate agent checks the response and records a disposition.
Recommended production method

Authenticate with ORCID, then require an exact arXiv match

MathAudit should accept the identity only when ORCID OAuth authenticates the claimant's iD and this exact arXiv paper appears in arXiv's public authority feed for that iD. A matching name alone is not sufficient.

ORCID OAuth and arXiv authority-record lookup are not connected in this local prototype.

Email fallback for papers without a linked ORCID

A production fallback could send a one-time link only when the submitted address matches an independently maintained author-contact allowlist for this paper. MathAudit must return the same message for every address so the form cannot reveal which contacts are on that list.

This demonstration does not send, store, or compare the address.

Structured response preview

This form remains unavailable until production identity verification succeeds. Nothing entered here is submitted.

This panel never establishes authorship in the local prototype. A production result should be described narrowly as an authenticated ORCID match or control of a separately allowlisted author-contact mailbox.