|
%I #21 May 02 2023 16:22:46
%S 3,43,349,433,439,443,449,499,3343,3433,3449,3499,3943,4339,4349,4493,
%T 4933,4943,4993,4999,9343,9349,9433,9439,9949,33343,33349,33493,34439,
%U 34499,34939,34949,39343,39439,39443,39499,43399,43499,43933,43943,44449,44939,49333,49339,49393,49433,49499,49939,49943,49993
%N Primes having only the (decimal) digits 3, 4 and 9.
%C A020461 and A020466 are subsequences. - _Vincenzo Librandi_, Jul 30 2015
%H Alois P. Heinz, <a href="/A199349/b199349.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Pri#PrimesWithDigits">Index to entries about primes with digits in a given set</a>.
%t Select[Prime[Range[2 10^4]], Complement[IntegerDigits[#], {3, 4, 9}]=={} &] (* _Vincenzo Librandi_, Jul 30 2015 *)
%t Select[Flatten[Table[FromDigits/@Tuples[{3,4,9},n],{n,5}]],PrimeQ] (* _Harvey P. Dale_, May 02 2023 *)
%o (PARI) a(n, list=0, L=[3,4,9], reqpal=0)={my(t); for(d=1, 1e9, u=vector(d, i, 10^(d-i))~; forvec(v=vector(d, i, [1+(i==1&!L[1]), #L]), isprime(t=vecextract(L,v)*u) || next; reqpal && !isprime(A004086(t)) && next; list && print1(t", "); n--||return(t)))}
%o (Magma) [p: p in PrimesUpTo(2*10^5) | Set(Intseq(p)) subset [3, 4, 9]]; // _Vincenzo Librandi_, Jul 30 2015
%Y Cf. Primes that contain only the digits (3,4,k): A199340 (k=0), A199341 (k=1), A199342 (k=2), A199345 (k=5), A199346 (k=6), A199347 (k=7), A199348 (k=8).
%Y Cf. A020449 - A020472, A199325 - A199329.
%K nonn,base
%O 1,1
%A _M. F. Hasler_, Nov 05 2011
|
