|
%I #17 Feb 20 2023 08:04:57
%S 7,14,15,19,21,22,23,30,40,42,46,47,48,56,59,61,62,67,75,78,85,86,87,
%T 93,94,96,98,104,105,106,110,111,112,115,117,119,120,124,131,138,139,
%U 143,147,157,159,163,166,167,173,176,180,181,183,184,187,194,195,199,200,202
%N Numbers that are not the sum of a square and 2 nonnegative cubes.
%H R. J. Mathar, <a href="/A022557/b022557.txt">Table of n, a(n) for n = 1..9477</a>
%p isA022557 := proc(n)
%p not isA022556(n) ;
%p end proc:
%p n := 1:
%p for c from 0 do
%p if isA022557(c) then
%p printf("%d %d\n",n,c);
%p n := n+1 ;
%p end if;
%p end do: # _R. J. Mathar_, Sep 02 2016
%t smax = 50000;
%t A022556 = Table[s = i^2 + j^3 + k^3; If[s <= smax, s, Nothing], {i, 0, Sqrt[smax] // Ceiling}, {j, 0, smax^(1/3) // Ceiling}, {k, j, smax^(1/3) // Ceiling}] // Flatten // Union;
%t Complement[Range[smax], A022556] (* _Jean-François Alcover_, Feb 20 2023 *)
%Y Cf. A022556 (complement).
%K nonn
%O 1,1
%A _N. J. A. Sloane_.
|
