OFFSET
0,3
LINKS
Neville Holmes, Integer Sequence Combinations
EXAMPLE
a(10) = 7(7+1)/2 = 28 because 7 is the sum of the prime factors of 10.
MATHEMATICA
f[n_]:=Module[{c=Total[Times@@@FactorInteger[n]]}, (c(c+1))/2]; Join[{0, 0}, Array[f, 60, 2]] (* Harvey P. Dale, Aug 21 2011 *)
PROG
(PARI) s(n)=sum(i=1, omega(n), component(component(factor(n), 1), i)*component(component(factor(n), 2), i))
a(n)=s(n)*(s(n)+1)/2
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
W. Neville Holmes, Aug 29 2002
EXTENSIONS
More terms from Benoit Cloitre, Sep 02 2002
STATUS
approved