cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A134739 Cubes of (positive numbers that are not the sum of three nonzero squares), that is, the terms of A004214, cubed.

This page as a plain text file.
%I A134739 #28 Jan 18 2016 10:26:45
%S A134739 1,8,64,125,343,512,1000,2197,3375,4096,8000,12167,15625,21952,29791,
%T A134739 32768,50653,59319,64000,103823,140608,166375,195112,216000,250047,
%U A134739 262144,357911,493039,512000,614125,658503,778688,857375,1000000
%N A134739 Cubes of (positive numbers that are not the sum of three nonzero squares), that is, the terms of A004214, cubed.
%C A134739 This sequence was inspired by e-mail from _Ray Chandler_, Nov 07 2007.
%C A134739 Original name was: Cubes which are not the sum of three nonzero squares. That definition would not include 125 = 5^2 + 6^2 + 8^2. - _Robert Israel_, Jan 12 2016
%C A134739 For "(cubes of positive numbers) that are not the sum of three nonzero squares", that is, the cubes in A004214, see A267189. - _N. J. A. Sloane_, Jan 18 2016
%H A134739 Robert Israel, <a href="/A134739/b134739.txt">Table of n, a(n) for n = 1..10000</a>
%F A134739 a(n) = A004214(n)^3. - _Ray Chandler_, Jan 29 2009
%e A134739 8 is in the sequence because it is not possible to express 2 as a sum of three nonzero squares and 2^3 = 8.
%e A134739 27 is not in the sequence because 3 = 1^2 + 1^2 + 1^2.
%p A134739 N:= 1000: # to get all terms <= N^3
%p A134739 A004214:= {$1..N} minus {seq(seq(seq(a^2 + b^2 + c^2, c = b .. floor(sqrt(N-a^2-b^2))), b = a .. floor(sqrt(N-a^2))),a=1..floor(sqrt(N/2)))}:
%p A134739 map(`^`,sort(convert(A004214,list)), 3); # _Robert Israel_, Jan 12 2016
%t A134739 searchMax = 16; Flatten[Position[Take[Rest[CoefficientList[Sum[x^(i^2), {i, searchMax}]^3, x]], searchMax^2], 0]]^3 (* Based on _Ray Chandler_'s program for A004214, _Alonso del Arte_, Jan 12 2016 *)
%o A134739 (PARI) is(n) = { my(a, b) ; a=1; while(a^2+1<n, b=1 ; while(b<=a && a^2+b^2<n, if(issquare(n-a^2-b^2), return(1) ) ; b++ ; ) ; a++ ; ) ; return(0) ; }
%o A134739 for(n=1, 1e3, if(!is(n), print1(n^3, ", "))); \\ _Altug Alkan_, Jan 13 2016
%Y A134739 Cf. A004214, A134738, A267189.
%K A134739 nonn
%O A134739 1,2
%A A134739 _Artur Jasinski_, Nov 07 2007
%E A134739 Definition corrected by _Robert Israel_, Jan 12 2016