A178251 Primes p such that p^3 - 2 is prime.

19, 31, 37, 67, 109, 151, 211, 241, 277, 367, 439, 457, 619, 691, 727, 787, 859, 967, 1087, 1171, 1471, 1489, 1531, 1579, 1951, 2131, 2287, 2791, 2851, 2971, 3061, 3319, 3511, 3547, 3559, 3739, 4129, 4357, 4447, 4507, 4591, 4651, 4789, 4801, 4831, 4951

Offset: 1

Author

Ulrich Krug (leuchtfeuer37(AT)gmx.de), May 24 2010

Keywords

• nonn
• easy

Examples

			6857 = prime(882) = 19^3 - 2, 19 = prime(8) is 1st term.
29789 = prime(3228) = 31^3 - 2, 31 = prime(11) is 2nd term.
		

Links

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

Crossrefs

Cf. A000040, A000578, A038599, A038600, A144953.

Programs

Magma

[p: p in PrimesUpTo(5000) | IsPrime(p^3-2)]; // Vincenzo Librandi, Nov 17 2010

Mathematica

Select[Prime[Range[10000]], PrimeQ[#^3 - 2] &] (* Vincenzo Librandi, Mar 20 2014 *)

PARI

list(lim)=my(v=List()); forprime(p=2,lim, if(isprime(p^3-2), listput(v, p))); Vec(v) \\ Charles R Greathouse IV, Feb 08 2016

Sage

a = list(p for p in primes(10000) if is_prime(p**3-2)) # D. S. McNeil, May 25 2010

Extensions

Base tag removed by D. S. McNeil, May 25 2010