OFFSET
0,2
COMMENTS
See Dershowitz (2017) for precise definition.
LINKS
Nachum Dershowitz, Touchard’s Drunkard, Journal of Integer Sequences, Vol. 20 (2017), #17.1.5.
FORMULA
MAPLE
b:= n-> binomial(n, floor(n/2))*binomial(n+1, floor((n+1)/2)):
C:= n-> binomial(2*n, n)/(n+1):
a:= n-> add(binomial(n, 2*k)*C(k)*b(n-2*k), k=0..n/2):
seq(a(n), n=0..25); # Alois P. Heinz, Dec 06 2024
# second Maple program:
a:= proc(n) option remember; `if`(n<4, [1, 2, 7, 24][n+1],
(8*(14*n^4+85*n^3+190*n^2+188*n+63)*a(n-1)+4*(n-1)*
(80*n^4+418*n^3+676*n^2+269*n-108)*a(n-2)-96*(n-1)*(n-2)*
(10*n^2+31*n+27)*a(n-3)-144*(n-1)*(n-2)*(n-3)*(8*n^2+33*n+36)*
a(n-4))/((n+4)*(n+3)*(n+2)*(8*n^2+17*n+11)))
end:
seq(a(n), n=0..25); # Alois P. Heinz, Dec 06 2024
MATHEMATICA
b[n_] := Binomial[n, Floor[n/2]]*Binomial[n+1, Floor[(n+1)/2]];
c[n_] := Binomial[2*n, n]/(n+1);
a[n_] := Sum[Binomial[n, 2*k]*c[k]*b[n - 2*k], {k, 0, n/2}];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Jan 28 2025, after Alois P. Heinz *)
PROG
(Python)
from math import comb as binomial
def C(n): return (binomial(2*n, n)//(n+1)) # Catalan numbers
def a(n):
return sum(binomial(n, k)*C((k+1)//2)*C(k//2)*(2*(k//2)+1)*binomial(n-k, (n-k)//2) for k in range(n+1))
print([a(n) for n in range(26)]) # Mélika Tebni, Dec 06 2024
CROSSREFS
KEYWORD
nonn,walk
AUTHOR
N. J. A. Sloane, Apr 09 2018
EXTENSIONS
a(13)-a(25) from Mélika Tebni, Dec 06 2024
STATUS
approved