a[0, 0] := 1; a[0, 1] := 1; a[0, n_] := a[0, n] = Sum[Binomial[2*n - k - 3, n - 2]*a[0, k], {k, 0, n - 1}]; a[1, 0] := 1; a[1, n_] := a[1, n] = a[1, n - 1] + a[0, n]; a[2, 0] := 1; a[2, n_] := a[2, n] = a[2, n - 1] + a[1, n]; Table[a[2, n], {n, 0, 50}] (*
G. C. Greubel, Sep 01 2017 *)