OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..332
Wikipedia, Permutation
FORMULA
a(n) = A344855(2n,n).
a(n) ~ c * (1 + exp(2))^n * (n-1)!, where c = sqrt((exp(2) + 1)/(exp(2) - 1))/(2*Pi) = 0.1823720711148962856100934464088354177502714116352616187167... - Vaclav Kotesovec, Jul 15 2021, updated Mar 17 2024
a(0) = 1; a(n) = Sum_{k=0..n} binomial(2*n, n + k + 1)*Stirling2(n + k + 1, k + 1). - Detlef Meya, Jan 18 2024
MAPLE
b:= proc(n) option remember; `if`(n=0, 1, add(expand(x*
b(n-j)*binomial(n-1, j-1)*ceil(2^(j-2))), j=1..n))
end:
a:= n-> coeff(b(2*n), x, n):
seq(a(n), n=0..18);
MATHEMATICA
b[n_] := b[n] = If[n == 0, 1, Sum[Expand[x b[n-j] Binomial[n-1, j-1]* Ceiling[2^(j-2)]], {j, n}]];
a[n_] := Coefficient[b[2n], x, n];
Table[a[n], {n, 0, 18}] (* Jean-François Alcover, Aug 25 2021, after Alois P. Heinz *)
a[0] := 1; a[n_] := Sum[Binomial[2*n, n + k + 1]*StirlingS2[n + k + 1, k + 1], {k, 0, n}]; Flatten[Table[a[n] , {n, 0, 17}]] (* Detlef Meya, Jan 18 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Jun 14 2021
STATUS
approved