OFFSET
1,4
COMMENTS
Obviously, all terms are odd.
LINKS
FORMULA
Multiplicative with a(p^e) = (p^(e+1)-1)/(p-1) if e is even and 1 otherwise.
Sum_{k=1..n} a(k) ~ c * n^(3/2), where c = (zeta(3)/3) * Product_{p prime} (1 + 1/p^(3/2) + 1/p^2 - 1/p^(5/2)) = 1.008259499413... . - Amiram Eldar, Nov 20 2022
Dirichlet g.f.: zeta(2*s) * zeta(2*s-2) * Product_{p prime} (1 + 1/p^s + 1/p^(2*s-1) - 1/p^(3*s-2)). - Amiram Eldar, Sep 03 2023
MATHEMATICA
f[p_, e_] := If[EvenQ[e], (p^(e + 1) - 1)/(p - 1), 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Feb 23 2022 *)
PROG
(PARI)
A350388(n) = { my(m=1, f=factor(n)); for(k=1, #f~, if(0==(f[k, 2]%2), m *= (f[k, 1]^f[k, 2]))); (m); };
(Python)
from math import prod
from sympy import factorint
def A351568(n): return prod(1 if e % 2 else (p**(e+1)-1)//(p-1) for p, e in factorint(n).items()) # Chai Wah Wu, Feb 24 2022
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Feb 23 2022
STATUS
approved