A023338 Primes that remain prime through 5 iterations of function f(x) = 3x + 10.
233, 709, 11699, 16417, 24481, 31663, 63317, 94999, 172279, 185243, 205081, 228713, 240347, 266681, 285539, 334511, 372059, 395719, 479327, 583501, 734197, 786613, 856627, 863711, 1146091, 1187167, 1188259, 1423067, 1534591, 1562269, 1594301
Offset: 1
Keywords
Links
- John Cerkan, Table of n, a(n) for n = 1..10000
Programs
-
Magma
[n: n in [1..25000000] | IsPrime(n) and IsPrime(3*n+10) and IsPrime(9*n+40) and IsPrime(27*n+130) and IsPrime(81*n+400) and IsPrime(243*n+1210)] // Vincenzo Librandi, Aug 05 2010
-
Mathematica
rp5Q[p_]:=AllTrue[Rest[NestList[3#+10&,p,5]],PrimeQ]; Select[Prime[ Range[ 150000]],rp5Q] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Mar 17 2020 *)
-
PARI
is(n)=isprime(n) && isprime(3*n+10) && isprime(9*n+40) && isprime(27*n+130) && isprime(81*n+400) && isprime(243*n+1210) \\ Charles R Greathouse IV, Oct 11 2016
Formula
a(n) == 9 (mod 14). - John Cerkan, Oct 11 2016
Extensions
Definition clarified by Harvey P. Dale, Mar 17 2020
Comments