A321239 a(n) = [x^(n^3)] Product_{k=1..n} Sum_{m>=0} x^(k^2*m).
1, 1, 3, 16, 141, 1534, 19111, 262103, 3853373, 59763670, 966945204, 16191250596, 278933800080, 4921604827876, 88627915588351, 1624349874930925, 30231112607904743, 570284342486800214, 10887435073866747752, 210086404047975194316, 4092940691144348506396, 80432925119259253535963
Offset: 0
Keywords
Examples
1^2* 0 + 2^2*0 + 3^2*3 = 27. 1^2* 1 + 2^2*2 + 3^2*2 = 27. 1^2* 2 + 2^2*4 + 3^2*1 = 27. 1^2* 3 + 2^2*6 + 3^2*0 = 27. 1^2* 5 + 2^2*1 + 3^2*2 = 27. 1^2* 6 + 2^2*3 + 3^2*1 = 27. 1^2* 7 + 2^2*5 + 3^2*0 = 27. 1^2* 9 + 2^2*0 + 3^2*2 = 27. 1^2*10 + 2^2*2 + 3^2*1 = 27. 1^2*11 + 2^2*4 + 3^2*0 = 27. 1^2*14 + 2^2*1 + 3^2*1 = 27. 1^2*15 + 2^2*3 + 3^2*0 = 27. 1^2*18 + 2^2*0 + 3^2*1 = 27. 1^2*19 + 2^2*2 + 3^2*0 = 27. 1^2*23 + 2^2*1 + 3^2*0 = 27. 1^2*27 + 2^2*0 + 3^2*0 = 27. So a(3) = 16.
Programs
-
PARI
{a(n) = polcoeff(prod(i=1, n, sum(j=0, n^3\i^2, x^(i^2*j)+x*O(x^(n^3)))), n^3)} -
PARI
{a(n) = polcoeff( 1/prod(k=1,n, 1 - x^(k^2) +x*O(x^(n^3)) ), n^3) } for(n=0,20, print1(a(n),", ")) \\ Paul D. Hanna, Feb 02 2024
Comments