OFFSET
1,2
COMMENTS
There is a family of multiplicative sequences based on trinomial numbers (A027907) and some fixed integer k >= 0. Let a(k,n), k >= 0, n > 0, be multiplicative with a(k,p^e) = Sum_{i=0..e} trinomial(k,i) for prime p and e >= 0, where trinomial(n,k) = 0 if 2*n < k. These sequences have the Dirichlet g.f.: Sum_{n>=1} a(k,n)/n^s = (zeta(s))^(k+1) / (zeta(3*s))^k. For several members of the family see A000012 (k=0), A073184 (k=1), and this sequence (k=2).
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(p^e) = Sum_{i=0..e} trinomial(2,i) for prime p and e >= 0, where trinomial(n,k) = 0 if 2*n < k.
Let b(n), n > 0, be the Dirichlet inverse of a(n); b(n) is multiplicative with b(p^(3*e)) = 1 for e >= 0, and b(p^(3*e-2)) = -3*e and b(p^(3*e-1)) = 3*e for e > 0 and prime p.
Sum_{k=1..n} a(k) ~ n * (log(n)^2/2 + (3*gamma - 6*zeta'(3)/zeta(3) - 1)*log(n) + 1 - 3*gamma + 3*gamma^2 + 6*(1 - 3*gamma)*zeta'(3)/zeta(3) + 27*zeta'(3)^2 / zeta(3)^2 - 9*zeta''(3)/zeta(3) - 3*sg1) / zeta(3)^2, where gamma is the Euler-Mascheroni constant A001620, zeta(3) = A002117, zeta'(3) = -A244115, zeta''(3) = A340442 and sg1 is the first Stieltjes constant (see A082633). - Vaclav Kotesovec, Nov 20 2021
PROG
(PARI) {T(n, k) = if( n<0, 0, polcoeff( (1 + x + x^2)^n, k))}; \\ A027907
a(n)={my(f=factor(n)); prod(k=1, #f[, 1], sum(i=0, f[k, 2], T(2, i)))};
for(j=1, 75, print1(a(j), ", ")) \\ Hugo Pfoertner, Feb 13 2021
(PARI) for(n=1, 100, print1(direuler(p=2, n, (1 - X^3)^2/(1 - X)^3)[n], ", ")) \\ Vaclav Kotesovec, Nov 20 2021
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Werner Schulte, Feb 09 2021
STATUS
approved