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.
For further motivation and cross-references, see sequence A235265 which is the main entry for this whole family of sequences.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
E.g., 89 = 1121_4 and 1121_3 = 43 both are prime.
MATHEMATICA
b4b3Q[n_]:=Module[{b4=IntegerDigits[n, 4]}, Max[b4]<3&&PrimeQ[ FromDigits[ b4, 3]]]; Select[Prime[Range[700]], b4b3Q] (* Harvey P. Dale, Dec 14 2021 *)
PROG
(PARI) is(p, b=3, c=4)=vecmax(d=digits(p, c))<b&&isprime(vector(#d, i, b^(#d-i))*d~)&&isprime(p)
(PARI) forprime(p=1, 1e3, is(p, 4, 3)&&print1(vector(#d=digits(p, 3), i, 4^(#d-i))*d~, ", ")) \\ To produce the terms, this is more efficient than to select them using straightforwardly is(.)=is(., 3, 4)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
M. F. Hasler, Jan 12 2014
STATUS
approved