%I #22 Sep 08 2022 08:44:47
%S 103,283,313,929,2549,2593,3931,4091,5279,6781,7759,8363,16229,16493,
%T 16703,19961,24071,25733,27059,32191,34961,35381,36341,44131,50093,
%U 52973,59671,66179,69833,69931,75853,90031,92993,93371,101879,103561,106949
%N Primes that remain prime through 4 iterations of function f(x) = 9x + 2.
%C Primes p such that 9*p+2, 81*p+20, 729*p+182 and 6561*p+1640 are also primes. - _Vincenzo Librandi_, Aug 04 2010
%H Harvey P. Dale, <a href="/A023324/b023324.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Prime[Range[11000]],AllTrue[Rest[NestList[9#+2&,#,4]],PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Feb 12 2015 *)
%o (Magma) [n: n in [1..5000000] | IsPrime(n) and IsPrime(9*n+2) and IsPrime(81*n+20) and IsPrime(729*n+182) and IsPrime(6561*n+1640)] // _Vincenzo Librandi_, Aug 04 2010
%K nonn
%O 1,1
%A _David W. Wilson_
|