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