A163570 Primes p such that p^6-6 is also prime.

5, 7, 17, 37, 113, 137, 157, 173, 223, 227, 337, 487, 503, 617, 743, 797, 857, 863, 877, 947, 1103, 1217, 1277, 1453, 1483, 1663, 1693, 1747, 1933, 2557, 2683, 2833, 2837, 2843, 3067, 3167, 3373, 3457, 3583, 3593, 3823, 4013, 4483, 4943, 4987, 5527, 5843

Offset

1,1

Links

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

Example

5^6-6=15619 = A000040(1821) puts p=5 into the sequence.
7^6-6=117643 = A000040(11098) puts p=7 into the sequence.

Mathematica

lst={}; Do[p=Prime[n]; If[PrimeQ[p^6-6], AppendTo[lst, p]], {n, 2*7!}]; lst
Select[Prime[Range[1000]], PrimeQ[(#^6 - 6)]&] (* Vincenzo Librandi, Apr 15 2013 *)

Prog

(Magma) [p: p in PrimesUpTo(6000) | IsPrime(p^6-6)]; // Vincenzo Librandi, Apr 15 2013
(PARI) for(n=1, 500, m=prime(n); if(isprime(m^6 -6), print1(m, ", "))) \\ G. C. Greubel, Jul 28 2017

Crossrefs

Sequence in context: A243400 A318568 A239414 * A113282 A323200 A096741

Adjacent sequences: A163567 A163568 A163569 * A163571 A163572 A163573

Keyword

nonn,easy

Author

Vladimir Joseph Stephan Orlovsky, Jul 31 2009

Extensions

Comment moved to the examples by R. J. Mathar, Oct 21 2009

Status

approved