|
%I #16 Nov 21 2019 15:02:47
%S 2,3,5,7,13,17,31,71,113,131,151,211,311,1117,1151,1171,1511,2111,
%T 11113,11117,11131,11171,11311,111121,111211,112111,113111,131111,
%U 311111,511111,1111151,1111211,1111711,1117111,1171111,11111117,11111131,11111171,11111311,11113111,11131111
%N Multiplicative primes: product of digits is a prime.
%C Primes with one prime digit and all other digits are 1. The linked table includes probable primes. - _Jens Kruse Andersen_, Jul 21 2014
%H Jens Kruse Andersen, <a href="/A046703/b046703.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Prime[Range[740000]],PrimeQ[Times@@IntegerDigits[#]]&] (* _Harvey P. Dale_, Oct 02 2011 *)
%t Select[FromDigits/@Flatten[Table[Permutations[PadRight[{p},n,1]],{n,8},{p,{2,3,5,7}}],2],PrimeQ]//Union (* _Harvey P. Dale_, Nov 21 2019 *)
%o (PARI) f(n,b,d) = if(d, f(10*n+1, b, d-1); if(!b, forprime(q=2, 9, f(10*n+q, 1, d-1))), if(b && isprime(n), print1(n", ")))
%o for(d=1, 8, f(0,0,d)) \\ f(0,0,d) prints d-digit terms. _Jens Kruse Andersen_, Jul 21 2014
%o (PARI) From _M. F. Hasler_, Apr 23 2019: (Start)
%o select( is_A046703(n)=isprime(vecprod(digits(n)))&&ispseudoprime(n), [0..9999]) \\ This defines is_A046703(). In older PARI versions, vecprod=factorback.
%o next_A046703(n)={if( n>1, until( ispseudoprime(n), my(d=digits(n)); n=fromdigits( apply( t->if(t>1, nextprime(t+1), 1), d))+(d[1]>5)); n, 2)}
%o A046703_vec(N=99)=vector(N, i, t=next_A046703(if(i>1, t))) \\ (End)
%Y Cf. A117835 ("noncomposite" variant), A007954 (product of digits), A028842 (product of digits is prime).
%K base,nonn
%O 1,1
%A _Felice Russo_
%E Corrected by _Harvey P. Dale_, Oct 02 2011
|
