problems / two-smooth-summands

Also in this folder: compute/ · figures/

How the two-smooth-summands search came together

0. What was actually missing

The missing object was not another table of the first nn that fails to be a sum of two yy-smooth numbers. That table is OEIS A062241, it is finite for each fixed yy, and Green and Bloom both warn that the uniform statement cannot be settled by a finite computation.

What would have counted is a template that lifts: a finite list of residue classes, or an explicit arithmetic construction depending on nn in a checkable way, that writes every large nn as a+ba+b with P+(ab)≤nεP^+(ab)\le n^\varepsilon for some ε<1/2\varepsilon<1/2. Isolated prefixes — even a prefix on which F(n)≤n2/5F(n)\le n^{2/5} has only sixteen exceptions — are an incomplete search unless the template continues.

The analytic gap is the same object Balog named in 1989: positivity of the binary smooth convolution

Ry(N)=∑a=1N−11P+(a)≤y 1P+(N−a)≤yR_y(N)=\sum_{a=1}^{N-1}\mathbf{1}_{P^+(a)\le y}\,\mathbf{1}_{P^+(N-a)\le y}

at y=Nεy=N^\varepsilon, uniformly in NN. Fouvry's large-modulus theorem gets him down to 4/(9e)4/(9\sqrt e). A covering search cannot replace that input.

1. Named false starts, with the specific obstruction

Square-plus-remainder as a new bound. Write n=m2+rn=m^2+r with 0≤r≤2m0\le r\le 2m, or (m−1)2+(2m−1)(m-1)^2+(2m-1) when nn is square. Then F(n)<2n+1F(n)<2\sqrt n+1 for every n≥2n\ge 2. The verifier in trivial_cover.py confirms this on [2,2⋅105][2,2\cdot 10^5] and on several 101210^{12}–101810^{18} spots. Balog's first page already calls the N1/2N^{1/2} statement “almost trivially true.” The obstruction is not mathematical; it is bibliographic. This is the covering we were told not to claim.

A fixed list of shifts, one per residue class. For each r(modM)r\pmod M try the first forty MM-smooth offsets aa and ask whether n−an-a is nαn^\alpha-smooth. Modulo 88 at α=0.4\alpha=0.4 this almost works on [2,20000][2,20000]: seven classes lift and class 55 dies at n=18181n=18181. At α=1/3\alpha=1/3 every class dies, the first deaths sitting between 937937 and 42824282. A fixed offset cannot follow a smoothness bound that grows with nn. That is the obstruction.

Square-adjustment. Replace m=⌊n⌋m=\lfloor\sqrt n\rfloor by m−km-k for k≤256k\le 256 and demand both m−km-k and n−(m−k)2n-(m-k)^2 be n0.4n^{0.4}-smooth. Three hundred and nine holes remain on [2,80000][2,80000], the last at 7992279922. Making the first summand a nearby square of a smooth integer does not buy an exponent below 1/21/2.

Even reduction. If m=a+bm=a+b then 2m=(2a)+(2b)2m=(2a)+(2b), so F(n)≤max⁡(2,F(odd part(n)))F(n)\le\max(2,F(\mathrm{odd\ part}(n))). The identity is correct and useless: the hard nn are odd, and the negative-pseudosquare construction lives in the class 7(mod8)7\pmod 8.

Short intervals of length nαn^\alpha as a proof of F(n)≤nαF(n)\le n^\alpha. This would lift if every interval (n−nα,n)(n-n^\alpha,n) contained an nαn^\alpha-smooth integer. For α=1/2\alpha=1/2 the only genuine holes on [2,106][2,10^6] are n=7,8,23,24n=7,8,23,24, and FF itself exceeds n\sqrt n only at 3,7,233,7,23. For α=0.4\alpha=0.4 the short-interval holes run up to 34963496; two large summands still save those nn as soon as one computes FF, but the short-interval template does not cover them. Unconditionally, the literature (Balog 1987, Harman, Matomäki) gives xεx^\varepsilon-smooth numbers in intervals of length x\sqrt x, which again returns F(n)≤nF(n)\le\sqrt n, not a new exponent.

2. The useful failure

The residue-class search dying at a single class modulo 88 looked, for a moment, like a near-cover. It taught the opposite lesson. The seven live classes were live because a short list of small smooth offsets happened to hit an n0.4n^{0.4}-smooth complement on a short prefix. Class 55 is not arithmetically special; it is the first place the same finite list ran out of luck. Lengthening the list, or the modulus, only postpones the first death. Smoothness is a multiplicative monoid. A finite set of additive residue templates cannot manufacture two large smooth summands for every nn.

The same failure, read against Balog's reduction, is sharper. Balog splits off the NαN^\alpha-smooth kernel of NN and reduces to NαN^\alpha-rough NN — primes and products of large primes. Those are exactly the inputs a small-shift cover never sees: the complement n−an-a is then comparable to nn, and asking it to be nαn^\alpha-smooth is the original problem.

3. The click

The click was not a cover. It was the realisation that the negative-pseudosquare obstruction, which is the only infinite local obstruction we can write down, is too weak to produce F(n)>n2/5F(n)>n^{2/5} past a few hundred, and too strong to let any cover we can certify go below the Burgess exponent.

If m≡7(mod8)m\equiv 7\pmod 8 and (−m/q)=1(-m/q)=1 for every odd prime q≤yq\le y, then mm is not a sum of two positive yy-smooth integers. The proof is Jacobi symbols plus quadratic reciprocity on the even and odd summands; it does not use Dirichlet or Linnik. Consequently G(y)≤A045535(π(y)−1)G(y)\le A045535(\pi(y)-1), and for primes p≡3(mod4)p\equiv 3\pmod 4 the least quadratic non-residue satisfies n(p)≤F(p)n(p)\le F(p). A uniform bound F(n)≤n0.1F(n)\le n^{0.1} would beat Burgess on that prime class. That is Green's / Wooley's wall, and it is real.

But the same construction, evaluated at the certified points, gives F(131486759)=83F(131486759)=83 while 1314867592/5≈1585131486759^{2/5}\approx 1585. The obstruction forbids y=79y=79, not y=n2/5y=n^{2/5}. Along the A062241 sequence the inequality F(G(y))>G(y)2/5F(G(y))>G(y)^{2/5} holds for the first few terms and then flips: F(1559)=17<15592/5≈18.93F(1559)=17<1559^{2/5}\approx 18.93. After that the obstruction no longer produces counterexamples to F(n)≤n2/5F(n)\le n^{2/5}. That is why a C sweep of [2,106][2,10^6] finds exactly sixteen exceptions to F(n)≤n2/5F(n)\le n^{2/5}, the last at n=479=G(11)n=479=G(11), and why the same sweep at exponent 0.270.27 still has exceptions at the right endpoint of [2,2⋅105][2,2\cdot 10^5].

So the cover we wanted — a finite residue recipe for some ε<1/2\varepsilon<1/2 — is blocked from below by the trivial square template and from above by Fouvry. The obstruction explains why ε=0.1\varepsilon=0.1 is hopeless and why ε=2/5\varepsilon=2/5 can look finite-exception on a long prefix without being a theorem we proved.

4. The argument, in the order it was found

Green #59 and Bloom #334 were fetched first. The live Bloom page (edited 03 April 2026) still states Balog's exponent and the no(1)n^{o(1)} conjecture; two comments add the A045535 comparison and the Vinogradov wall. Balog's five-page paper was fetched next. The first paragraph contains the two sentences that organised the night: N1/2N^{1/2} is almost trivial, and the method cannot reach NεN^\varepsilon.

The trivial covering was written down and certified so it could not be “discovered” later. Then the Jacobi obstruction was proved from scratch and run against every small negative pseudosquare, and against M=131486759M=131486759 with displayed roots and a two-pointer non-representation test. That gave the exact value F(M)=83F(M)=83.

Only then did the cover search run. Fixed shifts, square adjustments, short intervals, and the even reduction were each given a verifier and each produced holes. The useful picture was the exception sweep with integer tests F2>nF^2>n, F5>n2F^5>n^2, F3>nF^3>n: floating-point pow\mathtt{pow} had falsely flagged the cubes 8,27,1258,27,125 as n1/3n^{1/3} exceptions. The exact lists are finite on the computed prefixes and stop at values of G(y)G(y).

In parallel the C bitset independently reproduced A062241 through y=79y=79. The 73-smooth and 79-smooth counts through MM agreed with the Python generator to the integer. That is a finite certificate, not a new bound.

Stored under compute/certs/.

Figures: figures/g_of_y.png, figures/exceptions_two_fifths.png.

The bitset that certifies G(73)=131486759G(73)=131486759 ORs the 73-smooth characteristic vector, shifted by each of the 2217622176 73-smooth integers up to 200000200000, and finds no hole in [2,M−1][2,M-1]. The same MM has no 79-smooth splitting.

6. What is proved vs still open

Proved tonight, all checkable:

Still open, and not touched:

We did not beat the published record. A failed search with a verifier is the product.

2026-08-27 continuation

0. What was actually missing

The 17 August search killed finite shift lists. A dent still needed a formula a(n)a(n) that writes every large nn as two nεn^\varepsilon-smooth summands for some ε<1/2\varepsilon<1/2, or a two-factor shape that could be proved to work for all large nn. Isolated prefixes were already known not to count.

1. Named false starts

Floor-divisor as a new covering. Set u=⌊nε⌋u=\lfloor n^\varepsilon\rfloor and a=n−(n mod u)a=n-(n\bmod u). The remainder is nεn^\varepsilon-smooth by size. At ε=1/2\varepsilon=1/2 this is the square covering in different clothes. Below 1/21/2, ⌊n/u⌋\lfloor n/u\rfloor has size n1−ε>nεn^{1-\varepsilon}>n^\varepsilon and is nεn^\varepsilon-rough for a positive-density set of nn. First hole at n=3n=3 on every tested exponent; thousands of holes through 2000020000.

Largest power of two, triangular numbers, cubes. Each is a closed form. Each has a first hole at n=3n=3 and holes still appearing at n=20000n=20000. Cubes are worse than squares: the remainder window has exponent 2/32/3.

Two-factor with u≤n1/5u\le n^{1/5} as an effective 2/52/5. On [2,20000][2,20000] the holes are exactly the sixteen FF-exceptions, last 479479, and at 1/31/3 exactly the seventy-six exceptions, last 1819118191. The factor restriction is not losing splittings. It is also not N0N_0.

2. The useful failure

The floor-divisor rewrite taught that every “take a small remainder so it is automatically smooth” template is the short-interval template. The only closed-form construction that always puts a smooth number in (n−n1/2,n](n-n^{1/2},n] is a nearby square. Polynomials of degree d≥2d\ge 2 cannot do better by size: the gap is ≍n1−1/d≥n1/2\asymp n^{1-1/d}\ge n^{1/2}.

3. The click

For a prime PP and u∼Pε/(1−ε)u\sim P^{\varepsilon/(1-\varepsilon)}, the integer n=Pu+1n=Pu+1 makes floor-divisor choose a multiple of uu whose other factor is PP. The comparison P>nεP>n^\varepsilon is the same inequality as ε<1/2\varepsilon<1/2. That is an infinite failure family, not a prefix. The power-of-two formula has the same shape: n=2k+qn=2^k+q with qq the next prime after 2k−12^{k-1} fails even n1/2n^{1/2}.

So the cover we wanted is not hiding in a cleverer closed form of the same kind. A formula whose remainder is small cannot beat 1/21/2. A formula whose first summand is a sparse smooth set (powers of two, SS-units for a fixed SS, high powers) fails infinitely often. Making Balog–Sárközy effective still needs an analytic input we do not have.

4. The argument, in the order it was found

Green #59 and the Wayback copy of Erdős #334 were fetched first. Hildebrand–Tenenbaum pp. 468–469, the three arXiv neighbours, and Ki–Maier–Sankaranarayanan 2016 still name Balog 0.26950.2695 as the all-integers binary record. Then the parent verifier was started. Then the closed-form census, then the infinite families, then the polynomial size check, then two-factor. Two-factor matching FF exactly was a surprise and a disappointment: it looks like a theorem on a long prefix and is not one.

5. Computer search

Stored under compute/q1/certs/.

Replay: cd problems/two-smooth-summands/compute/q1 && ./run_all.sh.

6. What is proved vs still open

Proved this session, all checkable:

Still open, and not touched: the same four items as on 17 August. We did not beat the published record.