cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A065097 a(n) = ((2n+1) + (2n-1) - 1)!/((2n+1)!*(2n-1)!).

Original entry on oeis.org

1, 1, 7, 66, 715, 8398, 104006, 1337220, 17678835, 238819350, 3282060210, 45741281820, 644952073662, 9183676536076, 131873975875180, 1907493251046152, 27767032438524099, 406472021074865382, 5979899192930226746, 88366931393503350700, 1311063521138246054410
Offset: 0

Views

Author

Len Smiley, Nov 11 2001

Keywords

Comments

A Catalan-like formula using consecutive odd numbers. Recall that Catalan numbers (A000108) are given by ((n+1)+(n)-1)!/((n+1)!(n)!).
From David Callan, Jun 01 2006: (Start)
a(n) = number of Dyck (2n)-paths (i.e., semilength = 2n) all of whose interior returns to ground level (if any) occur at or before the (2n-2)-nd step, that is, they occur strictly before the midpoint of the path.
For example, a(2)=7 counts UUUUDDDD, UUUDUDDD, UUDUUDDD, UUDUDUDD, UUUDDUDD, UD.UUUDDD, UD.UUDUDD ("." denotes an interior return to ground level).
This result follows immediately from an involution on Dyck paths, due to Emeric Deutsch, defined by E->E, UPDQ -> UQDP (where E is the empty Dyck path; U=upstep, D=downstep and P,Q are arbitrary Dyck paths), because the involution is fixed-point-free on Dyck (2n)-paths and contains one path of the type being counted in each orbit.
a(n) = Sum_{k=0..n-1} C(2n-1-2k)*C(2k). This identity has the following combinatorial interpretation:
a(n) is the number of odd-GL-marked Dyck (2n-1)-paths. An odd-GL vertex is a vertex at location (2i,0) for some odd i >= 1 (path starts at origin). An odd-GL-marked Dyck path is a Dyck path with one of its odd-GL vertices marked. For example, a(2)=7 counts UUUDDD*, UUDUDD*, UD*UUDD, UDUUDD*, UD*UDUD, UDUDUD*, UUDDUD* (the * denotes the marked odd-GL vertex). (End)
a(n+1) = Sum_{k=0..n} C(k)*C(2*n+1-k), n >= 0, with C(n) = A000108(n), also gives the odd part of the bisection of the half-convolution of the Catalan sequence A000108 with itself. For the definition of the half-convolution of a sequence with itself see a comment on A201204. There one also finds the rule for the o.g.f. given below in the formula section. The even part of this bisection is found under A201205. - Wolfdieter Lang, Jan 05 2012
From Peter Bala, Dec 01 2015: (Start)
Let x = p/q be a positive rational in reduced form with p,q > 0. Define Cat(x) = (1/(2*p + q))*binomial(2*p + q, p). Then Cat(n) = Catalan(n). This sequence is Cat(n + 1/2) = (1/(4*n + 4))*binomial(4*n + 4, 2*n + 1). Cf. A265101 (Cat(n + 1/3)), A265102 (Cat(n + 1/4)) and A265103 (Cat(n + 1/5)).
Number of maximal faces of the rational associahedron Ass(2*n + 1, 2*n + 3). Number of lattice paths from (0, 0) to (2*n + 3, 2*n + 1) using steps of the form (1, 0) and (0, 1) and staying above the line y = (2*n + 1)/(2*n + 3)*x. See Armstrong et al. (End)
Also the number of ordered rooted trees with 2n nodes, most of which are leaves, i.e., the odd bisection of A358585. This follows from Callan's formula below. - Gus Wiseman, Nov 27 2022

Examples

			G.f.: 1 + x + 7*x^2 + 66*x^3 + 715*x^4 + 8398*x^5 + 104006*x^6 + ...
		

Crossrefs

Cf. A003150 (for analog with consecutive Fibonacci numbers).

Programs

  • Magma
    [Binomial(4*n-1, 2*n-1)/(2*n+1): n in [1..20]]; // Vincenzo Librandi, Dec 09 2015
  • Maple
    seq(binomial(4*n-1,2*n-1)/(2*n+1), n=0..30); # Robert Israel, Dec 08 2015
  • Mathematica
    a[ n_] := If[ n < 1, 0, Binomial[ 4 n - 1, 2 n - 1] / (2 n + 1)]; (* Michael Somos, Oct 25 2014 *)
  • MuPAD
    combinat::dyckWords::count(2*n)/2 $ n = 1..26 // Zerinvary Lajos, Apr 25 2007
    
  • PARI
    a(n) = { if(n==0, 1, (4*n - 1)!/((2*n + 1)!*(2*n - 1)!)) } \\ Harry J. Smith, Oct 07 2009
    
  • PARI
    vector(20, n, binomial(4*n-1, 2*n-1)/(2*n+1)) \\ Altug Alkan, Dec 08 2015
    
  • Sage
    A065097 = lambda n: hypergeometric([1-2*n,-2*n],[2],1)/2
    [Integer(A065097(n).n(500)) for n in (1..20)] # Peter Luschny, Sep 22 2014
    

Formula

a(n) = binomial(4*n-1, 2*n-1)/(2*n+1).
a(n) = C(2n)/2 where C(n) is the Catalan number A000108. - David Callan, Jun 01 2006
G.f.: 1/2 + (sqrt(2)/2)/sqrt(1+sqrt(1-16*x)). - Vladeta Jovovic, Sep 26 2003
G.f.: 1 + 3F2([1, 5/4, 7/4], [2, 5/2], 16*x). - Olivier Gérard, Feb 16 2011
O.g.f.: (1 + (cata(sqrt(x)) + cata(-sqrt(x)))/2)/2, with the o.g.f. cata(x) of the Catalan numbers. See the W. Lang comment above. - Wolfdieter Lang, Jan 05 2012
a(n) = hypergeometric([1-2*n,-2*n],[2],1)/2. - Peter Luschny, Sep 22 2014
a(n) = A001448(n) / (4*n + 2) if n>0. - Michael Somos, Oct 25 2014
n*(2*n+1)*a(n) - 2*(4*n-1)*(4*n-3)*a(n-1) = 0. - R. J. Mathar, Oct 31 2015
O.g.f. is 1 + Revert( x*(1 + x)/(1 + 2*x)^4 ). - Peter Bala, Dec 01 2015
Sum_{n>=0} 1/a(n) = 39/25 + 4*Pi/(9*sqrt(3)) - 24*log(phi)/(25*sqrt(5)), where phi is the golden ratio (A001622). - Amiram Eldar, Mar 02 2023
From Peter Bala, Apr 29 2024: (Start)
For n >= 1, a(n) = (1/8)*Sum_{k = 0..2*n-1} (-1)^k * 4^(2*n-k)*binomial(2*n-1, k)*Catalan(k+1).
For n >= 1, a(n) = (1/8)*(16^n)*hypergeom([1 - 2*n, 3/2], [3], 1). (End)

Extensions

a(0)=1 prepended by Alois P. Heinz, Nov 28 2021

A028364 Triangle T(n,m) = Sum_{k=0..m} Catalan(n-k)*Catalan(k).

Original entry on oeis.org

1, 1, 2, 2, 3, 5, 5, 7, 9, 14, 14, 19, 23, 28, 42, 42, 56, 66, 76, 90, 132, 132, 174, 202, 227, 255, 297, 429, 429, 561, 645, 715, 785, 869, 1001, 1430, 1430, 1859, 2123, 2333, 2529, 2739, 3003, 3432, 4862, 4862, 6292, 7150, 7810, 8398, 8986, 9646, 10504, 11934, 16796
Offset: 0

Views

Author

Keywords

Comments

There are several versions of a Catalan triangle: see A009766, A008315, A028364.
The subtriangle [1], [2, 3], [5, 7, 9], ..., namely T(N,M-1), for N >= 1, M=1..N, appears as one-point function in the totally asymmetric exclusion process for the parameters alpha=1=beta. See the Derrida et al. and Liggett references given under A067323, where these triangle entries are called T_{N,N+M-1} for the given alpha and beta values. See the row reversed triangle A067323.
Consider a Dyck path as a path with steps N=(0,1) and E=(1,0) from (0,0) to (n,n) that stays weakly above y=x. T(n,m) is the number of Dyck paths of semilength n+1 where the (m+1)st north step is followed by an east step. - Lara Pudwell, Apr 12 2023

Examples

			Triangle begins
   1;
   1,  2;
   2,  3,  5;
   5,  7,  9, 14;
  14, 19, 23, 28, 42;
		

Crossrefs

Cf. A000108 (column 0 and main diagonal), A001700 (row sums), A065097 (T(2*n-1, n-1)), A201205 (central terms).

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, add(
          expand(b(n-1, j)*`if`(i>n, x, 1)), j=1..i))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b((n+1)$2)):
    seq(T(n), n=0..10);  # Alois P. Heinz, Nov 28 2015
  • Mathematica
    t[n_, k_] = Sum[CatalanNumber[n-j]*CatalanNumber[j], {j, 0, k}]; Flatten[Table[t[n, k], {n, 0, 8}, {k, 0, n}]] (* Jean-François Alcover, Jul 22 2011 *)

Formula

T(n,k) = Sum_{j>=0} A039598(k,j)*A039599(n-k,j). - Philippe Deléham, Feb 18 2004
Sum_{k>=0} T(n,k) = A001700(n). T(n,k) = A067323(n,n-k), n >= k >= 0, otherwise 0. - Philippe Deléham, May 26 2005
G.f. for column sequences m >= 0: (-(c(m,x)-1)/x+c(m,x)*c(x))/x^m with the g.f. c(x) of A000108 (Catalan) and c(m,x):=sum(C(k)*x^k,k=0..m) with C(n):=A000108(n). - Wolfdieter Lang, Mar 24 2006
G.f. for column sequences m >= 0 (without leading zeros): c(x)*Sum_{k=0..m} C(m,k)*c(x)^k with the g.f. c(x) of A000108 (Catalan) and C(n,m) is the Catalan triangle A033184(n,m). - Wolfdieter Lang, Mar 24 2006
T(n,n) = T(n,k) + T(n,n-1-k) = A000108(n+1), n > 0, k = 0..floor((n+1)/2). - Yuchun Ji, Jan 09 2019
G.f. for triangle: Sum_{n>=0, m>=0} T(n, m)*x^n*y^m = (c(x)-c(xy))/(x(1-y)c(x)) with the g.f. c(x) of A000108 (Catalan). - Lara Pudwell, Apr 12 2023

A201204 Half-convolution of Catalan sequence A000108 with itself.

Original entry on oeis.org

1, 1, 3, 7, 23, 66, 227, 715, 2529, 8398, 30275, 104006, 380162, 1337220, 4939443, 17678835, 65844845, 238819350, 895451117, 3282060210, 12374186318, 45741281820, 173257703723, 644952073662, 2452607696798, 9183676536076, 35042725663002, 131873975875180, 504697422982484, 1907493251046152
Offset: 0

Views

Author

Wolfdieter Lang, Jan 02 2012

Keywords

Comments

In general the half-convolution of a sequence {b(n)}_0^infty with itself is defined by chat(n):=sum(b(k)*b(n-k), k=0..floor(n/2)), n>=0. The o.g.f. of the sequence {chat(n)} is obtained from the bisection 2*chat(2*k) - b(k)^2 = c(2*k), k>=0, with the ordinary convolution c(n):=sum(b(k)*b(n-k),k=0..n), n>=0, and 2*chat(2*k+1) = c(2*k+1), k>=0. This leads to the o.g.f.s for the corresponding even (e) and odd (o) parts:
2*Chate(x) - B2(x) = Ce(x) and 2*Chato(x) = Co(x), where Chate(x):= sum(chat(2*k)*x^k,k=0..infty), Chato(x):= sum(chat(2*k+1)*x^k,k=0..infty), B2(x) := sum(b(k)^2*x^k, k=0..infty), Ce(x) := sum(c(2*k)*x^k, k=0..infty) and Co(x) := sum(c(2*k+1)*x^k, k=0..infty). Thus Chate(x)=(Ce(x) + B2(x))/2 and Chato(x)=Co(x)/2. Expressing this in terms of C(x), the o.g.f. of {c(n)}, and B2(x) leads to the result: Chat(x)= (C(x) + B2(x^2))/2.
In the Catalan case b(n)=A000108(n), c(n)=b(n+1), C(x)= (cata(x)+1)/x, with the o.g.f. of A000108 cata(x)=(1-sqrt(1-4*x))/(2*x), and B2(x) is found under A001246 to be (-1 + hypergeom([-1/2,-1/2],[1],16*x))/(4*x). This produces the o.g.f. given in the formula section.
This computation was motivated by a question about the o.g.f. of A000992 ("half-Catalan numbers"). Note, however, that this sequence is not the half-convolution of the Catalan numbers presented here.
Apparently the number of hills to the left of or at the midpoint in all Dyck paths of semilength n+1. [David Scambler, Apr 30 2013]

Crossrefs

A000108, bisection: A201205 and A065097.

Programs

  • Maple
    C:= n -> binomial(2*n,n)/(n+1):
    A:= n -> add(C(k)*C(n-k),k=0..floor(n/2));
    seq(A(i),i=1..100); # Robert Israel, Jun 06 2014
  • Mathematica
    Table[Sum[CatalanNumber[k]CatalanNumber[n-k],{k,0,Floor[n/2]}],{n,0,30}] (* Harvey P. Dale, Jun 12 2012 *)
    Table[CatalanNumber[n + 1]/2 + 2^(2 n + 1) Binomial[1/2, n/2 + 1]^2, {n, 0, 30}] (* Vladimir Reshetnikov, Oct 03 2016 *)

Formula

a(n) = sum(Catalan(k)*Catalan(n-k),k=0..floor(n/2)), n>=0, with Catalan(n)=A000108(n).
O.g.f.: G(x)=(catalan(x)-1)/(2*x)+(-1+hypergeom([-1/2,-1/2],[1],16*x^2))/(8*x^2), with the o.g.f. catalan(x) of the Catalan numbers (see also the comment section).
a(n) ~ 2^(2*n+1) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Oct 15 2014
a(n) = A000108(n+1)/2 + 2^(2*n+1) * binomial(1/2, n/2+1)^2. - Vladimir Reshetnikov, Oct 03 2016
D-finite with recurrence: (n+1)*(n+2)^2*a(n) +6*(n-2)*(n+1)^2*a(n-1) +4*(-16*n^3+25*n^2+4*n-4)*a(n-2) +16*(-4*n^3+25*n^2-56*n+41)*a(n-3) +192*(4*n-7)*(n-3)^2*a(n-4) -256*(2*n-7)*(n-4)^2*a(n-5)=0. - R. J. Mathar, Feb 21 2020

A067323 Catalan triangle A028364 with row reversion.

Original entry on oeis.org

1, 2, 1, 5, 3, 2, 14, 9, 7, 5, 42, 28, 23, 19, 14, 132, 90, 76, 66, 56, 42, 429, 297, 255, 227, 202, 174, 132, 1430, 1001, 869, 785, 715, 645, 561, 429, 4862, 3432, 3003, 2739, 2529, 2333, 2123, 1859, 1430, 16796, 11934, 10504, 9646, 8986, 8398, 7810, 7150, 6292, 4862
Offset: 0

Views

Author

Wolfdieter Lang, Feb 05 2002

Keywords

Comments

a(N,p) equals X_{N}(N+1,p) := T_{N,p} for alpha= 1 =beta and N>=p>=1 in the Derrida et al. 1992 reference. The one-point correlation functions A000108(n)%20(Catalan)%20in%20this%20reference.%20See%20also%20the%20Derrida%20et%20al.%201993%20reference.%20In%20the%20Liggett%201999%20reference%20mu">{N} for alpha= 1 =beta equal a(N,K)/C(N+1) with C(n)=A000108(n) (Catalan) in this reference. See also the Derrida et al. 1993 reference. In the Liggett 1999 reference mu{N}{eta:eta(k)=1} of prop. 3.38, p. 275 is identical with _{N} and rho=0 and lambda=1.
Identity for each row n>=1: a(n,m)+a(n,n-m+1)= C(n+1), with C(n+1)=A000108(n+1)(Catalan) for every m=1..floor((n+1)/2). E.g., a(2k+1,k+1)=C(2*(k+1)).
The first column sequences (diagonals of A028364) are: A000108(n+1), A000245, A067324-6 for m=0..4.

Examples

			Triangle begins:
     1;
     2,    1;
     5,    3,    2;
    14,    9,    7,    5;
    42,   28,   23,   19,   14;
   132,   90,   76,   66,   56,   42;
   429,  297,  255,  227,  202,  174,  132;
  1430, 1001,  869,  785,  715,  645,  561,  429;
  4862, 3432, 3003, 2739, 2529, 2333, 2123, 1859, 1430;
  ...
		

References

  • B. Derrida, E. Domany and D. Mukamel, An exact solution of a one-dimensional asymmetric exclusion model with open boundaries, J. Stat. Phys. 69, 1992, 667-687; eqs. (19) - (23), p. 672.
  • B. Derrida, M. R. Evans, V. Hakim and V. Pasquier, Exact solution of a 1D asymmetric exclusion model using a matrix formulation, J. Phys. A 26, 1993, 1493-1517; eqs. (43), (44), pp. 1501-2 and eq.(81) with eqs.(80) and (81).
  • T. M. Liggett, Stochastic Interacting Systems: Contact, Voter and Exclusion Processes, Springer, 1999, pp. 269, 275.
  • G. Schuetz and E. Domany, Phase Transitions in an Exactly Soluble one-Dimensional Exclusion Process, J. Stat. Phys. 72 (1993) 277-295, eq. (2.18), p. 283, with eqs. (2.13)-(2.15).

Crossrefs

Cf. A001700 (row sums).
T(2n,n) gives A201205.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, add(
          expand(b(n-1, j)*`if`(i>n, x, 1)), j=1..i))
        end:
    T:= n-> (p-> seq(coeff(p, x, n-i), i=0..n))(b((n+1)$2)):
    seq(T(n), n=0..10);  # Alois P. Heinz, Nov 28 2015
  • Mathematica
    t[n_, k_] := Sum[ CatalanNumber[n - j]*CatalanNumber[j], {j, 0, k}]; Flatten[ Table[t[n, k], {n, 0, 9}, {k, n, 0, -1}]] (* Jean-François Alcover, Jul 17 2013 *)

Formula

a(n,m) = A028364(n,n-m), n>=m>=0, else 0.
G.f. for column m>=1 (without leading zeros): (c(x)^3)sum(C(m-1, k)*c(x)^k, k=0..m-1), with C(n, m) := (m+1)*binomial(2*n-m, n-m)/(n+1) (Catalan convolutions A033184); and for m=0: c^2(x), where c(x) is g.f. of A000108 (Catalan).
T(n,k) = Sum_{j>=0} A039598(n-k,j)*A039599(k,j). - Philippe Deléham, Feb 18 2004
G.f. for diagonal sequences: see g.f. for columns of A028364.
Showing 1-4 of 4 results.