%I #22 Sep 08 2022 08:44:47
%S 7,13,43,223,757,883,2143,2269,2521,2917,3253,3967,4297,4423,6229,
%T 6427,6679,7621,7741,10429,11353,12739,14593,14929,15259,15541,15937,
%U 16981,17449,17911,21187,23899,24697,25633,27481,28789,29611,29833,30253,30703
%N Primes that remain prime through 3 iterations of function f(x) = 5x + 8.
%C Primes p such that 5*p+8, 25*p+48 and 125*p+248 are also primes. - _Vincenzo Librandi_, Aug 04 2010
%H John Cerkan, <a href="/A023286/b023286.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) == 1 (mod 6). - _John Cerkan_, Sep 21 2016
%t it3Q[n_]:=AllTrue[Rest[NestList[5#+8&,n,3]],PrimeQ]; Select[ Prime[ Range[ 3500]],it3Q] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Jul 21 2015 *)
%o (Magma) [n: n in [1..150000] | IsPrime(n) and IsPrime(5*n+8) and IsPrime(25*n+48) and IsPrime(125*n+248)] // _Vincenzo Librandi_, Aug 04 2010
%Y Subsequence of A023220, A023255, and A111225.
%K nonn
%O 1,1
%A _David W. Wilson_
|