editing
approved
editing
approved
A034676 := proc(n)
a :=1 ;
for pe in ifactors(n)[2] do
p := op(1, pe) ;
e := op(2, pe) ;
a := a*(p^(2*e)+1) ;
end do:
a ;
end proc:
seq(A034676(n), n=1..40) ; # R. J. Mathar, Jul 12 2024
approved
editing
editing
approved
Sum_{k>=1} 1/a(k) = 1.5594563610641446770272272038182777336348840179730233519185104374159616326... - Vaclav Kotesovec, Sep 20 2020
approved
editing
proposed
approved
editing
proposed
Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/UnitaryDivisorFunction.html">Unitary Divisor Function</a>.
Wikipedia, <a href="http://en.wikipedia.org/wiki/Unitary_divisor">Unitary divisor</a>.
f[p_, e_] := p^(2*e)+1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 14 2020 *)
approved
editing
editing
approved
Sum_{k=1..n} a(k) ~ 30 * Zeta(3) * n^3 / Pi^4. - Vaclav Kotesovec, Jan 11 2019
approved
editing