OFFSET
0,20
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
Wikipedia, Unimodality, Unimodal function
FORMULA
G.f.: Sum_{n>=0} x^(9*n) / Product_{j=1..n-1} (1-x^j)^2.
a(n) ~ Pi^8 * 8! * exp(2*Pi*sqrt(n/3)) / (2^11 * 3^(19/4) * n^(21/4)). - Vaclav Kotesovec, Oct 24 2018
MAPLE
b:= proc(n, i) option remember; `if`(i>n, 0,
`if`(9*i=n, 1, 0)+add(b(n-i*j, i+1)*(j+1), j=0..n/i))
end:
a:= n-> `if`(n=0, 1, b(n, 1)):
seq(a(n), n=0..70);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 10 2018
STATUS
approved