cp's OEIS Frontend

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

Showing 1-9 of 9 results.

A190737 Diagonal sums of the Riordan matrix A104259.

Original entry on oeis.org

1, 2, 6, 19, 66, 244, 946, 3801, 15697, 66234, 284339, 1237983, 5453611, 24263355, 108865901, 492051006, 2238220336, 10238568080, 47070014643, 217363784060, 1007794226777, 4689545704246, 21893712581740, 102520882301832, 481393173378979
Offset: 0

Views

Author

Emanuele Munarini, May 18 2011

Keywords

Crossrefs

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); Coefficients(R!((1-x-2*x^2-Sqrt(1-6*x+5*x^2))/(2*x*(1-2*x+x^2+x^3)))); // G. C. Greubel, Apr 23 2018
  • Mathematica
    CoefficientList[Series[(1-x-2x^2-Sqrt[1-6x+5x^2])/(2x(1-2x+x^2+x^3)),{x,0,24}],x]
  • PARI
    x='x+O('x^30); Vec((1-x-2*x^2-sqrt(1-6*x+5*x^2))/(2*x*(1-2*x +x^2 +x^3))) \\ G. C. Greubel, Apr 23 2018
    

Formula

G.f.: (1-x-2*x^2-sqrt(1-6*x+5*x^2))/(2*x*(1-2*x+x^2+x^3)).
Recurrence: 0 = (n^2+17*n+72)*a(n+8) - (11*n^2+169*n+648)*a(n+7) + 3*(17*n^2+227*n+758)*a(n+6) - 2*(73*n^2+812*n+2259)*a(n+5)
+ 6*(41*n^2+375*n+846)*a(n+4) - (217*n^2+1559*n+2634)*a(n+3) + 3*(17*n^2+29*n-136)*a(n+2) + 10*(5*n^2+43*n+84)*a(n+1) - 25*(n^2+5*n+6)*a(n).
D-finite with recurrence: (n+1)*a(n) +(-8*n+1)*a(n-1) +3*(6*n-5)*a(n-2) +3*(-5*n+8)*a(n-3) +(-n-7)*a(n-4) +5*(n-2)*a(n-5)=0. - R. J. Mathar, Feb 24 2020

A190738 Central coefficients of the Riordan matrix A104259.

Original entry on oeis.org

1, 4, 27, 212, 1785, 15630, 140287, 1280592, 11833389, 110360150, 1036670272, 9794291556, 92972640761, 886023463122, 8471878678545, 81236546627920, 780898417097733, 7522708492892214, 72607180401922894, 701969331508141900, 6796919869909393140
Offset: 0

Views

Author

Emanuele Munarini, May 18 2011

Keywords

Crossrefs

Cf. A104259.

Programs

  • Mathematica
    Table[Sum[Binomial[2n,n+k]Binomial[n+2k,k](n+1)/(n+k+1),{k,0,n}],{n,0,20}]
  • Maxima
    makelist(sum(binomial(2*n,n+k)*binomial(n+2*k,k)*(n+1)/(n+k+1),k,0,n),n,0,20);

Formula

a(n) = T(2*n,n), where T(n,k) = A104259(n,k).
a(n) = sum(k=0..n, binomial(2*n,n+k)*binomial(n+2*k,k)*(n+1)/(n+k+1)).
Recurrence: 2*(n-1)*n*(2*n + 1)*(107*n^3 - 345*n^2 + 253*n + 12)*a(n) = (n-1)*(5029*n^5 - 16215*n^4 + 8284*n^3 + 13359*n^2 - 11675*n + 1974)*a(n-1) - 4*(2*n - 3)*(1498*n^5 - 4830*n^4 + 2899*n^3 + 3165*n^2 - 3254*n + 630)*a(n-2) + 100*(n-1)*(2*n - 5)*(2*n - 3)*(107*n^3 - 24*n^2 - 116*n + 27)*a(n-3). - Vaclav Kotesovec, Jul 05 2021
a(n) ~ sqrt(c) * d^n / sqrt(Pi*n), where d = 9.945658804810730213397409025621... is the real root of the equation -400 + 112*d - 47*d^2 + 4*d^3 = 0 and c = 0.3447849735035503206155951176700724872157... is the real root of the equation -125 - 173*c + 963*c^2 + 1712*c^3 = 0. - Vaclav Kotesovec, Jun 05 2022

A002212 Number of restricted hexagonal polyominoes with n cells.

Original entry on oeis.org

1, 1, 3, 10, 36, 137, 543, 2219, 9285, 39587, 171369, 751236, 3328218, 14878455, 67030785, 304036170, 1387247580, 6363044315, 29323149825, 135700543190, 630375241380, 2938391049395, 13739779184085, 64430797069375, 302934667061301, 1427763630578197
Offset: 0

Views

Author

N. J. A. Sloane, Ronald C. Read

Keywords

Comments

Number of Schroeder paths (i.e., consisting of steps U=(1,1), D=(1,-1), H=(2,0) and never going below the x-axis) from (0,0) to (2n,0) with no peaks at odd level. Example: a(2)=3 because we have UUDD, UHD and HH. - Emeric Deutsch, Dec 06 2003
Number of 3-Motzkin paths of length n-1 (i.e., lattice paths from (0,0) to (n-1,0) that do not go below the line y=0 and consist of steps U=(1,1), D=(1,-1) and three types of steps H=(1,0)). Example: a(4)=36 because we have 27 HHH paths, 3 HUD paths, 3 UHD paths and 3 UDH paths. - Emeric Deutsch, Jan 22 2004
Number of rooted, planar trees having edges weighted by strictly positive integers (multi-trees) with weight-sum n. - Roland Bacher, Feb 28 2005
Number of skew Dyck paths of semilength n. A skew Dyck path is a path in the first quadrant which begins at the origin, ends on the x-axis, consists of steps U=(1,1)(up), D=(1,-1)(down) and L=(-1,-1)(left) so that up and left steps do not overlap. The length of the path is defined to be the number of its steps. - Emeric Deutsch, May 10 2007
Equivalently, number of self-avoiding paths of semilength n in the first quadrant beginning at the origin, staying weakly above the diagonal, ending on the diagonal, and consisting of steps r=(+1,0) (right), U=(0,+1) (up), and D=(0,-1) (down). Self-avoidance implies that factors UD and DU and steps D reaching the diagonal before the end are forbidden. The a(3) = 10 such paths are UrUrUr, UrUUrD, UrUUrr, UUrrUr, UUrUrD, UUrUrr, UUUDrD, UUUrDD, UUUrrD, and UUUrrr. - Joerg Arndt, Jan 15 2024
Hankel transform of [1,3,10,36,137,543,...] is A000012 = [1,1,1,1,...]. - Philippe Deléham, Oct 24 2007
From Gary W. Adamson, May 17 2009: (Start)
Convolved with A026375, (1, 3, 11, 45, 195, ...) = A026378: (1, 4, 17, 75, ...)
(1, 3, 10, 36, 137, ...) convolved with A026375 = A026376: (1, 6, 30, 144, ...).
Starting (1, 3, 10, 36, ...) = INVERT transform of A007317: (1, 2, 5, 15, 51, ...). (End)
Binomial transform of A032357. - Philippe Deléham, Sep 17 2009
a(n) = number of rooted trees with n vertices in which each vertex has at most 2 children and in case a vertex has exactly one child, it is labeled left, middle or right. These are the hex trees of the Deutsch, Munarini, Rinaldi link. This interpretation yields the second MATHEMATICA recurrence below. - David Callan, Oct 14 2012
The left shift (1,3,10,36,...) of this sequence is the binomial transform of the left-shifted Catalan numbers (1,2,5,14,...). Example: 36 =1*14 + 3*5 + 3*2 + 1*1. - David Callan, Feb 01 2014
Number of Schroeder paths from (0,0) to (2n,0) with no level steps H=(2,0) at even level. Example: a(2)=3 because we have UUDD, UHD and UDUD. - José Luis Ramírez Ramírez, Apr 27 2015
This is the Riordan transform with the Riordan matrix A097805 (of the associated type) of the Catalan sequence A000108. See a Feb 17 2017 comment in A097805. - Wolfdieter Lang, Feb 17 2017
a(n) is the number of parking functions of size n avoiding the patterns 132 and 231. - Lara Pudwell, Apr 10 2023

Examples

			G.f. = 1 + x + 3*x^2 + 10*x^3 + 36*x^4 + 137*x^5 + 543*x^6 + 2219*x^7 + 9285*x^8 + ...
		

References

  • J. Brunvoll, B. N. Cyvin, and S. J. Cyvin, Studies of some chemically relevant polygonal systems: mono-q-polyhexes, ACH Models in Chem., 133 (3) (1996), 277-298, Eq 14.
  • S. J. Cyvin, J. Brunvoll, G. Xiaofeng, and Z. Fuji, Number of perifusenes with one internal vertex, Rev. Roumaine Chem., 38(1) (1993), 65-78.
  • 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

First differences of A007317.
Row sums of triangle A104259.

Programs

  • Magma
    I:= [1,3]; [1] cat [n le 2 select I[n]  else ((6*n-3)*Self(n-1)-5*(n-2)*Self(n-2)) div (n+1): n in [1..30]]; // Vincenzo Librandi, Jun 15 2015
  • Maple
    t1 := series(1+ (1-3*x-(1-x)^(1/2)*(1-5*x)^(1/2))/(2*x), x, 50):
    A002212_list := len -> seq(coeff(t1,x,n),n=0..len): A002212_list(40);
    a[0] := 1: a[1] := 1: for n from 2 to 50 do a[n] := (3*(2*n-1)*a[n-1]-5*(n-2)*a[n-2])/(n+1) od: print(convert(a,list)); # Zerinvary Lajos, Jan 01 2007
    a := n -> `if`(n=0,1,simplify(GegenbauerC(n-1, -n, -3/2)/n)):
    seq(a(n), n=0..23); # Peter Luschny, May 09 2016
  • Mathematica
    InverseSeries[Series[(y)/(1+3*y+y^2), {y, 0, 24}], x] (* then A(x)=1+y(x) *) (* Len Smiley, Apr 14 2000 *)
    (* faster *)
    a[0]=1;a[1]=1;
    a[n_]/;n>=2 := a[n] = a[n-1] +  Sum[a[i]a[n-1-i],{i,0,n-1}];
    Table[a[n],{n,0,14}] (* See COMMENTS above, [David Callan, Oct 14 2012] *)
    (* fastest *)
    s[0]=s[1]=1;
    s[n_]/;n>=2 := s[n] = (3(2n-1)s[n-1]-5(n-2)s[n-2])/(n+1);
    Table[s[n],{n,0,14 }] (* See Deutsch, Munarini, Rinaldi link, [David Callan, Oct 14 2012] *)
    (* 2nd fastest *)
    a[n_] := Hypergeometric2F1[3/2, 1-n, 3, -4]; a[0]=1; Table[a[n], {n, 0, 14}]  (* Jean-François Alcover, May 16 2013 *)
    CoefficientList[Series[(1 - x - Sqrt[1 - 6x + 5x^2])/(2x), {x, 0, 20}], x] (* Nikolaos Pantelidis, Jan 30 2023 *)
  • Maxima
    makelist(sum(binomial(n,k)*binomial(n-k,k)*3^(n-2*k)/(k+1),k,0,n/2),n,0,24); /* for a(n+1) */ /* Emanuele Munarini, May 18 2011 */
    
  • PARI
    {a(n) = polcoeff( (1 - x - sqrt(1 - 6*x + 5*x^2 + x^2 * O(x^n))) / 2, n+1)};
    
  • PARI
    {a(n) = if( n<1, n==0, polcoeff( serreverse( x / (1 + 3*x + x^2) + x * O(x^n)), n))}; /* Michael Somos */
    
  • PARI
    my(N=66,x='x+O('x^N)); Vec((1 - x - sqrt(1-6*x+5*x^2))/(2*x)) \\ Joerg Arndt, Jan 13 2024
    
  • Sage
    def A002212():
        x, y, n = 1, 1, 1
        while True:
            yield x
            n += 1
            x, y = y, ((6*n - 3)*y - (5*n - 10)*x) / (n + 1)
    a = A002212()
    [next(a) for i in range(24)]  # Peter Luschny, Oct 12 2013
    

Formula

a(0)=1, for n > 0: a(n) = Sum_{j=0..n-1} Sum_{i=0..j} a(i)*a(j-i). G.f.: A(x) = 1 + x*A(x)^2/(1-x). - Mario Catalani (mario.catalani(AT)unito.it), Jun 19 2003
a(n) = Sum_{i=ceiling((n-1)/2)..n-1} (3^(2i+1-n)*binomial(n, i)*binomial(i, n-i-1))/n. - Emeric Deutsch, Jul 23 2002
a(n) = Sum_{k=1..n} binomial(2k, k)*binomial(n-1, k-1)/(k+1), i.e., binomial transform of the Catalan numbers 1, 2, 5, 14, 42, ... (A000108). a(n) = Sum_{k=0..floor((n-1)/2)} 3^(n-1-2*k)*binomial(2k, k)*binomial(n-1, 2k)/(k+1). - Emeric Deutsch, Aug 05 2002
D-finite with recurrence: a(1)=1, a(n) = (3(2n-1)*a(n-1)-5(n-2)*a(n-2))/(n+1) for n > 1. - Emeric Deutsch, Dec 18 2002
a(n) is asymptotic to c*5^n/n^(3/2) with c=0.63.... - Benoit Cloitre, Jun 23 2003
In closed form, c = (1/2)*sqrt(5/Pi) = 0.63078313050504... - Vaclav Kotesovec, Oct 04 2012
Reversion of Sum_{n>0} a(n)x^n = -Sum_{n>0} A001906(n)(-x)^n.
G.f. A(x) satisfies xA(x)^2 + (1-x)(1-A(x)) = 0.
G.f.: (1 - x - sqrt(1 - 6x + 5x^2))/(2x). For n > 1, a(n) = 3*a(n-1) + Sum_{k=1..n-2} a(k)*a(n-k-1). - John W. Layman, Feb 22 2001
The Hankel transform of this sequence gives A001519 = 1, 2, 5, 13, 34, 89, ... E.g., Det([1, 1, 3, 10, 36; 1, 3, 10, 36, 137; 3, 10, 36, 137, 543; 10, 36, 137, 543, 2219; 36, 137, 543, 2219, 9285 ])= 34. - Philippe Deléham, Jan 25 2004
a(m+n+1) = Sum_{k>=0} A091965(m, k)*A091965(n, k) = A091965(m+n, 0). - Philippe Deléham, Sep 14 2005
a(n+1) = Sum_{k=0..n} 2^(n-k)*M(k)*binomial(n,k), where M(k) = A001006(k) is the k-th Motzkin number (from here it follows that a(n+1) and M(n) have the same parity). - Emeric Deutsch, May 10 2007
a(n+1) = Sum_{k=0..n} A097610(n,k)*3^k. - Philippe Deléham, Oct 02 2007
G.f.: 1/(1-x/(1-x-x/(1-x/(1-x-x/(1-x/(1-x-x/(1-... (continued fraction). - Paul Barry, May 16 2009
G.f.: (1-x)/(1-2x-x^2/(1-3x-x^2/(1-3x-x^2/(1-3x-x^2/(1-3x-x^2/(1-.... (continued fraction). - Paul Barry, Oct 17 2009
G.f.: 1/(1-z/(1-z/(1-z/(...)))) where z=x/(1-x) (continued fraction); more generally g.f. C(x/(1-x)) where C(x) is the g.f. for the Catalan numbers (A000108). - Joerg Arndt, Mar 18 2011
a(n) = -5^(1/2)/(10*(n+1)) * (5*hypergeom([1/2, n], [1], 4/5) -3*hypergeom([1/2, n+1], [1], 4/5)) (for n>0). - Mark van Hoeij, Nov 12 2009
For n >= 1, a(n) = (1/(2*Pi))*Integral_{x=1..5} x^(n-1)*sqrt((x-1)*(5-x)) dx. - Groux Roland, Mar 16 2011
a(n+1) = [x^n](1-x^2)(1+3*x+x^2)^n. - Emanuele Munarini, May 18 2011
From Gary W. Adamson, Jul 21 2011: (Start)
a(n) = upper left term in M^(n-1), M = an infinite square production matrix as follows (with 3,2,2,2,... as the main diagonal):
3, 1, 0, 0, 0, 0, ...
1, 2, 1, 0, 0, 0, ...
1, 1, 2, 1, 0, 0, ...
1, 1, 1, 2, 1, 0, ...
1, 1, 1, 1, 2, 0, ...
...
Alternatively, let M = the previous matrix but change the 3 to a 2. Then a(n) = sum of top row terms of M^(n-1). (End)
a(n) = hypergeometric([1-n,3/2],[3],-4), for n>0. - Peter Luschny, Aug 15 2012
a(n) = GegenbauerC(n-1, -n, -3/2)/n for n >= 1. - Peter Luschny, May 09 2016
E.g.f.: 1 + Integral (exp(3*x) * BesselI(1,2*x) / x) dx. - Ilya Gutkovskiy, Jun 01 2020
G.f.: 1 + x/G(0) with G(k) = (1 - 3*x - x^2/G(k+1)) (continued fraction). - Nikolaos Pantelidis, Dec 12 2022
From Peter Bala, Feb 03 2024: (Start)
G.f.: 1 + x/(1 - x) * c(x/(1 - x))^2 = 1 + x/(1 - 5*x) * c(-x/(1 - 5*x))^2, where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108.
a(n+1) = Sum_{k = 0..n} binomial(n, k)*Catalan(k+1).
a(n+1) = hypergeom([-n, 3/2], [3], -4).
a(n+1) = 5^n * Sum_{k = 0..n} (-5)^(-k)*binomial(n, k)*Catalan(k+1).
a(n+1) = 5^n * hypergeom([-n, 3/2], [3], 4/5). (End)

A007317 Binomial transform of Catalan numbers.

Original entry on oeis.org

1, 2, 5, 15, 51, 188, 731, 2950, 12235, 51822, 223191, 974427, 4302645, 19181100, 86211885, 390248055, 1777495635, 8140539950, 37463689775, 173164232965, 803539474345, 3741930523740, 17481709707825, 81912506777200, 384847173838501, 1812610804416698
Offset: 1

Views

Author

Keywords

Comments

Partial sums of A002212 (the restricted hexagonal polyominoes with n cells). Number of Schroeder paths (i.e., consisting of steps U=(1,1),D=(1,-1),H=(2,0) and never going below the x-axis) from (0,0) to (2n-2,0), with no peaks at even level. Example: a(3)=5 because among the six Schroeder paths from (0,0) to (4,0) only UUDD has a peak at an even level. - Emeric Deutsch, Dec 06 2003
Number of binary trees of weight n where leaves have positive integer weights. Non-commutative Non-associative version of partitions of n. - Michael Somos, May 23 2005
Appears also as the number of Euler trees with total weight n (associated with even switching class of matrices of order 2n). - David Garber, Sep 19 2005
Number of symmetric hex trees with 2n-1 edges; also number of symmetric hex trees with 2n-2 edges. A hex tree is a rooted tree where each vertex has 0, 1, or 2 children and, when only one child is present, it is either a left child, or a median child, or a right child (name due to an obvious bijection with certain tree-like polyhexes; see the Harary-Read reference). A hex tree is symmetric if it is identical with its reflection in a bisector through the root. - Emeric Deutsch, Dec 19 2006
The Hankel transform of [1, 2, 5, 15, 51, 188, ...] is [1, 1, 1, 1, 1, ...], see A000012 ; the Hankel transform of [2, 5, 15, 51, 188, 731, ...] is [2, 5, 13, 34, 89, ...], see A001519. - Philippe Deléham, Dec 19 2006
a(n) = number of 321-avoiding partitions of [n]. A partition is 321-avoiding if the permutation obtained from its canonical form (entries in each block listed in increasing order and blocks listed in increasing order of their first entries) is 321-avoiding. For example, the only partition of [5] that fails to be 321-avoiding is 15/24/3 because the entries 5,4,3 in the permutation 15243 form a 321 pattern. - David Callan, Jul 22 2008
The sequence 1,1,2,5,15,51,188,... has Hankel transform A001519. - Paul Barry, Jan 13 2009
From Gary W. Adamson, May 17 2009: (Start)
Equals INVERT transform of A033321: (1, 1, 2, 6, 21, 79, 311, ...).
Equals INVERTi transform of A002212: (1, 3, 10, 36, 137, ...).
Convolved with A026378, (1, 4, 17, 75, 339, ...) = A026376: (1, 6, 30, 144, ...)
(End)
a(n) is the number of vertices of the composihedron CK(n). The composihedra are a sequence of convex polytopes used to define maps of certain homotopy H-spaces. They are cellular quotients of the multiplihedra and cellular covers of the cubes. - Stefan Forcey (sforcey(AT)gmail.com), Dec 17 2009
a(n) is the number of Motzkin paths of length n-1 in which the (1,0)-steps at level 0 come in 2 colors and those at a higher level come in 3 colors. Example: a(4)=15 because we have 2^3 = 8 paths of shape UHD, 2 paths of shape HUD, 2 paths of shape UDH, and 3 paths of shape UHD; here U=(1,1), H=(1,0), and D=(1,-1). - Emeric Deutsch, May 02 2011
REVERT transform of (1, 2, -3, 5, -8, 13, -21, 34, ... ) where the entries are Fibonacci numbers, A000045. Equivalently, coefficients in the series reversion of x(1-x)/(1+x-x^2). This means that the substitution of the gf (1-x-(1-6x+5x^2)^(1/2))/(2(1-x)) for x in x(1-x)/(1+x-x^2) will simplify to x. - David Callan, Nov 11 2012
The number of plane trees with nodes that have positive integer weights and whose total weight is n. - Brad R. Jones, Jun 12 2014
From Tom Copeland, Nov 02 2014: (Start)
Let P(x) = x/(1+x) with comp. inverse Pinv(x) = x/(1-x) = -P[-x], and C(x)= [1-sqrt(1-4x)]/2, an o.g.f. for the shifted Catalan numbers A000108, with inverse Cinv(x) = x * (1-x).
Fin(x) = P[C(x)] = C(x)/[1 + C(x)] is an o.g.f. for the Fine numbers, A000957 with inverse Fin^(-1)(x) = Cinv[Pinv(x)] = Cinv[-P(-x)].
Mot(x) = C[P(x)] = C[-Pinv(-x)] gives an o.g.f. for shifted A005043, the Motzkin or Riordan numbers with comp. inverse Mot^(-1)(x) = Pinv[Cinv(x)] = (x - x^2) / (1 - x + x^2) (cf. A057078).
BTC(x) = C[Pinv(x)] gives A007317, a binomial transform of the Catalan numbers, with BTC^(-1)(x) = P[Cinv(x)] = (x-x^2) / (1 + x - x^2).
Fib(x) = -Fin[Cinv(Cinv(-x))] = -P[Cinv(-x)] = x + 2 x^2 + 3 x^3 + 5 x^4 + ... = (x+x^2)/[1-x-x^2] is an o.g.f. for the shifted Fibonacci sequence A000045, so the comp. inverse is Fib^(-1)(x) = -C[Pinv(-x)] = -BTC(-x) and Fib(x) = -BTC^(-1)(-x).
Generalizing to P(x,t) = x /(1 + t*x) and Pinv(x,t) = x /(1 - t*x) = -P(-x,t) gives other relations to lattice paths, such as the o.g.f. for A091867, C[P[x,1-t]], and that for A104597, Pinv[Cinv(x),t+1].
(End)
Starting with offset 0, a(n) is also the number of Schröder paths of semilength n avoiding UH (an up step directly followed by a long horizontal step). Example: a(2)=5 because among the six possible Schröder paths of semilength 2 only UHD contains UH. - Valerie Roitner, Jul 23 2020

Examples

			a(3)=5 since {3, (1+2), (1+(1+1)), (2+1), ((1+1)+1)} are the five weighted binary trees of weight 3.
G.f. = x + 2*x^2 + 5*x^3 + 15*x^4 + 51*x^5 + 188*x^6 + 731*x^7 + 2950*x^8 + 12235*x^9 + ... _Michael Somos_, Jan 17 2018
		

References

  • J. Brunvoll et al., Studies of some chemically relevant polygonal systems: mono-q-polyhexes, ACH Models in Chem., 133 (3) (1996), 277-298, Eq. 15.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

See A181768 for another version. - N. J. A. Sloane, Nov 12 2010
First column of triangle A104259. Row sums of absolute values of A091699.
Number of vertices of multiplihedron A121988.
m-th binomial transform of the Catalan numbers: A126930 (m = -2), A005043 (m = -1), A000108 (m = 0), A064613 (m = 2), A104455 (m = 3), A104498 (m = 4) and A154623 (m = 5).

Programs

  • Maple
    G := (1-sqrt(1-4*z/(1-z)))*1/2: Gser := series(G, z = 0, 30): seq(coeff(Gser, z, n), n = 1 .. 26); # Emeric Deutsch, Aug 12 2007
    seq(round(evalf(JacobiP(n-1,1,-n-1/2,9)/n,99)),n=1..25); # Peter Luschny, Sep 23 2014
  • Mathematica
    Rest@ CoefficientList[ InverseSeries[ Series[(y - y^2)/(1 + y - y^2), {y, 0, 26}], x], x] (* then A(x)=y(x); note that InverseSeries[Series[y-y^2, {y, 0, 24}], x] produces A000108(x) *) (* Len Smiley, Apr 10 2000 *)
    Range[0, 25]! CoefficientList[ Series[ Exp[ 3x] (BesselI[0, 2x] - BesselI[1, 2x]), {x, 0, 25}], x] (* Robert G. Wilson v, Apr 15 2011 *)
    a[n_] := Sum[ Binomial[n, k]*CatalanNumber[k], {k, 0, n}]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Aug 07 2012 *)
    Rest[CoefficientList[Series[3/2 - (1/2) Sqrt[(1 - 5 x)/(1 - x)], {x, 0, 40}], x]] (* Vincenzo Librandi, Nov 03 2014 *)
    Table[Hypergeometric2F1[1/2, -n+1, 2, -4], {n, 1, 30}] (* Vaclav Kotesovec, May 12 2022 *)
  • PARI
    {a(n) = my(A); if( n<2, n>0, A=vector(n); for(j=1,n, A[j] = 1 + sum(k=1,j-1, A[k]*A[j-k])); A[n])}; /* Michael Somos, May 23 2005 */
    
  • PARI
    {a(n) = if( n<1, 0, polcoeff( serreverse( (x - x^2) / (1 + x - x^2) + x * O(x^n)), n))}; /* Michael Somos, May 23 2005 */
    
  • PARI
    /* Offset = 0: */ {a(n)=local(A=1+x);for(i=1,n, A=sum(m=0,n, x^m*sum(k=0,m,A^k)+x*O(x^n))); polcoeff(A,n)} \\ Paul D. Hanna

Formula

(n+2)*a(n+2) = (6n+4)*a(n+1) - 5n*a(n).
G.f.: 3/2-(1/2)*sqrt((1-5*x)/(1-x)) [Gessel-Kim]. - N. J. A. Sloane, Jul 05 2014
G.f. for sequence doubled: (1/(2*x))*(1+x-(1-x)^(-1)*(1-x^2)^(1/2)*(1-5*x^2)^(1/2)).
a(n) = hypergeom([1/2, -n], [2], -4), n=0, 1, 2...; Integral representation as n-th moment of a positive function on a finite interval of the positive half-axis: a(n)=int(x^n*sqrt((5-x)/(x-1))/(2*Pi), x=1..5), n=0, 1, 2... This representation is unique. - Karol A. Penson, Sep 24 2001
a(1)=1, a(n)=1+sum(i=1, n-1, a(i)*a(n-i)). - Benoit Cloitre, Mar 16 2004
a(n) = Sum_{k=0..n} (-1)^k*3^(n-k)*binomial(n, k)*binomial(k, floor(k/2)) [offset 0]. - Paul Barry, Jan 27 2005
G.f. A(x) satisfies 0=f(x, A(x)) where f(x, y)=x-(1-x)(y-y^2). - Michael Somos, May 23 2005
G.f. A(x) satisfies 0=f(x, A(x), A(A(x))) where f(x, y, z)=x(z-z^2)+(x-1)y^2 . - Michael Somos, May 23 2005
G.f. (for offset 0): (-1+x+(1-6*x+5*x^2)^(1/2))/(2*(-x+x^2)).
G.f. =z*c(z/(1-z))/(1-z) = 1/2 - (1/2)sqrt(1-4z/(1-z)), where c(z)=(1-sqrt(1-4z))/(2z) is the Catalan function (follows from Michael Somos' first comment). - Emeric Deutsch, Aug 12 2007
G.f.: 1/(1-2x-x^2/(1-3x-x^2/(1-3x-x^2/(1-3x-x^2/(1-3x-x^2/(1-.... (continued fraction). - Paul Barry, Apr 19 2009
a(n) = Sum_{k, 0<=k<=n} A091965(n,k)*(-1)^k. - Philippe Deléham, Nov 28 2009
E.g.f.: exp(3x)*(I_0(2x)-I_1(2x)), where I_k(x) is a modified Bessel function of the first kind. - Emanuele Munarini, Apr 15 2011
If we prefix sequence with an additional term a(0)=1, g.f. is (3-3*x-sqrt(1-6*x+5*x^2))/(2*(1-x)). [See Kim, 2011] - N. J. A. Sloane, May 13 2011
From Gary W. Adamson, Jul 21 2011: (Start)
a(n) = upper left term in M^(n-1), M = an infinite square production matrix as follows:
2, 1, 0, 0, 0, 0, ...
1, 2, 1, 0, 0, 0, ...
1, 1, 2, 1, 0, 0, ...
1, 1, 1, 2, 1, 0, ...
1, 1, 1, 1, 2, 1, ...
1, 1, 1, 1, 1, 2, ...
... (End)
G.f. satisfies: A(x) = Sum_{n>=0} x^n * (1 - A(x)^(n+1))/(1 - A(x)); offset=0. - Paul D. Hanna, Nov 07 2011
G.f.: 1/x - 1/x/Q(0), 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
G.f.: (1-x - (1-5*x)*G(0))/(2*x*(1-x)), where G(k)= 1 + 4*x*(4*k+1)/( (4*k+2)*(1-x) - 2*x*(1-x)*(2*k+1)*(4*k+3)/(x*(4*k+3) + (1-x)*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 25 2013
Asymptotics (for offset 0): a(n) ~ 5^(n+3/2)/(8*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Jun 28 2013
G.f.: G(0)/(1-x), where G(k) = 1 + (4*k+1)*x/((k+1)*(1-x) - 2*x*(1-x)*(k+1)*(4*k+3)/(2*x*(4*k+3) + (2*k+3)*(1-x)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 29 2014
a(n) = JacobiP(n-1,1,-n-1/2,9)/n. - Peter Luschny, Sep 23 2014
0 = +a(n)*(+25*a(n+1) -50*a(n+2) +15*a(n+3)) +a(n+1)*(-10*a(n+1) +31*a(n+2) -14*a(n+3)) +a(n+2)*(+2*a(n+2) +a(n+3)) for all n in Z. - Michael Somos, Jan 17 2018
a(n+1) = (2/Pi) * Integral_{x = -1..1} (m + 4*x^2)^n*sqrt(1 - x^2) dx at m = 1. In general, the integral, qua sequence in n, gives the m-th binomial transform of the Catalan numbers. - Peter Bala, Jan 26 2020

A154930 Inverse of Fibonacci convolution array A154929.

Original entry on oeis.org

1, -2, 1, 5, -4, 1, -15, 14, -6, 1, 51, -50, 27, -8, 1, -188, 187, -113, 44, -10, 1, 731, -730, 468, -212, 65, -12, 1, -2950, 2949, -1956, 970, -355, 90, -14, 1, 12235, -12234, 8291, -4356, 1785, -550, 119, -16, 1, -51822, 51821, -35643, 19474, -8612, 3021
Offset: 0

Views

Author

Paul Barry, Jan 17 2009

Keywords

Comments

Alternating sign version of A104259. Row sums are (-1)^n*A033321. First column is (-1)^n*A007317.

Examples

			Triangle begins
1,
-2, 1,
5, -4, 1,
-15, 14, -6, 1,
51, -50, 27, -8, 1,
-188, 187, -113, 44, -10, 1,
731, -730, 468, -212, 65, -12, 1,
-2950, 2949, -1956, 970, -355, 90, -14, 1
Production array is
-2, 1,
1, -2, 1,
-1, 1, -2, 1,
1, -1, 1, -2, 1,
-1, 1, -1, 1, -2, 1,
1, -1, 1, -1, 1, -2, 1,
-1, 1, -1, 1, -1, 1, -2, 1
or ((1-x-x^2)/(1+x),x) beheaded.
		

Crossrefs

Formula

Riordan array ((1/(1+x))c(-x/(1+x)), (x/(1+x))c(x/(1+x))), c(x) the g.f. of A000108;
Riordan array ((sqrt(1+6x+5x^2)-x-1)/(2x(1+x)),(sqrt(1+6x+5x^2)-x-1)/ (2(1+x)));
Triangle T(n,k) = sum{j=0..n, (-1)^(n-k)*C(n,j)*C(2j-k,j-k)(k+1)/(j+1)}.
T(n,k) = T(n-1,k-1) -2*T(n-1,k) + Sum_{i, i>=0} T(n-1,k+1+i)*(-1)^i. - Philippe Deléham, Feb 23 2012

A171568 Riordan array (f(x), x*f(x)) where f(x) is the g.f. of A064613.

Original entry on oeis.org

1, 3, 1, 10, 6, 1, 37, 29, 9, 1, 150, 134, 57, 12, 1, 654, 622, 318, 94, 15, 1, 3012, 2948, 1686, 616, 140, 18, 1, 14445, 14317, 8781, 3693, 1055, 195, 21, 1, 71398, 71142, 45625, 21132, 7075, 1662, 259, 24, 1, 361114, 360602, 238170, 118042, 44303, 12345, 2464, 332, 27, 1
Offset: 0

Views

Author

Philippe Deléham, Dec 11 2009

Keywords

Comments

Equal to A171515*B = B*A104259, B = A007318.

Examples

			Triangle T(n,k) begins
[0]     1;
[1]     3,     1;
[2]    10,     6,     1;
[3]    37,    29,     9,     1;
[4]   150,   134,    57,    12,    1;
[5]   654,   622,   318,    94,   15,    1;
[6]  3012,  2948,  1686,   616,  140,   18,   1;
[7] 14445, 14317,  8781,  3693, 1055,  195,  21,  1;
[8] 71398, 71142, 45625, 21132, 7075, 1662, 259, 24, 1;
.
Production array begins
  3, 1
  1, 3, 1
  1, 1, 3, 1
  1, 1, 1, 3, 1
  1, 1, 1, 1, 3, 1
  1, 1, 1, 1, 1, 3, 1
- _Philippe Deléham_, Mar 05 2013
		

Crossrefs

Sum_{k=0..n} T(n,k)*x^k = A033543(n), A064613(n), A005572(n), A005573(n) for x = -1, 0, 1, 2 respectively.

Programs

  • Maple
    T := proc(n,k) option remember;
    if n < 0 or k < 0 then 0 elif n = k then 1 else
    T(n-1, k-1) + 3*T(n-1,k) + add(T(n-1, k+1+i), i=0..n) fi end:
    for n from 0 to 8 do seq(T(n,k), k = 0..n) od; # Peter Luschny, Oct 16 2022
  • Mathematica
    T[n_, k_] := T[n, k] = If[n < 0 || k < 0, 0, If[n == k, 1, T[n-1, k-1] + 3*T[n-1, k] + Sum[T[n-1, k+1+i], {i, 0, n}]]];
    Table[T[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 23 2024, after Peter Luschny *)

Formula

T(n, 0) - T(n, 1) = 2^n.
T(n, k) = T(n-1, k-1) + 3*T(n-1, k) + Sum_{i=0..n} T(n-1, k+1+i). - Philippe Deléham, Feb 23 2012

Extensions

Corrected and extended by Peter Luschny, Oct 16 2022

A171486 Riordan array (f(x), x*f(x)) where f(x) is the g.f. of A033321.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 6, 5, 3, 1, 21, 16, 9, 4, 1, 79, 58, 31, 14, 5, 1, 311, 224, 117, 52, 20, 6, 1, 1265, 900, 465, 205, 80, 27, 7, 1, 5275, 3720, 1910, 840, 330, 116, 35, 8, 1, 22431, 15713, 8034, 3532, 1396, 501, 161, 44, 9, 1, 96900, 67522, 34419, 15136, 6015, 2190
Offset: 0

Views

Author

Philippe Deléham, Dec 09 2009

Keywords

Comments

Equal to B*A065600 = A171224*B where B = A007318 ; equal to B*A039598*B^(-2).

Examples

			Triangle begins :
1
1, 1
2, 2, 1
6, 5, 3, 1
21, 16, 9, 4, 1
79, 58, 31, 14, 5, 1
311, 224, 117, 52, 20, 6, 1
		

Crossrefs

Formula

Sum_{k, 0<=k<=n} T(n,k)*x^k = A117641(n), A033321(n), A007317(n+1), A002212(n+1), A026378(n+1) for x = -1, 0, 1, 2, 3 respectively.
T(n,k) = T(n-1,k-1) + T(n-1,k) + sum_{i, i>=0} T(n-1,k+1+i)*2^i. - Philippe Deléham, Feb 23 2012

A155887 Riordan array (1, (1/(1-x))c(x/(1-x))), c(x) the g.f. of A000108.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 5, 4, 1, 0, 15, 14, 6, 1, 0, 51, 50, 27, 8, 1, 0, 188, 187, 113, 44, 10, 1, 0, 731, 730, 468, 212, 65, 12, 1, 0, 2950, 2949, 1956, 970, 355, 90, 14, 1, 0, 12235, 12234, 8291, 4356, 1785, 550, 119, 16, 1, 0, 51822, 51821, 35643, 19474, 8612, 3021
Offset: 0

Views

Author

Paul Barry, Jan 29 2009

Keywords

Comments

Inverse of (1, x(1-x)/(1+x-x^2)). Row sums are A002212. Augmented version of A104259.
Triangle T(n,k) given by [0,F(3)/F(1),F(1)/F(3),F(5)/F(3),F(3)/F(5),F(7)/F(5),F(5)/F(7),...] DELTA [1,0,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938, F(n)=A000045(n) (Fibonacci numbers). - Philippe Deléham, Jan 31 2009

Examples

			Triangle begins
  1;
  0,    1;
  0,    2,    1;
  0,    5,    4,    1;
  0,   15,   14,    6,   1;
  0,   51,   50,   27,   8,   1;
  0,  188,  187,  113,  44,  10,  1;
  0,  731,  730,  468, 212,  65, 12,  1;
  0, 2950, 2949, 1956, 970, 355, 90, 14, 1;
From _Paul Barry_, Sep 28 2009: (Start)
Production matrix is
  0, 1,
  0, 2, 1,
  0, 1, 2, 1,
  0, 1, 1, 2, 1,
  0, 1, 1, 1, 2, 1,
  0, 1, 1, 1, 1, 2, 1,
  0, 1, 1, 1, 1, 1, 2, 1,
  0, 1, 1, 1, 1, 1, 1, 2, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 2, 1 (End)
		

Formula

Riordan array (1, (1-sqrt((1-5x)/(1-x)))/(2x));
G.f.: 1/(1-xy/(1-x-x/(1-x/(1-x-x/(1-x/(1-x-x/(1-... (continued fraction).

A202847 Riordan array (f(x), x*f(x)) where f(x) is the g.f. of A126930.

Original entry on oeis.org

1, -1, 1, 2, -2, 1, -3, 5, -3, 1, 6, -10, 9, -4, 1, -10, 22, -22, 14, -5, 1, 20, -44, 54, -40, 20, -6, 1, -35, 93, -123, 109, -65, 27, -7, 1, 70, -186, 281, -276, 195, -98, 35, -8, 1, -126, 386, -618, 682, -541, 321, -140, 44, -9, 1
Offset: 0

Views

Author

Philippe Deléham, Mar 02 2013

Keywords

Examples

			Triangle begins
1
-1, 1
2, -2, 1
-3, 5, -3, 1
6, -10, 9, -4, 1
-10, 22, -22, 14, -5, 1
20, -44, 54, -40, 20, -6, 1
-35, 93, -123, 109, -65, 27, -7, 1
...
Production matrix begins
x, 1
1, x, 1
1, 1, x, 1
1, 1, 1, x, 1
1, 1, 1, 1, x, 1
1, 1, 1, 1, 1, x, 1
1, 1, 1, 1, 1, 1, x, 1
1, 1, 1, 1, 1, 1, 1, x, 1
1, 1, 1, 1, 1, 1, 1, 1, x, 1
1, 1, 1, 1, 1, 1, 1, 1, 1, x, 1
..., with x = -1.
		

Crossrefs

Cf. (sequences with similar production matrix) A097609 (x=0), A033184 (x=1), A104259 (x=2), A171568 (x=3), A171589 (x=4)

Formula

T(n,k) = (-1)^(n-k)*A054336(n,k).
Sum_{k, 0<=k<=n}T(n,k)*x^k = (-1)^n*A126931(n), (-1)^n*A054341(n), A126930(n), A126120(n), A001405(n), A054341(n), A126931(n) for x = -2, -1, 0, 1, 2, 3, 4 respectively.
Showing 1-9 of 9 results.