Preprint
Article

This version is not peer-reviewed.

On Problem 477 of Erdős and Graham: The Cubic Case

Submitted:

24 September 2026

Posted:

12 October 2026

You are already at the latest version

Abstract
Let \( C=\{m^3:m\in Z\} \). We prove that there is a set \( A\subseteq Z \) for which the addition map \( A\times C\to Z \) is a bijection. If \( k \) is a fixed noncube and \( [E_k(N)=\{n\in[N,2N)\cap Z:n^3-k\in C-C\}, \) we prove the quantitative estimate \( \#E_k(N)\ll_{k,\delta} N^{18/19+\delta} \qquad(\delta>0). \) The small-gap contribution is handled by a square sieve on the progression polynomial attached to the gap \( a=y-x \); all local branches, including sieve primes dividing \( a \) and the zero values of the polynomial, are treated explicitly. For larger gaps, a theorem of Heath–Brown on fixed levels of diagonal ternary forms gives the required power saving after a direct height calculation and an exclusion of every special solution. Density zero then gives simultaneous avoidance for every finite collection of noncubes, and a two-sided greedy construction produces \( A \). Independently of this source-level proof, a complete Lean 4 development checks the exact theorem and a self-contained determinant–conic–binary-sextic route to the coarser bound \( O_k(N^{99/100}) \). The formal sources, dependency lock, and axiom audit are supplied as a separate companion archive; they are not used as a substitute for the mathematical arguments in the paper.
Keywords: 
;  ;  

1. Introduction

1.1. The Problem and the Main Theorem

For subsets A , B ⊆ Z , write Z = A ⊕ B if every integer has exactly one representation a + b with ( a , b ) ∈ A × B . Erdos and Graham asked whether f ( Z ) can be a direct additive factor of Z for some integer-valued polynomial f of degree at least two; see [3,4]. For background on factorizations of the infinite cyclic group, see also Sekanina [8]. The quadratic case is impossible. Consequently, the cubic result below gives the least possible degree in that formulation.
The word unique is essential. Ordinary additive-complement results ask only for coverage, while a direct factorization requires both existence and uniqueness of every representation. Results on ordinary complements of positive powers, such as [2], therefore address a different problem. There are also recent non-refereed claims for higher powers, including a public manuscript treating thirteenth powers and an online claim covering degrees at least five; see [7,9]. They do not contain the cubic statement established here.
Our main theorem is as follows.
Theorem 1.
Let
C = { m 3 : m ∈ Z } .
There exists A ⊆ Z such that
Z = A ⊕ C .
The proof gives more than the qualitative conclusion. Put D = C − C and, for a fixed noncube k, count the cube indices n for which n 3 − k belongs to D. The power saving in the following theorem is the analytic input needed by the final combinatorial construction.
Theorem 2.
Fix k ∈ Z ∖ C . For N ≥ 1 , let
E k ( N ) = { n ∈ [ N , 2 N ) ∩ Z : n 3 − k ∈ D } .
For every δ > 0 there are constants C k , δ > 0 and N k , δ such that, for all N ≥ N k , δ ,
# E k ( N ) ≤ C k , δ N 18 / 19 + δ .
Theorem 2 is stated for distinct values of n, not for representations. Throughout the proof we are free to overcount by the auxiliary gap or by an integer triple: projecting any such collection to its n-coordinate can only decrease its cardinality.

1.2. Context and Scope

The result has two logically separate certificates. The proof printed here uses two published analytic estimates: Heath–Brown’s square sieve [5], together with the mixed character-sum estimate of Bonolis–Browning [1], for small gaps; and Heath–Brown’s fixed-level diagonal-form estimate [6] for large gaps. Their hypotheses and the substitutions used here are stated before either application.
The companion formal development follows a longer route because it does not postulate the fixed-level diagonal-form theorem as an axiom. It reconstructs the needed determinant method, common conics, binary sextics, height estimates, and factor counts in Lean. That development proves the exact statement of Theorem 1 and the independent quantitative bound O k ( N 99 / 100 ) . Machine verification establishes correctness relative to the Lean kernel and the stated foundations; it does not establish novelty, publication priority, or independent human peer review.

1.3. Proof Strategy

We now describe the proof in enough detail to display every change of scale.
1.
A finite simultaneous-escape property for a symmetric set S ⊆ Z is converted into a direct factorization Z = A ⊕ S by a two-sided greedy construction.
2.
For S = C , escape is reduced to showing that the bad indices E k ( N ) have density zero for every fixed noncube k.
3.
A representation n 3 − k = y 3 − x 3 is parametrized by the positive gap a = y − x and the sum v = x + y . Below G N , the resulting square condition is estimated by the square sieve after completing the relevant Jacobi-symbol correlations. With G N = N 17 / 19 this contributes O k ( N 18 / 19 log 4 N ) .
4.
Above G N , the gap identity gives | x | , | y | ≪ k N 20 / 19 . Thus every remaining representation gives a point of height B ≪ k N 20 / 19 on the fixed-level diagonal cubic x 3 + n 3 − y 3 = k .
5.
Heath–Brown’s theorem counts these points by O k , ε ( B 9 / 10 + ε ) . Its exceptional solutions would force k to be a cube, so none occur. Since ( 20 / 19 ) ( 9 / 10 ) = 18 / 19 , the large-gap count is O k , ε ( N 18 / 19 + ( 20 / 19 ) ε ) .
6.
Combining the two ranges proves Theorem 2; finite simultaneous escape and the greedy construction then prove Theorem 1. A squarefree-kernel reorganization is included as an arithmetic cross-check of the small-gap argument.
All asymptotic statements below are for sufficiently large N, with the threshold allowed to depend on the fixed parameters displayed in the subscript. The finitely many smaller values of N are absorbed into the implied constants.

1.4. Organization and Disclosure

Section 2 fixes the conventions and the three external estimates. Section 3 proves the finite-avoidance reduction and the greedy construction. Section 4 derives the gap algebra and the root-counting estimates. Section 5 proves the small-gap bound, including all local and zero-value branches. Section 6 proves the large-gap bound and checks every hypothesis of the fixed-level theorem. Section 7 completes the density and complement arguments. The appendices give an alternative squarefree-kernel organization, an exact source-notation dictionary, and the statement and audit boundary of the companion formal certificate.
Generative models were used to explore proof strategies, translate finite algebra into Lean, and search for counterexamples to intermediate claims. Model agreement was never treated as verification. The authors remain responsible for the mathematical statements; Appendix D records the role and limits of the computational assistance.

2. Preliminaries and External Inputs

All intervals of integers are half open. The symbols O k , ε ( · ) and ≪ k , ε allow constants depending only on the displayed subscripts. For an odd prime p, write
χ p ( u ) = u p , e p ( u ) = exp ( 2 π i u / p ) .
If a ≥ 1 , then ω ( a ) denotes its number of distinct prime factors and d 3 ( a ) its ternary divisor function.
The following three results are the only non-elementary analytic inputs in the source-level proof. We state the precise consequences used later so that the range, normalization, and exceptional sets can be checked locally.
Lemma 1
(Heath–Brown square sieve). Let w ( m ) ≥ 0 be finitely supported, let w ( 0 ) = 0 , and let P be a set of R odd primes. Suppose w ( m ) = 0 whenever | m | ≥ exp ( R ) . Then
∑ j ≥ 1 w ( j 2 ) ≪ 1 R ∑ m w ( m ) + 1 R 2 ∑ p , q ∈ P p ≠ q ∑ m w ( m ) m p q .
Proof. 
This is Theorem 1 of [5], with R denoting the number of sieve primes. The support and zero-value hypotheses are verified explicitly in Section 5.    □
Lemma 2
(Mixed cubic character sum). There is an absolute constant C such that, for every odd prime p, every b ∈ F p × , and every c ∈ F p ,
∑ u ∈ F p χ p ( u 3 + b ) e p ( c u ) ≤ C p .
If b = 0 , the same bound holds by the classical quadratic Gauss-sum formula.
Proof. 
For b ≠ 0 this is the degree-three specialization of [1]. When b = 0 , one has χ p ( u 3 ) = χ p ( u ) ; the zero frequency vanishes and every nonzero frequency is a quadratic Gauss sum of absolute value p .    □
Theorem 3
(Heath–Brown’s diagonal fixed-level estimate). Let d ≥ 3 , let
F ( X 1 , X 2 , X 3 ) = ± X 1 d ± X 2 d ± X 3 d
be nonsingular, and let R be a positive integer satisfying R ≪ F B 3 / 13 . For every ε > 0 , the number of integer points X in the height shell
F ( X ) = R , B / 2 < max i | X i | ≤ B ,
apart from the solutions for which one of the three displayed signed monomials already equals R, is
O F , ε B 9 / 10 + ε R 1 / 10 .
Proof. 
This is the final diagonal conclusion of [6]. That conclusion counts all integral points in the stated shell outside the special set; it does not impose a primitivity condition. Section 6 checks nonsingularity, the level range, the height substitution, the shell summation, and the absence of special points in the present application.    □

3. Finite Avoidance and the Greedy Construction

We first isolate the purely combinatorial part of the proof.
Definition 1.
A set S ⊆ Z has the finite simultaneous escape property if, for every finite set K ⊆ Z ∖ S , there is s ∈ S such that
k − s ∉ S − S ( k ∈ K ) .
Proposition 1.
If S = − S and S has the finite simultaneous escape property, then there is A ⊆ Z such that Z = A ⊕ S .
Proof. 
Choose a sequence z 0 , z 1 , … containing every integer. We construct finite sets B j such that distinct elements of B j have difference outside S − S . Start with B 0 = ⌀ . If z j ∈ B j + S , put B j + 1 = B j . Otherwise, for b ∈ B j put k b = z j − b . Every k b lies outside S, since otherwise z j ∈ B j + S . By finite simultaneous escape, choose s ∈ S such that k b − s ∉ S − S for every b ∈ B j , and put
b new = z j − s , B j + 1 = B j ∪ { b new } .
Then z j = b new + s , and b new − b = k b − s ∉ S − S . Since S − S is symmetric, the reverse difference is safe as well.
Let A = ⋃ j B j . Every integer is covered. If a + s = b + t with a , b ∈ A and s , t ∈ S , then a − b = t − s ∈ S − S . The construction forces a = b , and then s = t . Thus the representation is unique.    □
For cubes, C = − C and D = − D . The density theorem immediately supplies the escape hypothesis. 2
Corollary 1.
The set C has the finite simultaneous escape property.
Proof. 
Let K ⊆ Z ∖ C be finite. By 2, choose any 0 < δ < 1 / 19 . Then
# ⋃ k ∈ K E k ( N ) ≤ ∑ k ∈ K # E k ( N ) = O K , δ ( N 18 / 19 + δ ) = o K ( N ) .
For all sufficiently large N, some n ∈ [ N , 2 N ) lies in none of these sets. Hence n 3 − k ∉ D for every k ∈ K . Since D = − D , this is equivalent to k − n 3 ∉ D . Take s = n 3 .    □
Theorem 1 follows from Proposition 1 and Corollary 1. It remains to prove 2.
We also include the short quadratic obstruction, in order to make the least-degree assertion independent of the formulation on the problem page. 3
Proposition 2.
If f ∈ Q [ X ] is integer-valued on Z and has degree two, then f ( Z ) has no unique additive complement in Z .
Proof. 
There are integers A , B , C , with A ≠ 0 , such that
f ( n ) = A n 2 + B n + C .
If 2 B − A ≠ 0 , then
f ( n ) − f ( − n ) = ( 2 B − A ) n ,
so q Z ⊆ f ( Z ) − f ( Z ) for q = | 2 B − A | . If 2 B − A = 0 , then
f ( n + 1 ) − f ( n − 1 ) = 2 A n ,
and the same conclusion holds with q = 2 | A | .
In either case, suppose for contradiction that Z = T ⊕ f ( Z ) . If two distinct elements of T were congruent modulo q, their difference would belong to q Z ⊆ f ( Z ) − f ( Z ) , producing two representations of one integer. Thus T has at most q elements. But a finite translate set added to the one-sided bounded quadratic set f ( Z ) is still bounded on one side, and cannot equal Z .    □
Combining Theorem 1 and Proposition 2, degree three is the least possible degree in the Erdos–Graham problem.

4. Gap Algebra

Fix a nonzero noncube k. For sufficiently large N, every n ∈ [ N , 2 N ) satisfies n 3 − k > 0 . Suppose
n 3 − k = y 3 − x 3 .
Then y > x . Put
a = y − x > 0 , v = x + y .
A direct calculation gives
4 ( n 3 − k ) = a ( a 2 + 3 v 2 )
and
12 ( n 3 − k ) a − 3 a 2 = ( 3 v ) 2 .
In particular, a ∣ n 3 − k and
a 3 ≤ 4 ( n 3 − k ) ≪ k N 3 , a ≪ k N .
We count gaps in dyadic intervals H ≤ a < 2 H . Let Z k ( N , H ) be the number of distinct n ∈ [ N , 2 N ) admitting such a representation. We shall prove the following uniform estimate. 0
Proposition 3
(Small gaps). For fixed k ≠ 0 and η > 0 ,
Z k ( N , H ) ≪ k , η N H log 3 ( 2 N )
uniformly for 1 ≤ H ≤ N 1 − η .

5. The Small-Gap Square Sieve

5.1. Cubic Roots Modulo the Gap

Define
ρ k ( a ) = # { r ( mod a ) : r 3 ≡ k ( mod a ) } .
0
Lemma 3.
For fixed k ≠ 0 ,
ρ k ( a ) ≤ | k | 2 3 ω ( a ) .
Consequently,
∑ H ≤ a < 2 H ρ k ( a ) ≪ k H log 2 ( 2 H ) ,
∑ H ≤ a < 2 H ρ k ( a ) a ≪ k log 2 ( 2 H ) .
Proof. 
It is enough to work modulo prime powers. Fix p, put ν = v p ( k ) , and consider p e . If e ≤ ν , then the congruence is r 3 ≡ 0 ( mod p e ) , and hence
ρ k ( p e ) = p e − ⌈ e / 3 ⌉ = p ⌊ 2 e / 3 ⌋ .
Suppose e > ν . A root cannot vanish modulo p e ; writing r = p s u with u a unit, comparison of valuations below e gives 3 s = ν . There are no roots if 3 ∤ ν . If ν = 3 s , division by p ν gives a unit equation
u 3 ≡ k / p ν ( mod p e − ν ) .
The kernel of cubing on ( Z / p m Z ) × has size at most three, also for p = 2 , 3 . Each solution modulo p e − ν has exactly p 2 ν / 3 lifts for u modulo p e − s . Therefore
ρ k ( p e ) ≤ 3 p 2 ν / 3 .
Multiplicativity first gives the sharper real-valued estimate ρ k ( a ) ≤ | k | 2 / 3 3 ω ( a ) ; since | k | ≥ 1 , it implies the displayed integral bound | k | 2 3 ω ( a ) used by the formal proof. Finally, 3 ω ( a ) ≤ d 3 ( a ) , and the elementary mean-value estimate ∑ a ≤ X d 3 ( a ) ≪ X log 2 ( 2 X ) gives (4). The weighted estimate follows immediately on a dyadic interval.    □

5.2. Progressions and Local Fourier Coefficients

For every root class choose 0 ≤ r < a and write n = r + a t . Put
q 0 = r 3 − k a ∈ Z
and
f a , r ( t ) = 12 ( q 0 + 3 r 2 t + 3 r a t 2 + a 2 t 3 ) − 3 a 2 .
Then f a , r ( t ) is the left side of (2). Let I a , r be the integer interval defined by N ≤ r + a t < 2 N . It has length at most N / a + 1 . The total number of candidates before imposing the square condition satisfies
W ≤ ∑ H ≤ a < 2 H ρ k ( a ) ( N / a + 1 ) ≪ k N log 2 ( 2 N ) ,
where we used H ≤ N .
For a prime p ∤ 6 k , write χ p for the quadratic character and define
S p ( h ) = ∑ t mod p χ p ( f a , r ( t ) ) e p ( h t ) .
Lemma 4.
Uniformly in a , r , h and p ∤ 6 k ,
| S p ( h ) | ≪ p .
Proof. 
First suppose p ∤ a . The exact integer identity
a f a , r ( t ) = 12 ( r + a t ) 3 − 3 ( a 3 + 4 k )
gives, modulo p, after u = r + a t ,
S p ( h ) = χ p ( 12 a − 1 ) e p ( − h r a − 1 ) × ∑ u mod p χ p ( u 3 − k − a 3 / 4 ) e p ( h a − 1 u ) .
The result follows from 2. If a 3 + 4 k ≡ 0 ( mod p ) , the polynomial is a scalar multiple of u 3 ; equivalently, this special case reduces directly to the quadratic Gauss sum for χ p ( u ) .
Now suppose p ∣ a . Since r 3 ≡ k ( mod p ) and p ∤ k , we have r ¬ ≡ 0 ( mod p ) . Moreover,
f a , r ( t ) ≡ 12 q 0 + 36 r 2 t ( mod p ) ,
whose linear coefficient is nonzero. For a nonconstant linear polynomial α t + β ,
∑ t mod p χ p ( α t + β ) e p ( h t ) = 0 , h = 0 , e p ( − h β α − 1 ) χ p ( h α − 1 ) τ ( χ p ) , h ≠ 0 .
This again has absolute value at most p .    □

5.3. CRT and interval completion

Take distinct primes p , q ∈ ( Q , 4 Q ] , both coprime to 6 k , and put M = p q . Set
g ( t ) = χ p ( f a , r ( t ) ) χ q ( f a , r ( t ) ) , g ^ M ( h ) = ∑ b mod M g ( b ) e M ( h b ) .
Let q q ¯ ≡ 1 ( mod p ) and p p ¯ ≡ 1 ( mod q ) . The CRT gives the exact factorization
g ^ M ( h ) = S p ( h q ¯ mod p ) S q ( h p ¯ mod q ) .
Thus | g ^ M ( h ) | ≪ p q ≪ Q for every h. This includes all four possibilities for whether p and q divide a.
Lemma 5
(Interval completion). For every integer interval I of length L,
∑ t ∈ I g ( t ) ≪ L Q + Q log ( 2 Q ) .
Proof. 
Split I into complete periods modulo M and one remainder. Each complete period contributes g ^ M ( 0 ) = O ( Q ) , giving O ( L / Q ) . On the remainder use
g ( t ) = 1 M ∑ h mod M g ^ M ( h ) e M ( − h t )
and the standard bound
∑ h mod M ∑ j < ℓ e M ( h j ) ≪ M log ( 2 M ) .
Together with (8), this gives O ( Q log ( 2 M ) ) , as required.    □
Since
f a , r ′ ( t ) = 36 ( r + a t ) 2 > 0
on I a , r , the additional condition f a , r ( t ) > 0 cuts out at most one integer interval. Summing 5 over a and its root classes gives, for the off-diagonal correlation,
C p , q ≪ k N Q + H Q log 3 ( 2 N ) .
Indeed, the terms are bounded by
N Q ∑ a ≍ H ρ k ( a ) a + 1 Q + Q log ( 2 Q ) ∑ a ≍ H ρ k ( a ) ,
and 3 applies.

5.4. Application of the Square Sieve

Choose
Q = N / H , P Q = { p ∈ ( Q , 4 Q ] : p prime , p ∤ 6 k } , R = # P Q .
For H ≤ N 1 − η , we have Q ≥ N η / 2 . No prime number theorem is needed: Chebyshev’s explicit lower bound at 4 Q and upper bound at Q give
π ( 4 Q ) ≥ 4 Q log 2 − log ( 4 Q + 1 ) log ( 4 Q ) , π ( Q ) ≤ ( log 4 + ϵ ) Q log Q
for fixed ϵ > 0 and all sufficiently large Q. The leading constants leave a positive margin; the matching Chebyshev upper bound and deletion of finitely many primes therefore give
R ≍ Q / log Q .
Define w ( m ) to count triples ( a , r , t ) above for which f a , r ( t ) = m > 0 . It is a nonnegative, finitely supported weight with w ( 0 ) = 0 , total mass at most W, and
m ≤ 96 N 3 + 12 | k | .
Since R ≫ k , η N η / 2 / log N , the support condition m < exp ( R ) in 1 holds for all sufficiently large N. Furthermore,
∑ m w ( m ) m p q = C p , q .
Hence 1, (7), and (9) yield
∑ j ≥ 1 w ( j 2 ) ≪ W R + 1 R 2 ∑ p ≠ q | C p , q | ≪ k , η N Q log 3 ( 2 N ) + N Q + H Q log 3 ( 2 N ) ≪ k , η N H log 3 ( 2 N ) .
It remains to restore zero values. The equation f a , r ( t ) = 0 implies
4 ( n 3 − k ) = a 3 .
For fixed a there is at most one integer n, so all zero values contribute O ( H ) ≤ O ( N H ) . Every genuine representation is counted by its actual gap, root class, and progression index. This proves 3.

6. The Large-Gap Estimate

We now count the representations whose actual gap is larger than the small-gap cutoff. This part uses Theorem 3; no determinant-method assertion is left implicit.
For 0 < η < 1 , let L k ( N , N 1 − η ) be the number of distinct n ∈ [ N , 2 N ) for which there exist x , y ∈ Z satisfying
n 3 − k = y 3 − x 3 , a = y − x > N 1 − η .
Lemma 6
(Large-gap height reduction). For fixed k and all sufficiently large N, every triple counted by L k ( N , N 1 − η ) satisfies
| x | , | y | , n ≪ k N 1 + η / 2 .
Proof. 
Since n 3 − k > 0 , the difference y 3 − x 3 is positive and hence a > 0 . From (1),
3 a v 2 ≤ a ( a 2 + 3 v 2 ) = 4 ( n 3 − k ) ≤ 4 ( n 3 + | k | ) .
Consequently
| v | ≤ 4 ( n 3 + | k | ) 3 a 1 / 2 ≪ k N 3 / 2 a − 1 / 2 ≪ k N 1 + η / 2 .
Equation (3) gives a ≪ k N , which is also O k ( N 1 + η / 2 ) . Finally,
x = v − a 2 , y = v + a 2 ,
so the same bound holds for x and y; the bound for n is immediate.    □
Proposition 4
(Large gaps). Fix 0 < η < 1 . For every ε > 0 ,
L k ( N , N 1 − η ) ≪ k , ε N ( 1 + η / 2 ) ( 9 / 10 + ε ) .
Proof. 
Choose B = c k N 1 + η / 2 with c k large enough for Lemma 6. We count the associated integer triples; projecting them to their n-coordinate can only decrease the count.
Suppose first that k > 0 . Then
F ( x , n , y ) = k , F ( X 1 , X 2 , X 3 ) = X 1 3 + X 2 3 − X 3 3 .
The projective gradient 3 ( X 1 2 , X 2 2 , − X 3 2 ) has no common nonzero zero, so F is nonsingular. The level in Theorem 3 is the fixed positive integer R = k . Since B → ∞ , the condition R ≪ F B 3 / 13 holds for all sufficiently large N.
The source-special solutions would satisfy one of
x 3 = k , n 3 = k , − y 3 = k .
Each equality makes k an integer cube, the last because − y 3 = ( − y ) 3 . Thus the special set is empty.
If k < 0 , rewrite the equation as
y 3 − x 3 − n 3 = − k
and use
F ( X 1 , X 2 , X 3 ) = X 1 3 − X 2 3 − X 3 3 , ( X 1 , X 2 , X 3 ) = ( y , x , n ) , R = − k > 0 .
This form is nonsingular by the same gradient calculation. A special solution makes one of y 3 , − x 3 , − n 3 equal to − k , again forcing k to be a cube. Hence no special solution occurs in either sign case.
Partition the box max ( | x | , | n | , | y | ) ≤ B into dyadic shells T / 2 < max ( | x | , | n | , | y | ) ≤ T . Every nonempty shell has T ≥ N , because the second coordinate is n ≥ N . Hence the fixed-level condition | k | ≪ F T 3 / 13 holds on every nonempty shell once N is sufficiently large. Since 9 / 10 + ε > 0 , Theorem 3 and the geometric sum over T = B , B / 2 , B / 4 , … give
L k ( N , N 1 − η ) ≪ k , ε B 9 / 10 + ε | k | 1 / 10 ≪ k , ε N ( 1 + η / 2 ) ( 9 / 10 + ε ) .
There is no primitivity reduction, and the shell sum changes only the implied constant.    □

7. Completion of the Proof

Proof 
(Proof of Theorem 2). Set η = 2 / 19 and G = N 17 / 19 . Cover 1 ≤ a ≤ G by dyadic intervals H ≤ a < 2 H . Each lower endpoint satisfies H ≤ N 17 / 19 = N 1 − η , so Proposition 3 applies. There are O ( log ( 2 N ) ) intervals and
N G = N ( 1 + 17 / 19 ) / 2 = N 18 / 19 .
Thus the representations with a ≤ G contribute at most
O k N 18 / 19 log 4 ( 2 N ) .
For a > G , Proposition 4 gives
O k , ε N ( 20 / 19 ) ( 9 / 10 + ε ) = O k , ε N 18 / 19 + ( 20 / 19 ) ε .
Given δ > 0 , choose ε < 19 δ / 40 . The large-gap exponent is then at most 18 / 19 + δ / 2 , while log 4 ( 2 N ) ≪ δ N δ / 2 for sufficiently large N. Adding the two ranges proves
# E k ( N ) ≪ k , δ N 18 / 19 + δ .
The finitely many smaller values of N are absorbed into the constant.    □
Proof of Theorem 1.
The density theorem implies Corollary 1: for any finite K ⊆ Z ∖ C , a sufficiently large dyadic interval contains an n such that n 3 − k ∉ C − C for all k ∈ K . Since C = − C , this is the finite simultaneous escape property for C. Proposition 1 then gives A ⊆ Z with Z = A ⊕ C .    □

Appendix A. A Squarefree-Kernel Proof of the Small-Gap Estimate

For completeness, we give a second organization of the small-gap argument. It uses the same two analytic inputs, but it avoids the prime-power analysis of ρ k ( a ) .
Set
F a ( n ) = 4 n 3 − a 3 − 4 k .
For a genuine representation, (1) gives
F a ( n ) = 3 a v 2 , v = x + y .
Write uniquely
3 a = d s 2 , d = sf ( 3 a ) squarefree .
Then
F a ( n ) = d ( s v ) 2 .
For fixed d, let
A d = { a ∈ [ H , 2 H ) : sf ( 3 a ) = d } .
Since 3 a = d s 2 ,
# A d ≪ H / d + 1 , d ≤ 6 H .
If d ∣ F a ( r ) and n = r + d u , put
M a , r ( u ) = F a ( r ) d + 12 r 2 u + 12 r d u 2 + 4 d 2 u 3 .
This is exactly F a ( n ) / d , so every genuine representation produces a square value of (A3).
Because d is squarefree, CRT and the elementary root bound for a cubic give
# { r ( mod d ) : F a ( r ) ≡ 0 ( mod d ) } ≤ 3 ω ( d ) .
For clarity, if p ≥ 5 the reduction is a nonzero cubic and has at most three roots; modulo 3 it is the linear polynomial r − a − k ; modulo 2 there are at most two residue classes.
Fix ( d , a , r ) . The u-interval has length at most N / d + 1 . Put
Q d = N / d .
Since d ≤ 6 H ≤ 6 N 1 − η , we have Q d ≫ N η / 2 . Use sieve primes in ( Q d , 4 Q d ] not dividing 6 k .
If p ∤ d , the change t = r + d u gives
M a , r ( u ) ≡ d − 1 ( 4 t 3 − a 3 − 4 k ) ( mod p ) .
Lemma 2 gives the required O ( p ) bound for every additive twist; when a 3 + 4 k ≡ 0 ( mod p ) , the polynomial is pure cubic and the same statement is the elementary Gauss bound. If p ∣ d , then p ≥ 5 , p ∣ a , and d ∣ F a ( r ) gives r 3 ≡ k ( mod p ) . Thus r ¬ ≡ 0 ( mod p ) and
M a , r ( u ) ≡ F a ( r ) / d + 12 r 2 u ( mod p ) ,
a nonconstant linear polynomial. The linear Gauss bound applies.
On the interval under consideration, r + d u = n ≥ N , and hence
d d u M a , r ( u ) = 12 ( r + d u ) 2 > 0 .
Thus the additional condition M a , r ( u ) > 0 again cuts out a single integer interval.
The CRT completion of 5, followed by 1, now gives
# { u : M a , r ( u ) > 0 is a square } ≪ k , η Q d log ( 2 N ) .
Here are the details of this last use of the sieve. Let L ≤ N / d + 1 be the length of the u-interval, let P d be the allowed primes in ( Q d , 4 Q d ] , and put R d = # P d . Then
L ≪ Q d 2 + 1 , R d ≍ Q d / log Q d .
The total mass of the progression weight is at most L. For distinct p , q ∈ P d , the same CRT factorization and completion as before give
∑ u ∈ I M a , r ( u ) p q ≪ L Q d + Q d log ( 2 Q d ) .
Consequently the first square-sieve term is
L R d ≪ Q d log ( 2 N ) ,
and the average off-diagonal term has the same order. Moreover M a , r ( u ) = O k ( N 3 ) , so the support condition | M a , r ( u ) | < exp ( R d ) holds for sufficiently large N, uniformly in d ≤ 6 H . This proves (A5). Zero values contribute at most one n for each a.
Combining (A2), (A4), and (A5), the positive-square contribution at fixed d is
≪ k , η 3 ω ( d ) H / d + 1 N / d log ( 2 N ) .
Finally,
∑ d ≤ 6 H 3 ω ( d ) d ≪ log 3 ( 2 H ) , ∑ d ≤ 6 H 3 ω ( d ) d ≪ H log 2 ( 2 H ) .
Summing over d proves
Z k ( N , H ) ≪ k , η N H log 4 ( 2 N ) ,
again sufficient for the density theorem.

Appendix B. Normalization of the external inputs

The proof uses three published estimates with overlapping notation. We record the substitutions explicitly so that their hypotheses can be compared with the quantities defined in the body of the paper. The companion Lean certificate follows a separate route and expands the required finite, character-sum, determinant, height, and counting steps.

Appendix B.1. The Square Sieve

In Theorem 1 of [5], the letter P denotes the number of sieve primes. In this paper that quantity is R, while Q denotes their size. The full translation is
Heath-Brown (1984) Present paper
prime set P P Q = { p ∈ ( Q , 4 Q ] : p ∤ 6 k }
cardinality P R = # P Q ≍ Q / log Q
weight w ( n ) ≥ 0 multiplicity of positive progression values
w ( 0 ) = 0 zero values are separated before the sieve
| n | < exp ( P ) m ≤ 96 N 3 + 12 | k | < exp ( R )
square count ∑ j ≥ 1 w ( j 2 )
Jacobi correlation the completed sum C p , q
The condition involving exp ( P ) is the reason the small-gap range is kept a fixed power below N: it makes R ≫ N η / 2 / log N .

Appendix B.2. The Mixed Character Sum

In the branch p ∤ a , the local factor entering (8) is reduced by u = r + a t to
∑ u mod p χ p ( u 3 + b ) e p ( c u ) , b = − k − a 3 / 4 , c = h a − 1 .
Thus Lemma 2.3 of [1] is used with odd degree 3, linear coefficient 0, constant coefficient b, and additive frequency c. The factor
χ p ( 12 a − 1 ) e p ( − h r a − 1 )
has modulus one. When b = 0 , the cubic is a scalar multiple of u 3 and the estimate is the elementary quadratic Gauss bound. The branch p ∣ a does not use the cited cubic lemma; it is the nonconstant-linear calculation displayed in the proof of Lemma 4.

Appendix B.3. The Fixed-Level Ternary-Form Estimate

The variables in Theorem 1 of [6] are matched as follows.
Heath-Brown (2009) Present paper
degree d = 3
level N in the source R = | k | (fixed)
height B ≪ k N 1 + η / 2 , where N is our dyadic scale
range fixed R ≪ F B 3 / 13 for sufficiently large N
diagonal form X 1 3 + X 2 3 − X 3 3 or X 1 3 − X 2 3 − X 3 3
special solutions one signed cubic monomial already equals R
height-shell estimate B 9 / 10 + ε R 1 / 10
The source uses N for the level, whereas this paper uses N for the size of the cube index. The separate letter R prevents these roles from being confused. Since k is not a cube, every special solution is excluded in the application.

Appendix C. The Companion Formal Certificate

The companion Lean 4 development verifies the exact theorem
theorem erdos477_main :
    exists A : Set Int, IsUniqueAdditiveComplement A cubeValues.
Here IsUniqueAdditiveComplement says directly that every integer has exactly one pair of summands in A × C , and cubeValues is the image of m ↦ m 3 on all of Z .
The checked chain includes the greedy construction, the density-to-escape bridge, the small-gap square sieve, root averages, prime supply, determinant vanishing, common conics, conic charts, sextic factor routing, prime-power removal, height reconstruction, orbit-field transfer, smooth and singular active-cell counts, and the final density bound. This is an independent formal route: unlike Section 6, it does not introduce Theorem 3 as a formal axiom. Its quantitative endpoint is the coarser but still sublinear estimate O k ( N 99 / 100 ) .
The companion ZIP archive, rather than the mathematical body of this paper, contains the exact toolchain pin, dependency manifest, complete Lean sources, build configuration, source map, and checksums. Its default build checks the density theorem, the main theorem, and the axiom audit. The audit applies assert_no_sorry to the main theorem and reports the axiom closure
{ propext , Classical . choice , Quot . sound } .
There is no sorryAx, project-specific mathematical axiom, or synonymous analytic hypothesis in the final entrypoint. The README in the archive gives the one-command build and checksum procedures. Keeping these engineering details in the companion archive makes the paper readable while leaving the formal certificate fully reproducible.

Appendix D. Use of Generative AI and Review Status

The project used generative models in three roles. First, they proposed candidate reductions and searched for counterexamples to intermediate statements. Second, they translated elementary identities, finite sums, and asymptotic bookkeeping into Lean. Third, isolated model runs acted as adversarial readers of proposed arguments. These activities produced both useful lemmas and false starts.
Model agreement was never used as mathematical verification. A claim was retained only after it was supported by an explicit calculation, a theorem whose source and hypotheses were checked, a Lean kernel proof, or a finite computation used solely for falsification. In particular, provider failures, truncated responses, rejected routes, and counterexamples to earlier drafts do not count as evidence for the theorem.
The proof printed in Section 2–7 is designed for ordinary mathematical review and does not require access to model transcripts. The companion Lean development supplies a second, mechanically checked route, but kernel checking does not replace expert assessment of exposition, significance, or novelty. At the time of this version, the manuscript has not undergone independent human peer review. Final responsibility for every statement remains with the authors.

References

  1. Bonolis, D.; Browning, T. Uniform bounds for rational points on hyperelliptic fibrations. Ann. Sc. Norm. Super. Pisa Cl. Sci. 2023, arXiv:2007.1418224, 173–204. [Google Scholar] [CrossRef]
  2. Ding, Y.; Krause, B.; Sandor, C.; Sun, Y.-C.; Zhang, Z. Cross representations of additive complements of r-th powers. arXiv arXiv:2512.15407.
  3. Erdos Problems, Problem 477. Available online: https://www.erdosproblems.com/477.
  4. Erdos, P.; Graham, R. L. Old and New Problems and Results in Combinatorial Number Theory . In Monographies de L’Enseignement Mathematique; Universite de Geneve: Geneva, 1980; vol. 28, p. 95. [Google Scholar]
  5. Heath-Brown, D. R. The square sieve and consecutive square-free numbers. Math. Ann. 1984, 266, 251–259. [Google Scholar] [CrossRef]
  6. D. R. Heath-Brown, Sums and differences of three k-th powers. J. Number Theory 2009, arXiv:0806.4330129, 1579–1594. [CrossRef]
  7. L. Price, Proof claim for Erdos Problem 477, July 27, 2026. Available online: https://www.erdosproblems.com/forum/thread/477/proof-claims.
  8. Sekanina, M. Notes on the factorisation of infinite cyclic groups. Czechoslov. Math. J. 1959, 9, 485–495. [Google Scholar] [CrossRef]
  9. Anonymous, Erdos problem 477: a tiling complement for thirteenth powers, public manuscript. June 2026. Available online: https://github.com/Pengbinghui/pipeline-math/blob/main/papers/tiling-complement.pdf.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.