OFFSET
1,1
COMMENTS
f(x) = sum of prime factors with multiplicity, so that f(1500) = 2+2+3+5+5+5 = 22.
MATHEMATICA
spf[n_]:=Total[Flatten[Table[#[[1]], {#[[2]]}]&/@FactorInteger[n]]]; pf7Q[ n_]:= Module[ {lst=NestList[spf, n, 7]}, NoneTrue[Most[lst], PrimeQ] && !MemberQ[Most[lst], 4]&&(PrimeQ[Last[lst]]||Last[lst]==4)]; Select[Range[ 3600], pf7Q] (* Harvey P. Dale, Jul 11 2014 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved