OFFSET
0,1
COMMENTS
sopfr(k) = sum of the prime factors of k (with multiplicity).
LINKS
Paolo P. Lava, Table of n, a(n) for n = 0..200
EXAMPLE
sopfr(72) = 2+2+2+3+3 = 12 = 2+2+3+5 = sopfr(72 - sopfr(72)), so 72 is in the sequence.
MATHEMATICA
sopf[n_]:=Module[{sopfn=Total[Times@@@FactorInteger[n]], m}, m=n+sopfn; If[n==m-Total[Times@@@FactorInteger[m]], m, 0]]; DeleteCases[Table[ sopf[n], {n, 55000}], _?(#==0&)] (* Harvey P. Dale, Jun 15 2011 *)
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Patrick De Geest, Sep 15 1999
EXTENSIONS
Edited by Jon E. Schoenfield, Dec 25 2016
STATUS
approved