A118920 - OEIS

%I #25 Oct 27 2021 11:26:42

%S 2,4,2,10,8,2,28,28,12,2,84,96,54,16,2,264,330,220,88,20,2,858,1144,

%T 858,416,130,24,2,2860,4004,3276,1820,700,180,28,2,9724,14144,12376,

%U 7616,3400,1088,238,32,2,33592,50388,46512,31008,15504,5814,1596,304,36,2

%N Triangle read by rows: T(n,k) is the number of Grand Dyck paths of semilength n that cross the x-axis k times (n>=1, k>=0).

%C A Grand Dyck path of semilength n is a path in the half-plane x>=0, starting at (0,0), ending at (2n,0) and consisting of steps u=(1,1) and d=(1,-1).

%C Row sums are the central binomial coefficients (A000984). T(n,0)=2*A000108(n) (the Catalan numbers doubled). T(n,1)=2*A002057(n-2). Sum(k*T(n,k),k>=0)=2*A008549(n-1). For crossings of the x-axis in one direction, see A118919.

%C This triangle is related to paired pattern P_3 and P_4 defined in the Pan & Remmel link. - _Ran Pan_, Feb 01 2016

%H Ran Pan and Jeffrey B. Remmel, <a href="http://arxiv.org/abs/1601.07988">Paired patterns in lattice paths</a>, arXiv:1601.07988 [math.CO], 2016.

%F T(n,k) = 2*(k+1)*binomial(2*n,n-k-1)/n.

%F G.f.: G(t,z)=2*z*C^2/(1-t*z*C^2), where C=(1-sqrt(1-4*z))/(2*z) is the Catalan function.

%F More generally, the trivariate g.f. G=G(x,y,z), where x (y) marks number of downward (upward) crossings of the x-axis, is given by G = z*C^2*(2+(x+y)*z*C^2)/(1-x*y*z^2*C^4).

%F a(n) = 2 * A039598(n-1). - _Georg Fischer_, Oct 27 2021

%e T(3,1)=8 because we have ud|dudu,ud|dduu,udud|du,uudd|du,du|udud,du|uudd, dudu|ud and dduu|ud (the crossings of the x-axis are shown by |).

%e Triangle starts:

%e    2;

%e    4, 2;

%e   10, 8, 2;

%e   28,28,12, 2;

%e   84,96,54,16,2;

%p T:=(n,k)->2*(k+1)*binomial(2*n,n-k-1)/n: for n from 1 to 11 do seq(T(n,k),k=0..n-1) od; # yields sequence in triangular form

%t Table[2 (k + 1) Binomial[2 n, n - k - 1]/n, {n, 10}, {k, 0, n - 1}] // Flatten (* _Michael De Vlieger_, Feb 01 2016 *)

%o (Sage) # Algorithm of L. Seidel (1877)

%o # Prints the first n rows of the triangle.

%o def A118920_triangle(n) :

%o     D = [0]*(n+2); D[1] = 2

%o     b = True; h = 1

%o     for i in range(2*n) :

%o         if b :

%o             for k in range(h, 0, -1) : D[k] += D[k-1]

%o             h += 1

%o         else :

%o             for k in range(1, h, 1) : D[k] += D[k+1]

%o         b = not b

%o         if b : print([D[z] for z in (1..h-1)])

%o A118920_triangle(10)  # _Peter Luschny_, Oct 19 2012

%o (PARI) T(n,k) = 2*(k+1)*binomial(2*n,n-k-1)/n \\ _Charles R Greathouse IV_, Feb 01 2016

%Y Cf. A000984, A000108, A002057, A008549, A039598, A118919.

%K nonn,tabl

%O 1,1

%A _Emeric Deutsch_, May 06 2006