OFFSET
1,2
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
FORMULA
Dirichlet g.f.: zeta(s)^2 * Product_{primes p} (1 + 3/p^s - 1/p^(2*s)).
Dirichlet g.f.: zeta(s)^5 * Product_{primes p} (1 - 7/p^(2*s) + 11/p^(3*s) - 6/p^(4*s) + 1/p^(5*s)), (with a product that converges for s=1).
Sum_{k=1..n} a(k) ~ c * n * log(n)^4 / 24, where c = Product_{primes p} (1 - 7/p^2 + 11/p^3 - 6/p^4 + 1/p^5) = 0.091414252314317101861531055690354339957600046..., more precise (but very complicated) asymptotics can be obtained (in Mathematica notation) as Residue[Zeta[s]^5 * f[s] * n^s / s, {s, 1}], where f[s] = Product_{primes p} (1 - 7/p^(2*s) + 11/p^(3*s) - 6/p^(4*s) + 1/p^(5*s)).
MATHEMATICA
a[n_] := Times @@ ((3*Last[#] + 2) & /@ FactorInteger[n]); a[1] = 1; Array[a, 100] (* Amiram Eldar, Feb 25 2023 *)
PROG
(PARI) for(n=1, 100, print1(direuler(p=2, n, (1+3*X-X^2)/(1-X)^2)[n], ", "))
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Vaclav Kotesovec, Feb 25 2023
STATUS
approved