OFFSET
0,2
COMMENTS
Sequence A074374 is similar, based on the triangular numbers, giving s(s+1)/2 with s=sopfr(n). Here it is based on the square pyramidal numbers, giving s(s+1)(2s+1)/6 with s=sopfr(n).
MATHEMATICA
sopfr = Function[x, Plus @@ Map[Times @@ # &, FactorInteger[x]]]; Map[ #(# + 1)(2# + 1)/6 &, sopfr /@ Range[130]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Carlos Alves, Dec 16 2007
STATUS
approved