cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A023326 Primes that remain prime through 4 iterations of the function f(x) = 9x + 8.

This page as a plain text file.
%I A023326 #21 Sep 08 2022 08:44:47
%S A023326 397,467,907,1747,1901,4099,7237,8117,8581,9371,9587,23539,28081,
%T A023326 32611,36899,41729,54767,55207,57601,61991,64997,66449,73061,74821,
%U A023326 75527,86291,115021,118717,125659,126067,134287,140677,147011,148147,151531,171539
%N A023326 Primes that remain prime through 4 iterations of the function f(x) = 9x + 8.
%C A023326 Primes p such that 9*p+8, 81*p+80, 729*p+728 and 6561*p+6560 are also primes. - _Vincenzo Librandi_, Aug 04 2010
%H A023326 John Cerkan, <a href="/A023326/b023326.txt">Table of n, a(n) for n = 1..10000</a>
%t A023326 Select[Prime[Range[16000]],AllTrue[Rest[NestList[9#+8&,#,4]],PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Jun 25 2017 *)
%o A023326 (Magma) [n: n in [1..5000000] | IsPrime(n) and IsPrime(9*n+8) and IsPrime(81*n+80) and IsPrime(729*n+728) and IsPrime(6561*n+6560)] // _Vincenzo Librandi_, Aug 04 2010
%Y A023326 Subsequence of A023235, A023267, and A023298.
%K A023326 nonn
%O A023326 1,1
%A A023326 _David W. Wilson_