(MAGMAMagma) m:=40; R<q>:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1+q^k)^DivisorSigma(3, k): k in [1..m]]) )); // G. C. Greubel, Oct 30 2018
(MAGMAMagma) m:=40; R<q>:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1+q^k)^DivisorSigma(3, k): k in [1..m]]) )); // G. C. Greubel, Oct 30 2018
reviewed
approved
proposed
reviewed
editing
proposed
(PARI) m=5040; x='x+O('x^m); Vec(prod(k=1, m+2, , 1/(1+x^k)^sigma(k, 3))) \\ G. C. Greubel, Oct 30 2018
(MAGMA) m:=5040; R<q>:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1+q^k)^DivisorSigma(3, k): k in [1..(m+2)]]) )); // G. C. Greubel, Oct 30 2018
reviewed
editing
proposed
reviewed
editing
proposed
with(numtheory): seq(coeff(series(mul(1/(1+x^k)^(sigma[3](k)), k=1..n), x, n+1), x, n), n = 0 .. 30); # Muniru A Asiru, Oct 31 2018
proposed
editing
editing
proposed
(PARI) m=50; x='x+O('x^m); Vec(prod(k=1, m+2, 1/(1+x^k)^sigma(k, 3))) \\ G. C. Greubel, Oct 30 2018
(MAGMA) m:=50; R<q>:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1+q^k)^DivisorSigma(3, k): k in [1..(m+2)]]) )); // G. C. Greubel, Oct 30 2018
approved
editing