g:= proc(n, i, j) option remember; convert(series(`if`(j=0, 1,
`if`(i<0, 0, add(g(n, i-1, j-k)*x^(i*k)*binomial(
binomial(n, i), k), k=0..j))), x, 3), polynom)
end:
b:= proc(n, i) option remember; convert(series(`if`(n=0, 1, `if`(i<1, 0,
add(b(n-i*j, i-1)*g(i$2, j), j=0..n/i))), x, 3), polynom)
end:
a:= n-> coeff(b(n$2), x, 2):
seq(a(n), n=0..45);
