OFFSET
1,16
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A049419(n^2).
Multiplicative with a(p^e) = A183063(e), or equivalently, a(p^e) = 0 if e is odd, and A000005(e/2) if e is even.
a(n) >= 0, with equality if and only if n is not a square number (A000037).
a(n) <= A049419(n), with equality if and only if n = 1.
Sum_{k=1..n} a(k) ~ c * sqrt(n), where c = 1.602317... (A327837).
MATHEMATICA
f[p_, e_] := If[OddQ[e], 0, DivisorSigma[0, e/2]]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) a(n) = vecprod(apply(x -> if(x%2, 0, numdiv(x/2)), factor(n)[, 2]));
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Amiram Eldar, Jan 11 2024
STATUS
approved