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

1, 31, 527, 6510, 65162, 559271, 4259276, 29444575, 187761854, 1117650502, 6267006128, 33341859350, 169281687755, 824110662334, 3862233690794, 17483338870409, 76663232042998, 326438219863590

Offset: 0

Author

N. J. A. Sloane

Keywords

• nonn

Links

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

Programs

Magma

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

Mathematica

With[{nmax=50}, CoefficientList[Series[Product[(1+m*q^m)^31,{m,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)^31)) \\ G. C. Greubel, Jul 18 2018