%I #11 Mar 27 2024 21:43:26

%S 2,3,4,5,7,8,9,11,13,16,17,19,21,23,25,26,27,29,31,32,33,35,37,41,43,

%T 47,49,52,53,59,61,63,64,65,66,67,71,73,74,79,81,83,86,89,91,95,97,99,

%U 101,103,104,106,107,109,110,113,121,125,127,128,131,132,134,137,138,139

%N Numbers k divisible by A075860(k).

%p f := proc (n) option remember; if isprime(n) then n else procname(convert(numtheory:-factorset(n), `+`)) end if end proc:

%p g := proc (n) if `mod`(n, f(n)) = 0 then n end if end proc:

%p map(g, [$2 .. 100]);

%t 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 *)

%Y Cf. A075860.

%K nonn

%O 1,1

%A _Rafik Khalfi_, Mar 09 2024