A162737 G.f. is the polynomial (Product_{k=1..31} (1 - x^(3*k)))/(1-x)^31.

1, 31, 496, 5455, 46345, 324136, 1942335, 10249065, 48578364, 209961884, 837318680, 3110233844, 10844253964, 35718892168, 111747443068, 333597393565, 954069339889, 2623073605414, 6953911464088, 17823614247238, 44273067884650

Offset

0,2

Comments

This is a row of the triangle in A162499. Only finitely many terms are nonzero.

Links

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

Maple

m:=31: seq(coeff(series(mul((1-x^(3*i)), i=1..m)/(1-x)^m, x, n+1), x, n), n=0..21); # Muniru A Asiru, Jul 07 2018

Mathematica

CoefficientList[Series[Times@@(1-x^(3*Range[31]))/(1-x)^31, {x, 0, 50}], x] (* G. C. Greubel, Jul 07 2018 *)

Prog

(PARI) x='x+O('x^50); A = prod(k=1, 31, (1-x^(3*k)))/(1-x)^31; Vec(A) \\ G. C. Greubel, Jul 07 2018
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); F:=(&*[(1-x^(3*k)): k in [1..31]])/(1-x)^31; Coefficients(R!(F)); // G. C. Greubel, Jul 07 2018

Crossrefs

Sequence in context: A161636 A161977 A162378 * A010983 A022595 A125488

Adjacent sequences: A162734 A162735 A162736 * A162738 A162739 A162740

Keyword

nonn

Author

N. J. A. Sloane, Dec 02 2009

Status

approved