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