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.

Previous Showing 11-20 of 202 results. Next

A152250 Eigentriangle, row sums = A001850, the Delannoy numbers.

Original entry on oeis.org

1, 2, 1, 8, 2, 3, 36, 8, 6, 13, 172, 36, 24, 26, 63, 852, 172, 108, 104, 126, 321, 4324, 852, 516, 468, 504, 642, 1683, 22332, 4324, 2556, 2236, 2268, 2568, 3366, 8989, 116876, 22332, 12972, 11076, 10836, 11556, 13464, 17978, 48639
Offset: 0

Views

Author

Gary W. Adamson, Nov 30 2008

Keywords

Comments

Row sums = A001850, the Delannoy numbers: (1, 3, 13, 63, 321,...).
Sum of n-th row terms = rightmost term of next row.

Examples

			First few rows of the triangle =
  1;
  2, 1;
  8, 2, 3;
  36, 8, 6, 13;
  172, 36, 24, 26, 63;
  852, 172, 108, 104, 126, 321;
  4324, 852, 516, 468, 504, 642, 1683;
  22332, 4324, 2556, 2236, 2268, 2568, 3366, 8989;
  116876, 22332, 12972, 11076, 10836, 11556, 13464, 17978, 48639;
  ...
Row 3 = (36, 8, 6, 13) = termwise products of (36, 8, 2, 1) and (1, 1, 3, 13).
		

Crossrefs

Programs

  • Mathematica
    nmax = 8;
    T[0, 0] = 1;
    T[n_, 0] := SeriesCoefficient[1/(x + Sqrt[1 - 6x + x^2]), {x, 0, n}];
    T[n_, n_] :=  LegendreP[n - 1, 3];
    row[n_] := row[n] = Table[T[m, 0], {m, n, 0, -1}]*Table[T[m, m], {m, 0, n} ];
    T[n_, k_] /; 0 < k < n := row[n][[k + 1]];
    Table[T[n, k], {n, 0, nmax}, {k, 0, n}] // Flatten (* Jean-François Alcover, Aug 07 2018 *)

Formula

Triangle read by rows, M*Q. M = an infinite lower triangular matrix with A109980 in every column: (1, 2, 8, 36, 172,...); Q = a matrix with A001850 prefaced with a "1" as the main diagonal: (1, 1, 3, 13, 63, 321,...) and the rest zeros.

A330203 Composite numbers k such that D(k) == 3 (mod k), where D(k) is the k-th central Delannoy number (A001850).

Original entry on oeis.org

10, 15, 50, 370, 2418, 4371, 5341, 8430, 20535, 25338, 26958, 278674, 1194649, 4304445, 11984885, 12327121, 20746461, 27585010, 72363853, 79501818
Offset: 1

Views

Author

Amiram Eldar, Dec 05 2019

Keywords

Comments

Equivalently, composite numbers k such that P(k, 3) == 3 (mod k), where P(k, 3) = D(k) is the k-th Legendre polynomial evaluated at 3.
P(p, 3) == 3 (mod p) for all primes p. This is a special case of Schur congruences, named after Issai Schur, first published by his student Hildegard Ille in her Ph.D. thesis in 1924, and proven by Wahab in 1952. This sequence consists of the composite numbers for which the congruence holds.

Examples

			10 is in the sequence since it is composite and D(10) = 8097453 == 3 (mod 10).
		

References

  • Hildegard Ille, Zur Irreduzibilität der Kugelfunktionen, Jahrbuch der Dissertationen der Universität Berlin, (1924).
  • Peter S. Landweber, Elliptic Curves and Modular Forms in Algebraic Topology: Proceedings of a Conference held at the Institute for Advanced Study, Princeton, Sept. 15-17, 1986, Springer, 2006. See pp. 74-76.

Crossrefs

Programs

  • Mathematica
    Select[Range[2500], CompositeQ[#] && Divisible[LegendreP[#, 3] - 3, #] &]
  • Sage
    a, b = 1, 1
    for n in range(1, 10000):
        a, b = b, ((6*n-3)*b - (n-1)*a)//n
        if (b%n == 3) and (not Integer(n).is_prime()): print(n)  # Robin Visser, Aug 08 2023

Extensions

a(18) from Robin Visser, Aug 08 2023
a(19)-a(20) from Robin Visser, Sep 11 2023

A228192 a(n) = A001850(n^2), where A001850 forms the central Delannoy numbers.

Original entry on oeis.org

1, 3, 321, 1462563, 252055236609, 1569245074591690083, 345299757825442889707393857, 2653337188651000290233505189314055363, 706829476163540077094231781323762631545566527489, 6496844758902641761809431955916116052361210081093847336070467
Offset: 0

Views

Author

Paul D. Hanna, Aug 15 2013

Keywords

Examples

			L.g.f.: L(x) = 3*x + 321*x^2/2 + 1462563*x^3/3 + 252055236609*x^4/4 +...
where exponentiation yields the g.f. of A228193:
exp(L(x)) = 1 + 3*x + 165*x^2 + 488007*x^3 + 63015285321*x^4 + 313849204040245803*x^5 +...+ A228193(n)*X^n +...
		

Crossrefs

Programs

  • PARI
    {a(n)=sum(k=0,n^2,binomial(n^2,k)*binomial(n^2+k,k))}
    for(n=0,20,print1(a(n),", "))
    
  • PARI
    {A001850(n)=polcoeff(1/sqrt(1 - 6*x + x^2 + x*O(x^n)), n)}
    {a(n)=A001850(n^2)}
    for(n=0,20,print1(a(n),", "))

Formula

a(n) = Sum_{k=0..n^2} binomial(n^2, k) * binomial(n^2+k, k).
a(n) = [x^(n^2)] 1/sqrt(1 - 6*x + x^2).
Equals the logarithmic derivative of A228193, after ignoring the initial term.

A228193 G.f.: exp( Sum_{n>=1} A001850(n^2)*x^n/n ), where A001850 forms the central Delannoy numbers.

Original entry on oeis.org

1, 3, 165, 488007, 63015285321, 313849204040245803, 57549960579131376060801997, 379048169979935686476204047966170767, 88353684521579654155696728418892273040483607185, 721871639878336367921338532273490438662977816273231098545619
Offset: 0

Views

Author

Paul D. Hanna, Aug 15 2013

Keywords

Examples

			G.f.: A(x) = 1 + 3*x + 165*x^2 + 488007*x^3 + 63015285321*x^4 +...
where the logarithm of the g.f. begins:
log(A(x)) = 3*x + 321*x^2/2 + 1462563*x^3/3 + 252055236609*x^4/4 +...+ A001850(n^2)*x^n/n +...
		

Crossrefs

Programs

  • PARI
    {A228192(n)=sum(k=0,n^2,binomial(n^2,k)*binomial(n^2+k,k))}
    {a(n)=polcoeff(exp(sum(k=1,n+1,A228192(k)*x^k/k) +x*O(x^n)),n)}
    for(n=0,20,print1(a(n),", "))

Formula

Logarithmic derivative yields A228192.

A192674 Floor-Sqrt transform of large central Delannoy numbers (A001850).

Original entry on oeis.org

1, 1, 3, 7, 17, 41, 94, 220, 515, 1209, 2845, 6711, 15861, 37549, 89015, 211287, 502051, 1194083, 2842413, 6771220, 16141369, 38501639, 91888052, 219411132, 524154756, 1252695124, 2995034616, 7163349548, 17138661276, 41017945745, 98196952422, 235147807457, 563242204233
Offset: 0

Views

Author

Emanuele Munarini, Jul 07 2011

Keywords

Programs

  • Mathematica
    FSFromSeries[f_,x_,n_] := Map[Floor[Sqrt[#]]&,CoefficientList[Series[f,{x,0,n}],x]]
    FSFromSeries[1/Sqrt[1-6x+x^2],x,100]

Formula

a(n) = floor(sqrt(centralDelannoy(n))).

A358360 The 3-adic valuation of the central Delannoy numbers (sequence A001850).

Original entry on oeis.org

0, 1, 0, 2, 1, 2, 0, 1, 0, 3, 2, 3, 1, 2, 1, 3, 2, 3, 0, 1, 0, 2, 1, 2, 0, 1, 0, 4, 3, 4, 2, 3, 2, 4, 3, 4, 1, 2, 1, 3, 2, 3, 1, 2, 1, 4, 3, 4, 2, 3, 2, 4, 3, 4, 0, 1, 0, 2, 1, 2, 0, 1, 0, 3, 2, 3, 1, 2, 1, 3, 2, 3, 0, 1, 0, 2, 1, 2, 0, 1, 0, 5, 4, 5, 3, 4, 3
Offset: 0

Views

Author

Jeffrey Shallit, Nov 12 2022

Keywords

Comments

The 3-adic valuation of x is the exponent of the highest power of 3 dividing x.

References

  • J.-P. Allouche and J. Shallit, Automatic Sequences: Theory, Applications, Generalizations. Cambridge University Press (2003), p. 453.

Crossrefs

Programs

  • Mathematica
    Array[IntegerExponent[LegendreP[#, 3], 3] &, 105, 0] (* Michael De Vlieger, Nov 12 2022 *)

Formula

b(n) = b(floor(n/3)) + (floor(n/3) mod 2), if n == 0,2 (mod 3);
b(n) = b(floor(n/9)) + 1, if n == 1 (mod 3).
a(n) = A007949(A001850(n)).

A001006 Motzkin numbers: number of ways of drawing any number of nonintersecting chords joining n (labeled) points on a circle.

Original entry on oeis.org

1, 1, 2, 4, 9, 21, 51, 127, 323, 835, 2188, 5798, 15511, 41835, 113634, 310572, 853467, 2356779, 6536382, 18199284, 50852019, 142547559, 400763223, 1129760415, 3192727797, 9043402501, 25669818476, 73007772802, 208023278209, 593742784829, 1697385471211
Offset: 0

Views

Author

Keywords

Comments

Number of 4321-, (3412,2413)-, (3412,3142)- and 3412-avoiding involutions in S_n.
Number of sequences of length n-1 consisting of positive integers such that the first and last elements are 1 or 2 and the absolute difference between any 2 consecutive elements is 0 or 1. - Jon Perry, Sep 04 2003
From David Callan, Jul 15 2004: (Start)
Also number of Motzkin n-paths: paths from (0,0) to (n,0) in an n X n grid using only steps U = (1,1), F = (1,0) and D = (1,-1).
Number of Dyck n-paths with no UUU. (Given such a Dyck n-path, change each UUD to U, then change each remaining UD to F. This is a bijection to Motzkin n-paths. Example with n=5: U U D U D U U D D D -> U F U D D.)
Number of Dyck (n+1)-paths with no UDU. (Given such a Dyck (n+1)-path, mark each U that is followed by a D and each D that is not followed by a U. Then change each unmarked U whose matching D is marked to an F. Lastly, delete all the marked steps. This is a bijection to Motzkin n-paths. Example with n=6 and marked steps in small type: U U u d D U U u d d d D u d -> U U u d D F F u d d d D u d -> U U D F F D.) (End)
a(n) is the number of strings of length 2n+2 from the following recursively defined set: L contains the empty string and, for any strings a and b in L, we also find (ab) in L. The first few elements of L are e, (), (()), ((())), (()()), (((()))), ((()())), ((())()), (()(())) and so on. This proves that a(n) is less than or equal to C(n), the n-th Catalan number. See Orrick link (2024). - Saul Schleimer (saulsch(AT)math.rutgers.edu), Feb 23 2006 (Additional linked comment added by William P. Orrick, Jun 13 2024.)
a(n) = number of Dyck n-paths all of whose valleys have even x-coordinate (when path starts at origin). For example, T(4,2)=3 counts UDUDUUDD, UDUUDDUD, UUDDUDUD. Given such a path, split it into n subpaths of length 2 and transform UU->U, DD->D, UD->F (there will be no DUs for that would entail a valley with odd x-coordinate). This is a bijection to Motzkin n-paths. - David Callan, Jun 07 2006
Also the number of standard Young tableaux of height <= 3. - Mike Zabrocki, Mar 24 2007
a(n) is the number of RNA shapes of size 2n+2. RNA Shapes are essentially Dyck words without "directly nested" motifs of the form A[[B]]C, for A, B and C Dyck words. The first RNA Shapes are []; [][]; [][][], [[][]]; [][][][], [][[][]], [[][][]], [[][]][]; ... - Yann Ponty (ponty(AT)lri.fr), May 30 2007
The sequence is self-generated from top row A going to the left starting (1,1) and bottom row = B, the same sequence but starting (0,1) and going to the right. Take dot product of A and B and add the result to n-th term of A to get the (n+1)-th term of A. Example: a(5) = 21 as follows: Take dot product of A = (9, 4, 2, 1, 1) and (0, 1, 1, 2, 4) = (0, + 4 + 2 + 2 + 4) = 12; which is added to 9 = 21. - Gary W. Adamson, Oct 27 2008
Equals A005773 / A005773 shifted (i.e., (1,2,5,13,35,96,...) / (1,1,2,5,13,35,96,...)). - Gary W. Adamson, Dec 21 2008
Starting with offset 1 = iterates of M * [1,1,0,0,0,...], where M = a tridiagonal matrix with [0,1,1,1,...] in the main diagonal and [1,1,1,...] in the super and subdiagonals. - Gary W. Adamson, Jan 07 2009
a(n) is the number of involutions of {1,2,...,n} having genus 0. The genus g(p) of a permutation p of {1,2,...,n} is defined by g(p)=(1/2)[n+1-z(p)-z(cp')], where p' is the inverse permutation of p, c = 234...n1 = (1,2,...,n), and z(q) is the number of cycles of the permutation q. Example: a(4)=9; indeed, p=3412=(13)(24) is the only involution of {1,2,3,4} with genus > 0. This follows easily from the fact that a permutation p of {1,2,...,n} has genus 0 if and only if the cycle decomposition of p gives a noncrossing partition of {1,2,...,n} and each cycle of p is increasing (see Lemma 2.1 of the Dulucq-Simion reference). [Also, redundantly, for p=3412=(13)(24) we have cp'=2341*3412=4123=(1432) and so g(p)=(1/2)(4+1-2-1)=1.] - Emeric Deutsch, May 29 2010
Let w(i,j,n) denote walks in N^2 which satisfy the multivariate recurrence w(i,j,n) = w(i, j + 1, n - 1) + w(i - 1, j, n - 1) + w(i + 1, j - 1, n - 1) with boundary conditions w(0,0,0) = 1 and w(i,j,n) = 0 if i or j or n is < 0. Then a(n) = Sum_{i = 0..n, j = 0..n} w(i,j,n) is the number of such walks of length n. - Peter Luschny, May 21 2011
a(n)/a(n-1) tends to 3.0 as N->infinity: (1+2*cos(2*Pi/N)) relating to longest odd N regular polygon diagonals, by way of example, N=7: Using the tridiagonal generator [cf. comment of Jan 07 2009], for polygon N=7, we extract an (N-1)/2 = 3 X 3 matrix, [0,1,0; 1,1,1; 0,1,1] with an e-val of 2.24697...; the longest Heptagon diagonal with edge = 1. As N tends to infinity, the diagonal lengths tend to 3.0, the convergent of the sequence. - Gary W. Adamson, Jun 08 2011
Number of (n+1)-length permutations avoiding the pattern 132 and the dotted pattern 23\dot{1}. - Jean-Luc Baril, Mar 07 2012
Number of n-length words w over alphabet {a,b,c} such that for every prefix z of w we have #(z,a) >= #(z,b) >= #(z,c), where #(z,x) counts the letters x in word z. The a(4) = 9 words are: aaaa, aaab, aaba, abaa, aabb, abab, aabc, abac, abca. - Alois P. Heinz, May 26 2012
Number of length-n restricted growth strings (RGS) [r(1), r(2), ..., r(n)] such that r(1)=1, r(k)<=k, and r(k)!=r(k-1); for example, the 9 RGS for n=4 are 1010, 1012, 1201, 1210, 1212, 1230, 1231, 1232, 1234. - Joerg Arndt, Apr 16 2013
Number of length-n restricted growth strings (RGS) [r(1), r(2), ..., r(n)] such that r(1)=0, r(k)<=k and r(k)-r(k-1) != 1; for example, the 9 RGS for n=4 are 0000, 0002, 0003, 0004, 0022, 0024, 0033, 0222, 0224. - Joerg Arndt, Apr 17 2013
Number of (4231,5276143)-avoiding involutions in S_n. - Alexander Burstein, Mar 05 2014
a(n) is the number of increasing unary-binary trees with n nodes that have an associated permutation that avoids 132. For more information about unary-binary trees with associated permutations, see A245888. - Manda Riehl, Aug 07 2014
a(n) is the number of involutions on [n] avoiding the single pattern p, where p is any one of the 8 (classical) patterns 1234, 1243, 1432, 2134, 2143, 3214, 3412, 4321. Also, number of (3412,2413)-, (3412,3142)-, (3412,2413,3142)-avoiding involutions on [n] because each of these 3 sets actually coincides with the 3412-avoiding involutions on [n]. This is a complete list of the 8 singles, 2 pairs, and 1 triple of 4-letter classical patterns whose involution avoiders are counted by the Motzkin numbers. (See Barnabei et al. 2011 reference.) - David Callan, Aug 27 2014
From Tony Foster III, Jul 28 2016: (Start)
A series created using 2*a(n) + a(n+1) has Hankel transform of F(2n), offset 3, F being the Fibonacci bisection, A001906 (empirical observation).
A series created using 2*a(n) + 3*a(n+1) + a(n+2) gives the Hankel transform of Sum_{k=0..n} k*Fibonacci(2*k), offset 3, A197649 (empirical observation). (End)
Conjecture: (2/n)*Sum_{k=1..n} (2k+1)*a(k)^2 is an integer for each positive integer n. - Zhi-Wei Sun, Nov 16 2017
The Rubey and Stump reference proves a refinement of a conjecture of René Marczinzik, which they state as: "The number of 2-Gorenstein algebras which are Nakayama algebras with n simple modules and have an oriented line as associated quiver equals the number of Motzkin paths of length n." - Eric M. Schmidt, Dec 16 2017
Number of U_{k}-equivalence classes of Łukasiewicz paths. Łukasiewicz paths are P-equivalent iff the positions of pattern P are identical in these paths. - Sergey Kirgizov, Apr 08 2018
If tau_1 and tau_2 are two distinct permutation patterns chosen from the set {132,231,312}, then a(n) is the number of valid hook configurations of permutations of [n+1] that avoid the patterns tau_1 and tau_2. - Colin Defant, Apr 28 2019
Number of permutations of length n that are sorted to the identity by a consecutive-321-avoiding stack followed by a classical-21-avoiding stack. - Colin Defant, Aug 29 2020
From Helmut Prodinger, Dec 13 2020: (Start)
a(n) is the number of paths in the first quadrant starting at (0,0) and consisting of n steps from the infinite set {(1,1), (1,-1), (1,-2), (1,-3), ...}.
For example, denoting U=(1,1), D=(1,-1), D_ j=(1,-j) for j >= 2, a(4) counts UUUU, UUUD, UUUD_2, UUUD_3, UUDU, UUDD, UUD_2U, UDUU, UDUD.
This step set is inspired by {(1,1), (1,-1), (1,-3), (1,-5), ...}, suggested by Emeric Deutsch around 2000.
See Prodinger link that contains a bijection to Motzkin paths. (End)
Named by Donaghey (1977) after the Israeli-American mathematician Theodore Motzkin (1908-1970). In Sloane's "A Handbook of Integer Sequences" (1973) they were called "generalized ballot numbers". - Amiram Eldar, Apr 15 2021
Number of Motzkin n-paths a(n) is split into A107587(n), number of even Motzkin n-paths, and A343386(n), number of odd Motzkin n-paths. The value A107587(n) - A343386(n) can be called the "shadow" of a(n) (see A343773). - Gennady Eremin, May 17 2021
Conjecture: If p is a prime of the form 6m+1 (A002476), then a(p-2) is divisible by p. Currently, no counterexample exists for p < 10^7. Personal communication from Robert Gerbicz: mod such p this is equivalent to A066796 with comment: "Every A066796(n) from A066796((p-1)/2) to A066796(p-1) is divisible by prime p of form 6m+1". - Serge Batalov, Feb 08 2022
From Rob Burns, Nov 11 2024: (Start)
The conjecture is proved in the 2017 paper by Rob Burns in the Links below. The result is contained in Tables 4 and 5 of the paper, which show that a(p-2) == 0 (mod p) when p == 1 (mod 6) and a(p-2) == -1 (mod p) when p == -1 (mod 6).
In fact, the 2017 paper by Burns establishes more general congruences for a(p^k - 2) where k >= 1.
If p == 1 (mod 6) then a(p^k - 2) == 0 (mod p) for k >= 1.
If p == -1 (mod 6) then a(p^k - 2) == -1 (mod p) when k is odd and a(p^k - 2) == 0 (mod p) when k is even.
These are consequences of the transitions provided in Tables 4, 5 and 6 of the paper.
The 2024 paper by Nadav Kohen also proves the conjecture. Proposition 6 of the paper states that a prime p divides a(p-2) if and only if p = (1 mod 3). (End)
From Peter Bala, Feb 10 2022: (Start)
Conjectures:
(1) For prime p == 1 (mod 6) and n, r >= 1, a(n*p^r - 2) == -A005717(n-1) (mod p), where we take A005717(0) = 0 to match Batalov's conjecture above.
(2) For prime p == 5 (mod 6) and n >= 1, a(n*p - 2) == -A005773(n) (mod p).
(3) For prime p >= 3 and k >= 1, a(n + p^k) == a(n) (mod p) for 0 <= n <= (p^k - 3).
(4) For prime p >= 5 and k >= 2, a(n + p^k) == a(n) (mod p^2) for 0 <= n <= (p^(k-1) - 3). (End)
The Hankel transform of this sequence with a(0) omitted gives the period-6 sequence [1, 0, -1, -1, 0, 1, ...] which is A010892 with its first term omitted, while the Hankel transform of the current sequence is the all-ones sequence A000012, and also it is the unique sequence with this property which is similar to the unique Hankel transform property of the Catalan numbers. - Michael Somos, Apr 17 2022
The number of terms in which the exponent of any variable x_i is not greater than 2 in the expansion of Product_{j=1..n} Sum_{i=1..j} x_i. E.g.: a(4) = 9: 3*x1^2*x2^2, 4*x1^2*x2*x3, 2*x1^2*x2*x4, x1^2*x3^2, x1^2*x3*x4, 2*x1*x2^2*x3, x1*x2^2*x4, x1*x2*x3^2, x1*x2*x3*x4. - Elif Baser, Dec 20 2024

Examples

			G.f.: 1 + x + 2*x^2 + 4*x^3 + 9*x^4 + 21*x^5 + 51*x^6 + 127*x^7 + 323*x^8 + ...
.
The 21 Motzkin-paths of length 5: UUDDF, UUDFD, UUFDD, UDUDF, UDUFD, UDFUD, UDFFF, UFUDD, UFDUD, UFDFF, UFFDF, UFFFD, FUUDD, FUDUD, FUDFF, FUFDF, FUFFD, FFUDF, FFUFD, FFFUD, FFFFF.
		

References

  • F. Bergeron, L. Favreau, and D. Krob, Conjectures on the enumeration of tableaux of bounded height, Discrete Math, vol. 139, no. 1-3 (1995), 463-468.
  • F. R. Bernhart, Catalan, Motzkin, and Riordan numbers, Discr. Math., 204 (1999) 73-112.
  • R. Bojicic and M. D. Petkovic, Orthogonal Polynomials Approach to the Hankel Transform of Sequences Based on Motzkin Numbers, Bulletin of the Malaysian Mathematical Sciences, 2015, doi:10.1007/s40840-015-0249-3.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pp. 24, 298, 618, 912.
  • A. J. Bu, Automated counting of restricted Motzkin paths, Enumerative Combinatorics and Applications, ECA 1:2 (2021) Article S2R12.
  • Naiomi Cameron, JE McLeod, Returns and Hills on Generalized Dyck Paths, Journal of Integer Sequences, Vol. 19, 2016, #16.6.1.
  • L. Carlitz, Solution of certain recurrences, SIAM J. Appl. Math., 17 (1969), 251-259.
  • Michael Dairyko, Samantha Tyner, Lara Pudwell, and Casey Wynn, Non-contiguous pattern avoidance in binary trees. Electron. J. Combin. 19 (2012), no. 3, Paper 22, 21 pp. MR2967227.
  • D. E. Davenport, L. W. Shapiro, and L. C. Woodson, The Double Riordan Group, The Electronic Journal of Combinatorics, 18(2) (2012), #P33.
  • E. Deutsch and L. Shapiro, A survey of the Fine numbers, Discrete Math., 241 (2001), 241-265.
  • T. Doslic, D. Svrtan, and D. Veljan, Enumerative aspects of secondary structures, Discr. Math., 285 (2004), 67-82.
  • Tomislav Doslic and Darko Veljan, Logarithmic behavior of some combinatorial sequences. Discrete Math. 308 (2008), no. 11, 2182-2212. MR2404544 (2009j:05019).
  • S. Dulucq and R. Simion, Combinatorial statistics on alternating permutations, J. Algebraic Combinatorics, 8, 1998, 169-191.
  • M. Dziemianczuk, "Enumerations of plane trees with multiple edges and Raney lattice paths." Discrete Mathematics 337 (2014): 9-24.
  • Wenjie Fang, A partial order on Motzkin paths, Discrete Math., 343 (2020), #111802.
  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (5.2.10).
  • N. S. S. Gu, N. Y. Li, and T. Mansour, 2-Binary trees: bijections and related issues, Discr. Math., 308 (2008), 1209-1221.
  • Kris Hatch, Presentation of the Motzkin Monoid, Senior Thesis, Univ. Cal. Santa Barbara, 2012; http://ccs.math.ucsb.edu/senior-thesis/Kris-Hatch.pdf.
  • V. Jelinek, Toufik Mansour, and M. Shattuck, On multiple pattern avoiding set partitions, Advances in Applied Mathematics Volume 50, Issue 2, February 2013, pp. 292-326.
  • Hana Kim and R. P. Stanley, A refined enumeration of hex trees and related polynomials, http://www-math.mit.edu/~rstan/papers/hextrees.pdf, Preprint 2015.
  • S. Kitaev, Patterns in Permutations and Words, Springer-Verlag, 2011. See p. 399 Table A.7.
  • A. Kuznetsov et al., Trees associated with the Motzkin numbers, J. Combin. Theory, A 76 (1996), 145-147.
  • T. Lengyel, On divisibility properties of some differences of Motzkin numbers, Annales Mathematicae et Informaticae, 41 (2013) pp. 121-136.
  • W. A. Lorenz, Y. Ponty, and P. Clote, Asymptotics of RNA Shapes, Journal of Computational Biology. 2008, 15(1): 31-63. doi:10.1089/cmb.2006.0153.
  • Piera Manara and Claudio Perelli Cippo, The fine structure of 4321 avoiding involutions and 321 avoiding involutions, PU. M. A. Vol. 22 (2011), 227-238; http://www.mat.unisi.it/newsito/puma/public_html/22_2/manara_perelli-cippo.pdf.
  • Toufik Mansour, Restricted 1-3-2 permutations and generalized patterns, Annals of Combin., 6 (2002), 65-76.
  • Toufik Mansour, Matthias Schork, and Mark Shattuck, Catalan numbers and pattern restricted set partitions. Discrete Math. 312(2012), no. 20, 2979-2991. MR2956089.
  • T. S. Motzkin, Relations between hypersurface cross ratios and a combinatorial formula for partitions of a polygon, for permanent preponderance and for non-associative products, Bull. Amer. Math. Soc., 54 (1948), 352-360.
  • Jocelyn Quaintance and Harris Kwong, A combinatorial interpretation of the Catalan and Bell number difference tables, Integers, 13 (2013), #A29.
  • J. Riordan, Enumeration of plane trees by branches and endpoints, J. Combin. Theory, A 23 (1975), 214-222.
  • A. Sapounakis et al., Ordered trees and the inorder transversal, Disc. Math., 306 (2006), 1732-1741.
  • A. Sapounakis, I. Tasoulas, and P. Tsikouras, Counting strings in Dyck paths, Discrete Math., 307 (2007), 2909-2924.
  • E. Schroeder, Vier combinatorische Probleme, Z. f. Math. Phys., 15 (1870), 361-376.
  • L. W. Shapiro et al., The Riordan group, Discrete Applied Math., 34 (1991), 229-239.
  • Mark Shattuck, On the zeros of some polynomials with combinatorial coefficients, Annales Mathematicae et Informaticae, 42 (2013) pp. 93-101, http://ami.ektf.hu.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1.
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 6.37. Also Problem 7.16(b), y_3(n).
  • P. R. Stein and M. S. Waterman, On some new sequences generalizing the Catalan and Motzkin numbers, Discrete Math., 26 (1979), 261-272.
  • Z.-W. Sun, Conjectures involving arithmetical sequences, Number Theory: Arithmetic in Shangri-La (eds., S. Kanemitsu, H.-Z. Li and J.-Y. Liu), Proc. the 6th China-Japan Sem. Number Theory (Shanghai, August 15-17, 2011), World Sci., Singapore, 2013, pp. 244-258; http://math.nju.edu.cn/~zwsun/142p.pdf.
  • Chenying Wang, Piotr Miska, and István Mező, "The r-derangement numbers." Discrete Mathematics 340.7 (2017): 1681-1692.
  • Ying Wang and Guoce Xin, A Classification of Motzkin Numbers Modulo 8, Electron. J. Combin., 25(1) (2018), #P1.54.
  • Wen-Jin Woan, A combinatorial proof of a recursive relation of the Motzkin sequence by lattice paths. Fibonacci Quart. 40 (2002), no. 1, 3-8.
  • Wen-jin Woan, A Recursive Relation for Weighted Motzkin Sequences, Journal of Integer Sequences, Vol. 8 (2005), Article 05.1.6.
  • F. Yano and H. Yoshida, Some set partition statistics in non-crossing partitions and generating functions, Discr. Math., 307 (2007), 3147-3160.

Crossrefs

Bisections: A026945, A099250.
Sequences related to chords in a circle: A001006, A054726, A006533, A006561, A006600, A007569, A007678. See also entries for chord diagrams in Index file.
a(n) = A005043(n)+A005043(n+1).
A086246 is another version, although this is the main entry. Column k=3 of A182172.
Motzkin numbers A001006 read mod 2,3,4,5,6,7,8,11: A039963, A039964, A299919, A258712, A299920, A258711, A299918, A258710.
Cf. A004148, A004149, A023421, A023422, A023423, A290277 (inv. Euler Transf.).

Programs

  • Haskell
    a001006 n = a001006_list !! n
    a001006_list = zipWith (+) a005043_list $ tail a005043_list
    -- Reinhard Zumkeller, Jan 31 2012
    
  • Maple
    # Three different Maple scripts for this sequence:
    A001006 := proc(n)
        add(binomial(n,2*k)*A000108(k),k=0..floor(n/2)) ;
    end proc:
    A001006 := proc(n) option remember; local k; if n <= 1 then 1 else procname(n-1) + add(procname(k)*procname(n-k-2),k=0..n-2); end if; end proc:
    # n -> [a(0),a(1),..,a(n)]
    A001006_list := proc(n) local w, m, j, i; w := proc(i,j,n) option remember;
    if min(i,j,n) < 0 or max(i,j) > n then 0
    elif n = 0 then if i = 0 and j = 0 then 1 else 0 fi else
    w(i, j + 1, n - 1) + w(i - 1, j, n - 1) + w(i + 1, j - 1, n - 1) fi end:
    [seq( add( add( w(i, j, m), i = 0..m), j = 0..m), m = 0..n)] end:
    A001006_list(29); # Peter Luschny, May 21 2011
  • Mathematica
    a[0] = 1; a[n_Integer] := a[n] = a[n - 1] + Sum[a[k] * a[n - 2 - k], {k, 0, n - 2}]; Array[a, 30]
    (* Second program: *)
    CoefficientList[Series[(1 - x - (1 - 2x - 3x^2)^(1/2))/(2x^2), {x, 0, 29}], x] (* Jean-François Alcover, Nov 29 2011 *)
    Table[Hypergeometric2F1[(1-n)/2, -n/2, 2, 4], {n,0,29}] (* Peter Luschny, May 15 2016 *)
    Table[GegenbauerC[n,-n-1,-1/2]/(n+1),{n,0,100}] (* Emanuele Munarini, Oct 20 2016 *)
    MotzkinNumber = DifferenceRoot[Function[{y, n}, {(-3n-3)*y[n] + (-2n-5)*y[n+1] + (n+4)*y[n+2] == 0, y[0] == 1, y[1] == 1}]];
    Table[MotzkinNumber[n], {n, 0, 29}] (* Jean-François Alcover, Oct 27 2021 *)
  • Maxima
    a[0]:1$
    a[1]:1$
    a[n]:=((2*n+1)*a[n-1]+(3*n-3)*a[n-2])/(n+2)$
    makelist(a[n],n,0,12); /* Emanuele Munarini, Mar 02 2011 */
    
  • Maxima
    M(n) := coeff(expand((1+x+x^2)^(n+1)),x^n)/(n+1);
    makelist(M(n),n,0,60); /* Emanuele Munarini, Apr 04 2012 */
    
  • Maxima
    makelist(ultraspherical(n,-n-1,-1/2)/(n+1),n,0,12); /* Emanuele Munarini, Oct 20 2016 */
    
  • PARI
    {a(n) = polcoeff( ( 1 - x - sqrt((1 - x)^2 - 4 * x^2 + x^3 * O(x^n))) / (2 * x^2), n)}; /* Michael Somos, Sep 25 2003 */
    
  • PARI
    {a(n) = if( n<0, 0, n++; polcoeff( serreverse( x / (1 + x + x^2) + x * O(x^n)), n))}; /* Michael Somos, Sep 25 2003 */
    
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( exp(x + x * O(x^n)) * besseli(1, 2 * x + x * O(x^n)), n))}; /* Michael Somos, Sep 25 2003 */
    
  • Python
    from gmpy2 import divexact
    A001006 = [1, 1]
    for n in range(2, 10**3):
        A001006.append(divexact(A001006[-1]*(2*n+1)+(3*n-3)*A001006[-2],n+2))
    # Chai Wah Wu, Sep 01 2014
    
  • Python
    def mot():
        a, b, n = 0, 1, 1
        while True:
            yield b//n
            n += 1
            a, b = b, (3*(n-1)*n*a+(2*n-1)*n*b)//((n+1)*(n-1))
    A001006 = mot()
    print([next(A001006) for n in range(30)]) # Peter Luschny, May 16 2016
    
  • Python
    # A simple generator of Motzkin-paths (see the first comment of David Callan).
    C = str.count
    def aGen(n: int):
        a = [""]
        for w in a:
            if len(w) == n:
                if C(w, "U") == C(w, "D"): yield w
            else:
                for j in "UDF":
                    u = w + j
                    if C(u, "U") >= C(u, "D"): a += [u]
        return a
    for n in range(6):
        MP = [w for w in aGen(n)];
        print(len(MP), ":", MP)  # Peter Luschny, Dec 03 2024

Formula

G.f.: A(x) = ( 1 - x - (1-2*x-3*x^2)^(1/2) ) / (2*x^2).
G.f. A(x) satisfies A(x) = 1 + x*A(x) + x^2*A(x)^2.
G.f.: F(x)/x where F(x) is the reversion of x/(1+x+x^2). - Joerg Arndt, Oct 23 2012
a(n) = (-1/2) Sum_{i+j = n+2, i >= 0, j >= 0} (-3)^i*C(1/2, i)*C(1/2, j).
a(n) = (3/2)^(n+2) * Sum_{k >= 1} 3^(-k) * Catalan(k-1) * binomial(k, n+2-k). [Doslic et al.]
a(n) ~ 3^(n+1)*sqrt(3)*(1 + 1/(16*n))/((2*n+3)*sqrt((n+2)*Pi)). [Barcucci, Pinzani and Sprugnoli]
Limit_{n->infinity} a(n)/a(n-1) = 3. [Aigner]
a(n+2) - a(n+1) = a(0)*a(n) + a(1)*a(n-1) + ... + a(n)*a(0). [Bernhart]
a(n) = (1/(n+1)) * Sum_{i} (n+1)!/(i!*(i+1)!*(n-2*i)!). [Bernhart]
From Len Smiley: (Start)
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(n, k)*A000108(k+1), inv. Binomial Transform of A000108.
a(n) = (1/(n+1))*Sum_{k=0..ceiling((n+1)/2)} binomial(n+1, k)*binomial(n+1-k, k-1);
D-finite with recurrence: (n+2)*a(n) = (2*n+1)*a(n-1) + (3*n-3)*a(n-2). (End)
a(n) = Sum_{k=0..n} C(n, 2k)*A000108(k). - Paul Barry, Jul 18 2003
E.g.f.: exp(x)*BesselI(1, 2*x)/x. - Vladeta Jovovic, Aug 20 2003
a(n) = A005043(n) + A005043(n+1).
The Hankel transform of this sequence gives A000012 = [1, 1, 1, 1, 1, 1, ...]. E.g., Det([1, 1, 2, 4; 1, 2, 4, 9; 2, 4, 9, 21; 4, 9, 21, 51]) = 1. - Philippe Deléham, Feb 23 2004
a(m+n) = Sum_{k>=0} A064189(m, k)*A064189(n, k). - Philippe Deléham, Mar 05 2004
a(n) = (1/(n+1))*Sum_{j=0..floor(n/3)} (-1)^j*binomial(n+1, j)*binomial(2*n-3*j, n). - Emeric Deutsch, Mar 13 2004
a(n) = A086615(n) - A086615(n-1) (n >= 1). - Emeric Deutsch, Jul 12 2004
G.f.: A(x)=(1-y+y^2)/(1-y)^2 where (1+x)*(y^2-y)+x=0; A(x)=4*(1+x)/(1+x+sqrt(1-2*x-3*x^2))^2; a(n)=(3/4)*(1/2)^n*Sum_(k=0..2*n, 3^(n-k)*C(k)*C(k+1, n+1-k) ) + 0^n/4 [after Doslic et al.]. - Paul Barry, Feb 22 2005
G.f.: c(x^2/(1-x)^2)/(1-x), c(x) the g.f. of A000108. - Paul Barry, May 31 2006
Asymptotic formula: a(n) ~ sqrt(3/4/Pi)*3^(n+1)/n^(3/2). - Benoit Cloitre, Jan 25 2007
a(n) = A007971(n+2)/2. - Zerinvary Lajos, Feb 28 2007
a(n) = (1/(2*Pi))*Integral_{x=-1..3} x^n*sqrt((3-x)*(1+x)) is the moment representation. - Paul Barry, Sep 10 2007
Given an integer t >= 1 and initial values u = [a_0, a_1, ..., a_{t-1}], we may define an infinite sequence Phi(u) by setting a_n = a_{n-1} + a_0*a_{n-1} + a_1*a_{n-2} + ... + a_{n-2}*a_1 for n >= t. For example, Phi([1]) is the Catalan numbers A000108. The present sequence is Phi([0,1,1]), see the 6th formula. - Gary W. Adamson, Oct 27 2008
G.f.: 1/(1-x-x^2/(1-x-x^2/(1-x-x^2/(1-x-x^2/(1-x-x^2/.... (continued fraction). - Paul Barry, Dec 06 2008
G.f.: 1/(1-(x+x^2)/(1-x^2/(1-(x+x^2)/(1-x^2/(1-(x+x^2)/(1-x^2/(1-.... (continued fraction). - Paul Barry, Feb 08 2009
a(n) = (-3)^(1/2)/(6*(n+2)) * (-1)^n*(3*hypergeom([1/2, n+1],[1],4/3) - hypergeom([1/2, n+2],[1],4/3)). - Mark van Hoeij, Nov 12 2009
G.f.: 1/(1-x/(1-x/(1-x^2/(1-x/(1-x/(1-x^2/(1-x/(1-x/(1-x^2/(1-... (continued fraction). - Paul Barry, Mar 02 2010
G.f.: 1/(1-x/(1-x/(1+x-x/(1-x/(1+x-x/(1-x/(1+x-x/(1-x/(1+x-x/(1-... (continued fraction). - Paul Barry, Jan 26 2011 [Adds apparently a third '1' in front. - R. J. Mathar, Jan 29 2011]
Let A(x) be the g.f., then B(x)=1+x*A(x) = 1 + 1*x + 1*x^2 + 2*x^3 + 4*x^4 + 9*x^5 + ... = 1/(1-z/(1-z/(1-z/(...)))) where z=x/(1+x) (continued fraction); more generally B(x)=C(x/(1+x)) where C(x) is the g.f. for the Catalan numbers (A000108). - Joerg Arndt, Mar 18 2011
a(n) = (2/Pi)*Integral_{x=-1..1} (1+2*x)^n*sqrt(1-x^2). - Peter Luschny, Sep 11 2011
G.f.: (1-x-sqrt(1-2*x-3*(x^2)))/(2*(x^2)) = 1/2/(x^2)-1/2/x-1/2/(x^2)*G(0); G(k) = 1+(4*k-1)*x*(2+3*x)/(4*k+2-x*(2+3*x)*(4*k+1)*(4*k+2) /(x*(2+3*x)*(4*k+1)+(4*k+4)/G(k+1))), if -1 < x < 1/3; (continued fraction). - Sergei N. Gladkovskii, Dec 01 2011
G.f.: (1-x-sqrt(1-2*x-3*(x^2)))/(2*(x^2)) = (-1 + 1/G(0))/(2*x); G(k) = 1-2*x/(1+x/(1+x/(1-2*x/(1-x/(2-x/G(k+1)))))); (continued fraction). - Sergei N. Gladkovskii, Dec 11 2011
0 = a(n) * (9*a(n+1) + 15*a(n+2) - 12*a(n+3)) + a(n+1) * ( -3*a(n+1) + 10*a(n+2) - 5*a(n+3)) + a(n+2) * (a(n+2) + a(n+3)) unless n=-2. - Michael Somos, Mar 23 2012
a(n) = (-1)^n*hypergeometric([-n,3/2],[3],4). - Peter Luschny, Aug 15 2012
Representation in terms of special values of Jacobi polynomials P(n,alpha,beta,x), in Maple notation: a(n)= 2*(-1)^n*n!*JacobiP(n,2,-3/2-n,-7)/(n+2)!, n>=0. - Karol A. Penson, Jun 24 2013
G.f.: Q(0)/x - 1/x, where Q(k) = 1 + (4*k+1)*x/((1+x)*(k+1) - x*(1+x)*(2*k+2)*(4*k+3)/(x*(8*k+6)+(2*k+3)*(1+x)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 14 2013
Catalan(n+1) = Sum_{k=0..n} binomial(n,k)*a(k). E.g.: 42 = 1*1 + 4*1 + 6*2 + 4*4 + 1*9. - Doron Zeilberger, Mar 12 2015
G.f. A(x) with offset 1 satisfies: A(x)^2 = A( x^2/(1-2*x) ). - Paul D. Hanna, Nov 08 2015
a(n) = GegenbauerPoly(n,-n-1,-1/2)/(n+1). - Emanuele Munarini, Oct 20 2016
a(n) = a(n-1) + A002026(n-1). Number of Motzkin paths that start with an F step plus number of Motzkin paths that start with an U step. - R. J. Mathar, Jul 25 2017
G.f. A(x) satisfies A(x)*A(-x) = F(x^2), where F(x) is the g.f. of A168592. - Alexander Burstein, Oct 04 2017
G.f.: A(x) = exp(int((E(x)-1)/x dx)), where E(x) is the g.f. of A002426. Equivalently, E(x) = 1 + x*A'(x)/A(x). - Alexander Burstein, Oct 05 2017
G.f. A(x) satisfies: A(x) = Sum_{j>=0} x^j * Sum_{k=0..j} binomial(j,k)*x^k*A(x)^k. - Ilya Gutkovskiy, Apr 11 2019
From Gennady Eremin, May 08 2021: (Start)
G.f.: 2/(1 - x + sqrt(1-2*x-3*x^2)).
a(n) = A107587(n) + A343386(n) = 2*A107587(n) - A343773(n) = 2*A343386(n) + A343773(n). (End)
Revert transform of A049347 (after Michael Somos). - Gennady Eremin, Jun 11 2021
Sum_{n>=0} 1/a(n) = 2.941237337631025604300320152921013604885956025483079699366681494505960039781389... - Vaclav Kotesovec, Jun 17 2021
Let a(-1) = (1 - sqrt(-3))/2 and a(n) = a(-3-n)*(-3)^(n+3/2) for all n in Z. Then a(n) satisfies my previous formula relation from Mar 23 2012 now for all n in Z. - Michael Somos, Apr 17 2022
Let b(n) = 1 for n <= 1, otherwise b(n) = Sum_{k=2..n} b(k-1) * b(n-k), then a(n) = b(n+1) (conjecture). - Joerg Arndt, Jan 16 2023
From Peter Bala, Feb 03 2024: (Start)
G.f.: A(x) = 1/(1 + x)*c(x/(1 + x))^2 = 1 + x/(1 + x)*c(x/(1 + x))^3, where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108.
A(x) = 1/(1 - 3*x)*c(-x/(1 -3*x))^2.
a(n+1) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n, k)*A000245(k+1).
a(n) = 3^n * Sum_{k = 0..n} (-3)^(-k)*binomial(n, k)*Catalan(k+1).
a(n) = 3^n * hypergeom([3/2, -n], [3], 4/3). (End)
G.f. A(x) satisfies A(x) = exp( x*A(x) + Integral x*A(x)/(1 - x^2*A(x)) dx ). - Paul D. Hanna, Mar 04 2024
a(n) = hypergeom([-n/2,1/2-n/2],[2],4). - Karol A. Penson, May 18 2025

A001764 a(n) = binomial(3*n,n)/(2*n+1) (enumerates ternary trees and also noncrossing trees).

Original entry on oeis.org

1, 1, 3, 12, 55, 273, 1428, 7752, 43263, 246675, 1430715, 8414640, 50067108, 300830572, 1822766520, 11124755664, 68328754959, 422030545335, 2619631042665, 16332922290300, 102240109897695, 642312451217745, 4048514844039120, 25594403741131680, 162250238001816900
Offset: 0

Views

Author

Keywords

Comments

Smallest number of straight line crossing-free spanning trees on n points in the plane.
Number of dissections of some convex polygon by nonintersecting diagonals into polygons with an odd number of sides and having a total number of 2n+1 edges (sides and diagonals). - Emeric Deutsch, Mar 06 2002
Number of lattice paths of n East steps and 2n North steps from (0,0) to (n,2n) and lying weakly below the line y=2x. - David Callan, Mar 14 2004
With interpolated zeros, this has g.f. 2*sqrt(3)*sin(arcsin(3*sqrt(3)*x/2)/3)/(3*x) and a(n) = C(n+floor(n/2),floor(n/2))*C(floor(n/2),n-floor(n/2))/(n+1). This is the first column of the inverse of the Riordan array (1-x^2,x(1-x^2)) (essentially reversion of y-y^3). - Paul Barry, Feb 02 2005
Number of 12312-avoiding matchings on [2n].
Number of complete ternary trees with n internal nodes, or 3n edges.
Number of rooted plane trees with 2n edges, where every vertex has even outdegree ("even trees").
a(n) is the number of noncrossing partitions of [2n] with all blocks of even size. E.g.: a(2)=3 counts 12-34, 14-23, 1234. - David Callan, Mar 30 2007
Pfaff-Fuss-Catalan sequence C^{m}_n for m=3, see the Graham et al. reference, p. 347. eq. 7.66.
Also 3-Raney sequence, see the Graham et al. reference, p. 346-7.
The number of lattice paths from (0,0) to (2n,0) using an Up-step=(1,1) and a Down-step=(0,-2) and staying above the x-axis. E.g., a(2) = 3; UUUUDD, UUUDUD, UUDUUD. - Charles Moore (chamoore(AT)howard.edu), Jan 09 2008
a(n) is (conjecturally) the number of permutations of [n+1] that avoid the patterns 4-2-3-1 and 4-2-5-1-3 and end with an ascent. For example, a(4)=55 counts all 60 permutations of [5] that end with an ascent except 42315, 52314, 52413, 53412, all of which contain a 4-2-3-1 pattern and 42513. - David Callan, Jul 22 2008
Central terms of pendular triangle A167763. - Philippe Deléham, Nov 12 2009
With B(x,t)=x+t*x^3, the comp. inverse in x about 0 is A(x,t) = Sum_{j>=0} a(j) (-t)^j x^(2j+1). Let U(x,t)=(x-A(x,t))/t. Then DU(x,t)/Dt=dU/dt+U*dU/dx=0 and U(x,0)=x^3, i.e., U is a solution of the inviscid Burgers's, or Hopf, equation. Also U(x,t)=U(x-t*U(x,t),0) and dB(x,t)/dt = U(B(x,t),t) = x^3 = U(x,0). The characteristics for the Hopf equation are x(t) = x(0) + t*U(x(t),t) = x(0) + t*U(x(0),0) = x(0) + t*x(0)^3 = B(x(0),t). These results apply to all the Fuss-Catalan sequences with 3 replaced by n>0 and 2 by n-1 (e.g., A000108 with n=2 and A002293 with n=4), see also A086810, which can be generalized to A133437, for associahedra. - Tom Copeland, Feb 15 2014
Number of intervals (i.e., ordered pairs (x,y) such that x<=y) in the Kreweras lattice (noncrossing partitions ordered by refinement) of size n, see the Bernardi & Bonichon (2009) and Kreweras (1972) references. - Noam Zeilberger, Jun 01 2016
Number of sum-indecomposable (4231,42513)-avoiding permutations. Conjecturally, number of sum-indecomposable (2431,45231)-avoiding permutations. - Alexander Burstein, Oct 19 2017
a(n) is the number of topologically distinct endstates for the game Planted Brussels Sprouts on n vertices, see Ji and Propp link. - Caleb Ji, May 14 2018
Number of complete quadrillages of 2n+2-gons. See Baryshnikov p. 12. See also Nov 10 2014 comments in A134264. - Tom Copeland, Jun 04 2018
a(n) is the number of 2-regular words on the alphabet [n] that avoid the patterns 231 and 221. Equivalently, this is the number of 2-regular tortoise-sortable words on the alphabet [n] (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
a(n) is the number of Motzkin paths of length 3n with n steps of each type, with the condition that (1, 0) and (1, 1) steps alternate (starting with (1, 0)). - Helmut Prodinger, Apr 08 2019
a(n) is the number of uniquely sorted permutations of length 2n+1 that avoid the patterns 312 and 1342. - Colin Defant, Jun 08 2019
The compositional inverse o.g.f. pair in Copeland's comment above are related to a pair of quantum fields in Balduf's thesis by Theorem 4.2 on p. 92. - Tom Copeland, Dec 13 2019
The sequences of Fuss-Catalan numbers, of which this is the first after the Catalan numbers A000108 (the next is A002293), appear in articles on random matrices and quantum physics. See Banica et al., Collins et al., and Mlotkowski et al. Interpretations of these sequences in terms of the cardinality of specific sets of noncrossing partitions are provided by A134264. - Tom Copeland, Dec 21 2019
Call C(p, [alpha], g) the number of partitions of a cyclically ordered set with p elements, of cyclic type [alpha], and of genus g (the genus g Faa di Bruno coefficients of type [alpha]). This sequence counts the genus 0 partitions (non-crossing, or planar, partitions) of p = 3n into n parts of length 3: a(n) = C(3n, [3^n], 0). For genus 1 see A371250, for genus 2 see A371251. - Robert Coquereaux, Mar 16 2024
a(n) is the total number of down steps before the first up step in all 2_1-Dyck paths of length 3*n for n > 0. A 2_1-Dyck path is a lattice path with steps (1,2), (1,-1) that starts and ends at y = 0 and does not go below the line y = -1. - Sarah Selkirk, May 10 2020
a(n) is the number of pairs (A<=B) of noncrossing partitions of [n]. - Francesca Aicardi, May 28 2022
a(n) is the number of parking functions of size n avoiding the patterns 231 and 321. - Lara Pudwell, Apr 10 2023
Number of rooted polyominoes composed of n square cells of the hyperbolic regular tiling with Schläfli symbol {4,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {4,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 3 of the family {C(k, n)}A130564.%20-%20_Wolfdieter%20Lang">{n>=0} given in a comment in A130564. - _Wolfdieter Lang, Feb 05 2024
The number of Apollonian networks (planar 3-trees) with n+3 vertices with a given base triangle. - Allan Bickle, Feb 20 2024
Number of rooted polyominoes composed of n tetrahedral cells of the hyperbolic regular tiling with Schläfli symbol {3,3,oo}. A rooted polyomino has one external face identified, and chiral pairs are counted as two. a(n) = T(n) in the second Beineke and Pippert link. - Robert A. Russell, Mar 20 2024

Examples

			a(2) = 3 because the only dissections with 5 edges are given by a square dissected by any of the two diagonals and the pentagon with no dissecting diagonal.
G.f. = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + 1428*x^6 + 7752*x^7 + 43263*x^8 + ...
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 23.
  • I. M. H. Etherington, On non-associative combinations, Proc. Royal Soc. Edinburgh, 59 (Part 2, 1938-39), 153-162.
  • I. M. H. Etherington, Some problems of non-associative combinations (I), Edinburgh Math. Notes, 32 (1940), pp. i-vi. Part II is by A. Erdelyi and I. M. H. Etherington, and is on pages vii-xiv of the same issue.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347. See also the Pólya-Szegő reference.
  • W. Kuich, Languages and the enumeration of planted plane trees. Nederl. Akad. Wetensch. Proc. Ser. A 73 = Indag. Math. 32, (1970), 268-280.
  • T. V. Narayana, Lattice Path Combinatorics with Statistical Applications. Univ. Toronto Press, 1979, p. 98.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, Springer-Verlag, New York, Heidelberg, Berlin, 2 vols., 1972, Vol. 1, problem 211, p. 146 with solution on p. 348.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A001762, A001763, A002294 - A002296, A006013, A025174, A063548, A064017, A072247, A072248, A134264, A143603, A258708, A256311, A188687 (binomial transform), A346628 (inverse binomial transform).
A column of triangle A102537.
Bisection of A047749 and A047761.
Row sums of triangles A108410 and A108767.
Second column of triangle A062993.
Mod 3 = A113047.
2D Polyominoes: A005034 (oriented), A005036 (unoriented), A369315 (chiral), A047749 (achiral), A000108 {3,oo}, A002293 {5,oo}.
3D Polyominoes: A007173 (oriented), A027610 (unoriented), A371350 (chiral), A371351 (achiral).
Cf. A130564 (for C(k, n) cases).

Programs

  • GAP
    List([0..25],n->Binomial(3*n,n)/(2*n+1)); # Muniru A Asiru, Oct 31 2018
    
  • Haskell
    a001764 n = a001764_list !! n
    a001764_list = 1 : [a258708 (2 * n) n | n <- [1..]]
    -- Reinhard Zumkeller, Jun 23 2015
    
  • Magma
    [Binomial(3*n,n)/(2*n+1): n in [0..30]]; // Vincenzo Librandi, Sep 04 2014
    
  • Maple
    A001764 := n->binomial(3*n,n)/(2*n+1): seq(A001764(n), n=0..25);
    with(combstruct): BB:=[T,{T=Prod(Z,F),F=Sequence(B),B=Prod(F,Z,F)}, unlabeled]:seq(count(BB,size=i),i=0..22); # Zerinvary Lajos, Apr 22 2007
    with(combstruct):BB:=[S, {B = Prod(S,S,Z), S = Sequence(B)}, labelled]: seq(count(BB, size=n)/n!, n=0..21); # Zerinvary Lajos, Apr 25 2008
    n:=30:G:=series(RootOf(g = 1+x*g^3, g),x=0,n+1):seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 03 2015
    alias(PS=ListTools:-PartialSums): A001764List := proc(m) local A, P, n;
    A := [1,1]; P := [1]; for n from 1 to m - 2 do P := PS(PS([op(P), P[-1]]));
    A := [op(A), P[-1]] od; A end: A001764List(25); # Peter Luschny, Mar 26 2022
  • Mathematica
    InverseSeries[Series[y-y^3, {y, 0, 24}], x] (* then a(n)=y(2n+1)=ways to place non-crossing diagonals in convex (2n+4)-gon so as to create only quadrilateral tiles *) (* Len Smiley, Apr 08 2000 *)
    Table[Binomial[3n,n]/(2n+1),{n,0,25}] (* Harvey P. Dale, Jul 24 2011 *)
  • PARI
    {a(n) = if( n<0, 0, (3*n)! / n! / (2*n + 1)!)};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( serreverse( x - x^3 + O(x^(2*n + 2))), 2*n + 1))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = 1 + O(x); for( m=1, n, A = 1 + x * A^3); polcoeff(A, n))};
    
  • PARI
    b=vector(22);b[1]=1;for(n=2,22,for(i=1,n-1,for(j=1,n-1,for(k=1,n-1,if((i-1)+(j-1)+(k-1)-(n-2),NULL,b[n]=b[n]+b[i]*b[j]*b[k])))));a(n)=b[n+1]; print1(a(0));for(n=1,21,print1(", ",a(n))) \\ Gerald McGarvey, Oct 08 2008
    
  • PARI
    Vec(1 + serreverse(x / (1+x)^3 + O(x^30))) \\ Gheorghe Coserea, Aug 05 2015
    
  • Python
    from math import comb
    def A001764(n): return comb(3*n,n)//(2*n+1) # Chai Wah Wu, Nov 10 2022
  • Sage
    def A001764_list(n) :
        D = [0]*(n+1); D[1] = 1
        R = []; b = false; h = 1
        for i in range(2*n) :
            for k in (1..h) : D[k] += D[k-1]
            if not b : R.append(D[h])
            else : h += 1
            b = not b
        return R
    A001764_list(22) # Peter Luschny, May 03 2012
    

Formula

From Karol A. Penson, Nov 08 2001: (Start)
G.f.: (2/sqrt(3*x))*sin((1/3)*arcsin(sqrt(27*x/4))).
E.g.f.: hypergeom([1/3, 2/3], [1, 3/2], 27/4*x).
Integral representation as n-th moment of a positive function on [0, 27/4]: a(n) = Integral_{x=0..27/4} (x^n*((1/12) * 3^(1/2) * 2^(1/3) * (2^(1/3)*(27 + 3 * sqrt(81 - 12*x))^(2/3) - 6 * x^(1/3))/(Pi * x^(2/3)*(27 + 3 * sqrt(81 - 12*x))^(1/3)))), n >= 0. This representation is unique. (End)
G.f. A(x) satisfies A(x) = 1+x*A(x)^3 = 1/(1-x*A(x)^2) [Cyvin (1998)]. - Ralf Stephan, Jun 30 2003
a(n) = n-th coefficient in expansion of power series P(n), where P(0) = 1, P(k+1) = 1/(1 - x*P(k)^2).
G.f. Rev(x/c(x))/x, where c(x) is the g.f. of A000108 (Rev=reversion of). - Paul Barry, Mar 26 2010
From Gary W. Adamson, Jul 07 2011: (Start)
Let M = the production matrix:
1, 1
2, 2, 1
3, 3, 2, 1
4, 4, 3, 2, 1
5, 5, 4, 3, 2, 1
...
a(n) = upper left term in M^n. Top row terms of M^n = (n+1)-th row of triangle A143603, with top row sums generating A006013: (1, 2, 7, 30, 143, 728, ...). (End)
Recurrence: a(0)=1; a(n) = Sum_{i=0..n-1, j=0..n-1-i} a(i)a(j)a(n-1-i-j) for n >= 1 (counts ternary trees by subtrees of the root). - David Callan, Nov 21 2011
G.f.: 1 + 6*x/(Q(0) - 6*x); Q(k) = 3*x*(3*k + 1)*(3*k + 2) + 2*(2*(k^2) + 5*k +3) - 6*x*(2*(k^2) + 5*k + 3)*(3*k + 4)*(3*k + 5)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 27 2011
D-finite with recurrence: 2*n*(2n+1)*a(n) - 3*(3n-1)*(3n-2)*a(n-1) = 0. - R. J. Mathar, Dec 14 2011
REVERT transform of A115140. BINOMIAL transform is A188687. SUMADJ transform of A188678. HANKEL transform is A051255. INVERT transform of A023053. INVERT transform is A098746. - Michael Somos, Apr 07 2012
(n + 1) * a(n) = A174687(n).
G.f.: F([2/3,4/3], [3/2], 27/4*x) / F([2/3,1/3], [1/2], (27/4)*x) where F() is the hypergeometric function. - Joerg Arndt, Sep 01 2012
a(n) = binomial(3*n+1, n)/(3*n+1) = A062993(n+1,1). - Robert FERREOL, Apr 03 2015
a(n) = A258708(2*n,n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
0 = a(n)*(-3188646*a(n+2) + 20312856*a(n+3) - 11379609*a(n+4) + 1437501*a(n+5)) + a(n+1)*(177147*a(n+2) - 2247831*a(n+3) + 1638648*a(n+4) - 238604*a(n+5)) + a(n+2)*(243*a(n+2) + 31497*a(n+3) - 43732*a(n+4) + 8288*a(n+5)) for all integer n. - Michael Somos, Jun 03 2016
a(n) ~ 3^(3*n + 1/2)/(sqrt(Pi)*4^(n+1)*n^(3/2)). - Ilya Gutkovskiy, Nov 21 2016
Given g.f. A(x), then A(1/8) = -1 + sqrt(5), A(2/27) = (-1 + sqrt(3))*3/2, A(4/27) = 3/2, A(3/64) = -2 + 2*sqrt(7/3), A(5/64) = (-1 + sqrt(5))*2/sqrt(5), etc. A(n^2/(n+1)^3) = (n+1)/n if n > 1. - Michael Somos, Jul 17 2018
From Peter Bala, Sep 14 2021: (Start)
A(x) = exp( Sum_{n >= 1} (1/3)*binomial(3*n,n)*x^n/n ).
The sequence defined by b(n) := [x^n] A(x)^n = A224274(n) for n >= 1 and satisfies the congruence b(p) == b(1) (mod p^3) for prime p >= 3. Cf. A060941. (End)
G.f.: 1/sqrt(B(x)+(1-6*x)/(9*B(x))+1/3), with B(x):=((27*x^2-18*x+2)/54-(x*sqrt((-(4-27*x))*x))/(2*3^(3/2)))^(1/3). - Vladimir Kruchinin, Sep 28 2021
x*A'(x)/A(x) = (A(x) - 1)/(- 2*A(x) + 3) = x + 5*x^2 + 28*x^3 + 165*x^4 + ... is the o.g.f. of A025174. Cf. A002293 - A002296. - Peter Bala, Feb 04 2022
a(n) = hypergeom([1 - n, -2*n], [2], 1). Row sums of A108767. - Peter Bala, Aug 30 2023
G.f.: z*exp(3*z*hypergeom([1, 1, 4/3, 5/3], [3/2, 2, 2], (27*z)/4)) + 1.
- Karol A. Penson, Dec 19 2023
G.f.: hypergeometric([1/3, 2/3], [3/2], (3^3/2^2)*x). See the e.g.f. above. - Wolfdieter Lang, Feb 04 2024
a(n) = (3*n)! / (n!*(2*n+1)!). - Allan Bickle, Feb 20 2024
Sum_{n >= 0} a(n)*x^n/(1 + x)^(3*n+1) = 1. See A316371 and A346627. - Peter Bala, Jun 02 2024
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^5). - Seiichi Manyama, Jun 16 2025

A006318 Large Schröder numbers (or large Schroeder numbers, or big Schroeder numbers).

Original entry on oeis.org

1, 2, 6, 22, 90, 394, 1806, 8558, 41586, 206098, 1037718, 5293446, 27297738, 142078746, 745387038, 3937603038, 20927156706, 111818026018, 600318853926, 3236724317174, 17518619320890, 95149655201962, 518431875418926, 2832923350929742, 15521467648875090
Offset: 0

Views

Author

Keywords

Comments

For the little Schröder numbers (or little Schroeder numbers, or small Schroeder numbers) see A001003.
The number of perfect matchings in a triangular grid of n squares (n = 1, 4, 9, 16, 25, ...). - Roberto E. Martinez II, Nov 05 2001
a(n) is the number of subdiagonal paths from (0, 0) to (n, n) consisting of steps East (1, 0), North (0, 1) and Northeast (1, 1) (sometimes called royal paths). - David Callan, Mar 14 2004
Twice A001003 (except for the first term).
a(n) is the number of dissections of a regular (n+4)-gon by diagonals that do not touch the base. (A diagonal is a straight line joining two nonconsecutive vertices and dissection means the diagonals are noncrossing though they may share an endpoint. One side of the (n+4)-gon is designated the base.) Example: a(1)=2 because a pentagon has only 2 such dissections: the empty one and the one with a diagonal parallel to the base. - David Callan, Aug 02 2004
a(n) is the number of separable permutations, i.e., permutations avoiding 2413 and 3142 (see Shapiro and Stephens). - Vincent Vatter, Aug 16 2006
Eric W. Weisstein comments that the Schröder numbers bear the same relationship to the Delannoy numbers (A001850) as the Catalan numbers (A000108) do to the binomial coefficients. - Jonathan Vos Post, Dec 23 2004
a(n) is the number of lattice paths from (0, 0) to (n+1, n+1) consisting of unit steps north N = (0, 1) and variable-length steps east E = (k, 0), with k a positive integer, that stay strictly below the line y = x except at the endpoints. For example, a(2) = 6 counts 111NNN, 21NNN, 3NNN, 12NNN, 11N1NN, 2N1NN (east steps indicated by their length). If the word "strictly" is replaced by "weakly", the counting sequence becomes the little Schröder numbers, A001003 (offset). - David Callan, Jun 07 2006
a(n) is the number of dissections of a regular (n+3)-gon with base AB that do not contain a triangle of the form ABP with BP a diagonal. Example: a(1) = 2 because the square D-C | | A-B has only 2 such dissections: the empty one and the one with the single diagonal AC (although this dissection contains the triangle ABC, BC is not a diagonal). - David Callan, Jul 14 2006
a(n) is the number of (colored) Motzkin n-paths with each upstep and each flatstep at ground level getting one of 2 colors and each flatstep not at ground level getting one of 3 colors. Example: With their colors immediately following upsteps/flatsteps, a(2) = 6 counts U1D, U2D, F1F1, F1F2, F2F1, F2F2. - David Callan, Aug 16 2006
The Hankel transform of this sequence is A006125(n+1) = [1, 2, 8, 64, 1024, 32768, ...]; example: Det([1, 2, 6, 22; 2, 6, 22, 90; 6, 22, 90, 394; 22, 90, 394, 1806]) = 64. - Philippe Deléham, Sep 03 2006
Triangle A144156 has row sums equal to A006318 with left border A001003. - Gary W. Adamson, Sep 12 2008
a(n) is also the number of order-preserving and order-decreasing partial transformations (of an n-chain). Equivalently, it is the order of the Schröder monoid, PC sub n. - Abdullahi Umar, Oct 02 2008
Sum_{n >= 0} a(n)/10^n - 1 = (9 - sqrt(41))/2. - Mark Dols, Jun 22 2010
1/sqrt(41) = Sum_{n >= 0} Delannoy number(n)/10^n. - Mark Dols, Jun 22 2010
a(n) is also the dimension of the space Hoch(n) related to Hochschild two-cocycles. - Ph. Leroux (ph_ler_math(AT)yahoo.com), Aug 24 2010
Let W = (w(n, k)) denote the augmentation triangle (as at A193091) of A154325; then w(n, n) = A006318(n). - Clark Kimberling, Jul 30 2011
Conjecture: For each n > 2, the polynomial sum_{k = 0}^n a(k)*x^{n-k} is irreducible modulo some prime p < n*(n+1). - Zhi-Wei Sun, Apr 07 2013
From Jon Perry, May 24 2013: (Start)
Consider a Pascal triangle variant where T(n, k) = T(n, k-1) + T(n-1, k-1) + T(n-1, k), i.e., the order of performing the calculation must go from left to right (A033877). This sequence is the rightmost diagonal.
Triangle begins:
1;
1, 2;
1, 4, 6;
1, 6, 16, 22;
1, 8, 30, 68, 90;
... (End)
a(n) is the number of permutations avoiding 2143, 3142 and one of the patterns among 246135, 254613, 263514, 524361, 546132. - Alexander Burstein, Oct 05 2014
a(n) is the number of semi-standard Young tableaux of shape n x 2 with consecutive entries. That is, j in P and 1 <= i<= j imply i in P. - Graham H. Hawkes, Feb 15 2015
a(n) is the number of unary-rooted size n unary-binary trees (each node has either 1 or 2 degree out). - John Bodeen, May 29 2017
Conjecturally, a(n) is the number of permutations pi of length n such that s(pi) avoids the patterns 231 and 321, where s denotes West's stack-sorting map. - Colin Defant, Sep 17 2018
a(n) is the number of n X n permutation matrices which percolate under the 2-neighbor bootstrap percolation rule (see Shapiro and Stephens). The number of general n X n matrices of weight n which percolate is given in A146971. - Jonathan Noel, Oct 05 2018
a(n) is the number of permutations of length n+1 which avoid 3142 and 3241. The permutations are precisely the permutations that are sortable by a decreasing stack followed by an increasing stack in series. - Rebecca Smith, Jun 06 2019
a(n) is the number of permutations of length n+1 avoiding the partially ordered pattern (POP) {3>1, 4>1, 1>2} of length 4. That is, the number of length n+1 permutations having no subsequences of length 4 in which the second element is the smallest, and the first element is smaller than the third and fourth elements. - Sergey Kitaev, Dec 10 2020
Named after the German mathematician Ernst Schröder (1841-1902). - Amiram Eldar, Apr 15 2021
a(n) is the number of sequences of nonnegative integers (u_1, u_2, ..., u_n) such that (i) u_i <= i for all i, and (ii) the nonzero u_i are weakly increasing. For example, a(2) = 6 counts 00, 01, 02, 10, 11, 12. See link "Some bijections for lattice paths" at A001003. - David Callan, Dec 18 2021
a(n) is the number of separable elements of the Weyl group of type B_n/C_n (see Gaetz and Gao). - Fern Gossow, Jul 31 2023
The number of domino tilings of an Aztec triangle of order n. Dually, the number perfect matchings of the edges in the cellular graph formed by a triangular grid of n squares (n = 1, 4, 9, 16, 25, ...) as in Ciucu (1996). - Michael Somos, Sep 16 2024
a(n) is the number of dissections of a convex (n+3)-sided polygon by non-intersecting diagonals such that none of the dividing diagonals passes through a chosen vertex. - Muhammed Sefa Saydam, Mar 01 2025
a(n) is the number of dissections of a convex (n+m+1)-sided polygon by non-intersecting diagonals such that the selected m consecutive sides of the polygon will be in the same subpolygon. - Muhammed Sefa Saydam, Jul 02 2025

Examples

			a(3) = 22 since the top row of Q^n = (6, 6, 6, 4, 0, 0, 0, ...); where 22 = (6 + 6 + 6 + 4).
G.f. = 1 + 2*x + 6*x^2 + 22*x^3 + 90*x^4 + 394*x^5 + 1806*x^6 + 8858*x^7 + 41586*x^8 + ...
		

References

  • D. Andrica and E. J. Ionascu, On the number of polynomials with coefficients in [n], An. St. Univ. Ovidius Constanta, 2013, to appear.
  • Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.4.
  • Paul Barry, Riordan-Bernstein Polynomials, Hankel Transforms and Somos Sequences, Journal of Integer Sequences, Vol. 15 2012, #12.8.2.
  • Paul Barry, Riordan arrays, generalized Narayana triangles, and series reversion, Linear Algebra and its Applications, 491 (2016) 343-385.
  • Paul Barry, Jacobsthal Decompositions of Pascal's Triangle, Ternary Trees, and Alternating Sign Matrices, Journal of Integer Sequences, 19, 2016, #16.3.5.
  • O. Bodini, A. Genitrini, F. Peschanski, and N.Rolin, Associativity for binary parallel processes, CALDAM 2015.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pages 24, 618.
  • S. Brlek, E. Duchi, E. Pergola, and S. Rinaldi, On the equivalence problem for succession rules, Discr. Math., 298 (2005), 142-154.
  • Xiang-Ke Chang, XB Hu, H Lei, and YN Yeh, Combinatorial proofs of addition formulas, The Electronic Journal of Combinatorics, 23(1) (2016), #P1.8.
  • William Y. C. Chen and Carol J. Wang, Noncrossing Linked Partitions and Large (3, 2)-Motzkin Paths, Discrete Math., 312 (2012), 1918-1922.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81, #21, (4), q_n.
  • D. E. Davenport, L. W. Shapiro, and L. C. Woodson, The Double Riordan Group, The Electronic Journal of Combinatorics, 18(2) (2012), #P33.
  • Deng, Eva Y. P.; Dukes, Mark; Mansour, Toufik; and Wu, Susan Y. J.; Symmetric Schröder paths and restricted involutions. Discrete Math. 309 (2009), no. 12, 4108-4115. See p. 4109.
  • E. Deutsch, A bijective proof of an equation linking the Schroeder numbers, large and small, Discrete Math., 241 (2001), 235-240.
  • C. Domb and A. J. Barrett, Enumeration of ladder graphs, Discrete Math. 9 (1974), 341-358.
  • Doslic, Tomislav and Veljan, Darko. Logarithmic behavior of some combinatorial sequences. Discrete Math. 308 (2008), no. 11, 2182--2212. MR2404544 (2009j:05019) - From N. J. A. Sloane, May 01 2012
  • M. Dziemianczuk, Generalizing Delannoy numbers via counting weighted lattice paths, INTEGERS, 13 (2013), #A54.
  • Egge, Eric S., Restricted signed permutations counted by the Schröder numbers. Discrete Math. 306 (2006), 552-563. [Many applications of these numbers.]
  • S. Getu et al., How to guess a generating function, SIAM J. Discrete Math., 5 (1992), 497-499.
  • S. Gire, Arbres, permutations a motifs exclus et cartes planaire: quelques problemes algorithmiques et combinatoires, Ph.D. Thesis, Universite Bordeaux I, 1993.
  • N. S. S. Gu, N. Y. Li, and T. Mansour, 2-Binary trees: bijections and related issues, Discr. Math., 308 (2008), 1209-1221.
  • Guruswami, Venkatesan, Enumerative aspects of certain subclasses of perfect graphs. Discrete Math. 205 (1999), 97-117.
  • Silvia Heubach and Toufik Mansour, Combinatorics of Compositions and Words, CRC Press, 2010.
  • D. E. Knuth, The Art of Computer Programming, Vol. 1, Section 2.2.1, Problem 11.
  • D. Kremer, Permutations with forbidden subsequences and a generalized Schröder number, Discrete Math. 218 (2000) 121-130.
  • Kremer, Darla and Shiu, Wai Chee; Finite transition matrices for permutations avoiding pairs of length four patterns. Discrete Math. 268 (2003), 171-183. MR1983276 (2004b:05006). See Table 1.
  • Laradji, A. and Umar, A. Asymptotic results for semigroups of order-preserving partial transformations. Comm. Algebra 34 (2006), 1071-1075. - Abdullahi Umar, Oct 11 2008
  • L. Moser and W. Zayachkowski, Lattice paths with diagonal steps, Scripta Math., 26 (1961), 223-229.
  • L. Shapiro and A. B. Stephens, Bootstrap percolation, the Schröder numbers and the N-kings problem, SIAM J. Discrete Math., Vol. 4 (1991), pp. 275-280.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see page 178 and also Problems 6.39 and 6.40.
  • Lin Yang and S.-L. Yang, The parametric Pascal rhombus. Fib. Q., 57:4 (2019), 337-346.
  • Sheng-Liang Yang and Mei-yang Jiang, The m-Schröder paths and m-Schröder numbers, Disc. Math. (2021) Vol. 344, Issue 2, 112209. doi:10.1016/j.disc.2020.112209. See Table 1.

Crossrefs

Apart from leading term, twice A001003 (the small Schroeder numbers). Cf. A025240.
Sequences A085403, A086456, A103137, A112478 are essentially the same sequence.
Main diagonal of A033877.
Row sums of A104219. Bisections give A138462, A138463.
Row sums of A175124.
The sequences listed in Yang-Jiang's Table 1 appear to be A006318, A001003, A027307, A034015, A144097, A243675, A260332, A243676. - N. J. A. Sloane, Mar 28 2021

Programs

  • GAP
    Concatenation([1],List([1..25],n->(1/n)*Sum([0..n],k->2^k*Binomial(n,k)*Binomial(n,k-1)))); # Muniru A Asiru, Nov 29 2018
  • Haskell
    a006318 n = a004148_list !! n
    a006318_list = 1 : f [1] where
       f xs = y : f (y : xs) where
         y = head xs + sum (zipWith (*) xs $ reverse xs)
    -- Reinhard Zumkeller, Nov 13 2012
    
  • Maple
    Order := 24: solve(series((y-y^2)/(1+y),y)=x,y); # then A(x)=y(x)/x
    BB:=(-1-z-sqrt(1-6*z+z^2))/2: BBser:=series(BB, z=0, 24): seq(coeff(BBser, z, n), n=1..23); # Zerinvary Lajos, Apr 10 2007
    A006318_list := proc(n) local j, a, w; a := array(0..n); a[0] := 1;
    for w from 1 to n do a[w] := 2*a[w-1]+add(a[j]*a[w-j-1], j=1..w-1) od; convert(a,list)end: A006318_list(22); # Peter Luschny, May 19 2011
    A006318 := n-> add(binomial(n+k, n-k) * binomial(2*k, k)/(k+1), k=0..n): seq(A006318(n), n=0..22); # Johannes W. Meijer, Jul 14 2013
    seq(simplify(hypergeom([-n,n+1],[2],-1)), n=0..100); # Robert Israel, Mar 23 2015
  • Mathematica
    a[0] = 1; a[n_Integer] := a[n] = a[n - 1] + Sum[a[k]*a[n - 1 - k], {k, 0, n - 1}]; Array[a[#] &, 30]
    InverseSeries[Series[(y - y^2)/(1 + y), {y, 0, 24}], x] (* then A(x) = y(x)/x *) (* Len Smiley, Apr 11 2000 *)
    CoefficientList[Series[(1 - x - (1 - 6x + x^2)^(1/2))/(2x), {x, 0, 30}], x] (* Harvey P. Dale, May 01 2011 *)
    a[ n_] := 2 Hypergeometric2F1[ -n + 1, n + 2, 2, -1]; (* Michael Somos, Apr 03 2013 *)
    a[ n_] := With[{m = If[ n < 0, -1 - n, n]}, SeriesCoefficient[(1 - x - Sqrt[ 1 - 6 x + x^2])/(2 x), {x, 0, m}]]; (* Michael Somos, Jun 10 2015 *)
    Table[-(GegenbauerC[n+1, -1/2, 3] + KroneckerDelta[n])/2, {n, 0, 30}] (* Vladimir Reshetnikov, Nov 12 2016 *)
    CoefficientList[Nest[1+x(#+#^2)&, 1+O[x], 20], x] (* Oliver Seipel, Dec 21 2024 *)
  • PARI
    {a(n) = if( n<0, n = -1-n); polcoeff( (1 - x - sqrt( 1 - 6*x + x^2 + x^2 * O(x^n))) / 2, n+1)}; /* Michael Somos, Apr 03 2013 */
    
  • PARI
    {a(n) = if( n<1, 1, sum( k=0, n, 2^k * binomial( n, k) * binomial( n, k-1)) / n)};
    
  • Python
    from gmpy2 import divexact
    A006318 = [1, 2]
    for n in range(3,10**3):
        A006318.append(int(divexact(A006318[-1]*(6*n-9)-(n-3)*A006318[-2],n)))
    # Chai Wah Wu, Sep 01 2014
    
  • Sage
    # Generalized algorithm of L. Seidel
    def A006318_list(n) :
        D = [0]*(n+1); D[1] = 1
        b = True; h = 1; R = []
        for i in range(2*n) :
            if b :
                for k in range(h,0,-1) : D[k] += D[k-1]
                h += 1;
            else :
                for k in range(1,h, 1) : D[k] += D[k-1]
                R.append(D[h-1]);
            b = not b
        return R
    A006318_list(23) # Peter Luschny, Jun 02 2012
    

Formula

G.f.: (1 - x - (1 - 6*x + x^2)^(1/2))/(2*x).
a(n) = 2*hypergeom([-n+1, n+2], [2], -1). - Vladeta Jovovic, Apr 24 2003
For n > 0, a(n) = (1/n)*Sum_{k = 0..n} 2^k*C(n, k)*C(n, k-1). - Benoit Cloitre, May 10 2003
The g.f. satisfies (1 - x)*A(x) - x*A(x)^2 = 1. - Ralf Stephan, Jun 30 2003
For the asymptotic behavior, see A001003 (remembering that A006318 = 2*A001003). - N. J. A. Sloane, Apr 10 2011
From Philippe Deléham, Nov 28 2003: (Start)
Row sums of A088617 and A060693.
a(n) = Sum_{k = 0..n} C(n+k, n)*C(n, k)/(k+1). (End)
With offset 1: a(1) = 1, a(n) = a(n-1) + Sum_{i = 1..n-1} a(i)*a(n-i). - Benoit Cloitre, Mar 16 2004
a(n) = Sum_{k = 0..n} A000108(k)*binomial(n+k, n-k). - Benoit Cloitre, May 09 2004
a(n) = Sum_{k = 0..n} A011117(n, k). - Philippe Deléham, Jul 10 2004
a(n) = (CentralDelannoy(n+1) - 3 * CentralDelannoy(n))/(2*n) = (-CentralDelannoy(n+1) + 6 * CentralDelannoy(n) - CentralDelannoy(n-1))/2 for n >= 1, where CentralDelannoy is A001850. - David Callan, Aug 16 2006
From Abdullahi Umar, Oct 11 2008: (Start)
A123164(n+1) - A123164(n) = (2*n+1)*a(n) (n >= 0).
and 2*A123164(n) = (n+1)*a(n) - (n-1)*a(n-1) (n > 0). (End)
Define the general Delannoy numbers d(i, j) as in A001850. Then a(k) = d(2*k, k) - d(2*k, k-1) and a(0) = 1, Sum_{j=0..n} ((-1)^j * (d(n, j) + d(n-1, j-1)) * a(n-j)) = 0. - Peter E John, Oct 19 2006
Given an integer t >= 1 and initial values u = [a_0, a_1, ..., a_{t-1}], we may define an infinite sequence Phi(u) by setting a_n = a_{n-1} + a_0*a_{n-1} + a_1*a_{n-2} + ... + a_{n-2}*a_1 for n >= t. For example, Phi([1]) is the Catalan numbers A000108. The present sequence is (essentially) Phi([2]). - Gary W. Adamson, Oct 27 2008
G.f.: 1/(1-2x/(1-x/(1-2x/(1-x/(1-2x/(1-x/(1-2x/(1-x/(1-2x/(1-x.... (continued fraction). - Paul Barry, Dec 08 2008
G.f.: 1/(1 - x - x/(1 - x - x/(1 - x - x/(1 - x - x/(1 - x - x/(1 - ... (continued fraction). - Paul Barry, Jan 29 2009
a(n) ~ ((3 + 2*sqrt(2))^n)/(n*sqrt(2*Pi*n)*sqrt(3*sqrt(2) - 4))*(1-(9*sqrt(2) + 24)/(32*n) + ...). - G. Nemes (nemesgery(AT)gmail.com), Jan 25 2009
Logarithmic derivative yields A002003. - Paul D. Hanna, Oct 25 2010
a(n) = the upper left term in M^(n+1), M = the production matrix:
1, 1, 0, 0, 0, 0, ...
1, 1, 1, 0, 0, 0, ...
2, 2, 1, 1, 0, 0, ...
4, 4, 2, 1, 1, 0, ...
8, 8, 8, 2, 1, 1, ...
... - Gary W. Adamson, Jul 08 2011
a(n) is the sum of top row terms in Q^n, Q = an infinite square production matrix as follows:
1, 1, 0, 0, 0, 0, ...
1, 1, 2, 0, 0, 0, ...
1, 1, 1, 2, 0, 0, ...
1, 1, 1, 1, 2, 0, ...
1, 1, 1, 1, 1, 2, ...
... - Gary W. Adamson, Aug 23 2011
From Tom Copeland, Sep 21 2011: (Start)
With F(x) = (1 - 3*x - sqrt(1 - 6*x + x^2))/(2*x) an o.g.f. (nulling the n = 0 term) for A006318, G(x) = x/(2 + 3*x + x^2) is the compositional inverse.
Consequently, with H(x) = 1/ (dG(x)/dx) = (2 + 3*x + x^2)^2 / (2 - x^2),
a(n) = (1/n!)*[(H(x)*d/dx)^n] x evaluated at x = 0, i.e.,
F(x) = exp[x*H(u)*d/du] u, evaluated at u = 0. Also, dF(x)/dx = H(F(x)). (End)
a(n-1) = number of ordered complete binary trees with n leaves having k internal vertices colored black, the remaining n - 1 - k internal vertices colored white, and such that each vertex and its rightmost child have different colors ([Drake, Example 1.6.7]). For a refinement of this sequence see A175124. - Peter Bala, Sep 29 2011
D-finite with recurrence: (n-2)*a(n-2) - 3*(2*n-1)*a(n-1) + (n+1)*a(n) = 0. - Vaclav Kotesovec, Oct 05 2012
G.f.: A(x) = (1 - x - sqrt(1 - 6*x + x^2))/(2*x) = (1 - G(0))/x; G(k) = 1 + x - 2*x/G(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Jan 04 2012
G.f.: A(x) = (1 - x - sqrt(1 - 6*x + x^2))/(2*x) = (G(0) - 1)/x; G(k) = 1 - x/(1 - 2/G(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Jan 04 2012
a(n+1) = a(n) + Sum_{k=0..n} a(k)*(n-k). - Reinhard Zumkeller, Nov 13 2012
G.f.: 1/Q(0) where Q(k) = 1 + k*(1 - x) - x - x*(k+1)*(k+2)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Mar 14 2013
a(-1-n) = a(n). - Michael Somos, Apr 03 2013
G.f.: 1/x - 1 - U(0)/x, where U(k) = 1 - x - x/U(k+1); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013
G.f.: (2 - 2*x - G(0))/(4*x), where G(k) = 1 + 1/( 1 - x*(6 - x)*(2*k - 1)/(x*(6 - x)*(2*k - 1) + 2*(k + 1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013
a(n) = 1/(n + 1) * (Sum_{j=0..n} C(n+j, j)*C(n+j+1, j+1)*(Sum_{k=0..n-j} (-1)^k*C(n+j+k, k))). - Graham H. Hawkes, Feb 15 2015
a(n) = hypergeom([-n, n+1], [2], -1). - Peter Luschny, Mar 23 2015
a(n) = sqrt(2) * LegendreP(n, -1, 3) where LegendreP is the associated Legendre function of the first kind (in Maple's notation). - Robert Israel, Mar 23 2015
G.f. A(x) satisfies: A(x) = Sum_{j>=0} x^j * Sum_{k=0..j} binomial(j,k)*A(x)^k. - Ilya Gutkovskiy, Apr 11 2019
From Peter Bala, May 13 2024: (Start)
a(n) = 2 * Sum_{k = 0..floor(n/2)} binomial(n, 2*k)*binomial(2*n-2*k, n)/(n-2*k+1) for n >= 1.
a(n) = Integral_{x = 0..1} Legendre_P(n, 2*x+1) dx. (End)
G.f. A(x) = 1/(1 - x) * c(x/(1-x)^2), where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108. - Peter Bala, Aug 29 2024

Extensions

Edited by Charles R Greathouse IV, Apr 20 2010

A001109 a(n)^2 is a triangular number: a(n) = 6*a(n-1) - a(n-2) with a(0)=0, a(1)=1.

Original entry on oeis.org

0, 1, 6, 35, 204, 1189, 6930, 40391, 235416, 1372105, 7997214, 46611179, 271669860, 1583407981, 9228778026, 53789260175, 313506783024, 1827251437969, 10650001844790, 62072759630771, 361786555939836, 2108646576008245, 12290092900109634, 71631910824649559, 417501372047787720
Offset: 0

Views

Author

Keywords

Comments

8*a(n)^2 + 1 = 8*A001110(n) + 1 = A055792(n+1) is a perfect square. - Gregory V. Richardson, Oct 05 2002
For n >= 2, A001108(n) gives exactly the positive integers m such that 1,2,...,m has a perfect median. The sequence of associated perfect medians is the present sequence. Let a_1,...,a_m be an (ordered) sequence of real numbers, then a term a_k is a perfect median if Sum_{j=1..k-1} a_j = Sum_{j=k+1..m} a_j. See Puzzle 1 in MSRI Emissary, Fall 2005. - Asher Auel, Jan 12 2006
(a(n), b(n)) where b(n) = A082291(n) are the integer solutions of the equation 2*binomial(b,a) = binomial(b+2,a). - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de); comment revised by Michael Somos, Apr 07 2003
This sequence gives the values of y in solutions of the Diophantine equation x^2 - 8y^2 = 1. It also gives the values of the product xy where (x,y) satisfies x^2 - 2y^2 = +-1, i.e., a(n) = A001333(n)*A000129(n). a(n) also gives the inradius r of primitive Pythagorean triangles having legs whose lengths are consecutive integers, with corresponding semiperimeter s = a(n+1) = {A001652(n) + A046090(n) + A001653(n)}/2 and area rs = A029549(n) = 6*A029546(n). - Lekraj Beedassy, Apr 23 2003 [edited by Jon E. Schoenfield, May 04 2014]
n such that 8*n^2 = floor(sqrt(8)*n*ceiling(sqrt(8)*n)). - Benoit Cloitre, May 10 2003
For n > 0, ratios a(n+1)/a(n) may be obtained as convergents to continued fraction expansion of 3+sqrt(8): either successive convergents of [6;-6] or odd convergents of [5;1, 4]. - Lekraj Beedassy, Sep 09 2003
a(n+1) + A053141(n) = A001108(n+1). Generating floretion: - 2'i + 2'j - 'k + i' + j' - k' + 2'ii' - 'jj' - 2'kk' + 'ij' + 'ik' + 'ji' + 'jk' - 2'kj' + 2e ("jes" series). - Creighton Dement, Dec 16 2004
Kekulé numbers for certain benzenoids (see the Cyvin-Gutman reference). - Emeric Deutsch, Jun 19 2005
Number of D steps on the line y=x in all Delannoy paths of length n (a Delannoy path of length n is a path from (0,0) to (n,n), consisting of steps E=(1,0), N=(0,1) and D=(1,1)). Example: a(2)=6 because in the 13 (=A001850(2)) Delannoy paths of length 2, namely (DD), (D)NE, (D)EN, NE(D), NENE, NEEN, NDE, NNEE, EN(D), ENNE, ENEN, EDN and EENN, we have altogether six D steps on the line y=x (shown between parentheses). - Emeric Deutsch, Jul 07 2005
Define a T-circle to be a first-quadrant circle with integral radius that is tangent to the x- and y-axes. Such a circle has coordinates equal to its radius. Let C(0) be the T-circle with radius 1. Then for n > 0, define C(n) to be the smallest T-circle that does not intersect C(n-1). C(n) has radius a(n+1). Cf. A001653. - Charlie Marion, Sep 14 2005
Numbers such that there is an m with t(n+m)=2t(m), where t(n) are the triangular numbers A000217. For instance, t(20)=2*t(14)=210, so 6 is in the sequence. - Floor van Lamoen, Oct 13 2005
One half the bisection of the Pell numbers (A000129). - Franklin T. Adams-Watters, Jan 08 2006
Pell trapezoids: for n > 0, a(n) = (A000129(n-1)+A000129(n+1))*A000129(n)/2; see also A084158. - Charlie Marion, Apr 01 2006
Tested for 2 < p < 27: If and only if 2^p - 1 (the Mersenne number M(p)) is prime then M(p) divides a(2^(p-1)). - Kenneth J Ramsey, May 16 2006
If 2^p - 1 is prime then M(p) divides a(2^(p-1)-1). - Kenneth J Ramsey, Jun 08 2006; comment corrected by Robert Israel, Mar 18 2007
If 8*n+5 and 8*n+7 are twin primes then their product divides a(4*n+3). - Kenneth J Ramsey, Jun 08 2006
If p is an odd prime, then if p == 1 or 7 (mod 8), then a((p-1)/2) == 0 (mod p) and a((p+1)/2) == 1 (mod p); if p == 3 or 5 (mod 8), then a((p-1)/2) == 1 (mod p) and a((p+1)/2) == 0 (mod p). Kenneth J Ramsey's comment about twin primes follows from this. - Robert Israel, Mar 18 2007
a(n)*(a(n+b) - a(b-2)) = (a(n+1)+1)*(a(n+b-1) - a(b-1)). This identity also applies to any series a(0) = 0 a(1) = 1 a(n) = b*a(n-1) - a(n-2). - Kenneth J Ramsey, Oct 17 2007
For n < 0, let a(n) = -a(-n). Then (a(n+j) + a(k+j)) * (a(n+b+k+j) - a(b-j-2)) = (a(n+j+1) + a(k+j+1)) * (a(n+b+k+j-1) - a(b-j-1)). - Charlie Marion, Mar 04 2011
Sequence gives y values of the Diophantine equation: 0+1+2+...+x = y^2. If (a,b) and (c,d) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with aMohamed Bouhamida, Aug 29 2009
If (p,q) and (r,s) are two consecutive solutions of the Diophantine equation: 0+1+2+...+x = y^2 with p < r then r = 3*p+4*q+1 and s = 2*p+3*q+1. - Mohamed Bouhamida, Sep 02 2009
a(n)/A002315(n) converges to cos^2(Pi/8) (see A201488). - Gary Detlefs, Nov 25 2009
Binomial transform of A086347. - Johannes W. Meijer, Aug 01 2010
If x=a(n), y=A055997(n+1) and z = x^2+y, then x^4 + y^3 = z^2. - Bruno Berselli, Aug 24 2010
In general, if b(0)=1, b(1)=k and for n > 1, b(n) = 6*b(n-1) - b(n-2), then
for n > 0, b(n) = a(n)*k-a(n-1); e.g.,
for k=2, when b(n) = A038725(n), 2 = 1*2 - 0, 11 = 6*2 - 1, 64 = 35*2 - 6, 373 = 204*2 - 35;
for k=3, when b(n) = A001541(n), 3 = 1*3 - 0, 17 = 6*3 - 1; 99 = 35*3 - 6; 577 = 204*3 - 35;
for k=4, when b(n) = A038723(n), 4 = 1*4 - 0, 23 = 6*4 - 1; 134 = 35*4 - 6; 781 = 204*4 - 35;
for k=5, when b(n) = A001653(n), 5 = 1*5 - 0, 29 = 6*5 - 1; 169 = 35*5 - 6; 985 = 204*5 - 35.
- Charlie Marion, Dec 08 2010
See a Wolfdieter Lang comment on A001653 on a sequence of (u,v) values for Pythagorean triples (x,y,z) with x=|u^2-v^2|, y=2*u*v and z=u^2+v^2, with u odd and v even, generated from (u(0)=1,v(0)=2), the triple (3,4,5), by a substitution rule given there. The present a(n) appears there as b(n). The corresponding generated triangles have catheti differing by one length unit. - Wolfdieter Lang, Mar 06 2012
a(n)*a(n+2k) + a(k)^2 and a(n)*a(n+2k+1) + a(k)*a(k+1) are triangular numbers. Generalizes description of sequence. - Charlie Marion, Dec 03 2012
a(n)*a(n+2k) + a(k)^2 is the triangular square A001110(n+k). a(n)*a(n+2k+1) + a(k)*a(k+1) is the triangular oblong A029549(n+k). - Charlie Marion, Dec 05 2012
From Richard R. Forberg, Aug 30 2013: (Start)
The squares of a(n) are the result of applying triangular arithmetic to the squares, using A001333 as the "guide" on what integers to square, as follows:
a(2n)^2 = A001333(2n)^2 * (A001333(2n)^2 - 1)/2;
a(2n+1)^2 = A001333(2n+1)^2 * (A001333(2n+1)^2 + 1)/2. (End)
For n >= 1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,5}. - Milan Janjic, Jan 25 2015
Panda and Rout call these "balancing numbers" and note that the period of the sequence modulo a prime p is the same as that modulo p^2 when p = 13, 31, 1546463. But these are precisely the p in A238736 such that p^2 divides A000129(p - (2/p)), where (2/p) is a Jacobi symbol. In light of the above observation by Franklin T. Adams-Watters that the present sequence is one half the bisection of the Pell numbers, i.e., a(n) = A000129(2*n)/2, it follows immediately that modulo a fixed prime p, or any power thereof, the period of a(n) is half that of A000129(n). - John Blythe Dobson, Mar 06 2015
The triangular number = square number identity Tri((T(n, 3) - 1)/2) = S(n-1, 6)^2 with Tri, T, and S given in A000217, A053120 and A049310, is the special case k = 1 of the k-family of identities Tri((T(n, 2*k+1) - 1)/2) = Tri(k)*S(n-1, 2*(2*k+1))^2, k >= 0, n >= 0, with S(-1, x) = 0. For k=2 see A108741(n) for S(n-1, 10)^2. This identity boils down to the identities S(n-1, 2*x)^2 = (T(2*n, x) - 1)/(2*(x^2-1)) and 2*T(n, x)^2 - 1 = T(2*n, x) with x = 2*k+1. - Wolfdieter Lang, Feb 01 2016
a(2)=6 is perfect. For n=2*k, k > 0, k not equal to 1, a(n) is a multiple of a(2) and since every multiple (beyond 1) of a perfect number is abundant, then a(n) is abundant. sigma(a(4)) = 504 > 408 = 2*a(4). For n=2*k+1, k > 0, a(n) mod 10 = A000012(n), so a(n) is odd. If a(n) is a prime number, it is deficient; otherwise a(n) has one or two distinct prime factors and is therefore deficient again. So for n=2k+1, k > 0, a(n) is deficient. sigma(a(5)) = 1260 < 2378 = 2*a(5). - Muniru A Asiru, Apr 14 2016
Behera & Panda call these the balancing numbers, and A001541 are the balancers. - Michel Marcus, Nov 07 2017
In general, a second-order linear recurrence with constant coefficients having a signature of (c,d) will be duplicated by a third-order recurrence having a signature of (x,c^2-c*x+d,-d*x+c*d). The formulas of Olivares and Bouhamida in the formula section which have signatures of (7,-7,1) and (5,5,-1), respectively, are specific instances of this general rule for x = 7 and x = 5. - Gary Detlefs, Jan 29 2021
Note that 6 is the largest triangular number in the sequence, because it is proved that 8 and 9 are the largest perfect powers which are consecutive (Catalan's conjecture). 0 and 1 are also in the sequence because they are also perfect powers and 0*1/2 = 0^2 and 8*9/2 = (2*3)^2. - Metin Sariyar, Jul 15 2021

Examples

			G.f. = x + 6*x^2 + 35*x^3 + 204*x^4 + 1189*x^5 + 6930*x^6 + 40391*x^7 + ...
6 is in the sequence since 6^2 = 36 is a triangular number: 36 = 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8. - _Michael B. Porter_, Jul 02 2016
		

References

  • Julio R. Bastida, Quadratic properties of a linearly recurrent sequence. Proceedings of the Tenth Southeastern Conference on Combinatorics, Graph Theory and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1979), pp. 163--166, Congress. Numer., XXIII-XXIV, Utilitas Math., Winnipeg, Man., 1979. MR0561042 (81e:10009) - From N. J. A. Sloane, May 30 2012
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, pp. 193, 197.
  • D. M. Burton, The History of Mathematics, McGraw Hill, (1991), p. 213.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 10.
  • P. Franklin, E. F. Beckenbach, H. S. M Coxeter, N. H. McCoy, K. Menger, and J. L. Synge, Rings And Ideals, No 8, The Carus Mathematical Monographs, The Mathematical Association of America, (1967), pp. 144-146.
  • A. Patra, G. K. Panda, and T. Khemaratchatakumthorn. "Exact divisibility by powers of the balancing and Lucas-balancing numbers." Fibonacci Quart., 59:1 (2021), 57-64; see B(n).
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 257-258.
  • P.-F. Teilhet, Query 2376, L'Intermédiaire des Mathématiciens, 11 (1904), 138-139. - N. J. A. Sloane, Mar 08 2022

Crossrefs

Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), this sequence (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), A007655 (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.

Programs

  • GAP
    a:=[0,1];; for n in [3..25] do a[n]:=6*a[n-1]-a[n-2]; od; a; # Muniru A Asiru, Dec 18 2018
  • Haskell
    a001109 n = a001109_list !! n :: Integer
    a001109_list = 0 : 1 : zipWith (-)
       (map (* 6) $ tail a001109_list) a001109_list
    -- Reinhard Zumkeller, Dec 17 2011
    
  • Magma
    [n le 2 select n-1 else 6*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Jul 25 2015
    
  • Maple
    a[0]:=1: a[1]:=6: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n],n=0..26); # Emeric Deutsch
    with (combinat):seq(fibonacci(2*n,2)/2, n=0..20); # Zerinvary Lajos, Apr 20 2008
  • Mathematica
    Transpose[NestList[Flatten[{Rest[#],ListCorrelate[{-1,6},#]}]&, {0,1}, 30]][[1]]  (* Harvey P. Dale, Mar 23 2011 *)
    CoefficientList[Series[x/(1-6x+x^2),{x,0,30}],x]  (* Harvey P. Dale, Mar 23 2011 *)
    LinearRecurrence[{6, -1}, {0, 1}, 50] (* Vladimir Joseph Stephan Orlovsky, Feb 12 2012 *)
    a[ n_]:= ChebyshevU[n-1, 3]; (* Michael Somos, Sep 02 2012 *)
    Table[Fibonacci[2n, 2]/2, {n, 0, 20}] (* Vladimir Reshetnikov, Sep 16 2016 *)
    TrigExpand@Table[Sinh[2 n ArcCsch[1]]/(2 Sqrt[2]), {n, 0, 10}] (* Federico Provvedi, Feb 01 2021 *)
  • PARI
    {a(n) = imag((3 + quadgen(32))^n)}; /* Michael Somos, Apr 07 2003 */
    
  • PARI
    {a(n) = subst( poltchebi( abs(n+1)) - 3 * poltchebi( abs(n)), x, 3) / 8}; /* Michael Somos, Apr 07 2003 */
    
  • PARI
    {a(n) = polchebyshev( n-1, 2, 3)}; /* Michael Somos, Sep 02 2012 */
    
  • PARI
    is(n)=ispolygonal(n^2,3) \\ Charles R Greathouse IV, Nov 03 2016
    
  • Sage
    [lucas_number1(n,6,1) for n in range(27)] # Zerinvary Lajos, Jun 25 2008
    
  • Sage
    [chebyshev_U(n-1,3) for n in (0..20)] # G. C. Greubel, Dec 23 2019
    

Formula

G.f.: x / (1 - 6*x + x^2). - Simon Plouffe in his 1992 dissertation.
a(n) = S(n-1, 6) = U(n-1, 3) with U(n, x) Chebyshev's polynomials of the second kind. S(-1, x) := 0. Cf. triangle A049310 for S(n, x).
a(n) = sqrt(A001110(n)).
a(n) = A001542(n)/2.
a(n) = sqrt((A001541(n)^2-1)/8) (cf. Richardson comment).
a(n) = 3*a(n-1) + sqrt(8*a(n-1)^2+1). - R. J. Mathar, Oct 09 2000
a(n) = A000129(n)*A001333(n) = A000129(n)*(A000129(n)+A000129(n-1)) = ceiling(A001108(n)/sqrt(2)). - Henry Bottomley, Apr 19 2000
a(n) ~ (1/8)*sqrt(2)*(sqrt(2) + 1)^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
Limit_{n->oo} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 05 2002
a(n) = ((3 + 2*sqrt(2))^n - (3 - 2*sqrt(2))^n) / (4*sqrt(2)). - Gregory V. Richardson, Oct 13 2002. Corrected for offset 0, and rewritten. - Wolfdieter Lang, Feb 10 2015
a(2*n) = a(n)*A003499(n). 4*a(n) = A005319(n). - Mario Catalani (mario.catalani(AT)unito.it), Mar 21 2003
a(n) = floor((3+2*sqrt(2))^n/(4*sqrt(2))). - Lekraj Beedassy, Apr 23 2003
a(-n) = -a(n). - Michael Somos, Apr 07 2003
For n >= 1, a(n) = Sum_{k=0..n-1} A001653(k). - Charlie Marion, Jul 01 2003
For n > 0, 4*a(2*n) = A001653(n)^2 - A001653(n-1)^2. - Charlie Marion, Jul 16 2003
For n > 0, a(n) = Sum_{k = 0..n-1}((2*k+1)*A001652(n-1-k)) + A000217(n). - Charlie Marion, Jul 18 2003
a(2*n+1) = a(n+1)^2 - a(n)^2. - Charlie Marion, Jan 12 2004
a(k)*a(2*n+k) = a(n+k)^2 - a(n)^2; e.g., 204*7997214 = 40391^2 - 35^2. - Charlie Marion, Jan 15 2004
For j < n+1, a(k+j)*a(2*n+k-j) - Sum_{i = 0..j-1} a(2*n-(2*i+1)) = a(n+k)^2 - a(n)^2. - Charlie Marion, Jan 18 2004
From Paul Barry, Feb 06 2004: (Start)
a(n) = A000129(2*n)/2;
a(n) = ((1+sqrt(2))^(2*n) - (1-sqrt(2))^(2*n))*sqrt(2)/8;
a(n) = Sum_{i=0..n} Sum_{j=0..n} A000129(i+j)*n!/(i!*j!*(n-i-j)!)/2. (End)
E.g.f.: exp(3*x)*sinh(2*sqrt(2)*x)/(2*sqrt(2)). - Paul Barry, Apr 21 2004
A053141(n+1) + A055997(n+1) = A001541(n+1) + a(n+1). - Creighton Dement, Sep 16 2004
a(n) = Sum_{k=0..n} binomial(2*n, 2*k+1)*2^(k-1). - Paul Barry, Oct 01 2004
a(n) = A001653(n+1) - A038723(n); (a(n)) = chuseq[J]( 'ii' + 'jj' + .5'kk' + 'ij' - 'ji' + 2.5e ), apart from initial term. - Creighton Dement, Nov 19 2004, modified by Davide Colazingari, Jun 24 2016
a(n+1) = Sum_{k=0..n} A001850(k)*A001850(n-k), self convolution of central Delannoy numbers. - Benoit Cloitre, Sep 28 2005
a(n) = 7*(a(n-1) - a(n-2)) + a(n-3), a(1) = 0, a(2) = 1, a(3) = 6, n > 3. Also a(n) = ( (1 + sqrt(2) )^(2*n) - (1 - sqrt(2) )^(2*n) ) / (4*sqrt(2)). - Antonio Alberto Olivares, Oct 23 2003
a(n) = 5*(a(n-1) + a(n-2)) - a(n-3). - Mohamed Bouhamida, Sep 20 2006
Define f(x,s) = s*x + sqrt((s^2-1)*x^2+1); f(0,s)=0. a(n) = f(a(n-1),3), see second formula. - Marcos Carreira, Dec 27 2006
The perfect median m(n) can be expressed in terms of the Pell numbers P() = A000129() by m(n) = P(n + 2) * (P(n + 2) + P(n + 1)) for n >= 0. - Winston A. Richards (ugu(AT)psu.edu), Jun 11 2007
For k = 0..n, a(2*n-k) - a(k) = 2*a(n-k)*A001541(n). Also, a(2*n+1-k) - a(k) = A002315(n-k)*A001653(n). - Charlie Marion, Jul 18 2007
[A001653(n), a(n)] = [1,4; 1,5]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(n) = Sum_{k=0..n-1} 4^k*binomial(n+k,2*k+1). - Paul Barry, Apr 20 2009
a(n+1)^2 - 6*a(n+1)*a(n) + a(n)^2 = 1. - Charlie Marion, Dec 14 2010
a(n) = A002315(m)*A011900(n-m-1) + A001653(m)*A001652(n-m-1) - a(m) = A002315(m)*A053141(n-m-1) + A001653(m)*A046090(n-m-1) + a(m) with m < n; otherwise a(n) = A002315(m)*A053141(m-n) - A001653(m)*A011900(m-n) + a(m) = A002315(m)*A053141(m-n) - A001653(m)*A046090(m-n) - a(m) = (A002315(n) - A001653(n))/2. - Kenneth J Ramsey, Oct 12 2011
16*a(n)^2 + 1 = A056771(n). - James R. Buddenhagen, Dec 09 2011
A010054(A000290(a(n))) = 1. - Reinhard Zumkeller, Dec 17 2011
In general, a(n+k)^2 - A003499(k)*a(n+k)*a(n) + a(n)^2 = a(k)^2. - Charlie Marion, Jan 11 2012
a(n+1) = Sum_{k=0..n} A101950(n,k)*5^k. - Philippe Deléham, Feb 10 2012
PSUM transform of a(n+1) is A053142. PSUMSIGN transform of a(n+1) is A084158. BINOMIAL transform of a(n+1) is A164591. BINOMIAL transform of A086347 is a(n+1). BINOMIAL transform of A057087(n-1). - Michael Somos, May 11 2012
a(n+k) = A001541(k)*a(n) + sqrt(A132592(k)*a(n)^2 + a(k)^2). Generalizes formula dated Oct 09 2000. - Charlie Marion, Nov 27 2012
a(n) + a(n+2*k) = A003499(k)*a(n+k); a(n) + a(n+2*k+1) = A001653(k+1)*A002315(n+k). - Charlie Marion, Nov 29 2012
From Peter Bala, Dec 23 2012: (Start)
Product_{n >= 1} (1 + 1/a(n)) = 1 + sqrt(2).
Product_{n >= 2} (1 - 1/a(n)) = (1/3)*(1 + sqrt(2)). (End)
G.f.: G(0)*x/(2-6*x), where G(k) = 1 + 1/(1 - x*(8*k-9)/( x*(8*k-1) - 3/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 12 2013
G.f.: H(0)*x/2, where H(k) = 1 + 1/( 1 - x*(6-x)/(x*(6-x) + 1/H(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Feb 18 2014
a(n) = (a(n-1)^2 - a(n-3)^2)/a(n-2) + a(n-4) for n > 3. - Patrick J. McNab, Jul 24 2015
a(n-k)*a(n+k) + a(k)^2 = a(n)^2, a(n+k) + a(n-k) = A003499(k)*a(n), for n >= k >= 0. - Alexander Samokrutov, Sep 30 2015
Dirichlet g.f.: (PolyLog(s,3+2*sqrt(2)) - PolyLog(s,3-2*sqrt(2)))/(4*sqrt(2)). - Ilya Gutkovskiy, Jun 27 2016
4*a(n)^2 - 1 = A278310(n) for n > 0. - Bruno Berselli, Nov 24 2016
From Klaus Purath, Jan 18 2020: (Start)
a(n) = (a(n-3) + a(n+3))/198.
a(n) = Sum_{i=1..n} A001653(i), n>=1.
a(n) = sinh( 2 * n * arccsch(1) ) / ( 2 * sqrt(2) ). - Federico Provvedi, Feb 01 2021
(End)
a(n) = A002965(2*n)*A002965(2*n+1). - Jon E. Schoenfield, Jan 08 2022
a(n) = A002965(4*n)/2. - Gerry Martens, Jul 14 2023
a(n) = Sum_{k = 0..n-1} (-1)^(n+k+1)*binomial(n+k, 2*k+1)*8^k. - Peter Bala, Jul 17 2023

Extensions

Additional comments from Wolfdieter Lang, Feb 10 2000
Duplication of a formula removed by Wolfdieter Lang, Feb 10 2015
Previous Showing 11-20 of 202 results. Next