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