OFFSET
1,2
MATHEMATICA
Block[{k}, rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]];
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