OFFSET
1,3
COMMENTS
In general, if k>=1, m>=1 and a(n) = Sum_{j=1..n-1} sigma_k(j) * sigma_m(n-j), then Sum_{j=1..n} a(j) ~ Gamma(k+1) * Gamma(m+1) * zeta(k+1) * zeta(m+1) * n^(k+m+2) / Gamma(k+m+3).
FORMULA
Sum_{k=1..n} a(k) ~ Pi^4 * zeta(3) * n^7 / 37800.
MATHEMATICA
Table[Sum[DivisorSigma[2, k]*DivisorSigma[3, n-k], {k, n-1}], {n, 1, 50}]
PROG
(PARI) a(n) = sum(k=1, n-1, sigma(k, 2) * sigma(n-k, 3)); \\ Michel Marcus, Sep 19 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Sep 19 2024
STATUS
approved