OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..47
FORMULA
a(n) = [x^(9^n)] 1/Product_{j>=0}(1-x^(9^j)).
EXAMPLE
a(1) = 2, because there are 2 partitions of 9^1 into powers of 9: [1,1,1,1,1,1,1,1,1], [9].
MAPLE
g:= proc(b, n, k) option remember; local t; if b<0 then 0 elif b=0 or n=0 or k<=1 then 1 elif b>=n then add(g(b-t, n, k) *binomial(n+1, t) *(-1)^(t+1), t=1..n+1); else g(b-1, n, k) +g(b*k, n-1, k) fi end: a:= n-> g(1, n, 9): seq(a(n), n=0..13);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 11 2008
STATUS
approved