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

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

Author

Vincenzo Librandi, Dec 03 2010

Keywords

• nonn
• easy

Links

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

Crossrefs

Cf. A060891, A175170, A000040, A164989.

Subsequence of A002476.

Programs

Magma

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

Mathematica

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

PARI

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

Formula

A000040 INTERSECT A164989.