

A065097


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


12



1, 7, 66, 715, 8398, 104006, 1337220, 17678835, 238819350, 3282060210, 45741281820, 644952073662, 9183676536076, 131873975875180, 1907493251046152, 27767032438524099, 406472021074865382, 5979899192930226746, 88366931393503350700, 1311063521138246054410
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

A Catalanlike 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 (2n2)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 fixedpointfree on Dyck (2n)paths and contains one path of the type being counted in each orbit.
a(n) = Sum_{k=0..n1}C(2n12k)C(2k). This identity has the following combinatorial interpretation:
a(n) = number of oddGLmarked Dyck (2n1)paths. An oddGL vertex is a vertex at location (2i,0) for some odd i>=1 (path starts at origin). An oddGLmarked Dyck path is a Dyck path with one of its oddGL vertices marked. For example, a(2)=7 counts UUUDDD*, UUDUDD*, UD*UUDD, UDUUDD*, UD*UDUD, UDUDUD*, UUDDUD* (the * denotes the marked oddGL vertex). (End)
a(n+1) = Sum_{k=0..n}C(k)*C(2*n+1k), n>=0, with C(n) = A000108(n), gives also the odd part of the bisection of the halfconvolution of the Catalan sequence A000108 with itself. For the definition of the halfconvolution 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/4)).
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)


LINKS

Harry J. Smith, Table of n, a(n) for n=1..100
D. Armstrong, B. Rhoades, and N. Williams, Rational associahedra and noncrossing partitions arxiv:1305.7286v1 [math.CO], 2008.
E. Deutsch, An involution on Dyck paths and its consequences, Discrete Math., 204 (1999), no. 13, 163166.


FORMULA

a(n) = binomial(4*n1, 2*n1)/(2*n+1).
a(n) = C(2n)/2 where C(n) is the Catalan number A000108.  David Callan, Jun 01 2006
G.f.: (sqrt(2)/2)/sqrt(1+sqrt(116*x))1/2.  Vladeta Jovovic, Sep 26 2003
G.f.: 3F2( (1, 5/4, 7/4); (2, 5/2) )(16*x).  Olivier Gérard, Feb 16 2011
O.g.f.: ((cata(sqrt(x)) + cata(sqrt(x)))/2 1)/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([12*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*n1)*(4*n3)*a(n1)=0.  R. J. Mathar, Oct 31 2015
O.g.f. is Revert( x*(1 + x)/(1 + 2*x)^4 ).  Peter Bala, Dec 01 2015


EXAMPLE

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


MAPLE

seq(binomial(4*n1, 2*n1)/(2*n+1), n=1..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 *)


PROG

(MuPAD) combinat::dyckWords::count(2*n)/2 $ n = 1..26 // Zerinvary Lajos, Apr 25 2007
(PARI) { for (n=1, 100, a=(4*n  1)!/((2*n + 1)!*(2*n  1)!); write("b065097.txt", n, " ", a) ) } \\ Harry J. Smith, Oct 07 2009
(PARI) vector(100, n, binomial(4*n1, 2*n1)/(2*n+1)) \\ Altug Alkan, Dec 08 2015
(Sage)
A065097 = lambda n: hypergeometric([12*n, 2*n], [2], 1)/2
[Integer(A065097(n).n(500)) for n in (1..20)] # Peter Luschny, Sep 22 2014
(MAGMA) [Binomial(4*n1, 2*n1)/(2*n+1): n in [1..20]]; // Vincenzo Librandi, Dec 09 2015


CROSSREFS

Cf. A003150 (for analog with consecutive Fibonacci numbers).
Equals (1/2) A048990. Cf. A001795, A201204, A201205, A265101, A265102, A265103.
Sequence in context: A274822 A179880 A297310 * A300991 A122705 A185181
Adjacent sequences: A065094 A065095 A065096 * A065098 A065099 A065100


KEYWORD

nonn,easy


AUTHOR

Len Smiley, Nov 11 2001


STATUS

approved



