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