|
%I #9 Sep 08 2022 08:44:46
%S 1,31,527,6510,65162,559271,4259276,29444575,187761854,1117650502,
%T 6267006128,33341859350,169281687755,824110662334,3862233690794,
%U 17483338870409,76663232042998,326438219863590
%N Expansion of Product_{m>=1} (1+m*q^m)^31.
%H G. C. Greubel, <a href="/A022659/b022659.txt">Table of n, a(n) for n = 0..1000</a>
%t With[{nmax=50}, CoefficientList[Series[Product[(1+m*q^m)^31,{m,1,nmax}],{q,0,nmax}],q]] (* _G. C. Greubel_, Jul 18 2018 *)
%o (PARI) m=50; q='q+O('q^m); Vec(prod(n=1,m,(1+n*q^n)^31)) \\ _G. C. Greubel_, Jul 18 2018
%o (Magma) Coefficients(&*[(1+m*x^m)^31:m in [1..40]])[1..50] where x is PolynomialRing(Integers()).1; // _G. C. Greubel_, Jul 18 2018
%K nonn
%O 0,2
%A _N. J. A. Sloane_
|
