%I #16 Sep 08 2022 08:45:39
%S 11,19,41,61,89,149,181,191,199,419,449,461,491,499,619,641,661,691,
%T 811,881,911,919,941,991,1181,1481,1489,1499,1619,1669,1699,1811,1861,
%U 1889,1949,1999,4111,4441,4481,4649,4691,4861,4889,4919,4969,4999
%N Primes without 0's or primes in their decimal expansion.
%H Robert Israel, <a href="/A152313/b152313.txt">Table of n, a(n) for n = 1..10000</a>
%p F:= proc(d) local T, R, L, r;
%p R:= NULL;
%p T:= combinat:-cartprod([[1,4,6,8,9]$d]);
%p while not T[finished] do
%p L:= T[nextvalue]();
%p r:= add(L[i]*10^(d-i),i=1..d);
%p if isprime(r) then R:= R,r fi
%p od;
%p R
%p end proc:
%p seq(F(d),d=2..5); # _Robert Israel_, Dec 07 2017
%t Select[Prime[Range[800]], Complement[IntegerDigits[#], {1, 4, 6, 8, 9}] == {} &] (* _Vincenzo Librandi_, Oct 25 2016 *)
%o (Magma) [p: p in PrimesUpTo(5000) | Set(Intseq(p)) subset [1,4,6,8,9]]; // _Vincenzo Librandi_, Oct 25 2016
%Y Cf. A019546, A034844, A087363, A152312.
%K base,easy,nonn
%O 1,1
%A _Omar E. Pol_, Dec 02 2008
|