OFFSET
10,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 10..10000
G. E. Andrews and S. C. F. Rose, MacMahon's sum-of-divisors functions, Chebyshev polynomials, and Quasi-modular forms, arXiv:1010.5769 [math.NT], 2010.
FORMULA
G.f.: (1/9) * ( Sum_{k>=4} (-1)^k * (2*k+1) * binomial(k+4,8) * q^(k*(k+1)/2) ) / ( Sum_{k>=0} (-1)^k * (2*k+1) * q^(k*(k+1)/2) ).
a(n) = (5*sigma_7(n) - (126*n-441)*sigma_5(n) + (756*n^2-4410*n+4935)*sigma_3(n) - (840*n^3-5880*n^2+9870*n-3229)*sigma(n))/967680. - Seiichi Manyama, Jul 24 2024
PROG
(PARI) a(n) = (5*sigma(n, 7)-(126*n-441)*sigma(n, 5)+(756*n^2-4410*n+4935)*sigma(n, 3)-(840*n^3-5880*n^2+9870*n-3229)*sigma(n))/967680; \\ Seiichi Manyama, Jul 24 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Sep 15 2023
STATUS
approved