OFFSET
1,1
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 400 terms from Harry J. Smith)
MATHEMATICA
Select[Range[2, 10^5], Times @@ IntegerDigits @ # == Plus @@ Times @@@ FactorInteger @ # &] (* Giovanni Resta, Apr 23 2017 *)
PROG
(PARI) ProdD(x)= { local(p=1); while (x>9 && p>0, p*=(x-10*(x\10)); x\=10); return(p*x) }
sopfr(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]*f[i, 2]); return(s) }
{ n=0; for (m=1, 10^9, if (ProdD(m) == sopfr(m), write("b065774.txt", n++, " ", m); if (n==400, return)) ) } \\ Harry J. Smith, Oct 30 2009
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Jason Earls, Dec 04 2001
EXTENSIONS
Offset changed from 0 to 1 by Harry J. Smith, Oct 30 2009
STATUS
approved