A023319 Primes that remain prime through 4 iterations of function f(x) = 8x + 1.

2207, 126107, 230567, 366347, 590207, 788027, 831167, 919067, 985937, 1111427, 1154567, 1170857, 1367507, 1463597, 1624757, 1934417, 2034797, 2468027, 2502767, 2545307, 3019787, 3287027, 3385817, 3675197, 3692207, 4087757, 4565927, 4600577

Offset: 1

Author

David W. Wilson

Keywords

• nonn

Comments

Primes p such that 8*p+1, 64*p+9, 512*p+73 and 4096*p+585 are also primes. - Vincenzo Librandi, Aug 04 2010

Links

• Harvey P. Dale, Table of n, a(n) for n = 1..1000

Programs

Magma

[n: n in [1..5000000] | IsPrime(n) and IsPrime(8*n+1) and IsPrime(64*n+9) and IsPrime(512*n+73) and IsPrime(4096*n+585)] // Vincenzo Librandi, Aug 04 2010

Mathematica

rp4Q[n_]:=AllTrue[Rest[NestList[8#+1&,n,4]],PrimeQ]; Select[Prime[ Range[ 325000]],rp4Q] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Dec 23 2014 *)