%I #10 May 13 2017 23:29:51
%S 0,1,3,6,7,8,9,13,15,16,19,22,24,27,33,35,37,40,46,48,54,55,61,63,72,
%T 73,78,80,81,88,91,94,97,99,112,115,117,118,120,135,141,142,143,144,
%U 157,160,163,166,168,169,171,175,187,193,195,198,208,214,216,222,224
%N Numbers of the form x^2-3^y >= 0.
%t max = 300; Clear[f]; f[m_] := f[m] = Select[Table[x^2 - 3^y, {y, 0, m}, {x, Floor[3^(y/2)], Ceiling[Sqrt[3^y + max]]}] // Flatten // Union, 0 <= # <= max &]; f[1]; f[m = 2]; While[f[m] != f[m - 1], m++]; Print["m = ", m]; A051205 = f[m] (* _Jean-François Alcover_, May 13 2017 *)
%K nonn
%O 1,3
%A _David W. Wilson_
|