login
Expansion of Product_{m>=1} (1+m*q^m)^25.
2

%I #9 Sep 08 2022 08:44:46

%S 1,25,350,3625,30825,227005,1495225,8998625,50231225,262982425,

%T 1302361670,6141852925,27731605150,120415590250,504692324800,

%U 2048151994275,8069513499800,30937269482500,115647629802975

%N Expansion of Product_{m>=1} (1+m*q^m)^25.

%H G. C. Greubel, <a href="/A022653/b022653.txt">Table of n, a(n) for n = 0..1000</a>

%t With[{nmax=50}, CoefficientList[Series[Product[(1+m*q^m)^25,{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)^25)) \\ _G. C. Greubel_, Jul 18 2018

%o (Magma) Coefficients(&*[(1+m*x^m)^25: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_