OFFSET
1,1
COMMENTS
f(x) = sum of prime factors without multiplicity, so that f(1500) = 2+3+5 = 10.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..5000
MATHEMATICA
spf9Q[n_]:=Module[{spf=TakeDrop[NestList[Total[FactorInteger[#][[All, 1]]]&, n, 9], 9]}, AllTrue[spf[[1]], !PrimeQ[#]&]&&PrimeQ[spf][[2, 1]]]; Select[ Range[ 250000], spf9Q] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Apr 20 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved