login
A152313
Primes without 0's or primes in their decimal expansion.
8
11, 19, 41, 61, 89, 149, 181, 191, 199, 419, 449, 461, 491, 499, 619, 641, 661, 691, 811, 881, 911, 919, 941, 991, 1181, 1481, 1489, 1499, 1619, 1669, 1699, 1811, 1861, 1889, 1949, 1999, 4111, 4441, 4481, 4649, 4691, 4861, 4889, 4919, 4969, 4999
OFFSET
1,1
LINKS
MAPLE
F:= proc(d) local T, R, L, r;
R:= NULL;
T:= combinat:-cartprod([[1, 4, 6, 8, 9]$d]);
while not T[finished] do
L:= T[nextvalue]();
r:= add(L[i]*10^(d-i), i=1..d);
if isprime(r) then R:= R, r fi
od;
R
end proc:
seq(F(d), d=2..5); # Robert Israel, Dec 07 2017
MATHEMATICA
Select[Prime[Range[800]], Complement[IntegerDigits[#], {1, 4, 6, 8, 9}] == {} &] (* Vincenzo Librandi, Oct 25 2016 *)
PROG
(Magma) [p: p in PrimesUpTo(5000) | Set(Intseq(p)) subset [1, 4, 6, 8, 9]]; // Vincenzo Librandi, Oct 25 2016
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Omar E. Pol, Dec 02 2008
STATUS
approved