OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
134 is in the sequence because it is a semiprime and it is divisible by its multiplicative digital root, 2.
MATHEMATICA
spQ[n_] := Plus @@ Last /@ FactorInteger@n == 2; mdrQ[n_] := Mod[n, NestWhile[Times @@ IntegerDigits@# &, n, UnsameQ, All]] == 0; Select[ Range@9754, spQ@# && mdrQ@# &] (* Robert G. Wilson v, Aug 04 2006 *)
mdr[n_]:=Module[{c=NestWhile[Times@@IntegerDigits[#]&, n, #>9&]}, If[c>0, c, Pi]]; Select[ Range[ 8000], PrimeOmega[#]==2&&Divisible[#, mdr[#]]&] (* Harvey P. Dale, Feb 27 2024 *)
PROG
(PARI) A031347(n)= { local(resul, ncpy); if(n<10, return(n) ); ncpy=n; resul = ncpy % 10; ncpy = (ncpy - ncpy%10)/10; while( ncpy > 0, resul *= ncpy %10; ncpy = (ncpy - ncpy%10)/10; ); return(A031347(resul)); } { for(n=4, 5000, if( bigomega(n)==2, dr=A031347(n); if(dr !=0 && n % dr == 0, print1(n, ", "); ); ); ); } \\ R. J. Mathar, May 23 2006
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Luc Stevens (lms022(AT)yahoo.com), May 20 2006
EXTENSIONS
Corrected by R. J. Mathar, May 23 2006
More terms from Robert G. Wilson v, Aug 04 2006
STATUS
approved