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

1, 9, 54, 273, 1197, 4761, 17577, 60957, 200799, 633007, 1920510, 5633667, 16037700, 44439840, 120165858, 317762553, 823240341, 2092864401, 5228118701, 12848849154, 31100190048, 74208885351, 174708121455, 406132690635, 932871440739, 2118595079790, 4759875472491

Offset

0,2

Links

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

Formula

G.f.: exp(9*Sum_{j>=1} Sum_{k>=1} (-1)^(j+1)*k^j*x^(j*k)/j). - Ilya Gutkovskiy, Feb 08 2018

Mathematica

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

Prog

(PARI) m=50; q='q+O('q^m); Vec(prod(n=1, m, (1+n*q^n)^9)) \\ G. C. Greubel, Feb 17 2018
(Magma) Coefficients(&*[(1+m*x^m)^9:m in [1..40]])[1..40] where x is PolynomialRing(Integers()).1; // G. C. Greubel, Feb 17 2018

Crossrefs

Column k=9 of A297321.

Sequence in context: A169796 A359722 A027472 * A001392 A188428 A243415

Adjacent sequences: A022634 A022635 A022636 * A022638 A022639 A022640

Keyword

nonn

Author

N. J. A. Sloane

Status

approved