A365665 Expansion of Sum_{0
1, 3, 9, 22, 51, 108, 208, 390, 693, 1193, 1977, 3195, 4995, 7722, 11583, 17164, 24882, 35685, 50205, 70083, 96300, 131101, 176358, 235377, 310651, 407352, 529074, 682750, 874038, 1112085, 1405521, 1766259, 2206413, 2741431, 3389052, 4168089, 5103450, 6218469
Offset: 15
Keywords
Links
- Vaclav Kotesovec, Table of n, a(n) for n = 15..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.
Crossrefs
Programs
-
Mathematica
nmax = 60; Drop[CoefficientList[Series[-1/11 * Sum[(-1)^k*(2*k + 1)*Binomial[k + 5, 10]*x^(k*(k + 1)/2), {k, 5, nmax}]/Sum[(-1)^k*(2*k + 1)*x^(k*(k + 1)/2), {k, 0, nmax}], {x, 0, nmax}], x], 15] (* Vaclav Kotesovec, Jul 29 2025 *) (* or *) Table[(10679/17203200 - 1571*n/774144 + 133*n^2/92160 - n^3/3072 + n^4/46080) * DivisorSigma[1, n] + (1571/1548288 - 133*n/122880 + 3*n^2/10240 - n^3/46080) * DivisorSigma[3, n] + (133/1228800 - n/20480 + n^2/215040) * DivisorSigma[5, n] + (1/516096 - n/3096576) * DivisorSigma[7, n] + DivisorSigma[9, n]/154828800, {n, 15, 60}] (* Vaclav Kotesovec, Jul 29 2025 *)
Formula
G.f.: -(1/11) * ( Sum_{k>=5} (-1)^k * (2*k+1) * binomial(k+5,10) * q^(k*(k+1)/2) ) / ( Sum_{k>=0} (-1)^k * (2*k+1) * q^(k*(k+1)/2) ).
From Vaclav Kotesovec, Jul 29 2025: (Start)
a(n) = (10679/17203200 - 1571*n/774144 + 133*n^2/92160 - n^3/3072 + n^4/46080)*sigma(n) + (1571/1548288 - 133*n/122880 + 3*n^2/10240 - n^3/46080)*sigma_3(n) + (133/1228800 - n/20480 + n^2/215040)*sigma_5(n) + (1/516096 - n/3096576)*sigma_7(n) + sigma_9(n)/154828800.
Sum_{k=1..n} a(k) ~ Pi^10 * n^10 / 144850083840000.
(End)
Comments