A023281 Primes that remain prime through 3 iterations of function f(x) = 4x + 3.
2, 109, 179, 571, 677, 977, 1279, 1447, 1747, 1901, 2207, 2671, 3119, 3917, 5011, 5399, 5441, 5569, 5791, 6211, 6607, 7079, 7417, 8369, 8831, 9221, 9697, 9769, 11821, 11897, 12347, 13537, 13669, 13691, 13729, 13781, 13907, 14747, 14851, 15581, 17231, 17497
Offset: 1
Keywords
Links
- John Cerkan, Table of n, a(n) for n = 1..10000
Programs
-
Magma
[n: n in [1..150000] | IsPrime(n) and IsPrime(4*n+3) and IsPrime(16*n+15) and IsPrime(64*n+63)] // Vincenzo Librandi, Aug 04 2010
-
Maple
f:=proc(x) options operator, arrow: 4*x+3 end proc: a:=proc(n) if isprime(n)= true and isprime(f(n))=true and isprime(f(f(n)))=true and isprime(f(f(f(n)))) =true then n else end if end proc: seq(a(n),n=1..20000); # Emeric Deutsch, Jan 01 2008
-
Mathematica
Select[Prime@ Range@ 2100, Times @@ Boole@ PrimeQ@ Rest@ NestList[4 # + 3 &, #, 3] > 0 &] (* Michael De Vlieger, Sep 19 2016 *)
-
PARI
is(n)=isprime(n) && isprime(4*n+3) && isprime(16*n+15) && isprime(64*n+63) \\ Charles R Greathouse IV, Sep 20 2016
Comments