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 51-60 of 248 results. Next

A047781 a(n) = Sum_{k=0..n-1} binomial(n-1,k)*binomial(n+k,k). Also a(n) = T(n,n), array T as in A049600.

Original entry on oeis.org

0, 1, 4, 19, 96, 501, 2668, 14407, 78592, 432073, 2390004, 13286043, 74160672, 415382397, 2333445468, 13141557519, 74174404608, 419472490257, 2376287945572, 13482186743203, 76598310928096, 435730007006341, 2481447593848524, 14146164790774359
Offset: 0

Views

Author

Keywords

Comments

Also main diagonal of array: m(i,1)=1, m(1,j)=j, m(i,j)=m(i,j-1)+m(i-1,j-1)+m(i-1,j): 1 2 3 4 ... / 1 4 9 16 ... / 1 6 19 44 ... / 1 8 33 96 ... /. - Benoit Cloitre, Aug 05 2002
This array is now listed as A142978, where some conjectural congruences for the present sequence are given. - Peter Bala, Nov 13 2008
Convolution of central Delannoy numbers A001850 and little Schroeder numbers A001003. Hankel transform is 2^C(n+1,2)*A007052(n). - Paul Barry, Oct 07 2009
Define a finite triangle T(r,c) with T(r,0) = binomial(n,r) for 0 <= r <= n and the other terms recursively with T(r,c) = T(r-1,c-1) + 2*T(r,c-1). The sum of the last terms in the rows is Sum_{r=0..n} T(r,r) = a(n+1). Example: For n=4 the triangle has the rows 1; 4 9; 6 16 41; 4 14 44 129; 1 6 26 96 321 having sum of last terms 1 + 9 + 41 + 129 + 321 = 501 = a(5). - J. M. Bergot, Feb 15 2013
a(n) = A049600(2*n,n), when A049600 is seen as a triangle read by rows. - Reinhard Zumkeller, Apr 15 2014
a(n-1) for n > 1 is the number of assembly trees with the connected gluing rule for cycle graphs with n vertices. - Nick Mayers, Aug 16 2018

Crossrefs

Cf. A002003. Column 1 of A296129.

Programs

  • Haskell
    a047781 n = a049600 (2 * n) n  -- Reinhard Zumkeller, Apr 15 2014
    
  • Magma
    [n eq 0 select 0 else &+[Binomial(n-1, k)*Binomial(n+k, k): k in [0..n-1]]: n in [0..22]];  // Bruno Berselli, May 19 2011
    
  • Maple
    a := proc(n) local k; add(binomial(n-1,k)*binomial(n+k,k),k=0..n-1); end;
  • Mathematica
    Table[SeriesCoefficient[x*((1+x)-Sqrt[1-6*x+x^2])/(4*x*Sqrt[1-6*x+x^2]),{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 08 2012 *)
    a[n_] := Hypergeometric2F1[1-n, n+1, 1, -1]; a[0] = 0; Table[a[n], {n, 0, 23}] (* Jean-François Alcover, Feb 26 2013 *)
    a[n_] := Sum[ Binomial[n - 1, k] Binomial[n + k, k], {k, 0, n - 1}]; Array[a, 25] (* Robert G. Wilson v, Aug 08 2018 *)
  • Maxima
    makelist(if n=0 then 0 else sum(binomial(n-1, k)*binomial(n+k, k), k, 0, n-1), n, 0, 22); /* Bruno Berselli, May 19 2011 */
    
  • PARI
    A047781(n)=polcoeff((1+x)/sqrt(1+(O(x^n)-6)*x+x^2),n)\4  \\ M. F. Hasler, Oct 09 2012
    
  • Python
    from sympy import binomial
    def a(n):
        return sum(binomial(n - 1, k) * binomial(n + k, k) for k in range(n))
    print([a(n) for n in range(51)]) # Indranil Ghosh, Apr 18 2017
    
  • Python
    from math import comb
    def A047781(n): return sum(comb(n,k)**2*k<Chai Wah Wu, Mar 22 2023

Formula

D-finite with recurrence n*(2*n-3)*a(n) - (12*n^2-24*n+8)*a(n-1) + (2*n-1)*(n-2)*a(n-2) = 0. - Vladeta Jovovic, Aug 29 2004
a(n+1) = Sum_{k=0..n} binomial(n, k)*binomial(n+1, k+1)*2^k. - Paul Barry, Sep 20 2004
a(n) = Sum_{k=0..n} T(n, k), array T as in A008288.
If shifted one place left, the third binomial transform of A098660. - Paul Barry, Sep 20 2004
G.f.: ((1+x)/sqrt(1-6x+x^2)-1)/4. - Paul Barry, Sep 20 2004, simplified by M. F. Hasler, Oct 09 2012
E.g.f. for sequence shifted left: Sum_{n>=0} a(n+1)*x^n/n! = exp(3*x)*(BesselI(0, 2*sqrt(2)*x)+BesselI(1, 2*sqrt(2)*x)/sqrt(2)). - Paul Barry, Sep 20 2004
a(n) = Sum_{k=0..n-1} C(n,k)*C(n-1,k)*2^(n-k-1); a(n+1) = 2^n*Hypergeometric2F1(-n,-n-1;1;1/2). - Paul Barry, Feb 08 2011
a(n) ~ 2^(1/4)*(3+2*sqrt(2))^n/(4*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 08 2012
Recurrence (an alternative): n*a(n) = (6-n)*a(n-6) + 2*(5*n-27)*a(n-5) + (84-15*n)*a(n-4) + 52*(3-n)*a(n-3) + 3*(2-5*n)*a(n-2) + 2*(5*n-3)*a(n-1), n >= 7. - Fung Lam, Feb 05 2014
a(n) = A241023(n) / 4. - Reinhard Zumkeller, Apr 15 2014
a(n) = Hyper2F1([-n, n], [1], -1)/2 for n > 0. - Peter Luschny, Aug 02 2014
n^2*a(n) = Sum_{k=0..n-1} (2*k^2+2*k+1)*binomial(n-1,k)*binomial(n+k,k). By the Zeilberger algorithm, both sides of the equality satisfy the same recurrence. - Zhi-Wei Sun, Aug 30 2014
a(n) = [x^n] (1/2) * ((1+x)/(1-x))^n for n > 0. - Seiichi Manyama, Jun 07 2018

A144097 The 4-Schroeder numbers: a(n) = number of lattice paths (Schroeder paths) from (0,0) to (3n,n) with unit steps N=(0,1), E=(1,0) and D=(1,1) staying weakly above the line y = 3x.

Original entry on oeis.org

1, 2, 14, 134, 1482, 17818, 226214, 2984206, 40503890, 561957362, 7934063678, 113622696470, 1646501710362, 24098174350986, 355715715691350, 5289547733908510, 79163575684710818, 1191491384838325474
Offset: 0

Views

Author

Joachim Schroeder (schroderjd(AT)qwa.uovs.ac.za), Sep 10 2008

Keywords

Comments

a(n) is also the number of lattice path from (0,0) to (4n,0) with unit steps (1,3), (2,2) and (1,-1) staying weakly above the x-axis.
Also, the number of planar rooted trees with n non-leaf vertices such that each non-leaf vertex has either 3 or 4 children. - Cameron Marcott, Sep 18 2013
a(n) equals the number of ordered complete 4-ary trees with 3*n + 1 leaves, where the internal vertices come in two colors and such that each vertex and its rightmost child have different colors. See Drake, Example 1.6.9. - Peter Bala, Apr 30 2023

Examples

			a(2)=14, because
  01: NNNENNNE,
  02: NNDNNNE,
  03: NNNENND,
  04: NNDNND,
  05: NNNDNNE,
  06: NNNDND,
  07: NNNNENNE,
  08: NNNNEND,
  09: NNNNDNE,
  10: NNNNDD,
  11: NNNNNENE,
  12: NNNNNED,
  13: NNNNNDE,
  14: NNNNNNEE
are all the paths from (0,0) to (2,6) with steps N,E and D weakly above y=3x.
		

References

  • 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

Cf. A027307 (the case y=2x), A008288 (Delannoy numbers), A008412 (4-dimensional coordination numbers).
This appears to equal 2*A243675. - N. J. A. Sloane, Mar 28 2021
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

  • Maple
    Schr:=proc(n,m,l)(n-l*m+1)/m*sum(2^v*binomial(m,v)*binomial(n,v-1),v=1..m) end proc; where n=3m and l=3, also
    Schr:=proc(n,m,l)(n-l*m+1)/(n+1)*sum(2^v*binomial(m-1,v-1)*binomial(n+1,v),v=0..m) end proc; where n=3m and l=3, also
    Schr:=proc(n,m,l)(n-l*m+1)/m*sum(binomial(m,v)*binomial(n+v,m-1),v=0..m) end proc; where n=3m and l=3, also
    Schr:=proc(n,m,l)(n-l*m+1)/(n+1)*sum(binomial(n+1,v)*binomial(m-1+v,n),v=0..n+1) end proc; where n=3m and l=3.
    # alternative Maple program:
    a:= proc(n) option remember; `if`(n<2, n+1,
          ((15610*n^5 -67123*n^4 +106824*n^3 -77633*n^2
           +25514*n-3000)*a(n-1) -(3*(n-2))*(3*n-4)*
           (3*n-5)*(35*n^2-28*n+5)*a(n-2)) / ((3*(3*n-1))
           *(3*n+1)*n*(35*n^2-98*n+68)))
        end:
    seq(a(n), n=0..20);  # Alois P. Heinz, May 26 2015
  • Mathematica
    d[n_, k_] := Binomial[n+k, k] Hypergeometric2F1[-k, -n, -n-k, -1]; a[0] = 1; a[n_] = d[3n, n] - 3d[3n+1, n-1] - 2d[3n, n-1]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Feb 25 2017 *)
  • PARI
    {a(n) = sum(k=0, n, binomial(n, k) * binomial(3*n+k+1, n)/(3*n+k+1))} \\ Seiichi Manyama, Jul 25 2020
    
  • PARI
    {a(n) = if(n==0, 1, sum(k=1, n, 2^k*binomial(n, k) * binomial(3*n, k-1)/n))} \\ Seiichi Manyama, Jul 25 2020

Formula

G.f. A(z) satisfies A(z) = 1 + z(A(z)^3 + A(z)^4) a(n)= S_{3n+1}(n) - 3S_n(3n + 1), where S_a(b) are coordination numbers, i.e., the number of points in the a-dimensional cubic lattice Z^a having distance b in the L_1 norm.
Also a(n) = D(3n,n) - 3D(3n + 1,n-1) - 2D(3n,n-1), where D(a,b) are the Delannoy numbers, i.e., the number of paths with N, E and D steps from (0,0) to (a,b).
D-finite with recurrence 3*n*(3*n-1)*(3*n+1)*(35*n^2-98*n+68) *a(n) +(-15610*n^5+67123*n^4-106824*n^3+77633*n^2-25514*n+3000)*a(n-1) +3*(n-2) *(3*n-4) *(3*n-5) *(35*n^2-28*n+5) *a(n-2)=0. - R. J. Mathar, Sep 06 2016
From Seiichi Manyama, Jul 25 2020: (Start)
a(n) = Sum_{k=0..n} binomial(n,k) * binomial(3*n+k+1, n)/(3*n+k+1).
a(n) = (1/n) * Sum_{k=1..n} 2^k * binomial(n,k) * binomial(3*n,k-1) for n > 0. (End)
a(n) ~ sqrt(12160 + 3853*sqrt(10)) * 3^(3*n - 9/2) / (2*sqrt(5*Pi) * n^(3/2) * (223 - 70*sqrt(10))^(n - 1/2)). - Vaclav Kotesovec, Jul 31 2021
Series reversion of x*(1 - x^3)/(1 + x^3) = x + 2*x^4 + 14*x^7 + 134*x^10 + ... = Sum_{n >= 0} a(n)*x^(3*n+1). - Peter Bala, Apr 30 2023
From Peter Bala, Jun 16 2023: (Start)
The g.f. A(x) = 1 + 2*x + 14*x^2 + 134*x^3 + ... satisfies A(x)^3 = (1/x) * the series reversion of ((1 - x)/(1 + x))^3.
Define b(n) = [x^(3*n)] ( (1 + x)/(1 - x) )^n = (1/3) * [x^n] ((1 + x)/(1 - x))^(3*n) = A333715(n). Then A(x) = exp( Sum_{n >= 1} b(n)*x^n/n ).
a(n) = 2*hypergeom([1 - n, -3*n], [2], 2) for n >= 1. (End)
a(n) = (1/n) * Sum_{k=0..n-1} (-1)^k * 2^(n-k) * binomial(n,k) * binomial(4*n-k,n-1-k) for n > 0. - Seiichi Manyama, Aug 09 2023

A260332 Labelings of n diamond-shaped posets with 4 vertices per diamond where the labels follow the poset relations whose associated reading permutation avoids 231 in the classical sense.

Original entry on oeis.org

1, 2, 18, 226, 3298, 52450, 881970
Offset: 0

Views

Author

Manda Riehl, Jul 29 2015

Keywords

Comments

According to Yang-Jiang (2021) these are the 5-Schroeder numbers. If confirmed, this will prove Michael Weiner's conjectures and enable us to extend the sequence. Yang & Jiang (2021) give an explicit formula for the m-Schroeder numbers in Theorem 2.4. - N. J. A. Sloane, Mar 28 2021
By diamond-shaped poset with 4 vertices, we mean a poset on four elements with e_1 < e_2, e_1 < e_3, e_2 < e_4, e_3 < e_4, and no order relations between e_2 and e_3. In the Hasse diagram for such a poset, we have a least element, two elements in the level above, and one element in the top level, so the diagram resembles a diamond. The associated permutation is the permutation obtained by reading the labels of each poset by levels left to right, starting with the least element.
Also the number of labelings of n diamond-shaped posets with 4 vertices per diamond where the labels follow the poset relations whose associated reading permutation avoids 312 in the classical sense via reverse complement Wilf equivalence.
Conjecture: Also the number of lattice paths (Schroeder paths) from (0,0) to (n,4n) with unit steps N=(0,1), E=(1,0) and D=(1,1) staying weakly above the line y = 4x. - Michael D. Weiner, Jul 24 2019

Examples

			For a single diamond (n=1) poset with 4 vertices, we must label the least element 1 and the greatest element 4, and the two central elements can be labeled either 2, 3 or 3, 2 respectively. Thus the associated permutations are 1234 and 1324.
		

References

  • 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

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

Formula

There is a complicated recursive formula available in Paukner et al.
Yang & Jiang (2021) give an explicit formula for the 5-Schroeder numbers in Theorem 2.4. - N. J. A. Sloane, Mar 28 2021
Conjecture: a(n) = Sum_{k=1..n} binomial(n,k)*binomial(4*n,k-1)*2^k/n for n > 0. - Michael D. Weiner, Jul 23 2019
From Peter Bala, Jun 16 2023: (Start)
Conjectures: 1) the g.f. A(x) = 1 + 2*x + 18*x^2 + 226*x^3 + ... satisfies A(x)^4 = (1/x) * the series reversion of ((1 - x)/(1 + x))^4.
2) Define b(n) = (1/4) * [x^n] ((1 + x)/(1 - x))^(4*n). Then A(x) = exp( Sum_{n >= 1} b(n)*x^n/n ).
3) a(n) = 2 * hypergeom([1 - n, -4*n], [2], 2) for n >= 1 (equivalent to Weiner's conjecture above).
4) [x^n] A(x)^n = (2*n) * hypergeom([1 - n, 1 - 5*n], [2], 2) for n >= 1. (End)

A046736 Number of ways to place non-intersecting diagonals in convex n-gon so as to create no triangles.

Original entry on oeis.org

1, 0, 1, 1, 4, 8, 25, 64, 191, 540, 1616, 4785, 14512, 44084, 135545, 418609, 1302340, 4070124, 12785859, 40325828, 127689288, 405689020, 1293060464, 4133173256, 13246527139, 42557271268, 137032656700, 442158893833, 1429468244788
Offset: 2

Views

Author

Keywords

Examples

			a(4)=a(5)=1 because of null placement; a(6)=4 because in addition to not placing any, we might also place one between any of the 3 pairs of opposite vertices.
		

Crossrefs

Cf. A001003 (Schroeder), A001006 (Motzkin), A000108 (Catalan), A052524.

Programs

  • Magma
    A046736:= func< n | n eq 2 select 1 else (&+[Binomial(n+k-2,k)*Binomial(n-k-3, k-1)/(n-1): k in [0..Floor(n/2)-1]]) >;
    [A046736(n): n in [2..40]]; // G. C. Greubel, Jul 31 2024
    
  • Maple
    a := n->1/(n-1)*sum(binomial(n+k-2,k)*binomial(n-k-3,k-1),k=0..floor(n/2-1)); seq(a(i),i=2..30);
  • Mathematica
    (* Programs from Jean-François Alcover, Apr 14 2017: Start *)
    (* First program *)
    a[2]=1; a[n_] := Sum[Binomial[n+k-2, k]*Binomial[n-k-3, k-1], {k, 0, Floor[n/2]-1}]/(n-1);
    (* 2nd program: *)
    x*InverseSeries[Series[(y-y^2-y^3)/(1-y), {y, 0, 29}], x]
    (* 3rd program: *)
    a[2]=1; a[3]=0; a[n_] := HypergeometricPFQ[{2-n/2, 5/2-n/2, n}, {2, 4-n}, -4]; Table[a[n], {n, 2, 30}]
    (* End *)
  • PARI
    a(n)=if(n<2,0,polcoeff(serreverse((x-x^2-x^3)/(1-x)+x*O(x^n)),n-1))
    
  • SageMath
    def A046736(n): return 1 if n==2 else sum(binomial(n+k-2,k)*binomial(n-k-3, k-1)//(n-1) for k in range(n//2))
    [A046736(n) for n in range(2,41)] # G. C. Greubel, Jul 31 2024

Formula

G.f.: A(x) = Sum_{n>0} a(n)*x^(n-1) satisfies A(x) - A(x)^2 - A(x)^3 = x*(1 - A(x)).
a(n) = A052524(n-1)/(n-1)!, for n > 0.
Let g = (1-x)/(1-x-x^2) then a(m) = coeff. of x^(m-2) in g^(m-1)/(m-1).
D-finite with recurrence: 5*(n-1)*n*(37*n-95)*a(n) = 4*(n-1)*(74*n^2 -301*n +300)*a(n-1) + 8*(2*n-5)*(74*n^2 -301*n +297)*a(n-2) - 2*(n-3)*(2*n-7)*(37*n-58)*a(n-3). - Vaclav Kotesovec, Aug 10 2013
a(n) = A143363(n-3) + Sum_{k=0..n-6} ( A143363(k+1)*a(n-k-2) ). - Muhammed Sefa Saydam, Feb 14 2025 and Aug 05 2025

A086616 Partial sums of the large Schroeder numbers (A006318).

Original entry on oeis.org

1, 3, 9, 31, 121, 515, 2321, 10879, 52465, 258563, 1296281, 6589727, 33887465, 175966211, 921353249, 4858956287, 25786112993, 137604139011, 737922992937, 3974647310111, 21493266631001, 116642921832963, 635074797251889, 3467998148181631, 18989465797056721, 104239408386028035
Offset: 0

Views

Author

Paul D. Hanna, Jul 24 2003

Keywords

Comments

Row sums of triangle A086614. - Paul D. Hanna, Jul 24 2003
Hankel transform is A136577(n+1). - Paul Barry, Jun 03 2009

Examples

			a(1) = 2 + 1 = 3;
a(2) = 3 + 4 + 2 = 9;
a(3) = 4 + 10 + 12 + 5 = 31;
a(4) = 5 + 20 + 42 + 40 + 14 = 121.
		

Crossrefs

Cf. A086614 (triangle), A086615 (antidiagonal sums).
Cf. A006318.

Programs

  • Mathematica
    Table[SeriesCoefficient[(1-x-Sqrt[1-6*x+x^2])/(2*x*(1-x)),{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 14 2012 *)
  • PARI
    x='x+O('x^66); Vec((1-x-sqrt(1-6*x+x^2))/(2*x*(1-x))) \\ Joerg Arndt, May 10 2013
  • Sage
    # Generalized algorithm of L. Seidel
    def A086616_list(n) :
        D = [0]*(n+2); D[1] = 1
        b = True; h = 2; R = []
        for i in range(2*n) :
            if b :
                for k in range(h,0,-1) : D[k] += D[k-1]
            else :
                for k in range(1,h, 1) : D[k] += D[k-1]
                R.append(D[h-1]); h += 1;
            b = not b
        return R
    A086616_list(23) # Peter Luschny, Jun 02 2012
    

Formula

G.f.: A(x) = 1/(1 - x)^2 + x*A(x)^2.
a(1) = 1 and a(n) = n + Sum_{i=1..n-1} a(i)*a(n-i) for n >= 2. - Benoit Cloitre, Mar 16 2004
G.f.: (1 - x - sqrt(1 - 6*x + x^2))/(2*x*(1 - x)). Cf. A001003. - Ralf Stephan, Mar 23 2004
a(n) = Sum_{k=0..n} C(n+k+1, 2*k+1) * A000108(k). - Paul Barry, Jun 03 2009
Recurrence: (n+1)*a(n) = (7*n-2)*a(n-1) - (7*n-5)*a(n-2) + (n-2)*a(n-3). - Vaclav Kotesovec, Oct 14 2012
a(n) ~ sqrt(24 + 17*sqrt(2))*(3 + 2*sqrt(2))^n/(4*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 14 2012
A(x) = 1/(1 - x)^2 * 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

Name changed using a comment of Emeric Deutsch from Dec 20 2004. - Peter Luschny, Jun 03 2012

A118376 Number of all trees of weight n, where nodes have positive integer weights and the sum of the weights of the children of a node is equal to the weight of the node.

Original entry on oeis.org

1, 2, 6, 24, 112, 568, 3032, 16768, 95200, 551616, 3248704, 19389824, 117021824, 712934784, 4378663296, 27081760768, 168530142720, 1054464293888, 6629484729344, 41860283723776, 265346078982144, 1687918305128448, 10771600724946944, 68941213290561536
Offset: 1

Views

Author

Jeremy Johnson (jjohnson(AT)cs.drexel.edu), May 15 2006

Keywords

Comments

The number of trees with leaf nodes equal to 1 is counted by the sequence A001003 of super-Catalan numbers. The number of binary trees is counted by the sequence A007317 and the number of binary trees with leaf nodes equal to 1 is counted by the sequence A000108 of Catalan numbers.
Also the number of series-reduced enriched plane trees of weight n. A series-reduced enriched plane tree of weight n is either the number n itself or a finite sequence of at least two series-reduced enriched plane trees, one of each part of an integer composition of n. For example, the a(3) = 6 trees are: 3, (21), (12), (111), ((11)1), (1(11)). - Gus Wiseman, Sep 11 2018
Conjectured to be the number of permutations of length n avoiding the partially ordered pattern (POP) {1>2, 1>3, 3>4, 3>5} of length 5. That is, conjectured to be the number of length n permutations having no subsequences of length 5 in which the first element is the largest, and the third element is larger than the fourth and fifth elements. - Sergey Kitaev, Dec 13 2020
This conjecture has been proven. It can be restated as the number of size n permutations avoiding 51423, 51432, 52413, 52431, 53412, 53421, 54312, 54321. There are twelve sets of permutations avoiding eight size five permutations that are known to match this sequence. A further four are conjectured to match this sequence. - Christian Bean, Jul 24 2024

Examples

			T(3) = 6 because there are six trees
  3    3      3     3     3       3
      2 1    2 1   1 2   1 2    1 1 1
            1 1           1 1
From _Gus Wiseman_, Sep 11 2018: (Start)
The a(4) = 24 series-reduced enriched plane trees:
  4,
  (31), (13), (22), (211), (121), (112), (1111),
  ((21)1), ((12)1), (1(21)), (1(12)), (2(11)), ((11)2),
  ((111)1), (1(111)), ((11)(11)), ((11)11), (1(11)1), (11(11)),
  (((11)1)1), ((1(11))1), (1((11)1)), (1(1(11))).
(End)
		

Crossrefs

Programs

  • Maple
    T := proc(n) option remember; local C, s, p, tp, k, i; if n = 1 then return 1; else s := 1; for k from 2 to n do C := combinat[composition](n,k); for p in C do tp := map(T,p); s := s + mul(tp[i],i=1..nops(tp)); end do; end do; end if; return s; end;
  • Mathematica
    Rest[CoefficientList[Series[(Sqrt[1-8*x+8*x^2]-1)/(4*x-4), {x, 0, 20}], x]] (* Vaclav Kotesovec, Feb 03 2014 *)
    a[n_] := 1+Sum[Binomial[n-1, k-1]*Hypergeometric2F1[2-k, k+1, 2, -1], {k, 2, n}]; Table[a[n], {n, 1, 20}] (* Jean-François Alcover, Apr 03 2015, after Vladimir Kruchinin *)
    urp[n_]:=Prepend[Join@@Table[Tuples[urp/@ptn],{ptn,Join@@Permutations/@Select[IntegerPartitions[n],Length[#]>1&]}],n];
    Table[Length[urp[n]],{n,7}] (* Gus Wiseman, Sep 11 2018 *)
  • Maxima
    a(n):=sum((-1)^j*2^(n-j-1)*binomial(n,j)*binomial(2*n-2*j-2,n-2*j-1),j,0,(n-1)/2)/n; /* Vladimir Kruchinin, Sep 29 2020 */
  • PARI
    x='x+O('x^25); Vec((sqrt(1-8*x+8*x^2) - 1)/(4*x-4)) \\ G. C. Greubel, Feb 08 2017
    

Formula

Recurrence: T(1) = 1; For n > 1, T(n) = 1 + Sum_{n=n1+...+nt} T(n1)*...*T(nt).
G.f.: (-1+(1-8*z+8*z^2)^(1/2))/(-4+4*z).
From Vladimir Kruchinin, Sep 03 2010: (Start)
O.g.f.: A(x) = A001003(x/(1-x)).
a(n) = Sum_{k=1..n} binomial(n-1,k-1)*A001003(k), n>0. (End)
D-finite with recurrence: n*a(n) + 3*(-3*n+4)*a(n-1) + 4*(4*n-9)*a(n-2) + 8*(-n+3)*a(n-3) = 0. - R. J. Mathar, Sep 27 2013
a(n) ~ sqrt(sqrt(2)-1) * 2^(n-1/2) * (2+sqrt(2))^(n-1) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Feb 03 2014
From Peter Bala, Jun 17 2015: (Start)
With offset 0, binomial transform of A001003.
O.g.f. A(x) = series reversion of x*(2*x - 1)/(2*x^2 - 1); 2*(1 - x)*A^2(x) - A(x) + x = 0.
A(x) satisfies the differential equation (x - 9*x^2 + 16*x^3 - 8*x^4)*A'(x) + x*(3 - 4*x)*A(x) + x*(2*x - 1) = 0. Extracting coefficients gives Mathar's recurrence above. (End)
a(n) = Sum_{j=0..(n-1)/2} (-1)^j*2^(n-j-1)*C(n,j)*C(2*n-2*j-2,n-2*j-1)/n. - Vladimir Kruchinin, Sep 29 2020

A120803 Number of series-reduced balanced trees with n leaves.

Original entry on oeis.org

1, 1, 1, 2, 2, 4, 4, 8, 9, 16, 20, 37, 47, 80, 111, 183, 256, 413, 591, 940, 1373, 2159, 3214, 5067, 7649, 12054, 18488, 29203, 45237, 71566, 111658, 176710, 276870, 437820, 687354, 1085577, 1705080, 2688285, 4221333, 6644088, 10425748
Offset: 1

Views

Author

Keywords

Comments

In other words, rooted trees with all leaves at the same level and no node having exactly one child; the order of children is not significant.

Examples

			From _Gus Wiseman_, Oct 07 2018: (Start)
The a(10) = 16 series-reduced balanced rooted trees:
  (oooooooooo)
  ((ooooo)(ooooo))
  ((oooo)(oooooo))
  ((ooo)(ooooooo))
  ((oo)(oooooooo))
  ((ooo)(ooo)(oooo))
  ((oo)(oooo)(oooo))
  ((oo)(ooo)(ooooo))
  ((oo)(oo)(oooooo))
  ((oo)(oo)(ooo)(ooo))
  ((oo)(oo)(oo)(oooo))
  ((oo)(oo)(oo)(oo)(oo))
  (((oo)(ooo))((oo)(ooo)))
  (((oo)(oo))((ooo)(ooo)))
  (((oo)(oo))((oo)(oooo)))
  (((oo)(oo))((oo)(oo)(oo)))
(End)
		

Crossrefs

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    seq(n)={my(u=vector(n), v=vector(n)); u[1]=1; while(u, v+=u; u=EulerT(u)-u); v} \\ Andrew Howroyd, Oct 26 2018

Formula

Let s_0(n) = 1 if n = 1, 0 otherwise; s_{k+1}(n) = EULER(s_k)(n) - s_k(n), where EULER is the Euler transform. Then a_n = sum_k s_k(n). (s_k(n) is the number of such trees of height k.) Note that s_k(n) = 0 for n < 2^k.

A007853 Number of maximal antichains in rooted plane trees on n nodes.

Original entry on oeis.org

1, 2, 5, 15, 50, 178, 663, 2553, 10086, 40669, 166752, 693331, 2917088, 12398545, 53164201, 229729439, 999460624, 4374546305, 19250233408, 85120272755, 378021050306, 1685406494673, 7541226435054, 33852474532769, 152415463629568, 688099122024944
Offset: 1

Views

Author

Keywords

Comments

Also the number of initial subtrees (emanating from the root) of rooted plane trees on n vertices, where we require that an initial subtree contains either all or none of the branchings under any given node. The leaves of such a subtree comprise the roots of a corresponding antichain cover. Also, in the (non-commutative) multicategory of free pure multifunctions with one atom, a(n) is the number of composable pairs whose composite has n positions. - Gus Wiseman, Aug 13 2018
The g.f. is denoted by y_2 in Bacher 2004 Proposition 7.5 on page 20. - Michael Somos, Nov 07 2019

Examples

			G.f. = x + 2*x^2 + 5*x^3 + 15*x^4 + 50*x^5 + 178*x^6 + 663*x^7 + 2553*x^8 + ... - _Michael Somos_, Nov 07 2019
		

Crossrefs

Programs

  • Mathematica
    ie[t_]:=If[Length[t]==0,1,1+Product[ie[b],{b,t}]];
    allplane[n_]:=If[n==1,{{}},Join@@Function[c,Tuples[allplane/@c]]/@Join@@Permutations/@IntegerPartitions[n-1]];
    Table[Sum[ie[t],{t,allplane[n]}],{n,9}] (* Gus Wiseman, Aug 13 2018 *)
  • Maxima
    a(n):=1/(n+1)*binomial(2*n,n)+sum((k+2)/(n+1)*binomial(2*n-k-1,n-k-1)*(sum(((binomial(2*i,i))*(binomial(k+i,3*i)))/(i+1),i,0,floor(k/2))),k,0,n-1); /* Vladimir Kruchinin, Apr 05 2019 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = sqrt(1 - 4*x + x * O(x^n)); polcoeff( (3 - 2*x - A - sqrt(2 - 16*x + 4*x^2 + (2 + 4*x) * A)) / 4, n))}; /* Michael Somos, Nov 07 2019 */

Formula

G.f.: (1/4) * (3 - 2*x - sqrt(1-4*x) - sqrt(2) * sqrt((1+2*x) * sqrt(1-4*x) + 1 - 8*x + 2*x^2)) [from Klazar]. - Sean A. Irvine, Feb 06 2018
a(n) = (1/(n+1))*C(2*n,n) + Sum_{k=0..n-1} ((k+2)/(n+1))*C(2*n-k-1,n-k-1)*Sum_{i=0..floor(k/2)} C(2*i,i)*C(k+i,3*i)/(i+1). - Vladimir Kruchinin, Apr 05 2019
Given the g.f. A(x) and the g.f. of A213705 B(x), then -x = A(-B(x)). - Michael Somos, Nov 07 2019

Extensions

More terms from Sean A. Irvine, Feb 06 2018

A317658 Number of positions in the n-th free pure symmetric multifunction (with empty expressions allowed) with one atom.

Original entry on oeis.org

1, 2, 3, 3, 4, 4, 5, 4, 4, 5, 6, 5, 5, 6, 7, 4, 6, 6, 7, 8, 5, 7, 7, 8, 5, 9, 5, 6, 8, 8, 9, 5, 6, 10, 6, 5, 7, 9, 9, 10, 6, 7, 11, 7, 6, 8, 10, 10, 6, 11, 7, 8, 12, 8, 7, 9, 11, 11, 7, 12, 8, 9, 13, 5, 9, 8, 10, 12, 12, 8, 13, 9, 10, 14, 6, 10, 9, 11, 13, 13
Offset: 1

Views

Author

Gus Wiseman, Aug 03 2018

Keywords

Comments

Given a positive integer n > 1 we construct a unique free pure symmetric multifunction e(n) by expressing n as a power of a number that is not a perfect power to a product of prime numbers: n = rad(x)^(prime(y_1) * ... * prime(y_k)) where rad = A007916. Then e(n) = e(x)[e(y_1), ..., e(y_k)].
Also the number of positions in the orderless Mathematica expression with e-number n.

Examples

			The first twenty Mathematica expressions:
   1: o
   2: o[]
   3: o[][]
   4: o[o]
   5: o[][][]
   6: o[o][]
   7: o[][][][]
   8: o[o[]]
   9: o[][o]
  10: o[o][][]
  11: o[][][][][]
  12: o[o[]][]
  13: o[][o][]
  14: o[o][][][]
  15: o[][][][][][]
  16: o[o,o]
  17: o[o[]][][]
  18: o[][o][][]
  19: o[o][][][][]
  20: o[][][][][][][]
		

Crossrefs

First differs from A277615 at a(128) = 5, A277615(128) = 6.

Programs

  • Mathematica
    nn=100;
    radQ[n_]:=If[n===1,False,GCD@@FactorInteger[n][[All,2]]===1];
    rad[n_]:=rad[n]=If[n===0,1,NestWhile[#+1&,rad[n-1]+1,Not[radQ[#]]&]];
    Clear[radPi];Set@@@Array[radPi[rad[#]]==#&,nn];
    exp[n_]:=If[n===1,x,With[{g=GCD@@FactorInteger[n][[All,2]]},Apply[exp[radPi[Power[n,1/g]]],exp/@Flatten[Cases[FactorInteger[g],{p_?PrimeQ,k_}:>ConstantArray[PrimePi[p],k]]]]]];
    Table[exp[n],{n,1,nn}]

Formula

a(rad(x)^(prime(y_1) * ... * prime(y_k))) = a(x) + a(y_1) + ... + a(y_k).
e(2^(2^n)) = o[o,...,o].
e(2^prime(2^prime(2^...))) = o[o[...o[o]]].
e(rad(rad(rad(...)^2)^2)^2) = o[o][o]...[o].

A010683 Let S(x,y) = number of lattice paths from (0,0) to (x,y) that use the step set { (0,1), (1,0), (2,0), (3,0), ...} and never pass below y = x. Sequence gives S(n-1,n) = number of 'Schröder' trees with n+1 leaves and root of degree 2.

Original entry on oeis.org

1, 2, 7, 28, 121, 550, 2591, 12536, 61921, 310954, 1582791, 8147796, 42344121, 221866446, 1170747519, 6216189936, 33186295681, 178034219986, 959260792775, 5188835909516, 28167068630713, 153395382655222
Offset: 0

Views

Author

Robert Sulanke (sulanke(AT)diamond.idbsu.edu), N. J. A. Sloane

Keywords

Comments

a(n) is the number of compound propositions "on the negative side" that can be made from n simple propositions.
Convolution of A001003 (the little Schröder numbers) with itself. - Emeric Deutsch, Dec 27 2003
Number of dissections of a convex polygon with n+3 sides that have a triangle over a fixed side (the base) of the polygon. - Emeric Deutsch, Dec 27 2003
a(n-1) = number of royal paths from (0,0) to (n,n), A006318, with exactly one diagonal step on the line y=x. - David Callan, Mar 14 2004
Number of short bushes (i.e., ordered trees with no vertices of outdegree 1) with n+2 leaves and having root of degree 2. Example: a(2)=7 because, in addition to the five binary trees with 6 edges (they do have 4 leaves) we have (i) two edges rb, rc hanging from the root r with three edges hanging from vertex b and (ii) two edges rb, rc hanging from the root r with three edges hanging from vertex c. - Emeric Deutsch, Mar 16 2004
The a(n) equal the Fi2 sums, see A180662, of Schröder triangle A033877. - Johannes W. Meijer, Mar 26 2012
Row sums of A144944 and of A186826. - Reinhard Zumkeller, May 11 2013

Crossrefs

Second right-hand column of triangle A011117.
A177010 has a closely-related g.f..

Programs

  • Haskell
    a010683 = sum . a144944_row  -- Reinhard Zumkeller, May 11 2013
    
  • Magma
    [n le 2 select n else (6*(2*(n-1)^2-1)*Self(n-1) - (n-3)*(2*n-1)*Self(n-2))/((n+1)*(2*n-3)): n in [1..30]]; // G. C. Greubel, Mar 11 2023
  • Maple
    a := proc(n) local k: if n=0 then 1 else (2/n)*add(binomial(n, k)* binomial(n+k+1, k-1), k=1..n) fi: end:
    seq(a(n), n=0..21); # Johannes W. Meijer, Mar 26 2012, revised Mar 31 2015
  • Mathematica
    f[ x_, y_ ]:= f[ x, y ]= Module[ {return}, If[x==0, return =1, If[y==x-1, return =0, return= f[x,y-1] + Sum[f[k, y], {k,0,x-1} ]]]; return];
    (* Do[Print[Table[f[ k, j ], {k, 0, j}]], {j, 10, 0, -1}] *)
    Table[f[x, x+1], {x,0,21}]
    (* Second program: *)
    a[n_] := 2*Hypergeometric2F1[1-n, n+3, 2, -1]; a[0]=1;
    Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Dec 09 2014, after Wolfdieter Lang *)
  • PARI
    x='x+O('x^100); Vec(((1-x)^2-(1+x)*sqrt(1-6*x+x^2))/(8*x^2)) \\ Altug Alkan, Dec 19 2015
    
  • Sage
    a = lambda n: (n+1)*hypergeometric([1-n, -n], [3], 2)
    [simplify(a(n)) for n in range(22)] # Peter Luschny, Nov 19 2014
    

Formula

G.f.: ((1-t)^2-(1+t)*sqrt(1-6*t+t^2))/(8*t^2) = A(t)^2, with o.g.f. A(t) of A001003.
From Wolfdieter Lang, Sep 12 2005: (Start)
a(n) = (2/n)*Sum_{k=1..n} binomial(n, k)*binomial(n+k+1, k-1).
a(n) = 2*hypergeometric2F1([1-n, n+3], [2], -1), n>=1. a(0)=1. (End)
a(n) = ((2*n+1)*LegendreP(n+1,3) - (2*n+3)*LegendreP(n,3)) / (4*n*(n+2)) for n>0. - Mark van Hoeij, Jul 02 2010
From Gary W. Adamson, Jul 08 2011: (Start)
Let M = the production matrix:
1, 2, 0, 0, 0, 0, ...
1, 2, 1, 0, 0, 0, ...
1, 2, 1, 2, 0, 0, ...
1, 2, 1, 2, 1, 0, ...
1, 2, 1, 2, 1, 2, ...
...
a(n) is the upper entry in the vector (M(T))^n * [1,0,0,0,...]; where T is the transpose operation. (End)
D-finite with recurrence: (n+2)*(2*n-1)*a(n) = 6*(2*n^2-1)*a(n-1) - (n-2)*(2*n+1)*a(n-2). - Vaclav Kotesovec, Oct 07 2012
a(n) ~ sqrt(48+34*sqrt(2))*(3+2*sqrt(2))^n/(4*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 07 2012
Recurrence (an alternative): (n+2)*a(n) = (4-n)*a(n-4) + 2*(2*n-5)*a(n-3) + 10*(n-1)*a(n-2) + 2*(2*n+1)*a(n-1), n >= 4. - Fung Lam, Feb 18 2014
a(n) = (n+1)*hypergeometric2F1([1-n, -n], [3], 2). - Peter Luschny, Nov 19 2014
a(n) = (A001003(n) + A001003(n+1))/2 = sum(A001003(k) * A001003(n-k), k=0..n). - Johannes W. Meijer, Apr 29 2015

Extensions

Minor edits by Johannes W. Meijer, Mar 26 2012
Previous Showing 51-60 of 248 results. Next