OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
Dirichlet g.f.: Product_{primes p} (1 + p^3/(p^s - p)).
Dirichlet g.f.: zeta(s-3) * zeta(s-1) * Product_{primes p} (1 + p^(4-2*s) - p^(6-2*s) - p^(1-s)).
Sum_{k=1..n} a(k) ~ c * zeta(3) * n^4 / 4, where c = Product_{primes p} (1 - 1/p^2 - 1/p^3 + 1/p^4) = A330523 = 0.53589615382833799980850263131854595064822237...
From Amiram Eldar, Sep 01 2023: (Start)
Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + 1/(p^2*(p-1))) = A065483. (End)
MATHEMATICA
g[p_, e_] := p^(e+2); a[1] = 1; a[n_] := Times @@ g @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) for(n=1, 100, print1(direuler(p=2, n, 1 + p^3*X/(1 - p*X))[n], ", "))
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Vaclav Kotesovec, Mar 06 2023
STATUS
approved