%I #13 Sep 08 2022 08:46:00
%S 3,43,433,443,463,643,3343,3433,3463,3643,4363,4463,4643,4663,6343,
%T 33343,36343,36433,36643,43633,44633,46633,46643,46663,63443,63463,
%U 64333,64433,64633,64663,66343,66463,66643,333433,334333,334363,334643,336463,336643,343333,343433,344363,346433,363343,363463,364333,364433,364643
%N Primes having only the (decimal) digits 3, 4 and 6.
%C All terms end in 3 and have a number of digits '4' that is not divisible by 3.
%C A020461 is a subsequence. - _Vincenzo Librandi_, Jul 29 2015
%H Harvey P. Dale, <a href="/A199346/b199346.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Flatten[Table[FromDigits/@(Flatten[{#,3},1]&/@Tuples[{3,4,6},n]),{n,0,5}]],PrimeQ] (* _Harvey P. Dale_, Jan 01 2013 *)
%t Select[Prime[Range[10^5]], Complement[IntegerDigits[#], {3, 4, 6}]=={}&] (* _Vincenzo Librandi_, Jul 28 2015 *)
%o (PARI) a(n, list=0, L=[3, 4, 6], 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=vector(d, i, L[v[i]])*u) || next; reqpal & !isprime(A004086(t)) & next; list & print1(t", "); n--|return(t)))}
%o (Magma) [p: p in PrimesUpTo(4*10^5) | Set(Intseq(p)) subset [3, 4, 6]]; // _Vincenzo Librandi_, Jul 29 2015
%Y Cf. A020449 - A020472, A199325 - A199329.
%Y Cf. similar sequences listed in A199340.
%K nonn,easy,base
%O 1,1
%A _M. F. Hasler_, Nov 05 2011