OFFSET
1,2
COMMENTS
For terms > 10 the quotient (k / product of digits of k) is prime. - David A. Corneth, Mar 30 2021
LINKS
David A. Corneth, Table of n, a(n) for n = 1..12677 (terms <= 10^12; first 448 terms from Klaus Brockhaus)
EXAMPLE
21728 is in the sequence as 21728/(2*1*7*2*8) = 97 which is prime. - David A. Corneth, Mar 30 2021
MATHEMATICA
okQ[n_] := Block[{p = Times @@ IntegerDigits[n]}, n == p || PrimeQ[n/p]]; Select[ Range[36000], okQ]
Select[Range[40000], FreeQ[IntegerDigits[#], 0]&&!CompositeQ[#/Times@@IntegerDigits[#]]&] (* Harvey P. Dale, Mar 30 2024 *)
PROG
(Magma) IsA001103:=func< n | p ne 0 and n mod p eq 0 and (q eq 1 or IsPrime(q)) where q is (p eq 0 select 0 else n div p) where p is &*Intseq(n) >; [ n: n in [1..40000] | IsA001103(n) ]; // Klaus Brockhaus, Jan 24 2011
(Haskell)
a001103 n = a001103_list !! (n-1)
a001103_list = filter f a052382_list where
f x = m == 0 && (x' == 1 || a010051 x' == 1) where
(x', m) = divMod x $ a007954 x
-- Reinhard Zumkeller, Nov 02 2011
(PARI) is(n) = { my(vp = vecprod(digits(n))); if(vp > 0, c = n/vp; if(denominator(c) == 1, if(c == 1 || isprime(c), return(1)))); 0} \\ David A. Corneth, Mar 30 2021
CROSSREFS
KEYWORD
nonn,base,nice
AUTHOR
N. J. A. Sloane, Bill Moran (moran1(AT)llnl.gov)
STATUS
approved