A175891 - OEIS

%I #17 Nov 11 2020 08:44:13

%S 1,1,5,29,185,1226,8553,61642,455337,3429002,26229691,203237747,

%T 1591820564,12582288455,100241042348,804090987555,6488942266564,

%U 52644171729304,429123506792664,3512829202462126,28866426741057006,238031465396515626,1969001793889730276

%N Number of lattice paths from (0,0) to (n,n) using steps S={(k,0),(0,k)|0<k<=4} which never go above the line y=x.

%H Alois P. Heinz, <a href="/A175891/b175891.txt">Table of n, a(n) for n = 0..1000</a>

%F a(n) ~ c * d^n / n^(3/2), where d = 8.84734830841870961487278801886633962039798... is the real root of the equation 4 + 4*d - 8*d^2 - 8*d^3 + d^4 = 0 and c = 0.31736815701423989167651891084531024477617724724822148387263881713... - _Vaclav Kotesovec_, May 30 2017

%p b:= proc(x, y) option remember; `if`(y>x or y<0, 0,

%p       `if`(x=0, 1, add(b(x-j, y)+b(x, y-j), j=1..4)))

%p     end:

%p a:= n-> b(n$2):

%p seq(a(n), n=0..35);  # _Alois P. Heinz_, May 16 2017

%t b[x_, y_] := b[x, y] = If[y > x || y < 0, 0, If[x == 0, 1, Sum[b[x - j, y] + b[x, y - j], {j, 1, 4}]]];

%t a[n_] := b[n, n];

%t a /@ Range[0, 35] (* _Jean-François Alcover_, Nov 11 2020, after _Alois P. Heinz_ *)

%Y Cf. A000103, A122951, A175883, A286918.

%K nonn

%O 0,3

%A _Eric Werley_, Dec 05 2010