A301551 Expansion of Product_{k>=1} (1 + x^k)^(sigma_7(k)).

1, 1, 129, 2317, 26957, 385147, 5514889, 70250881, 866874825, 10634404922, 126906497939, 1470673175003, 16705788322140, 186487470519166, 2044203433733016, 22025647881901542, 233686866722213324, 2443978994099801452, 25211475391206919299, 256716054713570158748

Offset: 0

Author

Vaclav Kotesovec, Mar 23 2018

Keywords

• nonn

Links

• Seiichi Manyama, Table of n, a(n) for n = 0..1000

Crossrefs

Cf. A013955, A107742, A301545.

Programs

Mathematica

nmax = 30; CoefficientList[Series[Product[(1+x^k)^DivisorSigma[7, k], {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(9 * Pi^(8/9) * (17*Zeta(9))^(1/9) * n^(8/9) / 2^(29/9)) * (17*Zeta(9)/Pi)^(1/18) / (3 * 2^(883/1440) * n^(5/9)).

G.f.: exp(Sum_{k>=1} sigma_8(k)*x^k/(k*(1 - x^(2*k)))). - Ilya Gutkovskiy, Oct 26 2018