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

Author

Serge Batalov, Sep 03 2009

Keywords

• easy
• nonn

Links

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

Crossrefs

Cf. A153439, A175170.

Programs

Magma

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

Mathematica

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

PARI

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

Extensions

More terms from Vincenzo Librandi, Mar 25 2010