(MAGMAMagma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[(1 - x^(7*j))^56/(1 - x^j)^57: j in [1..m+2]]) )); // G. C. Greubel, Nov 18 2018
(MAGMAMagma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[(1 - x^(7*j))^56/(1 - x^j)^57: j in [1..m+2]]) )); // G. C. Greubel, Nov 18 2018
proposed
approved
editing
proposed
(PARI) my(N=30, x='x+O('x^30N)); Vec(prod(j=1, 30, N, (1 - x^(7*j))^56/(1 - x^j)^57)) \\ G. C. Greubel, Nov 18 2018
proposed
editing
editing
proposed
(MAGMA) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[(1 - x^(7*j))^56/(1 - x^j)^57: j in [1..30m+2]]) )); // G. C. Greubel, Nov 18 2018
proposed
editing
editing
proposed
(PARI) my(x='x+O('x^30)); Vec(prod(j=1, 5, 30, (1 - x^(7*j))^56/(1 - x^j)^57)) \\ G. C. Greubel, Nov 18 2018
(MAGMA) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[(1 - x^(7*j))^56/(1 - x^j)^57: j in [1..530]]) )); // G. C. Greubel, Nov 18 2018
proposed
editing
editing
proposed
prec = 30
x = R.gen().O(30prec)
s = prod((1 - x^(7*j))^56/(1 - x^j)^57 for j in (1..5prec))
print(s.coefficients() ) # G. C. Greubel, Nov 18 2018
approved
editing