cp's OEIS Frontend

A004999 Sums of two nonnegative cubes.

%I A004999 #41 Jul 02 2024 07:23:35
%S A004999 0,1,2,8,9,16,27,28,35,54,64,65,72,91,125,126,128,133,152,189,216,217,
%T A004999 224,243,250,280,341,343,344,351,370,407,432,468,512,513,520,539,559,
%U A004999 576,637,686,728,729,730,737,756,793,854,855,945,1000,1001
%N A004999 Sums of two nonnegative cubes.
%H A004999 T. D. Noe, <a href="/A004999/b004999.txt">Table of n, a(n) for n = 1..1000</a>
%H A004999 Kevin A. Broughan, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL6/Broughan/broughan25.html">Characterizing the sum of two cubes</a>, J. Integer Seqs., Vol. 6, 2003.
%H A004999 Samuel S. Wagstaff, Jr., <a href="https://cs.uwaterloo.ca/journals/JIS/VOL25/Wagstaff/wagstaff8.html">Equal Sums of Two Distinct Like Powers</a>, J. Int. Seq., Vol. 25 (2022), Article 22.3.1.
%H A004999 <a href="/index/Su#ssq">Index entries for sequences related to sums of cubes</a>
%t A004999 Union[(#[[1]]^3+#[[2]]^3)&/@Tuples[Range[0,20],{2}]] (* _Harvey P. Dale_, Dec 04 2010 *)
%o A004999 (PARI) is(n)=my(k1=ceil((n-1/2)^(1/3)), k2=floor((4*n+1/2)^(1/3)), L); fordiv(n,d,if(d>=k1 && d<=k2 && denominator(L=(d^2-n/d)/3)==1 && issquare(d^2-4*L), return(1))); 0
%o A004999 list(lim)=my(v=List());for(x=0,(lim+.5)^(1/3),for(y=0,min(x,(lim-x^3)^(1/3)),listput(v,x^3+y^3))); vecsort(Vec(v),,8) \\ _Charles R Greathouse IV_, Jun 12 2012
%o A004999 (PARI) is(n)=my(L=sqrtnint(n-1,3)+1,U=sqrtnint(4*n,3));fordiv(n,m,if(L<=m&m<=U,my(ell=(m^2-n/m)/3);if(denominator(ell)==1&&issquare(m^2-4*ell),return(1))));0 \\ _Charles R Greathouse IV_, Apr 16 2013
%o A004999 (PARI) T=thueinit('z^3+1);
%o A004999 is(n)=n==0 || #select(v->min(v[1],v[2])>=0, thue(T,n))>0 \\ _Charles R Greathouse IV_, Nov 29 2014
%o A004999 (Haskell)
%o A004999 a004999 n = a004999_list !! (n-1)
%o A004999 a004999_list = filter c2 [1..] where
%o A004999    c2 x = any (== 1) $ map (a010057 . fromInteger) $
%o A004999                        takeWhile (>= 0) $ map (x -) $ tail a000578_list
%o A004999 -- _Reinhard Zumkeller_, Dec 20 2013
%Y A004999 Subsequence of A045980; A003325 is a subsequence.
%Y A004999 Cf. A000578, A004825, A010057, A373972 (characteristic function).
%Y A004999 Indices of nonzero terms in A025446.
%K A004999 nonn,easy,nice
%O A004999 1,3
%A A004999 _N. J. A. Sloane_, _Steven Finch_