|
%I #16 Sep 08 2022 08:44:46
%S 1,13,104,663,3614,17576,78299,324766,1269242,4715204,16762551,
%T 57327556,189418658,606787572,1890046210,5738539729,17019191579,
%U 49394158541,140507716414,392299039821,1076369417474,2905414115877,7722941644821,20233362612424,52288914446548,133389316899462
%N Expansion of Product_{m>=1} (1 + m*q^m)^13.
%H G. C. Greubel, <a href="/A022641/b022641.txt">Table of n, a(n) for n = 0..1000</a>
%p [seq(coeff(series(mul((1+m*q^m)^(13), m=1..100),q,101),q,j),j=0..25)]; # _Muniru A Asiru_, Feb 18 2018
%t With[{nmax = 50}, CoefficientList[Series[Product[(1 + k*q^k)^13, {k, 1, nmax}], {q, 0, nmax}], q]] (* _G. C. Greubel_, Feb 17 2018 *)
%o (PARI) m=50; q='q+O('q^m); Vec(prod(n=1,m,(1+n*q^n)^13)) \\ _G. C. Greubel_, Feb 17 2018
%o (Magma) Coefficients(&*[(1+m*x^m)^13:m in [1..40]])[1..40] where x is PolynomialRing(Integers()).1; // _G. C. Greubel_, Feb 17 2018
%Y Column k=13 of A297321.
%K nonn
%O 0,2
%A _N. J. A. Sloane_
|
