OFFSET
1,1
LINKS
Vincenzo Librandi and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 250 terms from Librandi)
EXAMPLE
4614 = 2*3*769 and 4!+6!+1!+4! = 769.
MATHEMATICA
Select[Range[220000], FactorInteger[#][[-1, 1]]==Total[IntegerDigits[#]!]&] (* Harvey P. Dale, Mar 08 2014 *)
PROG
(PARI) gpf(n)=if(n<2, 1, my(f=factor(n)[, 1]); f[#f])
is(n)=my(d=digits(n), t=sum(i=1, #d, d[i]!)); n%t==0&&isprime(t)&&gpf(n/t)<=t \\ Charles R Greathouse IV, Mar 10 2014
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Jason Earls, Sep 20 2002
STATUS
approved