OFFSET
1,1
COMMENTS
Primes p such that A054054(p) is in the set {3, 5, 7}. - Felix Fröhlich, Jan 26 2017
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 4: # to get all terms with <+ N digits
for j in [3, 5, 7] do L[1, j]:= {$j..9}; C[1, j]:= {j} od;
for d from 2 to N do
for j in [3, 5, 7] do L[d, j]:= map(t -> $(10*t+j)..(10*t+9), L[d-1, j]) od:
C[d, 3]:= map(t -> (10*t+3), L[d-1, 3]) union map(t ->seq(10*t+j, j=4..9), C[d-1, 3]);
C[d, 5]:= map(t -> (10*t+5), L[d-1, 5]) union map(t -> seq(10*t+j, j=6..9), C[d-1, 5]);
C[d, 7]:= map(t -> (10*t+7), L[d-1, 7]) union map(t -> seq(10*t+j, j=8..9), C[d-1, 7]);
od:
sort(convert(`union`(seq(seq(select(isprime, C[d, j]), j=[3, 5, 7]), d=1..N)), list));
# Robert Israel, Jan 26 2017
PROG
(PARI) is(n) = ispseudoprime(n) && #setintersect(Set(vecmin(digits(n))), [3, 5, 7])==1 \\ Felix Fröhlich, Jan 26 2017
CROSSREFS
KEYWORD
nonn,base,less
AUTHOR
Juri-Stepan Gerasimov, Jan 19 2009
EXTENSIONS
Edited by N. J. A. Sloane, Jan 24 2009
Corrected by Robert Israel, Jan 26 2017
STATUS
approved