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