OFFSET
0,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..254
EXAMPLE
a(n) ~ c * (n-1)!^2, where c = A156648 = Product_{k>=1} (1 + 1/k^2) = sinh(Pi)/Pi = 3.67607791037497772...
MAPLE
b:= proc(n, i) option remember; `if`(i*(i+1)/2<n, 0, `if`(n=0, 1,
b(n, i-1)+b(n-i, min(n-i, i-1))*((i-1)!*binomial(n, i))^2))
end:
a:= n-> b(n$2):
seq(a(n), n=0..20); # Alois P. Heinz, Jul 27 2023
MATHEMATICA
nmax = 20; CoefficientList[Series[Product[(1+x^k/k^2), {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]!^2
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Jul 27 2019
STATUS
approved