OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
psi[n_] := n * Times @@ (1 + 1/FactorInteger[n][[;; , 1]]); psi[1] = 1; Join[{1}, psi /@ Select[Range[40000], AllTrue[Last /@ FactorInteger[#], #1 > 1 && OddQ[#1] &] &]]
PROG
(PARI) dedpsi(f) = prod(i = 1, #f~, (f[i, 1] + 1) * f[i, 1]^(f[i, 2]-1));
lista(max) = {my(f, ans); print1(1, ", "); for(k = 2, max, f = factor(k); ans = 1; for (i = 1, #f~, if (f[i, 2] == 1 || !(f[i, 2] % 2), ans = 0; break)); if(ans, print1(dedpsi(f), ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Mar 22 2024
STATUS
approved