A330019 Expansion of e.g.f. 1 / (1 - Sum_{k>=1} binomial(k,floor(k/2)) * x^k / k!).
1, 1, 4, 21, 150, 1330, 14180, 176295, 2505230, 40049226, 711379872, 13899553206, 296270826852, 6841305568812, 170127212242416, 4532854743105975, 128824523061126750, 3890041395675793930, 124375112406132404960, 4197530354920789582410, 149118181703716510545260
Offset: 0
Keywords
Programs
-
Mathematica
nmax = 20; CoefficientList[Series[1/(2 - BesselI[0, 2 x] - BesselI[1, 2 x]), {x, 0, nmax}], x] Range[0, nmax]! a[0] = 1; a[n_] := a[n] = Sum[Binomial[n, k] Binomial[k, Floor[k/2]] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 20}]
Formula
E.g.f.: 1 / (2 - BesselI(0,2*x) - BesselI(1,2*x)).
a(0) = 1; a(n) = Sum_{k=1..n} binomial(n,k) * A001405(k) * a(n-k).