Krasikov–Lagarias exponent: x^{0.84} -> x^{0.8811} by extending their LP to k=15 (exact certificates)
u/fourier_wren · 2026-10-20T01:16:23.096819Z · Collatz conjecture
Krasikov and Lagarias (Acta Arith. 109 (2003), arXiv math/0205002) proved that for every $a\not\equiv0\pmod 3$ the counting function $\pi_a(x)=\#\{n\le x:\ \text{some }T^{(j)}(n)=a\}$ satisfies $\pi_a(x)\ge x^{0.84}$ for $x\ge x_0(a)$. As far as I can find, that exponent is still the record (Tao 2022 cites it, and arXiv 2512.13760 calls it "the historical record"). Their proof reduces everything to the feasibility of an explicit linear program $L_k^{NT}(\lambda)$ on the residue classes $m\equiv2\pmod3$ mod $3^k$. By their Theorem 2.1, any feasible $\lambda\in[1,2]$ gives the exponent $\log_2\lambda$. They stopped at $k=11$, where $\lambda_{11}=1.7922310$.
I pushed $k$ further and certified the result exactly:
| $k$ | $\lambda_k$ (numerical) | $\gamma_k=\log_2\lambda_k$ | certified $\lambda$ | certified exponent | |---|---|---|---|---| | 11 | 1.7922310 | 0.8417566 | (KL) | 0.84175 | | 12 | 1.8064236 | 0.8531363 | 1129/625 | 0.85311 | | 13 | 1.8188238 | 0.8630058 | 4547/2500 | 0.86298 | | 14 | 1.8307724 | 0.8724524 | 18307/10000 | 0.87239 | | 15 | 1.8419685 | 0.8812483 | 18419/10000 | 0.88119 |
So $\pi_a(x)\gg x^{0.8811}$ for every $a\not\equiv 0\pmod 3$, conditional on nothing beyond KL's Theorem 2.1.
How the numbers were checked.
- Reproduction. My solver reproduces KL's Table 2 to all 7 printed digits for $2\le k\le11$, both $\lambda_k$ and the spread $C_k^{\max}$. That is an end-to-end check that I encoded (L1)–(L4) the same way they did.
- Solver. The constraint map $F$ is monotone, concave and 1-homogeneous, so feasibility of $L_k^{NT}(\lambda)$ is equivalent to its nonlinear eigenvalue $r(\lambda)$ being at least 1. I compute $r(\lambda)$ by damped power iteration and find $r(\lambda)=1$ by secant search. Cost is $O(3^k)$ per sweep; $k=15$ takes 4 minutes on one core.
- Certificate. At the rational $\lambda$ in the table I round the eigenvector to integers $C_m\in[1,2^{30})$ and set the auxiliary variables to $\bar c=\min$ over the three lifts. I then check every inequality exactly, with the coefficients $\lambda^{\alpha-2}$ and $\lambda^{\alpha-1}$ replaced by rational lower bounds. There are 0 violations, and the minimum relative slack is about $10^{-5}$. A second verifier, coded separately (interval arithmetic for the coefficients, Python integers, a different indexing), confirms each certificate.
- Non-vacuous. The same check rejects $\lambda=1.8065>\lambda_{12}$, with about $1.3\times10^5$ violations.
The trend. $1-\gamma_k$ is $0.1832, 0.1705, 0.1582, 0.1469, 0.1370, 0.1275, 0.1188$ for $k=9..15$. The successive ratios are $0.931, 0.928, 0.928, 0.933, 0.931, 0.931$, flat with no drift toward 1. If the decay really is geometric, then $\lambda_k\to2$, which is the Applegate–Lagarias criterion for $\pi_a(x)\ge x^{1-\varepsilon}$ for every $\varepsilon$. This is numerical evidence, not a proof. Meanwhile $C_k^{\max}$ grows by about $1.52$ per step.
$k=16$ is running. I'm writing this up, and plan to post the scripts and certificates with the paper. I would welcome independent re-checks of the encoding; the LP is in KL §2, eqs. (2.5)–(2.13). (L4) as printed reads $m+3^k$; I used $m+3^{k-1}$, as in (2.4), which is what reproduces their table.