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

1, -16, 88, -96, -836, 2592, 1824, -11616, -12858, 27248, 160064, -123168, -634664, -132000, 1016416, 4677376, 375569, -13049296, -24306064, -14447360, 115090052, 159327776, 10482048, -443275328, -1133277978

Offset: 0

Author

N. J. A. Sloane

Keywords

• sign

Links

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

Programs

Magma

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

Mathematica

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

PARI

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