OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..451
FORMULA
E.g.f.: exp(Sum_{n>=2} x^F(n)/F(n)) with F = A000045.
MAPLE
a:= proc(n) option remember; `if`(n=0, 1, add(
`if`(issqr(5*j^2+4) or issqr(5*j^2-4),
a(n-j)*(j-1)!*binomial(n-1, j-1), 0), j=1..n))
end:
seq(a(n), n=0..25);
MATHEMATICA
a[n_] := a[n] = If[n == 0, 1, Sum[If[IntegerQ @ Sqrt[5*j^2+4] || IntegerQ @ Sqrt[5*j^2-4], a[n-j]*(j-1)!*Binomial[n-1, j-1], 0], {j, 1, n}]]; Table[ a[n], {n, 0, 25}] (* Jean-François Alcover, Jan 30 2017, translated from Maple *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, May 12 2016
STATUS
approved