OFFSET
0,18
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
Wikipedia, Unimodality, Unimodal function
FORMULA
G.f.: Sum_{n>=0} x^(8*n) / Product_{j=1..n-1} (1-x^j)^2.
a(n) ~ Pi^7 * 7! * exp(2*Pi*sqrt(n/3)) / (2^10 * 3^(17/4) * n^(19/4)). - Vaclav Kotesovec, Oct 24 2018
MAPLE
b:= proc(n, i) option remember; `if`(i>n, 0,
`if`(8*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