%I #15 Nov 06 2012 12:38:05

%S 3,4,5,6,7,10,11,12,13,14,15,17,18,19,20,21,22,23,24,25,26,29,30,31,

%T 32,33,34,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,55,56,

%U 57,58,59,60,61,62,63,66,67,68,69,70,71,73,74,75,76,77

%N Positive integers that are not the sum of two nonnegative cubes.

%C Omits the positive cubes (A000578) since m^3 = 0^3 + m^3, so is different from A057903.

%H T. D. Noe, <a href="/A022555/b022555.txt">Table of n, a(n) for n=1..10000</a>

%H <a href="/index/Su#ssq">Index entries for sequences related to sums of cubes</a>

%p read transforms; cub := series(add(q^(n^3),n=0..100),q,1000001); t1 := series(cub^2,q,2000); t2 := POWERS(t1,q,2000); COMPl(t2);

%t r[n_] := Reduce[x >= 0 && y >= 0 && n == x^3 + y^3, {x, y}, Integers]; Select[ Range[80], r[#] === False &] (* _Jean-François Alcover_, Nov 06 2012 *)

%t Select[Range@100,PowersRepresentations[#,2,3]=={}&] (* A much faster solution given by _Giovanni Resta_, Nov 06 2012 *)

%Y Complement of A004999. Cf. A004825, A022561, A022566, A057903, A000578.

%K nonn

%O 1,1

%A _N. J. A. Sloane_.

%E Edited by _N. J. A. Sloane_, Sep 28 2007