OFFSET
1,1
EXAMPLE
38747603 is in the sequence because {1,769,50387} is the set of proper
divisors of 38747603 and 3!+8!+7!+4!+7!+6!+0!+3!=1+769+50387.
MATHEMATICA
Do[h = IntegerDigits[n]; l = Length[h]; If[DivisorSigma[1, n] - n == Sum[h[[k]]!, {k, l}], Print[n]], {n, 60000000}]
Select[Range[2413*10^5], Total[Most[Divisors[#]]]==Total[ IntegerDigits[ #]!]&] (* Harvey P. Dale, Oct 11 2018 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Farideh Firoozbakht, Dec 12 2004
EXTENSIONS
a(9)-a(22) from Donovan Johnson, Dec 07 2008
Comment deleted by Harvey P. Dale, Oct 11 2018
STATUS
approved