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-10 of 12 results. Next

A005043 Riordan numbers: a(n) = (n-1)*(2*a(n-1) + 3*a(n-2))/(n+1).

Original entry on oeis.org

1, 0, 1, 1, 3, 6, 15, 36, 91, 232, 603, 1585, 4213, 11298, 30537, 83097, 227475, 625992, 1730787, 4805595, 13393689, 37458330, 105089229, 295673994, 834086421, 2358641376, 6684761125, 18985057351, 54022715451, 154000562758, 439742222071, 1257643249140
Offset: 0

Views

Author

Keywords

Comments

Also called Motzkin summands or ring numbers.
The old name was "Motzkin sums", used in certain publications. The sequence has the property that Motzkin(n) = A001006(n) = a(n) + a(n+1), e.g., A001006(4) = 9 = 3 + 6 = a(4) + a(5).
Number of 'Catalan partitions', that is partitions of a set 1,2,3,...,n into parts that are not singletons and whose convex hulls are disjoint when the points are arranged on a circle (so when the parts are all pairs we get Catalan numbers). - Aart Blokhuis (aartb(AT)win.tue.nl), Jul 04 2000
Number of ordered trees with n edges and no vertices of outdegree 1. For n > 1, number of dissections of a convex polygon by nonintersecting diagonals with a total number of n+1 edges. - Emeric Deutsch, Mar 06 2002
Number of Motzkin paths of length n with no horizontal steps at level 0. - Emeric Deutsch, Nov 09 2003
Number of Dyck paths of semilength n with no peaks at odd level. Example: a(4)=3 because we have UUUUDDDD, UUDDUUDD and UUDUDUDD, where U=(1,1), D=(1,-1). Number of Dyck paths of semilength n with no ascents of length 1 (an ascent in a Dyck path is a maximal string of up steps). Example: a(4)=3 because we have UUUUDDDD, UUDDUUDD and UUDUUDDD. - Emeric Deutsch, Dec 05 2003
Arises in Schubert calculus as follows. Let P = complex projective space of dimension n+1. Take n projective subspaces of codimension 3 in P in general position. Then a(n) is the number of lines of P intersecting all these subspaces. - F. Hirzebruch, Feb 09 2004
Difference between central trinomial coefficient and its predecessor. Example: a(6) = 15 = 141 - 126 and (1 + x + x^2)^6 = ... + 126*x^5 + 141*x^6 + ... (Catalan number A000108(n) is the difference between central binomial coefficient and its predecessor.) - David Callan, Feb 07 2004
a(n) = number of 321-avoiding permutations on [n] in which each left-to-right maximum is a descent (i.e., is followed by a smaller number). For example, a(4) counts 4123, 3142, 2143. - David Callan, Jul 20 2005
The Hankel transform of this sequence give A000012 = [1, 1, 1, 1, 1, 1, 1, ...]; example: Det([1, 0, 1, 1; 0, 1, 1, 3; 1, 1, 3, 6; 1, 3, 6, 15]) = 1. - Philippe Deléham, May 28 2005
The number of projective invariants of degree 2 for n labeled points on the projective line. - Benjamin J. Howard (bhoward(AT)ima.umn.edu), Nov 24 2006
Define a random variable X=trA^2, where A is a 2 X 2 unitary symplectic matrix chosen from USp(2) with Haar measure. The n-th central moment of X is E[(X+1)^n] = a(n). - Andrew V. Sutherland, Dec 02 2007
Let V be the adjoint representation of the complex Lie algebra sl(2). The dimension of the invariant subspace of the n-th tensor power of V is a(n). - Samson Black (sblack1(AT)uoregon.edu), Aug 27 2008
Starting with offset 3 = iterates of M * [1,1,1,...], 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 08 2009
a(n) has the following standard-Young-tableaux (SYT) interpretation: binomial(n+1,k)*binomial(n-k-1,k-1)/(n+1)=f^(k,k,1^{n-2k}) where f^lambda equals the number of SYT of shape lambda. - Amitai Regev (amotai.regev(AT)weizmann.ac.il), Mar 02 2010
a(n) is also the sum of the numbers of standard Young tableaux of shapes (k,k,1^{n-2k}) for all 1 <= k <= floor(n/2). - Amitai Regev (amotai.regev(AT)weizmann.ac.il), Mar 10 2010
a(n) is the number of derangements 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(3)=1 because p=231=(123) is the only derangement of {1,2,3} with genus 0. Indeed, cp'=231*312=123=(1)(2)(3) and so g(p) = (1/2)(3+1-1-3)=0. - Emeric Deutsch, May 29 2010
Apparently: Number of Dyck 2n-paths with all ascents length 2 and no descent length 2. - David Scambler, Apr 17 2012
This is true. Proof: The mapping "insert a peak (UD) after each upstep (U)" is a bijection from all Dyck n-paths to those Dyck (2n)-paths in which each ascent is of length 2. It sends descents of length 1 in the n-path to descents of length 2 in the (2n)-path. But Dyck n-paths with no descents of length 1 are equinumerous with Riordan n-paths (Motzkin n-paths with no flatsteps at ground level) as follows. Given a Dyck n-path with no descents of length 1, split it into consecutive step pairs, then replace UU with U, DD with D, UD with a blue flatstep (F), DU with a red flatstep, and concatenate the new steps to get a colored Motzkin path. Each red F will be (immediately) preceded by a blue F or a D. In the latter case, transfer the red F so that it precedes the matching U of the D. Finally, erase colors to get the required Riordan path. For example, with lowercase f denoting a red flatstep, U^5 D^2 U D^4 U^4 D^3 U D^2 -> (U^2, U^2, UD, DU, D^2, D^2, U^2, U^2 D^2, DU, D^2) -> UUFfDDUUDfD -> UUFFDDUFUDD. - David Callan, Apr 25 2012
From Nolan Wallach, Aug 20 2014: (Start)
Let ch[part1, part2] be the value of the character of the symmetric group on n letters corresponding to the partition part1 of n on the conjucgacy class given by part2. Let A[n] be the set of (n+1) partitions of 2n with parts 1 or 2. Then deleting the first term of the sequence one has a(n) = Sum_{k=1..n+1} binomial(n,k-1)*ch[[n,n], A[n][[k]]])/2^n. This via the Frobenius Character Formula can be interpreted as the dimension of the SL(n,C) invariants in tensor^n (wedge^2 C^n).
Explanation: Let p_j denote sum (x_i)^j the sum in k variables. Then the Frobenius formula says then (p_1)^j_1 (p_2)^j_2 ... (p_r)^j_r is equal to sum(lambda, ch[lambda, 1^j_12^j_2 ... r^j_r] S_lambda) with S_lambda the Schur function corresponding to lambda. This formula implies that the coefficient of S([n,n]) in (((p_1)^1+p_2)/2)^n in its expansion in terms of Schur functions is the right hand side of our formula. If we specialize the number of variables to 2 then S[n,n](x,y)=(xy)^n. Which when restricted to y=x^(-1) is 1. That is it is 1 on SL(2).
On the other hand ((p_1)^2+p_2)/2 is the complete homogeneous symmetric function of degree 2 that is tr(S^2(X)). Thus our formula for a(n) is the same as that of Samson Black above since his V is the same as S^2(C^2) as a representation of SL(2). On the other hand, if we multiply ch(lambda) by sgn you get ch(Transpose(lambda)). So ch([n,n]) becomes ch([2,...,2]) (here there are n 2's). The formula for a(n) is now (1/2^n)*Sum_{j=0..n} ch([2,..,2], 1^(2n-2j) 2^j])*(-1)^j)*binomial(n,j), which calculates the coefficient of S_(2,...,2) in (((p_1)^2-p_2)/2)^n. But ((p_1)^2-p_2)/2 in n variables is the second elementary symmetric function which is the character of wedge^2 C^n and S_(2,...,2) is 1 on SL(n).
(End)
a(n) = number of noncrossing partitions (A000108) of [n] that contain no singletons, also number of nonnesting partitions (A000108) of [n] that contain no singletons. - David Callan, Aug 27 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)].
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).
Various relations among the o.g.f.s may be easily constructed, such as Fib[-Mot(-x)] = -P[P(-x)] = x/(1-2*x) a generating fct for 2^n.
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)
Consistent with David Callan's comment above, A249548, provides a refinement of the Motzkin sums into the individual numbers for the non-crossing partitions he describes. - Tom Copeland, Nov 09 2014
The number of lattice paths from (0,0) to (n,0) that do not cross below the x-axis and use up-step=(1,1) and down-steps=(1,-k) where k is a positive integer. For example, a(4) = 3: [(1,1)(1,1)(1,-1)(1,-1)], [(1,1)(1,-1)(1,1)(1,-1)] and [(1,1)(1,1)(1,1)(1,-3)]. - Nicholas Ham, Aug 19 2015
A series created using 2*(a(n) + a(n+1)) + (a(n+1) + a(n+2)) has Hankel transform of F(2n), offset 3, F being a Fibonacci number, A001906 (Empirical observation). - Tony Foster III, Jul 30 2016
The series a(n) + A001006(n) has Hankel transform F(2n+1), offset n=1, F being the Fibonacci bisection A001519 (empirical observation). - Tony Foster III, Sep 05 2016
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. Moreover, the number of such algebras having the double centraliser property with respect to a minimal faithful projective-injective module equals the number of Riordan paths, that is, Motzkin paths without level-steps at height zero, of length n." - Eric M. Schmidt, Dec 16 2017
A connection to the Thue-Morse sequence: (-1)^a(n) = (-1)^A010060(n) * (-1)^A010060(n+1) = A106400(n) * A106400(n+1). - Vladimir Reshetnikov, Jul 21 2019
Named by Bernhart (1999) after the American mathematician John Riordan (1903-1988). - Amiram Eldar, Apr 15 2021

Examples

			a(5)=6 because the only dissections of a polygon with a total number of 6 edges are: five pentagons with one of the five diagonals and the hexagon with no diagonals.
G.f. = 1 + x^2 + x^3 + 3*x^4 + 6*x^5 + 15*x^6 + 36*x^7 + 91*x^8 + 232*x^9 + ...
From _Gus Wiseman_, Nov 15 2022: (Start)
The a(0) = 1 through a(6) = 15 lone-child-avoiding (no vertices of outdegree 1) ordered rooted trees with n + 1 vertices (ranked by A358376):
  o  .  (oo)  (ooo)  (oooo)   (ooooo)   (oooooo)
                     ((oo)o)  ((oo)oo)  ((oo)ooo)
                     (o(oo))  ((ooo)o)  ((ooo)oo)
                              (o(oo)o)  ((oooo)o)
                              (o(ooo))  (o(oo)oo)
                              (oo(oo))  (o(ooo)o)
                                        (o(oooo))
                                        (oo(oo)o)
                                        (oo(ooo))
                                        (ooo(oo))
                                        (((oo)o)o)
                                        ((o(oo))o)
                                        ((oo)(oo))
                                        (o((oo)o))
                                        (o(o(oo)))
(End)
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Row sums of triangle A020474, first differences of A082395.
First diagonal of triangular array in A059346.
Binomial transform of A126930. - Philippe Deléham, Nov 26 2009
The Hankel transform of a(n+1) is A128834. The Hankel transform of a(n+2) is floor((2*n+4)/3) = A004523(n+2). - Paul Barry, Mar 08 2011
The Kn11 triangle sums of triangle A175136 lead to A005043(n+2), while the Kn12(n) = A005043(n+4)-2^(n+1), Kn13(n) = A005043(n+6)-(n^2+9*n+56)*2^(n-2) and the Kn4(n) = A005043(2*n+2) = A099251(n+1) triangle sums are related to the sequence given above. For the definitions of these triangle sums see A180662. - Johannes W. Meijer, May 06 2011
Cf. A187306 (self-convolution), A348210 (column 1).
Bisections: A099251, A099252.

Programs

  • Haskell
    a005043 n = a005043_list !! n
    a005043_list = 1 : 0 : zipWith div
       (zipWith (*) [1..] (zipWith (+)
           (map (* 2) $ tail a005043_list) (map (* 3) a005043_list))) [3..]
    -- Reinhard Zumkeller, Jan 31 2012
    
  • Maple
    A005043 := proc(n) option remember; if n <= 1 then 1-n else (n-1)*(2*A005043(n-1)+3*A005043(n-2))/(n+1); fi; end;
    Order := 20: solve(series((x-x^2)/(1-x+x^2),x)=y,x); # outputs g.f.
  • Mathematica
    a[0]=1; a[1]=0; a[n_]:= a[n] = (n-1)*(2*a[n-1] + 3*a[n-2])/(n+1); Table[ a[n], {n, 0, 30}] (* Robert G. Wilson v, Jun 14 2005 *)
    Table[(-3)^(1/2)/6 * (-1)^n*(3*Hypergeometric2F1[1/2,n+1,1,4/3]+ Hypergeometric2F1[1/2,n+2,1,4/3]), {n,0,32}] (* cf. Mark van Hoeij in A001006 *) (* Wouter Meeussen, Jan 23 2010 *)
    RecurrenceTable[{a[0]==1,a[1]==0,a[n]==(n-1) (2a[n-1]+3a[n-2])/(n+1)},a,{n,30}] (* Harvey P. Dale, Sep 27 2013 *)
    a[ n_]:= SeriesCoefficient[2/(1+x +Sqrt[1-2x-3x^2]), {x, 0, n}]; (* Michael Somos, Aug 21 2014 *)
    a[ n_]:= If[n<0, 0, 3^(n+3/2) Hypergeometric2F1[3/2, n+2, 2, 4]/I]; (* Michael Somos, Aug 21 2014 *)
    Table[3^(n+3/2) CatalanNumber[n] (4(5+2n)Hypergeometric2F1[3/2, 3/2, 1/2-n, 1/4] -9 Hypergeometric2F1[3/2, 5/2, 1/2 -n, 1/4])/(4^(n+3) (n+1)), {n, 0, 31}] (* Vladimir Reshetnikov, Jul 21 2019 *)
    Table[Sqrt[27]/8 (3/4)^n CatalanNumber[n] Hypergeometric2F1[1/2, 3/2, 1/2 - n, 1/4], {n, 0, 31}] (* Jan Mangaldan, Sep 12 2021 *)
  • Maxima
    a[0]:1$
    a[1]:0$
    a[n]:=(n-1)*(2*a[n-1]+3*a[n-2])/(n+1)$
    makelist(a[n],n,0,12); /* Emanuele Munarini, Mar 02 2011 */
    
  • PARI
    {a(n) = if( n<0, 0, n++; polcoeff( serreverse( (x - x^3) / (1 + x^3) + x * O(x^n)), n))}; /* Michael Somos, May 31 2005 */
    
  • PARI
    my(N=66); Vec(serreverse(x/(1+x*sum(k=1,N,x^k))+O(x^N))) \\ Joerg Arndt, Aug 19 2012
    
  • Python
    from functools import cache
    @cache
    def A005043(n: int) -> int:
        if n <= 1: return 1 - n
        return (n - 1) * (2 * A005043(n - 1) + 3 * A005043(n - 2)) // (n + 1)
    print([A005043(n) for n in range(32)]) # Peter Luschny, Nov 20 2022
  • Sage
    A005043 = lambda n: (-1)^n*jacobi_P(n,1,-n-3/2,-7)/(n+1)
    [simplify(A005043(n)) for n in (0..29)]
    # Peter Luschny, Sep 23 2014
    
  • Sage
    def ms():
        a, b, c, d, n = 0, 1, 1, -1, 1
        yield 1
        while True:
            yield -b + (-1)^n*d
            n += 1
            a, b = b, (3*(n-1)*n*a+(2*n-1)*n*b)/((n+1)*(n-1))
            c, d = d, (3*(n-1)*c-(2*n-1)*d)/n
    A005043 = ms()
    print([next(A005043) for  in range(32)]) # _Peter Luschny, May 16 2016
    

Formula

a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(n, k)*A000108(k). a(n) = (1/(n+1)) * Sum_{k=0..ceiling(n/2)} binomial(n+1, k)*binomial(n-k-1, k-1), for n > 1. - Len Smiley. [Comment from Amitai Regev (amitai.regev(AT)weizmann.ac.il), Mar 02 2010: the latter sum should be over the range k=1..floor(n/2).]
G.f.: (1 + x - sqrt(1-2*x-3*x^2))/(2*x*(1+x)).
G.f.: 2/(1+x+sqrt(1-2*x-3*x^2)). - Paul Peart (ppeart(AT)fac.howard.edu), May 27 2000
a(n+1) + (-1)^n = a(0)*a(n) + a(1)*a(n-1) + ... + a(n)*a(0). - Bernhart
a(n) = (1/(n+1)) * Sum_{i} (-1)^i*binomial(n+1, i)*binomial(2*n-2*i, n-i). - Bernhart
G.f. A(x) satisfies A = 1/(1+x) + x*A^2.
E.g.f.: exp(x)*(BesselI(0, 2*x) - BesselI(1, 2*x)). - Vladeta Jovovic, Apr 28 2003
a(n) = A001006(n-1) - a(n-1).
a(n+1) = Sum_{k=0..n} (-1)^k*A026300(n, k), where A026300 is the Motzkin triangle.
a(n) = Sum_{k=0..n} (-1)^k*binomial(n, k)*binomial(k, floor(k/2)). - Paul Barry, Jan 27 2005
a(n) = Sum_{k>=0} A086810(n-k, k). - Philippe Deléham, May 30 2005
a(n+2) = Sum_{k>=0} A064189(n-k, k). - Philippe Deléham, May 31 2005
Moment representation: a(n) = (1/(2*Pi))*Int(x^n*sqrt((1+x)(3-x))/(1+x),x,-1,3). - Paul Barry, Jul 09 2006
Inverse binomial transform of A000108 (Catalan numbers). - Philippe Deléham, Oct 20 2006
a(n) = (2/Pi)* Integral_{x=0..Pi} (4*cos(x)^2-1)^n*sin(x)^2 dx. - Andrew V. Sutherland, Dec 02 2007
G.f.: 1/(1-x^2/(1-x-x^2/(1-x-x^2/(1-x-x^2/(1-... (continued fraction). - Paul Barry, Jan 22 2009
G.f.: 1/(1+x-x/(1-x/(1+x-x/(1-x/(1+x-x/(1-... (continued fraction). - Paul Barry, May 16 2009
G.f.: 1/(1-x^2/(1-x/(1-x/(1-x^2/(1-x/(1-x/(1-x^2/(1-x/(1-... (continued fraction). - Paul Barry, Mar 02 2010
a(n) = -(-1)^n * hypergeom([1/2, n+2],[2],4/3) / sqrt(-3). - Mark van Hoeij, Jul 02 2010
a(n) = (-1)^n*hypergeometric([-n,1/2],[2],4). - Peter Luschny, Aug 15 2012
Let A(x) be the g.f., then x*A(x) is the reversion of x/(1 + x^2*Sum_{k>=0} x^k); see A215340 for the correspondence to Dyck paths without length-1 ascents. - Joerg Arndt, Aug 19 2012 and Apr 16 2013
a(n) ~ 3^(n+3/2)/(8*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 02 2012
G.f.: 2/(1+x+1/G(0)), where G(k) = 1 + x*(2+3*x)*(4*k+1)/( 4*k+2 - x*(2+3*x)*(4*k+2)*(4*k+3)/(x*(2+3*x)*(4*k+3) + 4*(k+1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 05 2013
D-finite (an alternative): (n+1)*a(n) = 3*(n-2)*a(n-3) + (5*n-7)*a(n-2) + (n-2)*a(n-1), n >= 3. - Fung Lam, Mar 22 2014
Asymptotics: a(n) = (3^(n+2)/(sqrt(3*n*Pi)*(8*n)))*(1-21/(16*n) + O(1/n^2)) (with contribution by Vaclav Kotesovec). - Fung Lam, Mar 22 2014
a(n) = T(2*n-1,n)/n, where T(n,k) = triangle of A180177. - Vladimir Kruchinin, Sep 23 2014
a(n) = (-1)^n*JacobiP(n,1,-n-3/2,-7)/(n+1). - Peter Luschny, Sep 23 2014
a(n) = Sum_{k=0..n} C(n,k)*(C(k,n-k)-C(k,n-k-1)). - Peter Luschny, Oct 01 2014
Conjecture: a(n) = A002426(n) - A005717(n), n > 0. - Mikhail Kurkov, Feb 24 2019 [The conjecture is true. - Amiram Eldar, May 17 2024]
a(n) = A309303(n) + A309303(n+1). - Vladimir Reshetnikov, Jul 22 2019
From Peter Bala, Feb 11 2022: (Start)
a(n) = A005773(n+1) - 2*A005717(n) for n >= 1.
Conjectures: for n >= 1, n divides a(2*n+1) and 2*n-1 divides a(2*n). (End)

Extensions

Thanks to Laura L. M. Yang (yanglm(AT)hotmail.com) for a correction, Aug 29 2004
Name changed to Riordan numbers following a suggestion from Ira M. Gessel. - N. J. A. Sloane, Jul 24 2020

A007477 Shifts 2 places left when convolved with itself.

Original entry on oeis.org

1, 1, 1, 2, 3, 6, 11, 22, 44, 90, 187, 392, 832, 1778, 3831, 8304, 18104, 39666, 87296, 192896, 427778, 951808, 2124135, 4753476, 10664458, 23981698, 54045448, 122041844, 276101386, 625725936, 1420386363, 3229171828, 7351869690, 16760603722, 38258956928, 87437436916, 200057233386, 458223768512, 1050614664580
Offset: 0

Views

Author

Keywords

Comments

Words of length n in language defined by L = 1 + a + (L)L: L(0)=1, L(1)=a, L(2)=(), L(3)=(a)+()a, L(4)=(())+(a)a+()(), ...
Series reversion of x*A(x) is x*A082582(-x). - Michael Somos, Jul 22 2003
a(n) = number of Motzkin n-paths (A001006) in which no flatstep (F) is immediately followed by either an upstep (U) or a flatstep, in other words, each flatstep is either followed by a downstep (D) or ends the path. For example, a(4)=3 counts UDUD, UFDF, UUDD. - David Callan, Jun 07 2006
a(n) = number of Dyck (n+1)-paths (A000108) containing no UDUs and no subpaths of the form UUPDD where P is a nonempty Dyck path. For example, a(4)=3 counts UUUDDUUDDD, UUDDUUDDUD, UUUDDUDDUD. - David Callan, Sep 25 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_0*a_{n-1} + a_1*a_{n-2} + ... + a_{n-1}*a_0 for n >= t. For example phi([1]) is the Catalan numbers A000108. The present sequence is (essentially) phi([0,1,1]). - Gary W. Adamson and R. J. Mathar, Oct 27 2008
The Kn21(n) triangle sums of A175136 lead to A007477(n+1), while the Kn22(n) = A007477(n+3)-1, Kn23(n) = A007477(n+5)-(4+n) and Kn3(n) = A007477(2*n+1) triangle sums of A175136 are related to the sequence given above. For the definition of these triangle sums see A180662. - Johannes W. Meijer, May 06 2011
For n>=2, a(n) gives number of possible, ways to parse an English sentence consisting of just n+1 copies of word "buffalo", with one particular "plausible" grammar. See the Wikipedia page and my Python source at OEIS Wiki. Whether these are really intelligible sentences is of course debatable. See A213705 for a more plausible example in the Finnish language. - Antti Karttunen, Sep 14 2012

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a007477 n = a007477_list !! n
    a007477_list = 1 : 1 : f [1,1] where
       f xs = y : f (y:xs) where y = sum $ zipWith (*) (tail xs) (reverse xs)
    -- Reinhard Zumkeller, Apr 09 2012
    
  • Maple
    A007477 := proc(n) option remember; local k; if n <= 1 then 1 else add(A007477(k)*A007477(n-k-2),k=0..n-2); fi; end;
    unprotect(phi);
    phi:=proc(t,u,M) local i,a;
    a:=Array(0..M); for i from 0 to t-1 do a[i]:=u[i+1]; od:
    for i from t to M do a[i]:=add(a[j]*a[i-1-j],j=0..i-1); od:
    [seq(a[i],i=0..M)]; end;
    phi(3,[0,1,1],30);
    # N. J. A. Sloane, Nov 02 2008
  • Mathematica
    f[x_] := (1 - Sqrt[1 - 4x^2 - 4x^3])/2; Drop[ CoefficientList[ Series[f[x], {x, 0, 32}], x], 2] (* Jean-François Alcover, Nov 22 2011, after Pari *)
    a[n_] := Sum[Binomial[2*k+2, n-k-2]*Binomial[n-k-2, k]/(k+1), {k, 0, n-2}]; a[0] = a[1] = 1; Array[a, 40, 0] (* Jean-François Alcover, Mar 04 2016, after Vladimir Kruchinin *)
  • Maxima
    a(n):=if n<2 then 1 else sum((binomial(2*k+2,n-k-2)*binomial(n-k-2,k))/(k+1),k,0,n-2); /* Vladimir Kruchinin, Nov 22 2014 */
  • PARI
    a(n)=polcoeff((1-sqrt(1-4*x^2-4*x^3+x^3*O(x^n)))/2,n+2)
    

Formula

a(n) = sum( a(k) * a(n-2-k) ), n>1.
G.f. A(x) satisfies the equation 0 = 1 + x - A(x) + (x*A(x))^2.
The g.f. satisfies A(x)-x^2*A(x)^2 = 1+x. - Ralf Stephan, Jun 30 2003
G.f.: (1-sqrt(1-4x^2-4x^3))/(2x^2).
G.f.: (1+x)c(x^2(1+x)) where c(x) is g.f. of A000108. - Paul Barry, May 31 2006
G.f.: 1/(1-x/(1-x^2/(1-x^2/(1-x/(1-x^2/(1-x^2/(1-x/(1-x^2/(1-x^2/(1-... (continued fraction). - Paul Barry, Jul 30 2010
D-finite with recurrence: (n+2)*a(n) +(n+1)*a(n-1) +4*(-n+1)*a(n-2) +2*(-4*n+9)*a(n-3) +2*(-2*n+7)*a(n-4)=0. - R. J. Mathar, Dec 02 2012
a(n) = Sum_{k=0..n-2} binomial(2*k+2,n-k-2)*binomial(n-k-2,k)/(k+1), n>1, a(0)=1, a(1)=1. - Vladimir Kruchinin, Nov 22 2014
a(n) = Sum_{k=0..n-1} (-1)^(n-1-k)*binomial(n-1,k)*A082582(k+2), for n>0. - Thomas Baruchel, Jan 22 2015
a(n) ~ sqrt(3 - 4*r^2) * (4*r)^n * (1+r)^(n+1) / (sqrt(Pi)*n^(3/2)), where r = 0.41964337760708056627592628232664330021208937304879612338939... is the root of the equation 4*r^2*(1+r) = 1. - Vaclav Kotesovec, Jul 03 2021

Extensions

Additional comments from Michael Somos, Aug 03 2000

A091894 Triangle read by rows: T(n,k) is the number of Dyck paths of semilength n, having k ddu's [here u = (1,1) and d = (1,-1)].

Original entry on oeis.org

1, 1, 2, 4, 1, 8, 6, 16, 24, 2, 32, 80, 20, 64, 240, 120, 5, 128, 672, 560, 70, 256, 1792, 2240, 560, 14, 512, 4608, 8064, 3360, 252, 1024, 11520, 26880, 16800, 2520, 42, 2048, 28160, 84480, 73920, 18480, 924, 4096, 67584, 253440, 295680, 110880, 11088, 132
Offset: 0

Views

Author

Emeric Deutsch, Mar 10 2004

Keywords

Comments

Number of Dyck paths of semilength n, having k uu's with midpoint at even height. Example: T(4,1) = 6 because we have u(uu)duddd, u(uu)udddd, udu(uu)ddd, u(uu)dddud, u(uu)ddudd and uud(uu)ddd [here u = (1,1), d = (1,-1) and the uu's with midpoint at even height are shown between parentheses]. Row sums are the Catalan numbers (A000108). T(2n+1,n) = A000108(n) (the Catalan numbers). Sum_{k>=0} k*T(n,k) = binomial(2n-2,n-3) = A002694(n-1).
Sometimes called the Touchard distribution (after Touchard's Catalan number identity). T(n,k) = number of full binary trees on 2n edges with k deep interior vertices (deep interior means you have to traverse at least 2 edges to reach a leaf) = number of binary trees on n-1 edges with k vertices having a full complement of 2 children. - David Callan, Jul 19 2004
From David Callan, Oct 25 2004: (Start)
T(n,k) = number of ordered trees on n edges with k prolific edges. A prolific edge is one whose child vertex has at least two children. For example with n=3, drawing ordered trees down from the root, /|\ has no prolific edge and the only tree with one prolific edge has the shape of an inverted Y, so T(3,1)=1.
Proof: Consider the following bijection, recorded by Emeric Deutsch, from ordered trees on n edges to Dyck n-paths. For a given ordered tree, traverse the tree in preorder (walk-around from root order). To each node of outdegree r there correspond r upsteps followed by 1 downstep; nothing corresponds to the last leaf. This bijection sends prolific edges to noninitial ascents of length >=2, that is, to DUU's. Then reverse the resulting Dyck n-path so that prolific edges correspond to DDU's. (End)
T(n,k) is the number of Łukasiewicz paths of length n having k fall steps (1,-1) that start at an even level. A Łukasiewicz path of length n is a path in the first quadrant from (0,0) to (n,0) using rise steps (1,k) for any positive integer k, level steps (1,0) and fall steps (1,-1) (see R. P. Stanley, Enumerative Combinatorics, Vol. 2, Cambridge Univ. Press, Cambridge, 1999, p. 223, Exercise 6.19w; the integers are the slopes of the steps). Example: T(3,1) = 1 because we have U(2)(D)D, where U(2) = (1,2), D = (1,-1) and the fall steps that start at an even level are shown between parentheses. Row n contains ceiling(n/2) terms (n >= 1). - Emeric Deutsch, Jan 06 2005
Number of binary trees with n-1 edges and k+1 leaves (a binary tree is a rooted tree in which each vertex has at most two children and each child of a vertex is designated as its left or right child). - Emeric Deutsch, Jul 31 2006
Number of full binary trees with 2n edges and k+1 vertices both children of which are leaves (n >= 1; a full binary tree is a rooted tree in which each vertex has either 0 or two children). - Emeric Deutsch, Dec 26 2006
Number of ordered trees with n edges and k jumps. In the preorder traversal of an ordered tree, any transition from a node at a deeper level to a node on a strictly higher level is called a jump. - Emeric Deutsch, Jan 18 2007
It is remarkable that we can generate the coefficients of the right hand columns of triangle A175136 with the aid of the coefficients in the rows of the triangle given above. See A175136 for more information. - Johannes W. Meijer, May 06 2011
The antidiagonal sums equal A152225. - Johannes W. Meijer, Sep 13 2012
This array also counts 231-avoiding permutations according to the number of peaks, i.e., positions w[i-1] < w[i] > w[i+1]. For example, 123, 213, 312, and 321 have no peaks, while 132 has one peak. Note also T(n,k) = 2^(n - 1 - 2*k)*A055151(n-1,k). - Kyle Petersen, Aug 02 2013

Examples

			T(4,1) = 6 because we have uduu(ddu)d, uu(ddu)dud, uuu(ddu)dd, uu(ddu)udd, uudu(ddu)d and uuud(ddu)d [here u = (1,1), d = (1,-1) and the ddu's are shown between parentheses].
Triangle begins:
    1;
    1;
    2;
    4,    1;
    8,    6;
   16,   24,    2;
   32,   80,   20;
   64,  240,  120,   5;
  128,  672,  560,  70;
  256, 1792, 2240, 560, 14;
  ...
		

References

  • T. K. Petersen, Eulerian Numbers, Birkhauser, 2015, Section 4.3.

Crossrefs

The first few columns equal A011782, A001788, 2*A003472, 5*A002409, 14*A140325 and 42*A172242. - Johannes W. Meijer, Sep 13 2012

Programs

  • GAP
    T:=Concatenation([1],Flat(List([1..13],n->List([0..Int((n-1)/2)],k->2^(n-2*k-1)*Binomial(n-1,2*k)*Binomial(2*k,k)/(k+1))))); # Muniru A Asiru, Nov 29 2018
    
  • Maple
    a := proc(n,k) if n=0 and k=0 then 1 elif n=0 then 0 else 2^(n-2*k-1)*binomial(n-1,2*k)*binomial(2*k,k)/(k+1) fi end: 1,seq(seq(a(n,k),k=0..(n-1)/2),n=1..15);
  • Mathematica
    A091894[n_] := Prepend[Table[ CoefficientList[ 2^i (1 - z)^((2 i + 3)/2) Hypergeometric2F1[(i + 3)/2, (i + 4)/2, 2, z], z], {i, 0, n}], {1}] (* computes a table of the first n rows. Stumbled accidentally on it. Perhaps someone can find a relationship here? Thies Heidecke (theidecke(AT)astrophysik.uni-kiel.de), Sep 23 2008 *)
    Join[{1},Select[Flatten[Table[2^(n-2k-1) Binomial[n-1,2k] Binomial[2k,k]/ (k+1), {n,20},{k,0,n}]],#!=0&]] (* Harvey P. Dale, Mar 05 2012 *)
    p[n_] := 2^n Hypergeometric2F1[(1 - n)/2, -n/2, 2, x]; Flatten[Join[{{1}}, Table[CoefficientList[p[n], x], {n, 0, 12}]]] (* Peter Luschny, Jan 23 2018 *)
  • PARI
    {T(n, k) = if( n<1, n==0 && k==0, polcoeff( polcoeff( serreverse( x / (1 + 2*x*y + x^2) + x * O(x^n)), n), n-1 - 2*k))} /* Michael Somos, Sep 25 2006 */
    
  • Sage
    [1] + [[2^(n-2*k-1)*binomial(n-1,2*k)*binomial(2*k,k)/(k+1) for k in (0..floor((n-1)/2))] for n in (1..12)] # G. C. Greubel, Nov 30 2018

Formula

T(n,k) = 2^(n - 2*k - 1)*binomial(n-1,2*k)*binomial(2*k,k)/(k + 1), T(0,0) = 1, for 0 <= k <= floor((n-1)/2).
G.f.: G = G(t,z) satisfies: t*z*G^2 - (1 - 2*z + 2*t*z)*G + 1 - z + t*z = 0.
With first row zero, the o.g.f. is g(x,t) = (1 - 2*x - sqrt((1 - 2*x)^2 - 4*t*x^2)) / (2*t*x) with the inverse ginv(x,t) = x / (1 + 2*x + t*x^2), an o.g.f. for shifted A207538 and A133156 mod signs, so A134264 and A125181 can be used to interpret the polynomials of this entry. Cf. A097610. - Tom Copeland, Feb 08 2016
If we delete the first 1 from the data these are the coefficients of the polynomials p(n) = 2^n*hypergeom([(1 - n)/2, - n/2], [2], x). - Peter Luschny, Jan 23 2018

A099251 Bisection of Motzkin sums (A005043).

Original entry on oeis.org

1, 1, 3, 15, 91, 603, 4213, 30537, 227475, 1730787, 13393689, 105089229, 834086421, 6684761125, 54022715451, 439742222071, 3602118427251, 29671013856627, 245613376802185, 2042162142208813, 17047255430494497, 142816973618414817
Offset: 0

Views

Author

N. J. A. Sloane, Nov 16 2004

Keywords

Comments

The Kn4 triangle sums of A175136 lead to the sequence given above (n >= 1). For the definition of the Kn4 and other triangle sums see A180662. - Johannes W. Meijer, May 06 2011
Equals the expected value of trace(O)^(2n), where O is a 3 X 3 orthogonal matrix randomly selected according to Haar measure (see MathOverflow link). - Nathaniel Johnston, Sep 05 2014
From Petros Hadjicostas, Jul 23 2020: (Start)
In Smith (1985), we apparently have a(n) = P(2*n), where P(n) is the number of linearly independent three-dimensional n-th order isotropic tensors. In the paper, he refers to Smith (1968) for more details. It is not clear why he does not list the values of P(2*n+1). See also the 1978 letter of D. L. Andrews to N. J. A. Sloane.
Eric Weisstein gives some details on how the material in Smith (1968) about isotropic tensors is related to Motzkin sums. (End)

References

  • G. F. Smith, On isotropic tensors and rotation tensors of dimension m and order n, Tensor (N.S.), Vol. 19 (1968), 79-88 (MR0224008).

Crossrefs

Programs

  • Maple
    G := (1+x-sqrt(1-2*x-3*x^2))/(2*x*(1+x)): Gser := series(G,x=0,60):
    1, seq(coeff(Gser, x^(2*n)), n=1..25); # Emeric Deutsch
    a := n -> hypergeom([1/2, -2*n], [2], 4):
    seq(simplify(a(n)), n=0..21); # Peter Luschny, Jul 25 2020
  • Mathematica
    Take[CoefficientList[Series[(1 + x - Sqrt[1 - 2 * x - 3 * x^2])/(2 * x * (1 + x)), {x, 0, 60}], x], {1, -1, 2}] (* Vaclav Kotesovec, Oct 17 2012 *)
  • Maxima
    a(n):=sum(binomial(2*j,j)*(-1)^(j)*binomial(2*n+1,j+1),j,0,2*n+1)/(2*n+1); /*Vladimir Kruchinin, Apr 02 2017*/
  • PARI
    x='x+O('x^66); v=Vec((1+x-sqrt(1-2*x-3*x^2))/(2*x*(1+x))); vector(#v\2,n,v[2*n-1]) \\ Joerg Arndt, May 12 2013
    

Formula

Recurrence: n*(2*n + 1)*a(n) = (2*n - 1)*(13*n - 10)*a(n-1) - 3*(26*n^2 - 87*n + 76)*a(n-2) + 27*(n - 2)*(2*n - 5)*a(n-3). - Vaclav Kotesovec, Oct 17 2012
a(n) ~ 3^(2*n + 3/2)/(16*sqrt(2*Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 17 2012
Conjecture: a(n) = (2/Pi)*Integral_{t=0..1} sqrt((1 - t)/t)*(1 - 8*t + 16*t^2)^n. - Benedict W. J. Irwin, Oct 05 2016
a(n) = Sum_{j=0..2*n+1} (C(2*j,j)*(-1)^(j)*C(2*n+1,j+1))/(2*n+1). - Vladimir Kruchinin, Apr 02 2017
a(n) = hypergeom([1/2, -2*n], [2], 4). - Peter Luschny, Jul 25 2020

Extensions

More terms from Emeric Deutsch, Nov 18 2004

A162148 a(n) = n*(n+1)*(5*n+7)/6.

Original entry on oeis.org

0, 4, 17, 44, 90, 160, 259, 392, 564, 780, 1045, 1364, 1742, 2184, 2695, 3280, 3944, 4692, 5529, 6460, 7490, 8624, 9867, 11224, 12700, 14300, 16029, 17892, 19894, 22040, 24335, 26784, 29392, 32164, 35105, 38220, 41514, 44992, 48659, 52520, 56580
Offset: 0

Views

Author

Keywords

Comments

Partial sums of A147875.
Equals the fourth right hand column of A175136 for n>=1. - Johannes W. Meijer, May 06 2011
a(n) is the number of triples (w,x,y) havingt all terms in {0,...,n} and x+y>w. - Clark Kimberling, Jun 14 2012

Crossrefs

Programs

Formula

a(n) = A162147(n) + A000217(n).
From Johannes W. Meijer, May 06 2011: (Start)
G.f.: x*(4+x)/(1-x)^4.
a(n) = 4*binomial(n+2,3) + binomial(n+1,3).
a(n) = A091894(3,0)*binomial(n+2,3) + A091894(3,1)*binomial(n+1,3). (End)
a(n) = (n+1)*A000290(n+1) - Sum_{i=1..n+1} A000217(i).
a(n) = 4*a(n-1) -6*a(n-2) +4*a(n-3) -a(n-4), a(0)=0, a(1)=4, a(2)=17, a(3)=44. - Harvey P. Dale, May 20 2014
E.g.f.: x*(24 +27*x +5*x^2)*exp(x)/6. - G. C. Greubel, Mar 31 2021

Extensions

Definition rephrased by R. J. Mathar, Jun 27 2009

A166300 Number of Dyck paths of semilength n, having no ascents and no descents of length 1, and having no UUDD's starting at level 0.

Original entry on oeis.org

1, 0, 0, 1, 1, 2, 5, 10, 22, 50, 113, 260, 605, 1418, 3350, 7967, 19055, 45810, 110637, 268301, 653066, 1594980, 3907395, 9599326, 23643751, 58374972, 144442170, 358136905, 889671937, 2214015802, 5518884019, 13778312440, 34448765740
Offset: 0

Views

Author

Emeric Deutsch, Nov 07 2009

Keywords

Comments

a(n) = A166299(n,0).
a(n) is the number of peakless Motzkin paths of length n with no (1,0)-steps at level 0. Example: a(5)=2 because, denoting U=(1,1), H=(1,0), and D=(1,-1), we have UHHHD and UUHDD. - Emeric Deutsch, May 03 2011
From Paul Barry, Mar 31 2011: (Start)
The Hankel transform of a(n+3) is A188444(n+1).
a(n+3) gives the diagonal sums of the triangle A100754.
a(n+3) has g.f. 1/(1-x-x^2/(1-2x+3x^2/(1+2x+x^2/(1-2x-(1/3)x^2/(1-x-(2/3)x^2/(1-2x+(5/2)x^2/(1+2x+(3/2)x^2/(1-...)))))))) (continued fraction) where the coefficients of x^2 have denominators A188442 and numerators A188443. (End)
The Ca1 triangle sums of triangle A175136 lead to this sequence (n>=3). For the definitions of the Ca1 and other triangle sums see A180662. - Johannes W. Meijer, May 06 2011
a(n) is the number of closed Deutsch paths of n steps with all peaks at even height. A Deutsch path is a lattice path of up-steps (1,1) and down-steps (1,-j), j>=1, starting at the origin that never goes below the x-axis, and it is closed if it ends on the x-axis. For example a(5) = 2 counts UUUU4, UU1U2, where U denotes an up-step and a down-step is denoted by its length, and a(6) = 5 counts UUUU13, UUUU22, UUUU31, UU1U11, UU2UU2. - David Callan, Dec 08 2021

Examples

			a(5)=2 because we have UUUDDUUDDD and UUUUUDDDDD.
G.f. = 1 + x^3 + x^4 + 2*x^5 + 5*x^6 + 10*x^7 + 22*x^8 + 50*x^9 + 113*x^10 + ...
		

Crossrefs

Programs

  • Magma
    m:=40; R:=PowerSeriesRing(Rationals(), m); Coefficients(R!(2/(1+x+x^2+Sqrt((1+x+x^2)*(1-3*x+x^2))))); // G. C. Greubel, Sep 22 2018
  • Maple
    G := 2/(1+z+z^2+sqrt((1+z+z^2)*(1-3*z+z^2))): Gser := series(G, z = 0, 35): seq(coeff(Gser, z, n), n = 0 .. 32);
  • Mathematica
    CoefficientList[Series[2/(1+x+x^2+Sqrt[(1+x+x^2)*(1-3*x+x^2)]), {x, 0, 20}], x] (* Vaclav Kotesovec, Feb 12 2014 *)
  • Maxima
    a(n):=sum((j+1)*sum((binomial(j+2*i+1,i)*sum(binomial(k,n-k-j-2*i)*binomial(k+j+2*i,k)*(-1)^(n-k),k,0,n-j-2*i))/(j+2*i+1),i,0,n-j),j,0,n); /*  Vladimir Kruchinin, Mar 07 2016 */
    
  • PARI
    {a(n) = local(A); A = 1 + O(x); for( k=1, ceil(n / 5), A = 1 / (1 - x^3 / (1 - x / (1 - x * A)))); polcoeff( A, n)}; /* Michael Somos, May 12 2012 */
    
  • PARI
    x='x+O('x^40); Vec(2/(1+x+x^2+((1+x+x^2)*(1-3*x+x^2))^(1/2))) \\ Altug Alkan, Sep 23 2018
    

Formula

G.f. = G(z)=2/(1 + z + z^2 + sqrt((1 + z + z^2)*(1 - 3*z + z^2))).
G.f.: 1 / (1 - x^3 / (1 - x / (1 - x / (1 - x^3 / (1 - x / (1 - x / ...)))))). - Michael Somos, May 12 2012
G.f. A(x) satisfies A(x) = 1 / (1 - x^3 / (1 - x / (1 - x *A(x)))). - Michael Somos, May 12 2012
Conjecture: (n+1)*a(n) +2*(-n+1)*a(n-1) +(-n+1)*a(n-2) +2*(-n+1)*a(n-3) +(n-3)*a(n-4)=0. - R. J. Mathar, Nov 24 2012
a(n) ~ (3+sqrt(5))^(n+2) * sqrt(7*sqrt(5)-15) / (2 * sqrt(Pi) * n^(3/2) * 2^(n+9/2)). - Vaclav Kotesovec, Feb 12 2014. Equivalently, a(n) ~ 5^(1/4) * phi^(2*n + 2) / (8 * sqrt(Pi) * n^(3/2)), where phi = A001622 is the golden ratio. - Vaclav Kotesovec, Dec 08 2021
a(n) = Sum_{j=0..n}((j+1)*Sum_{i=0..n-j}((binomial(j+2*i+1,i)*Sum_{k=0..n-j-2*i}(binomial(k,n-k-j-2*i)*binomial(k+j+2*i,k)*(-1)^(n-k)))/(j+2*i+1))). - Vladimir Kruchinin, Mar 07 2016

A190048 Expansion of (8+6*x)/(1-x)^5.

Original entry on oeis.org

8, 46, 150, 370, 770, 1428, 2436, 3900, 5940, 8690, 12298, 16926, 22750, 29960, 38760, 49368, 62016, 76950, 94430, 114730, 138138, 164956, 195500, 230100, 269100, 312858, 361746, 416150, 476470, 543120, 616528, 697136, 785400, 881790, 986790, 1100898
Offset: 0

Views

Author

Johannes W. Meijer, May 06 2011

Keywords

Comments

Equals the fifth right hand column of A175136.

Crossrefs

Related to A000332 and A091894.

Programs

  • Magma
    [(7*n^4+58*n^3+173*n^2+218*n+96)/12: n in [0..50]]; // Vincenzo Librandi, May 07 2011
    
  • Maple
    A190048 := proc(n) option remember; a(n):=(7*n^4+58*n^3+173*n^2+218*n+96)/12 end: seq(A190048(n),n=0..35);
  • Mathematica
    LinearRecurrence[{5,-10,10,-5,1}, {8,46,150,370,770}, 30] (* or *) CoefficientList[Series[(8+6*x)/(1-x)^5, {x, 0, 50}], x] (* G. C. Greubel, Jan 10 2018 *)
  • PARI
    x='x+O('x^30); Vec((8+6*x)/(1-x)^5) \\ G. C. Greubel, Jan 10 2018
    
  • PARI
    for(n=0,50, print1((7*n^4 +58*n^3 +173*n^2 +218*n +96)/12, ", ")) \\ G. C. Greubel, Jan 10 2018

Formula

G.f.: (8+6*x)/(1-x)^5.
a(n) = 8*binomial(n+4,4) + 6*binomial(n+3,4).
a(n) = A091894(4,0)*binomial(n+4,4) + A091894(4,1)*binomial(n+3,4).
a(n) = (7*n^4 +58*n^3 +173*n^2 +218*n +96)/12.

A190049 Expansion of (16+24*x+2*x^2)/(x-1)^6.

Original entry on oeis.org

16, 120, 482, 1412, 3402, 7168, 13692, 24264, 40524, 64504, 98670, 145964, 209846, 294336, 404056, 544272, 720936, 940728, 1211098, 1540308, 1937474, 2412608, 2976660, 3641560, 4420260, 5326776, 6376230, 7584892
Offset: 0

Views

Author

Johannes W. Meijer, May 06 2011

Keywords

Comments

Equals the sixth right hand column of A175136.

Crossrefs

Related to A000389 and A091894.

Programs

  • Magma
    [(21*n^5+245*n^4+1105*n^3+2395*n^2+2474*n+960)/60: n in [0..50]]; // Vincenzo Librandi, May 07 2011
    
  • Maple
    A190049 := proc(n) option remember; a(n):=(21*n^5 +245*n^4 +1105*n^3 +2395*n^2 +2474*n +960)/60 end: seq(A190049(n),n=0..27);
  • Mathematica
    LinearRecurrence[{6,-15,20,-15,6,-1}, {16,120,482,1412,3402,7168}, 30] (* or *) CoefficientList[Series[(16 +24*x +2*x^2)/(1-x)^6, {x, 0, 50}], x] (* G. C. Greubel, Jan 10 2018 *)
  • PARI
    x='x+O('x^30); Vec((16 +24*x +2*x^2)/(1-x)^6) \\ G. C. Greubel, Jan 10 2018
    
  • PARI
    for(n=0,30, print1((21*n^5 +245*n^4 +1105*n^3 +2395*n^2 +2474*n +960)/60, ", ")) \\ G. C. Greubel, Jan 10 2018

Formula

G.f.: (16 +24*x +2*x^2)/(1-x)^6.
a(n) = 16*binomial(n+5,5) +24*binomial(n+4,5) +2*binomial(n+3,5).
a(n) = A091894(5,0)*binomial(n+5,5) + A091894(5,1)*binomial(n+4,5) + A091894(5,2)*binomial(n+3,5).
a(n) = (21*n^5 +245*n^4 +1105*n^3 +2395*n^2 +2474*n +960)/60.

A190050 Expansion of ((1-x)*(3*x^2-3*x+1))/(1-2*x)^3.

Original entry on oeis.org

1, 2, 6, 17, 46, 120, 304, 752, 1824, 4352, 10240, 23808, 54784, 124928, 282624, 634880, 1417216, 3145728, 6946816, 15269888, 33423360, 72876032, 158334976, 342884352, 740294656, 1593835520, 3422552064
Offset: 0

Views

Author

Johannes W. Meijer, May 06 2011

Keywords

Comments

The second left hand column of triangle A175136.

Crossrefs

Related to A001788.

Programs

  • Magma
    [1] cat [(n^2 + 5*n + 10)*2^(n-4): n in [1..30]]; // G. C. Greubel, Jan 10 2018
  • Maple
    A190050:= proc(n) option remember; if n=0 then A190050(n):=1: else A190050(n):=(n^2+5*n+10)*2^(n-4) fi: end: seq (A190050(n), n=0..26);
  • Mathematica
    Join[{1}, LinearRecurrence[{6,-12,8}, {2,6,17}, 30]] (* or *) CoefficientList[Series[((1-x)*(3*x^2-3*x+1))/(1-2*x)^3, {x, 0, 50}], x] (* G. C. Greubel, Jan 10 2018 *)
  • PARI
    x='x+O('x^30); Vec(((1-x)*(3*x^2-3*x+1))/(1-2*x)^3) \\ G. C. Greubel, Jan 10 2018
    
  • PARI
    for(n=0,30, print1(if(n==0,1,(n^2 + 5*n + 10)*2^(n-4)), ", ")) \\ G. C. Greubel, Jan 10 2018
    

Formula

G.f.: ((1-x)*(3*x^2-3*x+1))/(1-2*x)^3.
a(n) = (n^2 + 5*n + 10)*2^(n-4) for n >=1 with a(0)=1.
a(n) = A001788(n+1) -4*A001788(n) +6*A001788(n-1) -3*A001788(n-2) for n >=1 with a(0)=1.

A190051 Expansion of (1-x)*(10*x^4-20*x^3+16*x^2-6*x+1)/(1-2*x)^5.

Original entry on oeis.org

1, 3, 12, 44, 150, 482, 1476, 4344, 12368, 34240, 92544, 244992, 636928, 1629696, 4111360, 10242048, 25227264, 61505536, 148570112, 355860480, 845807616, 1996095488, 4680056832, 10906763264, 25275924480, 58271465472
Offset: 0

Views

Author

Johannes W. Meijer, May 06 2011

Keywords

Comments

The third left hand column of triangle A175136.

Crossrefs

Related to A003472.

Programs

  • Magma
    [1] cat [(264 + 214*n + 14*n^3 + 83*n^2 + n^4)*2^(n-7)/3: n in [1..30]]; // G. C. Greubel, Jan 10 2018
  • Maple
    A190051:= proc(n) option remember; if n=0 then A190051(n):=1 else A190051(n):= (264+214*n+14*n^3+83*n^2+n^4)*2^(n-7)/3 fi: end: seq (A190051(n), n=0..25);
  • Mathematica
    Join[{1}, LinearRecurrence[{10,-40,80,-80,32}, {3,12,44,150,482}, 30]] (* or *) CoefficientList[Series[(1 - x)*(10*x^4 -20*x^3 +16*x^2 -6*x + 1)/(1 -2*x)^5, {x, 0, 50}], x] (* G. C. Greubel, Jan 10 2018 *)
  • PARI
    x='x+O('x^30); Vec((1-x)*(10*x^4-20*x^3+16*x^2-6*x+1)/(1-2*x)^5) \\ G. C. Greubel, Jan 10 2018
    
  • PARI
    for(n=0,30, print1(if(n==0,1,(264 + 214*n + 14*n^3 + 83*n^2 + n^4)*2^(n-7)/3), ", ")) \\ G. C. Greubel, Jan 10 2018
    

Formula

G.f.: (1-x)*(10*x^4-20*x^3+16*x^2-6*x+1)/(1-2*x)^5.
a(n) = (264 + 214*n + 14*n^3 + 83*n^2 + n^4)*2^(n-7)/3 for n >=1 with a(0)=1.
a(n-4) = A003472(n) -7*A003472(n-1) +22*A003472(n-2) -36*A003472(n-3) +30*A003472(n-4) -10*A003472(n-5) for n>=5 with a(0) = 1.
Showing 1-10 of 12 results. Next