OFFSET
9,4
LINKS
G. C. Greubel, Table of n, a(n) for n = 9..1000
FORMULA
G.f.: ( Sum_{k>=1} x^binomial(k+2,3) )^9.
MATHEMATICA
nmax = 66; CoefficientList[Series[Sum[x^Binomial[k + 2, 3], {k, 1, nmax}]^9, {x, 0, nmax}], x] // Drop[#, 9] &
PROG
(Magma)
R<x>:=PowerSeriesRing(Integers(), 70);
Coefficients(R!( (&+[x^Binomial(j+2, 3): j in [1..70]])^9 )); // G. C. Greubel, Jul 18 2022
(SageMath)
def f(m, x): return ( sum( x^(binomial(j+2, 3)) for j in (1..8) ) )^m
def A341808_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( f(9, x) ).list()
a=A341808_list(100); a[9:71] # G. C. Greubel, Jul 18 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 20 2021
STATUS
approved