OFFSET
1,2
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
MATHEMATICA
rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]];
Block[{k}, Table[k = PrimeOmega[n];
Total@ Select[Range[n^PrimeNu[n]],
Divisible[n, rad[#]] && PrimeOmega[#] <= k &], {n, 60}]]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Nov 14 2024
STATUS
approved