A175170 -id:A175170 - cp's OEIS frontend

A164989 Numbers n such that n^6 - n^3 + 1 is prime.

6, 7, 9, 42, 49, 54, 66, 73, 81, 88, 96, 100, 106, 118, 120, 121, 126, 175, 190, 196, 207, 208, 234, 235, 240, 271, 301, 312, 322, 342, 343, 367, 378, 379, 381, 384, 397, 399, 415, 418, 423, 442, 444, 480, 505, 513, 520, 556, 574, 609, 628, 664, 666, 700, 702

Offset: 1

Serge Batalov, Sep 03 2009

G. C. Greubel, Table of n, a(n) for n = 1..10000

Cf. A153439, A175170.

[n: n in [0..1000] |IsPrime(n^3*(n^3-1)+1)] // Vincenzo Librandi, Dec 03 2010

Select[Range[10^3], PrimeQ[#^6 - #^3 + 1] &] (* Michael De Vlieger, Dec 02 2017 *)

for(k=1,300,if(isprime(k^6-k^3+1),print1(k, ", ")))

More terms from Vincenzo Librandi, Mar 25 2010

A175205 Primes p such that p^6-p^3+1 is also prime.

Original entry on oeis.org

7, 73, 271, 367, 379, 397, 739, 823, 859, 1291, 1297, 1627, 1747, 2083, 2203, 2221, 2269, 2311, 3067, 3079, 3109, 3163, 3463, 3607, 3733, 3907, 4093, 4339, 4603, 4933, 5077, 5431, 5527, 5821, 5869, 6091, 6301, 6373, 6421, 6529, 6883, 7927, 7951, 8011, 8467, 8689, 8737

Offset: 1

Vincenzo Librandi, Dec 03 2010

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Cf. A060891, A175170, A000040, A164989.

Subsequence of A002476.

[p: p in PrimesUpTo(9000) | IsPrime(p^6-p^3+1)]; // Vincenzo Librandi, Apr 16 2013

Select[Prime[Range[2000]], PrimeQ[#^6 -  #^3 + 1]&] (* Vincenzo Librandi, Apr 16 2013 *)

forprime(p=2,10000,if(isprime(p^6-p^3+1),print1(p,", ")))

A000040 INTERSECT A164989.

cp's OEIS Frontend

A164989 Numbers n such that n^6 - n^3 + 1 is prime.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Extensions