A022644 Expansion of Product_{m>=1} (1 + m*q^m)^16.

1, 16, 152, 1120, 6972, 38368, 191968, 889184, 3862214, 15881616, 62275840, 234205472, 848652120, 2974133152, 10112507808, 33448941824, 107874784017, 339879773648, 1047953793136, 3166817754880, 9391718326404, 27366626142688, 78435144301696, 221322772710464, 615375631077094

Offset: 0

Author

N. J. A. Sloane

Keywords

• nonn

Links

• G. C. Greubel, Table of n, a(n) for n = 0..1000

Crossrefs

Column k=16 of A297321.

Programs

Magma

Coefficients(&*[(1+m*x^m)^16:m in [1..40]])[1..40] where x is PolynomialRing(Integers()).1; // G. C. Greubel, Feb 17 2018

Maple

[seq(coeff(series(mul((1+m*q^m)^(16), m=1..100),q,101),q,j),j=0..25)]; # Muniru A Asiru, Feb 18 2018

Mathematica

With[{nmax = 50}, CoefficientList[Series[Product[(1 + k*q^k)^16, {k, 1, nmax}], {q, 0, nmax}], q]] (* G. C. Greubel, Feb 17 2018 *)

PARI

m=50; q='q+O('q^m); Vec(prod(n=1,m,(1+n*q^n)^16)) \\ G. C. Greubel, Feb 17 2018