OFFSET
1,1
MATHEMATICA
a[n_]:=NestWhile[#+Total[Times@@@FactorInteger[#]]&, n, !PrimeQ[#]&]; t={}; Do[If[!PrimeQ[n], AppendTo[t, a[n]]], {n, 4, nn=678}]; Complement[Prime[Range[PrimePi[nn]]], Select[Union[t], #<nn&]] (* Jayanta Basu, Jun 01 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Patrick De Geest, Sep 15 1999
STATUS
approved