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.

A283019 Primes which are the sum of three nonzero 8th powers.

This page as a plain text file.
%I A283019 #11 Nov 05 2017 21:01:28
%S A283019 3,6563,72353,137633,787811,1745153,7444673,44726593,49202147,
%T A283019 61503553,86093443,91858243,100006817,100072097,101686177,107444417,
%U A283019 143046977,200006561,214756067,257412163,300452323,430372577,431661313,435812033,447149537,452523713,489805633,530372321,744340577
%N A283019 Primes which are the sum of three nonzero 8th powers.
%C A283019 Primes of form x^8 + y^8 + z^8 where x, y, z > 0.
%H A283019 Chai Wah Wu, <a href="/A283019/b283019.txt">Table of n, a(n) for n = 1..10000</a>
%e A283019 3 = 1^8 + 1^8 + 1^8;
%e A283019 6563 = 1^8 + 1^8 + 3^8;
%e A283019 72353 = 2^8 + 3^8 + 4^8, etc.
%t A283019 nn = 13; Select[Union[Plus @@@ (Tuples[Range[nn], {3}]^8)], # <= nn^8 && PrimeQ[#] &]
%o A283019 (PARI) list(lim)=my(v=List(),A,B,t); lim\=1; for(a=1,sqrtnint(lim-2,8), A=a^8; for(b=1,min(sqrtnint(lim-A-1,8),a), B=A+b^8; forstep(c=if(B%2,2,1),sqrtnint(lim-B,8),2, if(isprime(t=B+c^8), listput(v,t))))); Set(v) \\ _Charles R Greathouse IV_, Nov 05 2017
%Y A283019 Cf. A001016, A003381, A006686, A007490, A085317, A085318, A085319, A283017, A283018.
%K A283019 nonn
%O A283019 1,1
%A A283019 _Ilya Gutkovskiy_, Feb 26 2017