OFFSET
1,1
MAPLE
f := proc (n) option remember; if isprime(n) then n else procname(convert(numtheory:-factorset(n), `+`)) end if end proc:
g := proc (n) if `mod`(n, f(n)) = 0 then n end if end proc:
map(g, [$2 .. 100]);
MATHEMATICA
a[n_] := If[n == 1, 0, FixedPoint[Total[FactorInteger[#][[All, 1]]]&, n]]; r=140; Select[Range[2, r], Divisible[#, Part[Array[a, r], #]]&] (* James C. McMahon, Mar 10 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Rafik Khalfi, Mar 09 2024
STATUS
approved