A010736 Let S(x,y) = number of lattice paths from (0,0) to (x,y) that use the step set { (0,1), (1,0), (2,0), (3,0), ....} and never pass below y = x. Sequence gives S(n-2,n).

1, 3, 12, 52, 237, 1119, 5424, 26832, 134913, 687443, 3541932, 18421524, 96585597, 509960223, 2709067968, 14469453632, 77655751329, 418567792899, 2264867271852, 12298297439892, 66993811842477, 366009125766463

Offset

0,2

Comments

Threefold convolution of A001003 with itself. Number of dissections of a convex polygon with n+4 sides that have a quadrilateral over a fixed side (the base) of the polygon. Example: a(1)=3 because the only dissections of the convex pentagon ABCDE (AB being the base), that have a quadrilateral over AB are the dissections made by the diagonals EC, AD and BD, respectively. - Emeric Deutsch, Dec 27 2003

a(n-1) = number of royal paths (A006318) from (0,0) to (n,n) with exactly 2 diagonal steps on the line y=x. - David Callan, Jul 15 2004

Links

Vincenzo Librandi, Table of n, a(n) for n = 0..1000

Formula

G.f.: (1+z-sqrt(1-6*z+z^2))^3/(64*z^3). - Emeric Deutsch, Dec 27 2003

a(n) = (3/n)*Sum_{k = 1..n} binomial(n, k)*binomial(n+k+2, k-1) = 3*hypergeom([1-n, n+4], [2], -1), n>=1, a(0)=1.

Recurrence: (n+3)*(2*n-1)*a(n) = (12*n^2+11*n-11)*a(n-1) - (n-3)*(2*n-1)*a(n-2) + (3-n)*a(n-3). - Vaclav Kotesovec, Oct 07 2012

a(n) ~ 3 * (1 + sqrt(2))^(2*n+3) / (2^(11/4) * sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Oct 07 2012, simplified Dec 24 2017

From Peter Bala, Jul 31 2024: (Start)

a(n) = (3/4) * Sum_{k = 0..n+1} binomial(n+1, k)*binomial(n+k+1, k)/(k+2) for n >= 1.

Second-order recurrence: (n + 3)*n^2*a(n) = (2*n + 1)*(3*n^2 + 3*n - 2)*a(n-1) - (n - 2)*(n + 1)^2*a(n-2), with a(0) = 1 and a(1) = 3.

a(n) = (3/8) * hypergeom([2, n + 2, - n - 1], [1, 3], -1) for n >= 1.

a(n) = (3/4) * Integral_{x = 0..1} x*Legendre_P(n+1, 2*x+1) for n >= 1. Note that A006318(n) = Integral_{x = 0..1} Legendre_P(n, 2*x+1). (End)

Mathematica

f[ x_, y_ ] := f[ x, y ] = Module[ {return}, If[ x == 0, return = 1, If[ y == x-1, return = 0, return = f[ x, y-1 ] + Sum[ f[ k, y ], {k, 0, x-1} ] ] ]; return ]; Do[ Print[ Table[ f[ k, j ], {k, 0, j} ] ], {j, 10, 0, -1} ]
CoefficientList[Series[(1+x-Sqrt[1-6x+x^2])^3/(64x^3), {x, 0, 30}], x] (* Harvey P. Dale, Apr 18 2012 *)

Prog

(PARI) my(x='x+O('x^66)); Vec((1+x-sqrt(1-6*x+x^2))^3/(64*x^3)) \\ Joerg Arndt, May 04 2013

Crossrefs

Cf. A001003, A006318.

Right-hand column 3 of triangle A011117.

Third column of convolution triangle A011117.

Sequence in context: A151195 A151196 A363090 * A389436 A151197 A007198

Adjacent sequences: A010733 A010734 A010735 * A010737 A010738 A010739

Keyword

nonn,easy

Author

Robert Sulanke (sulanke(AT)diamond.idbsu.edu)

Extensions

More terms from Emeric Deutsch, Dec 27 2003

Status

approved