cp's OEIS Frontend

A120730 Another version of Catalan triangle A009766.

%I A120730 #42 Sep 14 2024 06:52:06
%S A120730 1,0,1,0,1,1,0,0,2,1,0,0,2,3,1,0,0,0,5,4,1,0,0,0,5,9,5,1,0,0,0,0,14,
%T A120730 14,6,1,0,0,0,0,14,28,20,7,1,0,0,0,0,0,42,48,27,8,1,0,0,0,0,0,42,90,
%U A120730 75,35,9,1,0,0,0,0,0,0,132,165,110,44,10,1
%N A120730 Another version of Catalan triangle A009766.
%C A120730 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, -1, 0, 0, 1, -1, 0, 0, 1, -1, 0, 0, ...] DELTA [1, 0, 0, -1, 1, 0, 0, -1, 1, 0, 0, -1, 1, ...] where DELTA is the operator defined in A084938.
%C A120730 Aerated version gives A165408. - _Philippe Deléham_, Sep 22 2009
%C A120730 T(n,k) is the number of length n left factors of Dyck paths having k up steps. Example: T(5,4)=4 because we have UDUUU, UUDUU, UUUDU, and UUUUD, where U=(1,1) and D=(1,-1). - _Emeric Deutsch_, Jun 19 2011
%C A120730 With zeros omitted: 1,1,1,1,2,1,2,3,1,5,4,1,... = A008313. - _Philippe Deléham_, Nov 02 2011
%H A120730 Alois P. Heinz, <a href="/A120730/b120730.txt">Rows n = 0..200, flattened</a>
%F A120730 G.f.: G(t,z) = 4*z/((2*z-1+sqrt(1-4*t*z^2))*(1+sqrt(1-4*t*z^2))). - _Emeric Deutsch_, Jun 19 2011
%F A120730 Sum_{k=0..n} x^k*T(n,n-k) = A001405(n), A126087(n), A128386(n), A121724(n), A128387(n), A132373(n), A132374(n), A132375(n), A121725(n) for x=1,2,3,4,5,6,7,8,9 respectively. [corrected by _Philippe Deléham_, Oct 16 2008]
%F A120730 T(2*n,n) = A000108(n); A000108: Catalan numbers.
%F A120730 From _Philippe Deléham_, Oct 18 2008: (Start)
%F A120730 Sum_{k=0..n} T(n,k)^2 = A000108(n) and Sum_{n>=k} T(n,k) = A000108(k+1).
%F A120730 Sum_{k=0..n} T(n,k)^3 = A003161(n).
%F A120730 Sum_{k=0..n} T(n,k)^4 = A129123(n). (End)
%F A120730 Sum_{k=0..n}, T(n,k)*x^k = A000007(n), A001405(n), A151281(n), A151162(n), A151254(n), A156195(n), A156361(n), A156362(n), A156566(n), A156577(n) for x=0,1,2,3,4,5,6,7,8,9 respectively. - _Philippe Deléham_, Feb 10 2009
%F A120730 From _G. C. Greubel_, Nov 07 2022: (Start)
%F A120730 T(n, k) = 0 if n > 2*k, otherwise binomial(n, k)*(2*k-n+1)/(k+1).
%F A120730 Sum_{k=0..n} (-1)^k*T(n,k) = A105523(n).
%F A120730 Sum_{k=0..n} (-1)^k*T(n,k)^2 = -A132889(n), n >= 1.
%F A120730 Sum_{k=0..floor(n/2)} T(n-k, k) = A357654(n).
%F A120730 T(n, n-1) = A001477(n).
%F A120730 T(n, n-2) = [n=2] + A000096(n-3), n >= 2.
%F A120730 T(n, n-3) = 2*[n<5] + A005586(n-5), n >= 3.
%F A120730 T(n, n-4) = 5*[n<7] - 2*[n=4] + A005587(n-7), n >= 4.
%F A120730 T(2*n+1, n+1) = A000108(n+1), n >= 0.
%F A120730 T(2*n-1, n+1) = A099376(n-1), n >= 1. (End)
%e A120730 As a triangle, this begins:
%e A120730   1;
%e A120730   0,  1;
%e A120730   0,  1,  1;
%e A120730   0,  0,  2,  1;
%e A120730   0,  0,  2,  3,  1;
%e A120730   0,  0,  0,  5,  4,  1;
%e A120730   0,  0,  0,  5,  9,  5,  1;
%e A120730   0,  0,  0,  0, 14, 14,  6,  1;
%e A120730   ...
%p A120730 G := 4*z/((2*z-1+sqrt(1-4*z^2*t))*(1+sqrt(1-4*z^2*t))): Gser := simplify(series(G, z = 0, 13)): for n from 0 to 12 do P[n] := sort(coeff(Gser, z, n)) end do: for n from 0 to 12 do seq(coeff(P[n], t, k), k = 0 .. n) end do; # yields sequence in triangular form  # _Emeric Deutsch_, Jun 19 2011
%p A120730 # second Maple program:
%p A120730 b:= proc(x, y) option remember; `if`(y<0 or y>x, 0,
%p A120730      `if`(x=0, 1, add(b(x-1, y+j), j=[-1, 1])))
%p A120730     end:
%p A120730 T:= (n, k)-> b(n, 2*k-n):
%p A120730 seq(seq(T(n, k), k=0..n), n=0..14);  # _Alois P. Heinz_, Oct 13 2022
%t A120730 b[x_, y_]:= b[x, y]= If[y<0 || y>x, 0, If[x==0, 1, Sum[b[x-1, y+j], {j, {-1, 1}}] ]];
%t A120730 T[n_, k_] := b[n, 2 k - n];
%t A120730 Table[Table[T[n, k], {k, 0, n}], {n, 0, 14}] // Flatten (* _Jean-François Alcover_, Oct 21 2022, after _Alois P. Heinz_ *)
%t A120730 T[n_, k_]:= If[n>2*k, 0, Binomial[n, k]*(2*k-n+1)/(k+1)];
%t A120730 Table[T[n, k], {n,0,13}, {k,0,n}]//Flatten (* _G. C. Greubel_, Nov 07 2022 *)
%o A120730 (Magma)
%o A120730 A120730:= func< n,k | n gt 2*k select 0 else Binomial(n, k)*(2*k-n+1)/(k+1) >;
%o A120730 [A120730(n,k): k in [0..n], n in [0..13]]; // _G. C. Greubel_, Nov 07 2022
%o A120730 (SageMath)
%o A120730 def A120730(n,k): return 0 if (n>2*k) else binomial(n, k)*(2*k-n+1)/(k+1)
%o A120730 flatten([[A120730(n,k) for k in range(n+1)] for n in range(14)]) # _G. C. Greubel_, Nov 07 2022
%Y A120730 Cf. A000007, A000096, A000108, A001405, A001477, A003161, A005586, A005587.
%Y A120730 Cf. A008313, A009766, A039598, A039599, A084938, A099376, A105523, A121725.
%Y A120730 Cf. A126087, A128386, A121724, A128387, A129123, A132373, A132374, A132375.
%Y A120730 Cf. A132889, A151162, A151254, A151281, A156195, A156361, A156362, A156566.
%Y A120730 Cf. A156577, A165408, A357654.
%K A120730 nonn,tabl
%O A120730 0,9
%A A120730 _Philippe Deléham_, Aug 17 2006, corrected Sep 15 2006