|
%I #20 Sep 08 2022 08:44:47
%S 11,4931,21991,23021,37691,44531,46021,57601,62081,95621,124171,
%T 135731,142111,189271,191251,195511,216571,222841,238291,251611,
%U 269431,319211,368021,391031,413411,427711,468581,476591,625831,638431,682951,701671
%N Primes that remain prime through 5 iterations of function f(x) = 2x + 9.
%C Primes p such that 2*p+9, 4*p+27, 8*p+63, 16*p+135 and 32*p+279 are also primes. - _Vincenzo Librandi_, Aug 04 2010
%H Harvey P. Dale, <a href="/A023334/b023334.txt">Table of n, a(n) for n = 1..500</a>
%t it5Q[n_]:=AllTrue[Rest[NestList[2#+9&,n,5]],PrimeQ]; Select[Prime[ Range[ 60000]],it5Q] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Dec 01 2015 *)
%o (Magma) [n: n in [1..5000000] | IsPrime(n) and IsPrime(2*n+9) and IsPrime(4*n+27) and IsPrime(8*n+63) and IsPrime(16*n+135) and IsPrime(32*n+279)] // _Vincenzo Librandi_, Aug 04 2010
%K nonn
%O 1,1
%A _David W. Wilson_
|
