OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..347
FORMULA
a(n) = A001100(2n+1,n).
MAPLE
S:= proc(n) option remember; `if`(n<4, [1, 1, 2*t, 4*t+2*t^2]
[n+1], expand((n+1-t)*S(n-1) -(1-t)*(n-2+3*t)*S(n-2)
-(1-t)^2*(n-5+t)*S(n-3) +(1-t)^3*(n-3)*S(n-4)))
end:
a:= n-> coeff(S(2*n+1), t, n):
seq(a(n), n=0..20);
MATHEMATICA
S[n_] := S[n] = If[n < 4, {1, 1, 2*t, 4*t + 2*t^2}[[n + 1]], Expand[
(n + 1 - t)*S[n - 1] -
(1 - t)*(n - 2 + 3*t)*S[n - 2] -
(1 - t)^2*(n - 5 + t)*S[n - 3] +
(1 - t)^3*(n - 3)*S[n - 4]]];
a[n_] := Coefficient[S[2*n + 1], t, n];
Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Apr 21 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Dec 02 2018
STATUS
approved