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.

Showing 1-2 of 2 results.

A283017 Primes which are the sum of three nonzero 6th powers.

Original entry on oeis.org

3, 857, 1459, 4889, 50753, 51481, 66377, 119107, 210961, 262937, 308801, 525017, 531569, 539633, 562691, 766739, 797681, 840241, 1000793, 1046657, 1078507, 1772291, 1864873, 2303003, 2834443, 2986777, 3032641, 3107729, 3365777, 4757609, 4804201, 5135609, 5987593, 7530329, 7534361, 7743529, 8061041
Offset: 1

Views

Author

Ilya Gutkovskiy, Feb 26 2017

Keywords

Comments

Primes of form x^6 + y^6 + z^6 where x, y, z > 0.

Examples

			3 = 1^6 + 1^6 + 1^6;
857 = 2^6 + 2^6 + 3^6;
1459 = 1^6 + 3^6 + 3^6, etc.
		

Crossrefs

Programs

  • Maple
    N:= 10^8: # to get all terms <= N
    S:= [seq(i^6, i=1..floor(N^(1/6)))]:
    S3:= {seq(seq(seq(S[i]+S[j]+S[k],k=1..j),j=1..i),i=1..nops(S))}:
    sort(convert(select(t -> t <= N and isprime(t), S3), list)); # Robert Israel, Mar 09 2017
  • Mathematica
    nn = 15; Select[Union[Plus @@@ (Tuples[Range[nn], {3}]^6)], # <= nn^6 && PrimeQ[#] &]
  • PARI
    list(lim)=my(v=List(),a6,a6b6,t); lim\=1; for(a=1,sqrtnint(lim-2,6), a6=a^6; for(b=1,min(sqrtnint(lim-a6-1,6),a), a6b6=a6+b^6; forstep(c=if(a6b6%2,2,1),min(sqrtnint(lim-a6b6,6),b),2, if(isprime(t=a6b6+c^6), listput(v,t))))); Set(v) \\ Charles R Greathouse IV, Mar 09 2017

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

Original entry on oeis.org

3, 6563, 72353, 137633, 787811, 1745153, 7444673, 44726593, 49202147, 61503553, 86093443, 91858243, 100006817, 100072097, 101686177, 107444417, 143046977, 200006561, 214756067, 257412163, 300452323, 430372577, 431661313, 435812033, 447149537, 452523713, 489805633, 530372321, 744340577
Offset: 1

Views

Author

Ilya Gutkovskiy, Feb 26 2017

Keywords

Comments

Primes of form x^8 + y^8 + z^8 where x, y, z > 0.

Examples

			3 = 1^8 + 1^8 + 1^8;
6563 = 1^8 + 1^8 + 3^8;
72353 = 2^8 + 3^8 + 4^8, etc.
		

Crossrefs

Programs

  • Mathematica
    nn = 13; Select[Union[Plus @@@ (Tuples[Range[nn], {3}]^8)], # <= nn^8 && PrimeQ[#] &]
  • 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
Showing 1-2 of 2 results.