OFFSET
1,3
MATHEMATICA
max = 2000; Clear[f]; f[m_] := f[m] = Select[Table[10^x - y^2, {x, 0, m}, {y, 0, Ceiling[10^(x/2)]}] // Flatten // Union, 0 <= # <= max &]; f[1]; f[m = 2]; While[f[m] != f[m - 1], m++]; Print["m = ", m]; A051221 = f[m] (* Jean-François Alcover, May 13 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved