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-7 of 7 results.

A053221 Row sums of triangle A053218.

Original entry on oeis.org

1, 5, 16, 43, 106, 249, 568, 1271, 2806, 6133, 13300, 28659, 61426, 131057, 278512, 589807, 1245166, 2621421, 5505004, 11534315, 24117226, 50331625, 104857576, 218103783, 452984806, 939524069, 1946157028, 4026531811, 8321499106
Offset: 1

Views

Author

Asher Auel, Jan 01 2000

Keywords

Comments

Considered as a vector, the sequence = A074909 * [1, 2, 3, ...], where A074909 is the beheaded Pascal's triangle as a matrix. - Gary W. Adamson, Mar 06 2012
a(n) is the sum of the upper left n X n subarray of A052509 (viewed as an infinite square array). For example (1+1+1) + (1+2+2) + (1+3+4) = 16. - J. M. Bergot, Nov 06 2012
Number of ternary strings of length n that contain at least one 2 and at most one 0. For example, a(3) = 16 since the strings are the 6 permutations of 201, the 3 permutations of 211, the 3 permutations of 220, the 3 permutations of 221, and 222. - Enrique Navarrete, Jul 25 2021

Examples

			a(4) = 4 + 7 + 12 + 20 = 43.
		

Crossrefs

Programs

  • Magma
    [(n+2)*2^(n-1)-n-1: n in [1..50]]; // G. C. Greubel, Sep 03 2018
  • Maple
    A053221 := proc(n) (n+2)*2^(n-1)-n-1 ; end proc: # R. J. Mathar, Sep 02 2011
  • Mathematica
    Table[(n + 2)*2^(n - 1) - n - 1, {n, 29}] (* or *)
    Rest@ CoefficientList[Series[-x (-1 + x + x^2)/((2 x - 1)^2*(x - 1)^2), {x, 0, 29}], x] (* Michael De Vlieger, Sep 22 2017 *)
    LinearRecurrence[{6,-13,12,-4},{1,5,16,43},30] (* Harvey P. Dale, Jun 28 2021 *)
  • PARI
    vector(50,n, (n+2)*2^(n-1)-n-1) \\ G. C. Greubel, Sep 03 2018
    

Formula

a(n) = (n+2)*2^(n-1)-n-1. - Vladeta Jovovic, Feb 28 2003
G.f.: -x*(-1+x+x^2) / ( (2*x-1)^2*(x-1)^2 ). - R. J. Mathar, Sep 02 2011
a(n) = (1/2) * Sum_{k=1..n} Sum_{i=1..n} C(k,i) + C(n,k). - Wesley Ivan Hurt, Sep 22 2017
E.g.f.: exp(x)*(exp(x)-1)*(1+x). - Enrique Navarrete, Jul 25 2021
a(n+1) = 2*a(n) + A006127(n). - Ya-Ping Lu, Jan 01 2024

A053219 Reverse of triangle A053218, read by rows.

Original entry on oeis.org

1, 3, 2, 8, 5, 3, 20, 12, 7, 4, 48, 28, 16, 9, 5, 112, 64, 36, 20, 11, 6, 256, 144, 80, 44, 24, 13, 7, 576, 320, 176, 96, 52, 28, 15, 8, 1280, 704, 384, 208, 112, 60, 32, 17, 9, 2816, 1536, 832, 448, 240, 128, 68, 36, 19, 10, 6144, 3328, 1792, 960, 512, 272, 144, 76, 40
Offset: 1

Views

Author

Asher Auel, Jan 01 2000

Keywords

Comments

First element in each row gives A001792. Difference between center element of row 2n-1 and row sum of row n (A053220(n+4) - A053221(n+4)) gives A045618(n).
Subtriangle of triangle in A062111. - Philippe Deléham, Nov 21 2011
Can be seen as the transform of 1, 2, 3, 4, 5, ... by a variant of the boustrophedon algorithm (see the Sage implementation). - Peter Luschny, Oct 30 2014

Examples

			Triangle begins:
1
3, 2
8, 5, 3
20, 12, 7, 4
48, 28, 16, 9, 5 ...
		

Crossrefs

Cf. A053218 (reverse of this triangle), A053220 (center elements), A053221 (row sums), A001792, A045618, A062111.

Programs

  • Mathematica
    Map[Reverse,NestList[FoldList[Plus,#[[1]]+1,#]&,{1},10]]//Grid (* Geoffrey Critzer, Jun 27 2013 *)
  • Sage
    def u():
        for n in PositiveIntegers():
            yield n
    def bous_variant(f):
        k = 0
        am = next(f)
        a = [am]
        while True:
            yield list(a)
            am = next(f)
            a.append(am)
            for m in range(k,-1,-1):
                am += a[m]
                a[m] = am
            k += 1
    b = bous_variant(u())
    [next(b) for  in range(8)] # _Peter Luschny, Oct 30 2014

A001792 a(n) = (n+2)*2^(n-1).

Original entry on oeis.org

1, 3, 8, 20, 48, 112, 256, 576, 1280, 2816, 6144, 13312, 28672, 61440, 131072, 278528, 589824, 1245184, 2621440, 5505024, 11534336, 24117248, 50331648, 104857600, 218103808, 452984832, 939524096, 1946157056, 4026531840, 8321499136, 17179869184, 35433480192
Offset: 0

Views

Author

Keywords

Comments

Number of parts in all compositions (ordered partitions) of n + 1. For example, a(2) = 8 because in 3 = 2 + 1 = 1 + 2 = 1 + 1 + 1 we have 8 parts. Also number of compositions (ordered partitions) of 2n + 1 with exactly 1 odd part. For example, a(2) = 8 because the only compositions of 5 with exactly 1 odd part are 5 = 1 + 4 = 2 + 3 = 3 + 2 = 4 + 1 = 1 + 2 + 2 = 2 + 1 + 2 = 2 + 2 + 1. - Emeric Deutsch, May 10 2001
Binomial transform of natural numbers [1, 2, 3, 4, ...].
For n >= 1 a(n) is also the determinant of the n X n matrix with 3's on the diagonal and 1's elsewhere. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 06 2001
The arithmetic mean of first n terms of the sequence is 2^(n-1). - Amarnath Murthy, Dec 25 2001, corrected by M. F. Hasler, Dec 17 2016
Also the number of "winning paths" of length n across an n X n Hex board. Satisfies the recursion a(n) = 2a(n-1) + 2^(n-2). - David Molnar (molnar(AT)stolaf.edu), Apr 10 2002
Diagonal in A053218. - Benoit Cloitre, May 08 2002
Let M_n be the n X n matrix m_(i, j) = 1 + abs(i-j) then det(M_n) = (-1)^(n-1)*a(n-1). - Benoit Cloitre, May 28 2002
Absolute value of determinant of n X n matrix of form: [1 2 3 4 5 / 2 1 2 3 4 / 3 2 1 2 3 / 4 3 2 1 2 / 5 4 3 2 1]. - Benoit Cloitre, Jun 20 2002
Number of ones in all (n+1)-bit integers (cf. A000120). - Ralf Stephan, Aug 02 2003
This sequence also emerges as a floretion force transform of powers of 2 (see program code). Define a(-1) = 0 (as the sequence is returned by FAMP). Then a(n-1) + A098156(n+1) = 2*a(n) (conjecture). - Creighton Dement, Mar 14 2005
This sequence gives the absolute value of the determinant of the Toeplitz matrix with first row containing the first n integers. - Paul Max Payton, May 23 2006
Equals sums of rows right of left edge of A102363 divided by three, + 2^K. - David G. Williams (davidwilliams(AT)paxway.com), Oct 08 2007
If X_1, X_2, ..., X_n are 2-blocks of a (2n+1)-set X then, for n >= 1, a(n) is the number of (n+1)-subsets of X intersecting each X_i, (i = 1, 2, ..., n). - Milan Janjic, Nov 18 2007
Also, a(n-1) is the determinant of the n X n matrix with A[i, j] = n - |i-j|. - M. F. Hasler, Dec 17 2008
1/2 the number of permutations of 1 .. (n+2) arranged in a circle with exactly one local maximum. - R. H. Hardin, Apr 19 2009
The first corrector line for transforming 2^n offset 0 with a leading 1 into the Fibonacci sequence. - Al Hakanson (hawkuu(AT)gmail.com), Jun 01 2009
a(n) is the number of runs of consecutive 1's in all binary sequences of length (n+1). - Geoffrey Critzer, Jul 02 2009
Let X_j (0 < j <= 2^n) all the subsets of N_n; m(i, j) := if {i} in X_j then 1 else 0. Let A = transpose(M).M; Then a(i, j) = (number of elements)|X_i intersect X_j|. Determinant(X*I-A) = (X-(n+1)*2^(n-2))*(X-2^(n-2))^(n-1)*X^(2^n-n).
Eigenvector for (n+1)*2^(n-2) is V_i=|X_i|.
Sum_{k=1..2^n} |X_i intersect X_k|*|X_k| = (n+1)*2^(n-2)*|X_i|.
Eigenvectors for 2^(n-2) are {line(M)[i] - line(M)[j], 1 <= i, j <= n}. - CLARISSE Philippe (clarissephilippe(AT)yahoo.fr), Mar 24 2010
The sequence b(n) = 2*A001792(n), for n >= 1 with b(0) = 1, is an elephant sequence, see A175655. For the central square four A[5] vectors, with decimal values 187, 190, 250 and 442, lead to the b(n) sequence. For the corner squares these vectors lead to the companion sequence A134401. - Johannes W. Meijer, Aug 15 2010
Equals partial sums of A045623: (1, 2, 5, 12, 28, ...); where A045623 = the convolution square of (1, 1, 2, 4, 8, 16, 32, ...). - Gary W. Adamson, Oct 26 2010
Equals (1, 2, 4, 8, 16, ...) convolved with (1, 1, 2, 4, 8, 16, ...); e.g., a(3) = 20 = (1, 1, 2, 4) dot (8, 4, 2, 1) = (8 + 4 + 4 + 4). - Gary W. Adamson, Oct 26 2010
This sequence seems to give the first x+1 nonzero terms in the sequence derived by subtracting the m-th term in the x_binacci sequence (where the first term is one and the y-th term is the sum of x terms immediately preceding it) from 2^(m-2). - Dylan Hamilton, Nov 03 2010
Recursive formulas for a(n) are in many cases derivable from its property wherein delta^k(a(n)) - a(n) = k*2^n where delta^k(a(n)) represents the k-th forward difference of a(n). Provable with a difference table and a little induction. - Ethan Beihl, May 02 2011
Let f(n,k) be the sum of numbers in the subsets of size k of {1, 2, ..., n}. Then a(n-1) is the average of the numbers f(n, 0), ... f(n, n). Example: (f(3, 1), f(3, 2), f(3, 3)) = (6, 12, 6), with average (6+12+6)/3. - Clark Kimberling, Feb 24 2012
a(n) is the number of length-2n binary sequences that contain a subsequence of ones with length n or more. To derive this result, note that there are 2^n sequences where the initial one of the subsequence occurs at entry one. If the initial one of the subsequence occurs at entry 2, 3, ..., or n + 1, there are 2^(n-1) sequences since a zero must precede the initial one. Hence a(n) = 2^n + n*2^(n-1)=(n+2)2^(n-1). An example is given in the example section below. - Dennis P. Walsh, Oct 25 2012
As the total number of parts in all compositions of n+1 (see the first line in Comments) the equivalent sequence for partitions is A006128. On the other hand, as the first differences of A001787 (see the first line in Crossrefs) the equivalent sequence for partitions is A138879. - Omar E. Pol, Aug 28 2013
a(n) is the number of spanning trees of the complete tripartite graph K_{n,1,1}. - James Mahoney, Oct 24 2013
a(n-1) = denominator of the mean (2n/(n+1), after reduction), of the compositions of n; numerator is given by A022998(n). - Clark Kimberling, Mar 11 2014
From Tom Copeland, Nov 09 2014: (Start)
The shifted array belongs to an interpolated family of arrays associated to the Catalan A000108 (t=1), and Riordan, or Motzkin sums A005043 (t=0), with the interpolating o.g.f. (1-sqrt(1-4x/(1+(1-t)x)))/2 and inverse x(1-x)/(1+(t-1)x(1-x)). See A091867 for more info on this family. Here the interpolation is t=-3 (mod signs in the results).
Let C(x) = (1 - sqrt(1-4x))/2, an o.g.f. for the Catalan numbers A000108, with inverse Cinv(x) = x*(1-x) and P(x,t) = x/(1+t*x) with inverse P(x,-t).
Shifted o.g.f: G(x) = x*(1-x)/(1 - 4x*(1-x)) = P[Cinv(x),-4].
Inverse o.g.f: Ginv(x) = [1 - sqrt(1 - 4*x/(1+4x))]/2 = C[P(x, 4)] (signed shifted A001700). Cf. A030528. (End)
For n > 0, element a(n) of the sequence is equal to the gradients of the (n-1)-th row of Pascal triangle multiplied with the square of the integers from n+1,...,1. I.e., row 3 of Pascal's triangle 1,3,3,1 has gradients 1,2,0,-2,-1, so a(4) = 1*(5^2) + 2*(4^2) + 0*(3^2) - 2*(2^2) - 1*(1^2) = 48. - Jens Martin Carlsson, May 18 2017
Number of self-avoiding paths connecting all the vertices of a convex (n+2)-gon. - Ivaylo Kortezov, Jan 19 2020
a(n-1) is the total number of elements of subsets of {1,2,..,n} that contain n. For example, for n = 3, a(2) = 8, and the subsets of {1,2,3} that contain 3 are {3}, {1,3}, {2,3}, {1,2,3}, with a total of 8 elements. - Enrique Navarrete, Aug 01 2020

Examples

			a(0) = 1, a(1) = 2*1 + 1 = 3, a(2) = 2*3 + 2 = 8, a(3) = 2*8 + 4 = 20, a(4) = 2*20 + 8 = 48, a(5) = 2*48 + 16 = 112, a(6) = 2*112 + 32 = 256, ... - _Philippe Deléham_, Apr 19 2009
a(2) = 8 since there are 8 length-4 binary sequences with a subsequence of ones of length 2 or more, namely, 1111, 1110, 1101, 1011, 0111, 1100, 0110, and 0011. - _Dennis P. Walsh_, Oct 25 2012
G.f. = 1 + 3*x + 8*x^2 + 20*x^3 + 48*x^4 + 112*x^5 + 256*x^6 + 576*x^7 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • 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).
  • A. M. Stepin and A. T. Tagi-Zade, Words with restrictions, pp. 67-74 of Kvant Selecta: Combinatorics I, Amer. Math. Soc., 2001 (G_n on p. 70).

Crossrefs

First differences of A001787.
a(n) = A049600(n, 1), a(n) = A030523(n + 1, 1).
Cf. A053113.
Row sums of triangles A008949 and A055248.
a(n) = -A039991(n+2, 2).
If the exponent E in a(n) = Sum_{m=0..n} (Sum_{k=0..m} C(n,k))^E is 1, 2, 3, 4, 5 we get A001792, A003583, A007403, A294435, A294436 respectively.

Programs

  • GAP
    List([0..35],n->(n+2)*2^(n-1)); # Muniru A Asiru, Sep 25 2018
    
  • Haskell
    a001792 n = a001792_list !! n
    a001792_list = scanl1 (+) a045623_list
    -- Reinhard Zumkeller, Jul 21 2013
    
  • Magma
    [(n+2)*2^(n-1): n in [0..40]]; // Vincenzo Librandi, Nov 10 2014
    
  • Maple
    A001792 := n-> (n+2)*2^(n-1);
    spec := [S, {B=Set(Z, 0 <= card), S=Prod(Z, B, B)}, labeled]: seq(combstruct[count](spec, size=n)/4, n=2..30); # Zerinvary Lajos, Oct 09 2006
    A001792:=-(-3+4*z)/(2*z-1)^2; # Simon Plouffe in his 1992 dissertation, which gives the sequence without the initial 1
    G(x):=1/exp(2*x)*(1-x): f[0]:=G(x): for n from 1 to 54 do f[n]:=diff(f[n-1],x) od: x:=0: seq(abs(f[n]),n=0..28 ); # Zerinvary Lajos, Apr 17 2009
    a := n -> hypergeom([-n, 2], [1], -1);
    seq(round(evalf(a(n),32)), n=0..31); # Peter Luschny, Aug 02 2014
  • Mathematica
    matrix[n_Integer /; n >= 1] := Table[Abs[p - q] + 1, {q, n}, {p, n}]; a[n_Integer /; n >= 1] := Abs[Det[matrix[n]]] (* Josh Locker (joshlocker(AT)macfora.com), Apr 29 2004 *)
    g[n_,m_,r_] := Binomial[n - 1, r - 1] Binomial[m + 1, r] r; Table[1 + Sum[g[n, k - n, r], {r, 1, k}, {n, 1, k - 1}], {k, 1, 29}] (* Geoffrey Critzer, Jul 02 2009 *)
    a[n_] := (n + 2)*2^(n - 1); a[Range[0, 40]] (* Vladimir Joseph Stephan Orlovsky, Feb 09 2011 *)
    LinearRecurrence[{4, -4}, {1, 3}, 40] (* Harvey P. Dale, Aug 29 2011 *)
    CoefficientList[Series[(1 - x) / (1 - 2 x)^2, {x, 0, 40}], x] (* Vincenzo Librandi, Nov 10 2014 *)
    b[i_]:=i; a[n_]:=Abs[Det[ToeplitzMatrix[Array[b, n], Array[b, n]]]]; Array[a, 40] (* Stefano Spezia, Sep 25 2018 *)
    a[n_]:=Hypergeometric2F1[2,-n+1,1,-1];Array[a,32] (* Giorgos Kalogeropoulos, Jan 04 2022 *)
  • PARI
    A001792(n)=(n+2)<<(n-1) \\ M. F. Hasler, Dec 17 2008
    
  • Python
    for n in range(0,40): print(int((n+2)*2**(n-1)), end=' ') # Stefano Spezia, Oct 16 2018

Formula

a(n) = (n+2)*2^(n-1).
G.f.: (1 - x)/(1 - 2*x)^2 = 2F1(1,3;2;2x).
a(n) = 4*a(n-1) - 4*a(n-2).
G.f. (-1 + (1-2*x)^(-2))/(x*2^2). - Wolfdieter Lang
a(n) = A018804(2^n). - Matthew Vandermast, Mar 01 2003
a(n) = Sum_{k=0..n+2} binomial(n+2, 2k)*k. - Paul Barry, Mar 06 2003
a(n) = (1/4)*A001787(n+2). - Emeric Deutsch, May 24 2003
With a leading 0, this is ((n+1)2^n - 0^n)/4 = Sum_{m=0..n} binomial(n - 1, m - 1)*m, the binomial transform of A004526(n+1). - Paul Barry, Jun 05 2003
a(n) = Sum_{k=0..n} binomial(n, k)*(k + 1). - Lekraj Beedassy, Jun 24 2004
a(n) = A000244(n) - A066810(n). - Ross La Haye, Apr 29 2006
Row sums of triangle A130585. - Gary W. Adamson, Jun 06 2007
Equals A125092 * [1/1, 1/2, 1/3, ...]. - Gary W. Adamson, Nov 16 2007
a(n) = (n+1)*2^n - n*2^(n-1). Equals A128064 * A000079. - Gary W. Adamson, Dec 28 2007
G.f.: F(3, 1; 2; 2x). - Paul Barry, Sep 03 2008
a(n) = 1 + Sum_{k=1..n} (n - k + 4)2^(n - k - 1). This follows from the result that the number of parts equal to k in all compositions of n is (n - k + 3)2^(n - k - 2) for 0 < k < n. - Geoffrey Critzer, Sep 21 2008
a(n) = 2^(n-1) + 2 a(n-1) ; a(n-1) = det(n - |i - j|){i, j = 1..n}. - _M. F. Hasler, Dec 17 2008
a(n) = 2*a(n-1) + 2^(n-1). - Philippe Deléham, Apr 19 2009
a(n) = A164910(2^n). - Gary W. Adamson, Aug 30 2009
a(n) = Sum_{i=1..2^n} gcd(i, 2^n) = A018804(2^n). So we have: 2^0 * phi(2^n) + ... + 2^n * phi(2^0) = (n + 2)*2^(n-1), where phi is the Euler totient function. - Jeffrey R. Goodwin, Nov 11 2011
a(n) = Sum_{j=0..n} Sum_{i=0..n} binomial(n, i + j). - Yalcin Aktar, Jan 17 2012
Eigensequence of an infinite lower triangular matrix with 2^n as the left border and the rest 1's. - Gary W. Adamson, Jan 30 2012
G.f.: 1 + 2*x*U(0) where U(k) = 1 + (k + 1)/(2 - 8*x/(4*x + (k + 1)/U(k + 1))); (continued fraction, 3 - step). - Sergei N. Gladkovskii, Oct 19 2012
a(n) = Sum_{k=0..n} Sum_{j=0..k} binomial(n,j). - Peter Luschny, Dec 03 2013
a(n) = Hyper2F1([-n, 2], [1], -1). - Peter Luschny, Aug 02 2014
G.f.: 1 / (1 - 3*x / (1 + x / (3 - 4*x))). - Michael Somos, Aug 26 2015
a(n) = -A053120(2+n, n), n >= 0, the negative of the third (sub)diagonal of the triangle of Chebyshev's T polynomials. - Wolfdieter Lang, Nov 26 2019
From Amiram Eldar, Jan 12 2021: (Start)
Sum_{n>=0} 1/a(n) = 8*log(2) - 4.
Sum_{n>=0} (-1)^n/a(n) = 4 - 8*log(3/2). (End)
E.g.f.: exp(2*x)*(1 + x). - Stefano Spezia, Jun 11 2021

A053220 a(n) = (3*n-1) * 2^(n-2).

Original entry on oeis.org

1, 5, 16, 44, 112, 272, 640, 1472, 3328, 7424, 16384, 35840, 77824, 167936, 360448, 770048, 1638400, 3473408, 7340032, 15466496, 32505856, 68157440, 142606336, 297795584, 620756992, 1291845632, 2684354560, 5570035712, 11542724608, 23890755584, 49392123904
Offset: 1

Views

Author

Asher Auel, Jan 01 2000

Keywords

Comments

Coefficients in the hypergeometric series identity 1 - 5*x/(x + 4) + 16*x*(x - 1)/((x + 4)*(x + 6)) - 44*x*(x - 1)*(x - 2)/((x + 4)*(x + 6)*(x + 8)) + ... = 0, valid in the half-plane Re(x) > 0. Cf. A276289. - Peter Bala, May 30 2019
For n>=2, a(n) is the total number of ones in runs of ones of length >=5 over all binary strings of length n+3. - Félix Balado, Aug 06 2025

Crossrefs

Center elements from triangle A053218. Also a diagonal of triangle A056242.

Programs

  • Haskell
    a053220 n = a056242 (n + 1) n  -- Reinhard Zumkeller, May 08 2014
  • Magma
    [(3*n-1)*2^(n-2): n in [1..50]]; // Vincenzo Librandi, May 09 2011
    
  • Mathematica
    ListCorrelate[{1, 1}, Table[n 2^(n - 1), {n, 0, 28}]] (* or *) ListConvolve[{1, 1}, Table[n 2^(n - 1), {n, 0, 28}]] (* Ross La Haye, Feb 24 2007 *)
    LinearRecurrence[{4, -4}, {1, 5}, 35] (* Vladimir Joseph Stephan Orlovsky, Jan 29 2012 *)
    Array[(3# - 1) 2^(# - 2) &, 35] (* Alonso del Arte, Sep 04 2018 *)
    CoefficientList[Series[(1 + x)/(1 - 2 * x)^2, {x, 0, 50}], x] (* Stefano Spezia, Sep 04 2018 *)
  • PARI
    a(n)=if(n<1,0,(3*n-1)*2^(n-2))
    
  • PARI
    a(n)=(3*n-1)<<(n-2) \\ Charles R Greathouse IV, Apr 17 2012
    

Formula

G.f.: x*(1+x)/(1-2*x)^2.
a(n) = (3*n-1) * 2^(n-2).
E.g.f.: exp(2*x)*(1+3*x). The sequence 0, 1, 5, 16, ... has a(n) = ((3n-1)*2^n + 0^n)/4 (offset 0). It is the binomial transform of A032766. The sequence 1, 5, 16, ... has a(n) = (2+3n)*2^(n-1) (offset 0). It is the binomial transform of A016777. - Paul Barry, Jul 23 2003
Row sums of A132776(n-1). - Gary W. Adamson, Aug 29 2007

A064319 Triangle with a(n,1) = n and a(n,k) = a(n,k-1) * a(n-1,k-1).

Original entry on oeis.org

1, 2, 2, 3, 6, 12, 4, 12, 72, 864, 5, 20, 240, 17280, 14929920, 6, 30, 600, 144000, 2488320000, 37150418534400000, 7, 42, 1260, 756000, 108864000000, 270888468480000000000, 10063619980174622195712000000000000000, 8, 56
Offset: 1

Views

Author

Henry Bottomley, Sep 10 2001

Keywords

Examples

			Rows start
  1;
  2,  2;
  3,  6, 12;
  4, 12, 72, 864;
  ...
		

Crossrefs

Columns include A000027, A002378, A047928, A064321. Right-hand side is A064320. Cf. A053218 which uses addition rather than multiplication to produce binomial transform.

Extensions

Name corrected by Sean A. Irvine, Jun 30 2023

A003469 Number of minimal covers of an (n + 1)-set by a collection of n nonempty subsets, a(n) = A035348(n,n-1).

Original entry on oeis.org

1, 6, 22, 65, 171, 420, 988, 2259, 5065, 11198, 24498, 53157, 114583, 245640, 524152, 1113959, 2359125, 4980546, 10485550, 22019865, 46137091, 96468716, 201326292, 419430075, 872414881, 1811938950, 3758095978, 7784627789, 16106126895, 33285996048
Offset: 1

Views

Author

Keywords

Comments

A cover of a set S is a collection of nonempty subsets of S whose union is S. A cover of S is called minimal if none of its proper subsets covers S. [from the Hearne/Wagner reference]
Partial sums of A053221.
Construct an inverted triangle table with n rows as follows: the first row are numbers from 1 to n; for the other rows, each number is the sum of the two numbers above it. Then a(n) is the sum of all numbers in the table. See examples below. - Jianing Song, Sep 04 2018

Examples

			From _Jianing Song_, Sep 04 2018: (Start)
For n = 4 the inverted triangle table is:
1     2     3     4
   3     5     7
      8    12
        20
So a(4) = 1 + 2 + 3 + 4 + 3 + 5 + 7 + 8 + 12 + 20 = 65.
For n = 5 the inverted triangle table is:
1     2     3     4     5
   3     5     7     9
      8    12    16
        20    28
           48
So a(5) = 1 + 2 + 3 + 4 + 5 + 3 + 5 + 7 + 9 + 8 + 12 + 16 + 20 + 28 + 48 = 171. (End)
		

References

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

Crossrefs

Cf. A053218, A053221 (first differences).

Programs

  • Magma
    [2^n*(n+1)-(n^2+3*n+2)/2: n in [1..30]]; // Vincenzo Librandi, Aug 19 2011
  • Maple
    a := n -> add((n+1)*binomial(n+1, k+1)/2, k=1..n):
    seq(a(n), n=1..30); # Zerinvary Lajos, May 08 2007
    A003469:=(-1+z+z**2)/(2*z-1)**2/(z-1)**3; # conjectured by Simon Plouffe in his 1992 dissertation
  • Mathematica
    Table[(n+1)2^n-(n+1)(n+2)/2, {n, 200}] (* Vladimir Joseph Stephan Orlovsky, Jun 30 2011 *)
    CoefficientList[Series[((2*x + 1)*Exp[2*x] - (x^2/2 + 2*x + 1)*Exp[x])/x, {x, 0, 200}], x]*Table[(k+1)!, {k, 0, 200}] (* Stefano Spezia, Sep 04 2018 *)
  • PARI
    a(n) = (n+1)*2^n-(n+1)*(n+2)/2;
    

Formula

G.f.: x*(1 - x - x^2)/((1 - x)^3*(1 - 2*x)^2).
a(n) = (n + 1)*2^n - (n + 1)*(n + 2)/2. - Paul Barry, Jan 27 2003
E.g.f.: (2*x + 1)*exp(2*x) - (x^2/2 + 2*x + 1)*exp(x). - Jianing Song, Sep 04 2018

Extensions

Offset changed from 2 to 1 by Vincenzo Librandi, Aug 19 2011
Title corrected by Geoffrey Critzer, Jun 29 2013

A132379 Array read by columns down to the main diagonal where the first element in row n is n and element k+1 in row n is the sum of element k in row n and element k-1 in row n-1.

Original entry on oeis.org

1, 1, 4, 1, 5, 12, 1, 6, 17, 1, 7, 23, 49, 80, 1, 8, 30, 72, 129, 192, 1, 9, 38, 102, 201, 321, 448, 1, 10, 47, 140, 303, 522, 769, 1024, 1, 11, 57, 187, 443, 825, 1291, 1785, 2304, 1, 12, 68, 244, 630, 1268, 2116, 3076, 4089, 5120
Offset: 1

Views

Author

Jonathan R. Love (japanada11(AT)yahoo.ca), Nov 20 2007

Keywords

Examples

			The array begins:
(1) 1 1 1 1 1 1 1...
(2 3) 4 5 6 7 8 9...
(3 5 8) 12 17 23 30 38...
(4 7 12 20)32 49 72 102...
(5 9 16 28 48)80 129 201...,
each term being the sum of the number to its left and the number to its top right. Terms in parentheses are given in sequence A053218, the rest included in this sequence.
		

Crossrefs

Cf. A053218.
Showing 1-7 of 7 results.