%I #23 Sep 08 2022 08:44:47
%S 3181,61981,134291,342821,459091,882451,984931,1016011,1028471,
%T 1181701,1391561,1897801,2009311,2272301,2476421,2769791,3048041,
%U 3085421,3128821,3207221,3545111,4092931,4690591,5015321,5863651,6027941,6361351,6796541
%N Primes that remain prime through 5 iterations of function f(x) = 3x + 8.
%C Primes p such that 3*p+8, 9*p+32, 27*p+104, 81*p+320, and 243*p+968 are also primes. - _Vincenzo Librandi_, Aug 05 2010
%H John Cerkan, <a href="/A023337/b023337.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) == 31 (mod 70). - _John Cerkan_, Oct 11 2016
%t Select[Prime[Range[500000]],And@@PrimeQ[Rest[NestList[3#+8&,#,5]]]&] (* _Harvey P. Dale_, Apr 07 2014 *)
%o (Magma) [n: n in [1..25000000] | IsPrime(n) and IsPrime(3*n+8) and IsPrime(9*n+32) and IsPrime(27*n+104) and IsPrime(81*n+320) and IsPrime(243*n+968)] // _Vincenzo Librandi_, Aug 05 2010
%o (PARI) is(n)=isprime(n) && isprime(3*n+8) && isprime(9*n+32) && isprime(27*n+104) && isprime(81*n+320) && isprime(243*n+968) \\ _Charles R Greathouse IV_, Oct 11 2016
%Y Subsequence of A023210, A023248, A023279, and A023309.
%K nonn
%O 1,1
%A _David W. Wilson_