OFFSET
1,1
COMMENTS
This sequence is part of the two-dimensional array of sequences based on this same idea for any two different bases b, c > 1. Sequence A235265 and A235266 are the most elementary ones in this list. Sequences A089971, A089981 and A090707 through A090721, and sequences A065720 - A065727, follow the same idea with one base equal to 10.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
EXAMPLE
19 is a term: 19 = 21_9 and 21_8 = 17, also a prime.
79 is not a term: 79 = 87_9 and 87 is not a valid base-8 representation.
MATHEMATICA
b9b8pQ[n_]:=Module[{id=IntegerDigits[n, 9]}, Max[id]<8&&PrimeQ[FromDigits[ id, 8]]]; Select[Prime[Range[500]], b9b8pQ] (* Harvey P. Dale, Mar 12 2018 *)
PROG
(PARI) is(p, b=8, c=9)=vecmax(d=digits(p, c))<b&&isprime(vector(#d, i, b^(#d-i))*d~)&&isprime(p)
(PARI) forprime(p=1, 3e3, is(p, 9, 8)&&print1(vector(#d=digits(p, 8), i, 9^(#d-i))*d~, ", ")) \\ To produce the terms, this is more efficient than to select them using straightforwardly is(.)=is(., 8, 9)
(PARI) isok(p) = isprime(p) && (q = digits(p, 9)) && (vecmax(q) < 8) && isprime(fromdigits(q, 8)); \\ Michel Marcus, Mar 12 2018
CROSSREFS
KEYWORD
nonn,base
AUTHOR
M. F. Hasler, Jan 13 2014
STATUS
approved