OFFSET
2,1
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 2..100
Emeric Deutsch, Problem 10658: Another Type of Lattice Path, American Math. Monthly, 107, 2000, 368-370.
FORMULA
a(n) = Sum_{k=1..2n-1} k*A108443(n,k). Example: a(3) = 1*24 + 2*15 + 3*3 + 4*1 = 73.
G.f.: zA(2A^2-2zA^2-zA-2)/(1-2zA-3zA^2), where A=1+zA^2+zA^3 or, equivalently, A=(2/3)*sqrt((z+3)/z)*sin((1/3)*arcsin(sqrt(z)*(z+18)/(z+3)^(3/2)))-1/3 (the g.f. of A027307).
Recurrence: n*(2*n+1)*(40*n^5 - 100*n^4 - 758*n^3 + 3649*n^2 - 5474*n + 2727)*a(n) = (880*n^7 - 2200*n^6 - 15316*n^5 + 79354*n^4 - 145332*n^3 + 125379*n^2 - 48111*n + 5220)*a(n-1) + (n-3)*(2*n - 5)*(40*n^5 + 100*n^4 - 758*n^3 + 1175*n^2 - 650*n + 84)*a(n-2). - Vaclav Kotesovec, Mar 18 2014
a(n) ~ 5^(3/4) * ((11+5*sqrt(5))/2)^n / (10*sqrt(Pi*n)). - Vaclav Kotesovec, Mar 18 2014
EXAMPLE
a(2)=5 because in the ten paths udud, udUdd, uudd, uU(ddd), Uddud, UddUdd, Ududd, UdU(ddd), Uu(ddd) and UU(d[dd)d] (see A027307) we have 5 ddd's (shown between parentheses).
MAPLE
A:=(2/3)*sqrt((z+3)/z)*sin((1/3)*arcsin(sqrt(z)*(z+18)/(z+3)^(3/2)))-1/3: G:=z*A*(-z*A-2*z*A^2-2+2*A^2)/(1-3*z*A^2-2*z*A): Gser:=series(G, z=0, 26): seq(coeff(Gser, z^n), n=2..21);
CROSSREFS
KEYWORD
nonn
AUTHOR
Emeric Deutsch, Jun 10 2005
STATUS
approved