OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
f[p_, e_] := If[e > 1, p^e, 1]; powPart[n_] := Times @@ f @@@ FactorInteger[n]; Select[Array[powPart, 200], # > 1 &]
PROG
(PARI) lista(kmax) = {my(p, f); for(k = 1, kmax, f = factor(k); p = prod(i=1, #f~, if(f[i, 2] > 1, f[i, 1]^f[i, 2], 1)); if(p > 1, print1(p, ", "))); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Dec 09 2023
STATUS
approved