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 11 results. Next

A128254 Erroneous version of A103450.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 5, 1, 1, 7, 12, 7, 1, 1, 9, 22, 22, 9, 1, 1, 11, 35, 50, 35, 11, 11, 13, 51, 95, 95, 51, 13, 1, 1, 15, 70, 161, 210, 161, 70, 15, 1, 1, 17, 92, 252, 406, 406, 252, 92, 17, 1
Offset: 1

Views

Author

Keywords

A000326 Pentagonal numbers: a(n) = n*(3*n-1)/2.

Original entry on oeis.org

0, 1, 5, 12, 22, 35, 51, 70, 92, 117, 145, 176, 210, 247, 287, 330, 376, 425, 477, 532, 590, 651, 715, 782, 852, 925, 1001, 1080, 1162, 1247, 1335, 1426, 1520, 1617, 1717, 1820, 1926, 2035, 2147, 2262, 2380, 2501, 2625, 2752, 2882, 3015, 3151
Offset: 0

Views

Author

Keywords

Comments

The average of the first n (n > 0) pentagonal numbers is the n-th triangular number. - Mario Catalani (mario.catalani(AT)unito.it), Apr 10 2003
a(n) is the sum of n integers starting from n, i.e., 1, 2 + 3, 3 + 4 + 5, 4 + 5 + 6 + 7, etc. - Jon Perry, Jan 15 2004
Partial sums of 1, 4, 7, 10, 13, 16, ... (1 mod 3), a(2k) = k(6k-1), a(2k-1) = (2k-1)(3k-2). - Jon Perry, Sep 10 2004
Starting with offset 1 = binomial transform of [1, 4, 3, 0, 0, 0, ...]. Also, A004736 * [1, 3, 3, 3, ...]. - Gary W. Adamson, Oct 25 2007
If Y is a 3-subset of an n-set X then, for n >= 4, a(n-3) is the number of 4-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
Solutions to the duplication formula 2*a(n) = a(k) are given by the index pairs (n, k) = (5,7), (5577, 7887), (6435661, 9101399), etc. The indices are integer solutions to the pair of equations 2(6n-1)^2 = 1 + y^2, k = (1+y)/6, so these n can be generated from the subset of numbers [1+A001653(i)]/6, any i, where these are integers, confined to the cases where the associated k=[1+A002315(i)]/6 are also integers. - R. J. Mathar, Feb 01 2008
a(n) is a binomial coefficient C(n,4) (A000332) if and only if n is a generalized pentagonal number (A001318). Also see A145920. - Matthew Vandermast, Oct 28 2008
Even octagonal numbers divided by 8. - Omar E. Pol, Aug 18 2011
Sequence found by reading the line from 0, in the direction 0, 5, ... and the line from 1, in the direction 1, 12, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Sep 08 2011
The hyper-Wiener index of the star-tree with n edges (see A196060, example). - Emeric Deutsch, Sep 30 2011
More generally the n-th k-gonal number is equal to n + (k-2)*A000217(n-1), n >= 1, k >= 3. In this case k = 5. - Omar E. Pol, Apr 06 2013
Note that both Euler's pentagonal theorem for the partition numbers and Euler's pentagonal theorem for the sum of divisors refer more exactly to the generalized pentagonal numbers, not this sequence. For more information see A001318, A175003, A238442. - Omar E. Pol, Mar 01 2014
The Fuss-Catalan numbers are Cat(d,k)= [1/(k*(d-1)+1)]*binomial(k*d,k) and enumerate the number of (d+1)-gon partitions of a (k*(d-1)+2)-gon (cf. Schuetz and Whieldon link). a(n)= Cat(n,3), so enumerates the number of (n+1)-gon partitions of a (3*(n-1)+2)-gon. Analogous sequences are A100157 (k=4) and A234043 (k=5). - Tom Copeland, Oct 05 2014
Binomial transform of (0, 1, 3, 0, 0, 0, ...) (A169585 with offset 1) and second partial sum of (0, 1, 3, 3, 3, ...). - Gary W. Adamson, Oct 05 2015
For n > 0, a(n) is the number of compositions of n+8 into n parts avoiding parts 2 and 3. - Milan Janjic, Jan 07 2016
a(n) is also the number of edges in the Mycielskian of the complete graph K[n]. Indeed, K[n] has n vertices and n(n-1)/2 edges. Then its Mycielskian has n + 3n(n-1)/2 = n(3n-1)/2. See p. 205 of the West reference. - Emeric Deutsch, Nov 04 2016
Sum of the numbers from n to 2n-1. - Wesley Ivan Hurt, Dec 03 2016
Also the number of maximal cliques in the n-Andrásfai graph. - Eric W. Weisstein, Dec 01 2017
Coefficients in the hypergeometric series identity 1 - 5*(x - 1)/(2*x + 1) + 12*(x - 1)*(x - 2)/((2*x + 1)*(2*x + 2)) - 22*(x - 1)*(x - 2)*(x - 3)/((2*x + 1)*(2*x + 2)*(2*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A002412 and A002418. Column 2 of A103450. - Peter Bala, Mar 14 2019
A generalization of the Comment dated Apr 10 2003 follows. (k-3)*A000292(n-2) plus the average of the first n (2k-1)-gonal numbers is the n-th k-gonal number. - Charlie Marion, Nov 01 2020
a(n+1) is the number of Dyck paths of size (3,3n+1); i.e., the number of NE lattice paths from (0,0) to (3,3n+1) which stay above the line connecting these points. - Harry Richman, Jul 13 2021
a(n) is the largest sum of n positive integers x_1, ..., x_n such that x_i | x_(i+1)+1 for each 1 <= i <= n, where x_(n+1) = x_1. - Yifan Xie, Feb 21 2025

Examples

			Illustration of initial terms:
.
.                                       o
.                                     o o
.                          o        o o o
.                        o o      o o o o
.                o     o o o    o o o o o
.              o o   o o o o    o o o o o
.        o   o o o   o o o o    o o o o o
.      o o   o o o   o o o o    o o o o o
.  o   o o   o o o   o o o o    o o o o o
.
.  1    5     12       22           35
- _Philippe Deléham_, Mar 30 2013
		

References

  • Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, pages 2 and 311.
  • Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; p. 129.
  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 189.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 38, 40.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 6.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 1.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 291.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 284.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 64.
  • Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 52-53, 129-130, 132.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 7-10.
  • André Weil, Number theory: an approach through history; from Hammurapi to Legendre, Birkhäuser, Boston, 1984; see p. 186.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, 1987, pp. 98-100.
  • Douglas B. West, Introduction to Graph Theory, 2nd ed., Prentice-Hall, NJ, 2001.

Crossrefs

The generalized pentagonal numbers b*n+3*n*(n-1)/2, for b = 1 through 12, form sequences A000326, A005449, A045943, A115067, A140090, A140091, A059845, A140672, A140673, A140674, A140675, A151542.
Cf. A001318 (generalized pentagonal numbers), A049452, A033570, A010815, A034856, A051340, A004736, A033568, A049453, A002411 (partial sums), A033579.
See A220083 for a list of numbers of the form n*P(s,n)-(n-1)*P(s,n-1), where P(s,n) is the n-th polygonal number with s sides.
Cf. A240137: sum of n consecutive cubes starting from n^3.
Cf. similar sequences listed in A022288.
Partial sums of A016777.

Programs

  • GAP
    List([0..50],n->n*(3*n-1)/2); # Muniru A Asiru, Mar 18 2019
    
  • Haskell
    a000326 n = n * (3 * n - 1) `div` 2  -- Reinhard Zumkeller, Jul 07 2012
    
  • Magma
    [n*(3*n-1)/2 : n in [0..100]]; // Wesley Ivan Hurt, Oct 15 2015
    
  • Maple
    A000326 := n->n*(3*n-1)/2: seq(A000326(n), n=0..100);
    A000326:=-(1+2*z)/(z-1)**3; # Simon Plouffe in his 1992 dissertation
    a[0]:=0:a[1]:=1:for n from 2 to 50 do a[n]:=2*a[n-1]-a[n-2]+3 od: seq(a[n], n=0..50); # Miklos Kristof, Zerinvary Lajos, Feb 18 2008
  • Mathematica
    Table[n (3 n - 1)/2, {n, 0, 60}] (* Stefan Steinerberger, Apr 01 2006 *)
    Array[# (3 # - 1)/2 &, 47, 0] (* Zerinvary Lajos, Jul 10 2009 *)
    LinearRecurrence[{3, -3, 1}, {0, 1, 5}, 61] (* Harvey P. Dale, Dec 27 2011 *)
    pentQ[n_] := IntegerQ[(1 + Sqrt[24 n + 1])/6]; pentQ[0] = True; Select[Range[0, 3200], pentQ@# &] (* Robert G. Wilson v, Mar 31 2014 *)
    Join[{0}, Accumulate[Range[1, 312, 3]]] (* Harvey P. Dale, Mar 26 2016 *)
    (* For Mathematica 10.4+ *) Table[PolygonalNumber[RegularPolygon[5], n], {n, 0, 46}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
    CoefficientList[Series[x (-1 - 2 x)/(-1 + x)^3, {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)
    PolygonalNumber[5, Range[0, 20]] (* Eric W. Weisstein, Dec 01 2017 *)
  • PARI
    a(n)=n*(3*n-1)/2
    
  • PARI
    vector(100, n, n--; binomial(3*n, 2)/3) \\ Altug Alkan, Oct 06 2015
    
  • PARI
    is_a000326(n) = my(s); n==0 || (issquare (24*n+1, &s) && s%6==5); \\ Hugo Pfoertner, Aug 03 2023
    
  • Python
    # Intended to compute the initial segment of the sequence, not isolated terms.
    def aList():
         x, y = 1, 1
         yield 0
         while True:
             yield x
             x, y = x + y + 3, y + 3
    A000326 = aList()
    print([next(A000326) for i in range(47)]) # Peter Luschny, Aug 04 2019

Formula

Product_{m > 0} (1 - q^m) = Sum_{k} (-1)^k*x^a(k). - Paul Barry, Jul 20 2003
G.f.: x*(1+2*x)/(1-x)^3.
E.g.f.: exp(x)*(x+3*x^2/2).
a(n) = n*(3*n-1)/2.
a(-n) = A005449(n).
a(n) = binomial(3*n, 2)/3. - Paul Barry, Jul 20 2003
a(n) = A000290(n) + A000217(n-1). - Lekraj Beedassy, Jun 07 2004
a(0) = 0, a(1) = 1; for n >= 2, a(n) = 2*a(n-1) - a(n-2) + 3. - Miklos Kristof, Mar 09 2005
a(n) = Sum_{k=1..n} (2*n - k). - Paul Barry, Aug 19 2005
a(n) = 3*A000217(n) - 2*n. - Lekraj Beedassy, Sep 26 2006
a(n) = A126890(n, n-1) for n > 0. - Reinhard Zumkeller, Dec 30 2006
a(n) = A049452(n) - A022266(n) = A033991(n) - A005476(n). - Zerinvary Lajos, Jun 12 2007
Equals A034856(n) + (n - 1)^2. Also equals A051340 * [1,2,3,...]. - Gary W. Adamson, Jul 27 2007
a(n) = binomial(n+1, 2) + 2*binomial(n, 2).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0) = 0, a(1) = 1, a(2) = 5. - Jaume Oliver Lafont, Dec 02 2008
a(n) = a(n-1) + 3*n-2 with n > 0, a(0)=0. - Vincenzo Librandi, Nov 20 2010
a(n) = A000217(n) + 2*A000217(n-1). - Vincenzo Librandi, Nov 20 2010
a(n) = A014642(n)/8. - Omar E. Pol, Aug 18 2011
a(n) = A142150(n) + A191967(n). - Reinhard Zumkeller, Jul 07 2012
a(n) = (A000290(n) + A000384(n))/2 = (A000217(n) + A000566(n))/2 = A049450(n)/2. - Omar E. Pol, Jan 11 2013
a(n) = n*A000217(n) - (n-1)*A000217(n-1). - Bruno Berselli, Jan 18 2013
a(n) = A005449(n) - n. - Philippe Deléham, Mar 30 2013
From Oskar Wieland, Apr 10 2013: (Start)
a(n) = a(n+1) - A016777(n),
a(n) = a(n+2) - A016969(n),
a(n) = a(n+3) - A016777(n)*3 = a(n+3) - A017197(n),
a(n) = a(n+4) - A016969(n)*2 = a(n+4) - A017641(n),
a(n) = a(n+5) - A016777(n)*5,
a(n) = a(n+6) - A016969(n)*3,
a(n) = a(n+7) - A016777(n)*7,
a(n) = a(n+8) - A016969(n)*4,
a(n) = a(n+9) - A016777(n)*9. (End)
a(n) = A000217(2n-1) - A000217(n-1), for n > 0. - Ivan N. Ianakiev, Apr 17 2013
a(n) = A002411(n) - A002411(n-1). - J. M. Bergot, Jun 12 2013
Sum_{n>=1} a(n)/n! = 2.5*exp(1). - Richard R. Forberg, Jul 15 2013
a(n) = floor(n/(exp(2/(3*n)) - 1)), for n > 0. - Richard R. Forberg, Jul 27 2013
From Vladimir Shevelev, Jan 24 2014: (Start)
a(3*a(n) + 4*n + 1) = a(3*a(n) + 4*n) + a(3*n+1).
A generalization. Let {G_k(n)}_(n >= 0) be sequence of k-gonal numbers (k >= 3). Then the following identity holds: G_k((k-2)*G_k(n) + c(k-3)*n + 1) = G_k((k-2)*G_k(n) + c(k-3)*n) + G_k((k-2)*n + 1), where c = A000124. (End)
A242357(a(n)) = 1 for n > 0. - Reinhard Zumkeller, May 11 2014
Sum_{n>=1} 1/a(n)= (1/3)*(9*log(3) - sqrt(3)*Pi). - Enrique Pérez Herrero, Dec 02 2014. See the decimal expansion A244641.
a(n) = (A000292(6*n+k-1)-A000292(k))/(6*n-1)-A000217(3*n+k), for any k >= 0. - Manfred Arens, Apr 26 2015 [minor edits from Wolfdieter Lang, May 10 2015]
a(n) = A258708(3*n-1,1) for n > 0. - Reinhard Zumkeller, Jun 23 2015
a(n) = A007584(n) - A245301(n-1), for n > 0. - Manfred Arens, Jan 31 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*(sqrt(3)*Pi - 6*log(2))/3 = 0.85501000622865446... - Ilya Gutkovskiy, Jul 28 2016
a(m+n) = a(m) + a(n) + 3*m*n. - Etienne Dupuis, Feb 16 2017
In general, let P(k,n) be the n-th k-gonal number. Then P(k,m+n) = P(k,m) + (k-2)mn + P(k,n). - Charlie Marion, Apr 16 2017
a(n) = A023855(2*n-1) - A023855(2*n-2). - Luc Rousseau, Feb 24 2018
a(n) = binomial(n,2) + n^2. - Pedro Caceres, Jul 28 2019
Product_{n>=2} (1 - 1/a(n)) = 3/5. - Amiram Eldar, Jan 21 2021
(n+1)*(a(n^2) + a(n^2+1) + ... + a(n^2+n)) = n*(a(n^2+n+1) + ... + a(n^2+2n)). - Charlie Marion, Apr 28 2024
a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * binomial(k, 2)*binomial(3*n+k-1, 2*k). - Peter Bala, Nov 04 2024

Extensions

Incorrect example removed by Joerg Arndt, Mar 11 2010

A002412 Hexagonal pyramidal numbers, or greengrocer's numbers.

Original entry on oeis.org

0, 1, 7, 22, 50, 95, 161, 252, 372, 525, 715, 946, 1222, 1547, 1925, 2360, 2856, 3417, 4047, 4750, 5530, 6391, 7337, 8372, 9500, 10725, 12051, 13482, 15022, 16675, 18445, 20336, 22352, 24497, 26775, 29190, 31746, 34447, 37297, 40300
Offset: 0

Views

Author

Keywords

Comments

Binomial transform of (1, 6, 9, 4, 0, 0, 0, ...). - Gary W. Adamson, Oct 16 2007
a(n) is the sum of the maximum(m,n) over {(m,n): m,n in positive integers, m<=n}. - Geoffrey Critzer, Oct 11 2009
We obtain these numbers for d=2 in the identity n*(n*(d*n-d+2)/2)-sum(k*(d*k-d+2)/2, k=0..n-1) = n*(n+1)*(2*d*n-2*d+3)/6 (see Klaus Strassburger in Formula lines). - Bruno Berselli, Apr 21 2010, Nov 16 2010
q^a(n) is the Hankel transform of the q-Catalan numbers. - Paul Barry, Dec 15 2010
Row 1 of the convolution array A213835. - Clark Kimberling, Jul 04 2012
From Ant King, Oct 24 2012: (Start)
For n>0, the digital roots of this sequence A010888(A002412(n)) form the purely periodic 27-cycle {1,7,4,5,5,8,9,3,3,4,1,7,8,8,2,3,6,6,7,4,1,2,2,5,6,9,9}.
For n>0, the units' digits of this sequence A010879(A002412(n)) form the purely periodic 20-cycle {1,7,2,0,5,1,2,2,5,5,6,2,7,5,0,6,7,7,0,0}.
(End)
Partial sums of A000384. - Omar E. Pol, Jan 12 2013
Row sums of A094728. - J. M. Bergot, Jun 14 2013
Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 40320. - Philippe A.J.G. Chevalier, Dec 28 2015
Coefficients in the hypergeometric series identity 1 - 7*(x - 1)/(3*x + 1) + 22*(x - 1)*(x - 2)/((3*x + 1)*(3*x + 2)) - 50*(x - 1)*(x - 2)*(x - 3)/((3*x + 1)*(3*x + 2)*(3*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A000326 and A002418. Column 3 of A103450. - Peter Bala, Mar 14 2019

Examples

			Let n=5, 2*n=10. Since 10 = 1 + 9 = 2 + 8 = 3 + 7 = 4 + 6 = 5 + 5, a(5) = 1*9 + 2*8 + 3*7 + 4*6 + 5*5 = 95. - _Vladimir Shevelev_, May 11 2012
		

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
  • E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 2.
  • T. A. Gulliver, Sequences from Arrays of Integers, Int. Math. Journal, Vol. 1, No. 4, pp. 323-332, 2002.
  • I. Siap, Linear codes over F_2 + u*F_2 and their complete weight enumerators, in Codes and Designs (Ohio State, May 18, 2000), pp. 259-271. De Gruyter, 2002.
  • 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

Bisection of A002623. Equals A000578(n) - A000330(n-1).
a(n) = A093561(n+2, 3), (4, 1)-Pascal column.
Cf. A220084 for a list of numbers of the form n*P(k,n)-(n-1)*P(k,n-1), where P(k,n) is the n-th k-gonal pyramidal number (see Adamson's formula).
Cf. similar sequences listed in A237616.
Orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008585, A005843, A001477, A000217.

Programs

  • GAP
    List([0..40],n->n*(n+1)*(4*n-1)/6); # Muniru A Asiru, Mar 18 2019
    
  • Magma
    [n*(n+1)*(4*n-1)/6: n in [0..40]]; // Vincenzo Librandi, Nov 28 2015
    
  • Maple
    seq(sum(i*(2*k-i), i=1..k), k=0..100); # Wesley Ivan Hurt, Sep 25 2013
  • Mathematica
    Figurate[ ngon_, rank_, dim_] := Binomial[rank + dim - 2, dim - 1] ((rank - 1)*(ngon - 2) + dim)/dim; Table[ Figurate[6, r, 3], {r, 0, 40}] (* Robert G. Wilson v, Aug 22 2010 *)
    Table[n(n+1)(4n-1)/6, {n,0,40}] (* or *) LinearRecurrence[{4,-6,4,-1}, {0,1,7,22}, 40] (* Harvey P. Dale, Jul 16 2011 *)
  • Maxima
    A002412(n):=n*(n+1)*(4*n-1)/6$ makelist(A002412(n),n,0,20); /* Martin Ettl, Dec 12 2012 */
    
  • PARI
    v=vector(40,i,(i*(i+1))\2); s=0; print1(s","); forstep(i=1,40,2,s+=v[i]; print1(s","))
    
  • Python
    print([n*(n+1)*(4*n-1)//6 for n in range(40)]) # Michael S. Branicky, Mar 28 2022

Formula

a(n) = n(n + 1)(4n - 1)/6.
G.f.: x*(1+3*x)/(1-x)^4. - Simon Plouffe in his 1992 dissertation.
a(n) = n^3 - Sum_{i=1..n-1} i^2. - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de)
Partial sums of n odd-indexed triangular numbers, e.g., a(3) = t(1)+t(3)+t(5) = 1+6+15 = 22. - Jon Perry, Jul 23 2003
a(n) = Sum_{i=0..n-1} (n - i)*(n + i). - Jon Perry, Sep 26 2004
a(n) = n*A000292(n) - (n-1)*A000292(n-1) = n*binomial((n+2),3) - (n-1)*binomial((n+1),3); e.g., a(5) = 95 = 5*35 - 4*20. - Gary W. Adamson, Dec 28 2007
a(n) = Sum_{i=0..n} (2i^2 + 3i + 1), for n >= 0 (Omits the leading 0). - William A. Tedeschi, Aug 25 2010
a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4), with a(0)=0, a(1)=1, a(2)=7, a(3)=22. - Harvey P. Dale, Jul 16 2011
a(n) = sum a*b, where the summing is over all unordered partitions 2*n = a+b. - Vladimir Shevelev, May 11 2012
From Ant King, Oct 24 2012: (Start)
a(n) = a(n-1) + n*(2*n-1).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 4.
a(n) = (n+1)*(2*A000384(n) + n)/6 = (4*n-1)*A000217(n)/3.
a(n) = A000292(n) + 3*A000292(n-1) = A002411(n) + A000292(n-1).
a(n) = binomial(n+2,3) + 3*binomial(n+1,3) = (4*n-1)*binomial(n+1,2)/3.
Sum_{n>=1} 1/a(n) = 6*(12*log(2)-2*Pi-1)/5 = 1.2414...
(End)
a(n) = Sum_{i=1..n} Sum_{j=1..n} max(i,j) = Sum_{i=1..n} i*(2*n-i). - Enrique Pérez Herrero, Jan 15 2013
a(n) = A005900(n+1) - A000326(n+1) = Octahedral - Pentagonal Numbers. - Richard R. Forberg, Aug 07 2013
a(n) = n*A000217(n) + Sum_{i=0..n-1} A000217(i). - Bruno Berselli, Dec 18 2013
a(n) = 2n * A000217(n) - A000330(n). - J. M. Bergot, Apr 05 2014
a(n) = A080851(4,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: x*(6 + 15*x + 4*x^2)*exp(x)/6. - Ilya Gutkovskiy, May 12 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = 6*(1 + 2*sqrt(2)*Pi - 2*(3+sqrt(2))*log(2) + 4*sqrt(2)*log(2-sqrt(2)))/5. - Amiram Eldar, Jan 04 2022

A045623 Number of 1's in all compositions of n+1.

Original entry on oeis.org

1, 2, 5, 12, 28, 64, 144, 320, 704, 1536, 3328, 7168, 15360, 32768, 69632, 147456, 311296, 655360, 1376256, 2883584, 6029312, 12582912, 26214400, 54525952, 113246208, 234881024, 486539264, 1006632960, 2080374784, 4294967296, 8858370048, 18253611008, 37580963840
Offset: 0

Views

Author

Keywords

Comments

Let M_n be the n X n matrix m_(i,j) = 2 + abs(i-j) then det(M_n) = (-1)^(n-1)*a(n-1). - Benoit Cloitre, May 28 2002
a(n) is the number of triangulations of a regular (n+3)-gon in which every triangle shares at least one side with the polygon itself. - David Callan, Mar 25 2004
Number of compositions of j+n, j>n and j the maximum part. E.g. a(4) is derived from the number of compositions of, for example: 54(2), 531(6), 522(3), 5211(12) and 51111(5) giving 2+6+3+12+5=28. - Jon Perry, Sep 13 2005
If X_1,X_2,...,X_n are 2-blocks of a (2n+2)-set X then, for n>=1, a(n+1) is the number of (n+1)-subsets of X intersecting each X_i, (i=1,2,...,n). - Milan Janjic, Nov 18 2007
Generated from iterates of M * [1,1,1,...], where M = an infinite triadiagonal matrix with (1,1,1,...) in the main and superdiagonals and (1,0,0,0,...) in the subdiagonal. - Gary W. Adamson, Jan 04 2009
a(n) is the number of weak compositions of n with exactly 1 part equal to 0. - Milan Janjic, Jun 27 2010
An elephant sequence, see A175654. For the corner squares 16 A[5] vectors, with decimal values between 19 and 400, lead to this sequence. For the central square these vectors lead to the companion sequence A045891 (without the first leading 1). - Johannes W. Meijer, Aug 15 2010
Equals first finite difference row of A001792: (1, 3, 8, 20, 48, 112, ...). - Gary W. Adamson, Oct 26 2010
With alternating signs the g.f. is: (1 + x)^2/(1 + 2*x)^2.
Number of 132-avoiding permutations of [n+2] containing exactly one 213 pattern. - David Scambler, Nov 07 2011
a(n) is the number of 1's in all compositions of n+1 = the number of 2's in all compositions of n+2 = the number of 3's in all compositions of n+3 = ... So the partial sums = A001792. - Geoffrey Critzer, Feb 12 2012
Also number of compositions of n into 2 sorts of parts where all parts of the first sort precede all parts of the second sort; see example. - Joerg Arndt, Apr 28 2013
a(n) is also the difference of the total number of parts between all compositions of n+1 and all compositions of n. The equivalent sequence for partitions is A138137. - Omar E. Pol, Aug 28 2013
Except for an initial 1, this is the p-INVERT of (1,1,1,1,1,...) for p(S) = (1 - S)^2; see A291000. - Clark Kimberling, Aug 24 2017
For a composition of n, the total number of runs of parts of size k is a(n-k) - a(n-2k). - Gregory L. Simay, Feb 17 2018
a(n) is the number of binary trees on n+1 nodes that are isomorphic to a path graph. The ratio of a(n)/A000108(n+1) gives the probability that a random Catalan tree on n+1 nodes is isomorphic to a path graph. - Marcel K. Goh, May 09 2020
a(n) is the number of words of length n over the alphabet {0,1,2} such that the first letter is not 2 and the last 1 occurs before the first 0. - Henri Mühle, Mar 08 2021
Also the number of "special permutations" in the Weng and Zagier reference. - F. Chapoton, Sep 30 2022
a(n-k) is the total number of runs of 1s of length k over all binary n-strings. - Félix Balado, Dec 11 2022

Examples

			E.g. a(2)=5 because in the compositions of 3, namely 3,2+1,1+2,1+1+1, we have five 1's altogether.
There are a(3)=12 compositions of 3 into 2 sorts of parts where all parts of the first sort precede all parts of the second sort. Here p:s stands for "part p of sort s":
01:  [ 1:0  1:0  1:0  ]
02:  [ 1:0  1:0  1:1  ]
03:  [ 1:0  1:1  1:1  ]
04:  [ 1:0  2:0  ]
05:  [ 1:0  2:1  ]
06:  [ 1:1  1:1  1:1  ]
07:  [ 1:1  2:1  ]
08:  [ 2:0  1:0  ]
09:  [ 2:0  1:1  ]
10:  [ 2:1  1:1  ]
11:  [ 3:0  ]
12:  [ 3:1  ]
- _Joerg Arndt_, Apr 28 2013
For the compositions of 6, the total number of runs of parts of size 2 is a(6-2) - a(6-2*2) = 28 - 5 = 23, enumerated as follows (with the runs of 2 enclosed in []): 4,[2]; [2],4; [2],3,1; [2],1,3; 3,[2],1; 1,[2],3; 3,1,[2]; 1,3,[2]; [2,2,2]; [2,2],1,1; 1,[2,2],1; 1,1,[2,2]; [2],1,[2],1; 1,[2],1,[2]; [2],1,1,[2]; [2],1,1,1,1; 1,[2],1,1,1; 1,1,[2],1,1; 1,1,1,[2],1; and 1,1,1,1[2]. - _Gregory L. Simay_, Feb 17 2018
There are a(3)=12 triwords of length 3: (0,0,0), (0,0,2), (0,2,0), (0,2,2), (1,0,0), (1,0,2), (1,1,0), (1,1,1), (1,1,2), (1,2,0), (1,2,1), (1,2,2). - _Henri Mühle_, Mar 08 2021
		

References

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

Crossrefs

Convolution of A011782.
Row sums of A103450, A152195, A177992, A198069.
Cf. A001792.

Programs

  • GAP
    a:=[2,5];; for n in [3..40] do a[n]:=4*a[n-1]-4*a[n-2]; od; Concatenation([1],a); # Muniru A Asiru, Oct 16 2018
    
  • Haskell
    a045623 n = a045623_list !! n
    a045623_list = tail $ f a011782_list [] where
       f (u:us) vs = sum (zipWith (*) vs $ reverse ws) : f us ws
         where ws = u : vs
    -- Reinhard Zumkeller, Jul 21 2013
    
  • Maple
    seq(ceil(1/4*2^n*(n+3)),n=0..50);
  • Mathematica
    Table[If[n==0, 1, 2^(n-2)(n+3)], {n, 0, 29}] (* Robert G. Wilson v, Jun 27 2005 *)
    CoefficientList[Series[(1 -2x +x^2)/(1-2x)^2, {x, 0, 30}], x] (* or *)
    LinearRecurrence[{4, -4}, {1, 2, 5}, 31] (* Robert G. Wilson v, Feb 18 2018 *)
  • Maxima
    a(n):=sum(((2*m+2)*n-2*m^2+1)*binomial(2*n+2,2*m+1),m,0,n)/((4*n+2)*2^n); /* Vladimir Kruchinin, Nov 01 2020 */
  • PARI
    a(n)=if(n<1,n==0,(n+3)*2^(n-2))
    

Formula

Sum_{k = 0..n} (k+2)*binomial(n,k) gives the sequence but with a different offset: 2, 5, 12, 28, 64, 144, 320, 704, 1536, ... - N. J. A. Sloane, Jan 30 2008 - formula corrected by Robert G. Wilson v, Feb 26 2018
Binomial transform of 1,1,2,2,3,3,... . - Paul Barry, Mar 06 2003
a(0)=1, a(n) = (n+3)*2^(n-2), n >= 1.
a(n+1) = 2*a(n) + 2^(n-1), n>0.
G.f.: (1-x)^2/(1-2*x)^2. - Detlef Pauly (dettodet(AT)yahoo.de), Mar 03 2003
G.f.: 1/(1-x-x^2-x^3-...)^2. - Jon Perry, Jul 04 2004
a(n) = Sum_{0 <= j <= k <= n} binomial(n, j+k). - Benoit Cloitre, Oct 14 2004
a(n) = Sum_{k=0..n} C(n, k)*floor((k+2)/2). - Paul Barry, Mar 06 2003
a(n+1) - 2*a(n) = A131577(n). - Paul Curtz, May 18 2008
G.f.: 1/(1-x) + Q(0)*x/(1-x)^3, where Q(k)= 1 + (k+1)*x/(1 - x - x*(1-x)/(x + (k+1)*(1-x)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 25 2013
a(n) = Sum_{k=0..n} (k+1)*C(n-1,n-k). - Peter Luschny, Apr 20 2015
a(n) = Sum_{k=0..n-1} a(k) + 2^(n-1) = A001787(n-1) + 2^n, a(0)=1. - Yuchun Ji, May 22 2020
a(n) = Sum_{m=0..n}((2*m+2)*n-2*m^2+1)*C(2*n+2,2*m+1)/((4*n+2)*2^n). - Vladimir Kruchinin, Nov 01 2020
E.g.f.: (1 + exp(2*x)*(3 + 2*x))/4. - Stefano Spezia, Dec 19 2021
From Amiram Eldar, Jan 05 2022: (Start)
Sum_{n>=0} 1/a(n) = 32*log(2) - 61/3.
Sum_{n>=0} (-1)^n/a(n) = 32*log(3/2) - 37/3. (End)

A002418 4-dimensional figurate numbers: a(n) = (5*n-1)*binomial(n+2,3)/4.

Original entry on oeis.org

0, 1, 9, 35, 95, 210, 406, 714, 1170, 1815, 2695, 3861, 5369, 7280, 9660, 12580, 16116, 20349, 25365, 31255, 38115, 46046, 55154, 65550, 77350, 90675, 105651, 122409, 141085, 161820, 184760, 210056, 237864, 268345, 301665, 337995
Offset: 0

Views

Author

Keywords

Comments

Partial sums of A002413.
Principal diagonal of the convolution array A213550, for n>0. - Clark Kimberling, Jun 17 2012
Convolution of A000027 with A000566. - Bruno Berselli, Dec 06 2012
Coefficients in the hypergeometric series identity 1 - 9*(x - 1)/(4*x + 1) + 35*(x - 1)*(x - 2)/((4*x + 1)*(4*x + 2)) - 95*(x - 1)*(x - 2)*(x - 3)/((4*x + 1)*(4*x + 2)*(4*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A000326 and A002412. Column 4 of A103450. - Peter Bala, Mar 14 2019

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A093562 ((5, 1) Pascal, column m=4).
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.

Programs

  • GAP
    List([0..40],n->(5*n-1)*Binomial(n+2,3)/4); # Muniru A Asiru, Mar 18 2019
    
  • Magma
    [(5*n - 1)*Binomial(n + 2, 3)/4: n in [0..40]]; // Vincenzo Librandi, Oct 17 2012
    
  • Magma
    /* A000027 convolved with A000566: */ A000566:=func; [&+[(n-i+1)*A000566(i): i in [0..n]]: n in [0..35]]; // Bruno Berselli, Dec 06 2012
    
  • Mathematica
    Table[(5n-1) Binomial[n+2,3]/4,{n,0,40}] (* or *) LinearRecurrence[ {5,-10,10,-5,1},{0,1,9,35,95},40] (* Harvey P. Dale, Oct 16 2012 *)
    CoefficientList[Series[x*(1 + 4*x)/(1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Oct 17 2012 *)
  • PARI
    a(n)=(5*n-1)*binomial(n+2,3)/4 \\ Charles R Greathouse IV, Sep 24 2015
    
  • Sage
    [(5*n-1)*binomial(n+2,3)/4 for n in (0..40)] # G. C. Greubel, Jul 03 2019

Formula

G.f.: x*(1+4*x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation.
Starting (1, 9, 35, 95, ...), = A128064 * A000332, (A000332 starting 1, 5, 15, 35, 70, ...), such that a(n) = n*C(n+3,4) - (n-1)*C(n+2,4). E.g., a(5) = 210 = 5*C(8,4) - 4*C(7,4) = 5*70 - 4*35. - Gary W. Adamson, Dec 28 2007
Unit digit, A010879(a(n)), is one of {0,1,9,5,6,4} [Eric Desbiaux] because a(n) mod 5 = 0,1,4,0,0, periodic with period 5. [Proof: A002413(n) mod 5 = 1,3,1,0,0 with period 5 and a(n) are the partial sums of A002413.] - R. J. Mathar, Mar 19 2008
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Harvey P. Dale, Oct 16 2012
a(n) = A080852(5,n-1). - R. J. Mathar, Jul 28 2016
a(n) = Sum_{i=0..n} (n-i) * Sum_{j=i..n} j. - J. M. Bergot, May 30 2017
E.g.f.: x*(24 + 84*x + 44*x^2 + 5*x^3)*exp(x)/4!. - G. C. Greubel, Jul 03 2019
Sum_{n>=1} 1/a(n) = (50*sqrt(5)*log(phi) + 125*log(5) - 50*sqrt(1+2/sqrt(5))*Pi - 26)/11, where phi is the golden ratio (A001622). - Amiram Eldar, Feb 11 2022

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

Original entry on oeis.org

1, 5, 22, 95, 406, 1722, 7260, 30459, 127270, 529958, 2200276, 9111830, 37650172, 155266100, 639191160, 2627302995, 10784089350, 44208873390, 181025067300, 740483276610, 3026059513620, 12355464845100
Offset: 0

Views

Author

Paul Barry, Jun 14 2008

Keywords

Comments

Apply Riordan array (1/sqrt(1-4x), xc(x)) to A131056, c(x) the g.f. of A000108.
Apply Riordan array (c(x)/sqrt(1-4*x), x*c(x)^2) to A131055.
Hankel transform appears to be (-1)^n*A085046(n).
Coefficients T(2*n+1,n) of triangle A103450. [Emanuele Munarini, Jun 01 2012, corrected by Werner Schulte, Nov 27 2021]

Programs

  • Mathematica
    Table[((1+3*n+n^2)*Binomial[2*n, n])/(n+1),{n,0,20}] (* Vaclav Kotesovec, Feb 13 2014 *)
    CoefficientList[Series[-1/(2*x)+(2*x-1)^2/(2*x*(1-4x)^(3/2)),{x,0,20}],x] (* Vaclav Kotesovec, Feb 13 2014 *)
    a[n_] := (1 + 3 n + n^2) CatalanNumber[n];
    Table[a[n], {n, 0, 21}] (* Peter Luschny, Nov 28 2021 *)
  • Maxima
    a(n):=sum(binomial(2*n,k)*binomial(n+1,2*n-k),k,0,n); makelist(a(n),n,0,40); /* Emanuele Munarini, Jun 01 2012 */

Formula

a(n) = Sum_{k=0..n} (1 + (k+1)*2^(k-1) - 0^k/2)*C(2n-k,n-k); a(n) = Sum_{k=0..n} C(2n,k)*C(n+1,2n-k).
Equals the Narayana transform (A001263) of integer squares. - Gary W. Adamson, Jul 29 2011
Conjecture: (n+1)*a(n) + 2*(-3*n-1)*a(n-1) + 4*(2*n-3)*a(n-2) = 0. - R. J. Mathar, Nov 24 2012
From Vaclav Kotesovec, Feb 13 2014: (Start)
G.f.: -1/(2*x) + (2*x-1)^2/(2*x*(1-4x)^(3/2)).
a(n) = (1 + 3*n + n^2) * C(2*n,n) / (n+1).
Recurrence: (n+1)*(n^2 + n - 1)*a(n) = 2*(2*n-1)*(n^2 + 3*n + 1)*a(n-1).
(End)

Extensions

Name of the sequence corrected by Vaclav Kotesovec, Feb 13 2014

A382436 Triangle read by rows, defined by the two-variable g.f. 1/(1 - (y + 1)*x - y*x^2 - (y^2 + y)*x^3).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 6, 1, 1, 9, 17, 9, 1, 1, 12, 36, 36, 12, 1, 1, 15, 64, 101, 64, 15, 1, 1, 18, 101, 227, 227, 101, 18, 1, 1, 21, 147, 440, 627, 440, 147, 21, 1, 1, 24, 202, 767, 1459, 1459, 767, 202, 24, 1, 1, 27, 266, 1235, 2994, 3999, 2994, 1235, 266, 27, 1
Offset: 0

Views

Author

F. Chapoton, Mar 25 2025

Keywords

Comments

The original definition was "Decomposition of A077938".
Every row is symmetric.

Examples

			Triangle begins:
  1;
  1,  1;
  1,  3,   1;
  1,  6,   6,    1;
  1,  9,  17,    9,    1;
  1, 12,  36,   36,   12,    1;
  1, 15,  64,  101,   64,   15,    1;
  1, 18, 101,  227,  227,  101,   18,    1;
  1, 21, 147,  440,  627,  440,  147,   21,   1;
  1, 24, 202,  767, 1459, 1459,  767,  202,  24,  1;
  1, 27, 266, 1235, 2994, 3999, 2994, 1235, 266, 27, 1;
  ...
		

Crossrefs

Similar to A008288, A103450, and A382444.
Row sums are A077938.
T(2n, n) gives A339565.
Cf. A056594.

Programs

  • Sage
    y = polygen(QQ, 'y')
    x = y.parent()[['x']].gen()
    inverse = 1 + (-y - 1)*x - y*x^2 + (-y^2 - y)*x^3
    gf = 1 / inverse
    [list(u) for u in list(gf.O(11))]

Formula

G.f. 1/(1 - (y + 1)*x - y*x^2 - (y^2 + y)*x^3).
Sum_{k=0..n} (-1)^k * T(n,k) = A056594(n). - Alois P. Heinz, Mar 25 2025

A128255 A114219(signed) * A007318.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 2, 6, 5, 1, 3, 10, 13, 7, 1, 3, 15, 27, 23, 9, 1, 4, 21, 48, 57, 36, 1, 4, 28, 78, 118, 104, 52, 13, 1, 5, 36, 118, 218, 246, 172, 71, 15, 1, 5, 45, 170, 370, 510, 458, 265, 93, 17, 1
Offset: 1

Views

Author

Gary W. Adamson, Feb 20 2007

Keywords

Comments

Row sums = A059570: (1, 2, 6, 14, 34, 78, 178,...).

Examples

			First few rows of the triangle are:
1;
1, 1;
2, 3, 1;
2, 6, 5, 1;
3, 10, 13, 7, 1;
3, 15, 27, 23, 9, 1;
4, 21, 48, 57, 36, 11, 1;
...
		

Crossrefs

Formula

Let the signed version of A114219 {1; 0,1; 0,1,1; 0,-1,2,1; 0,1,-2,-3,4;...} = M; and P = Pascal's triangle, A007318. Then A128255 = A114219(signed) * A007318.

A382444 Triangle read by rows, defined by the two-variable g.f. (1 + y*x^2 + (y^2 + y)*x^3)/(1-(1+y)*x-y*x^2).

Original entry on oeis.org

1, 1, 1, 1, 4, 1, 1, 7, 7, 1, 1, 9, 18, 9, 1, 1, 11, 34, 34, 11, 1, 1, 13, 54, 86, 54, 13, 1, 1, 15, 78, 174, 174, 78, 15, 1, 1, 17, 106, 306, 434, 306, 106, 17, 1, 1, 19, 138, 490, 914, 914, 490, 138, 19, 1, 1, 21, 174, 734, 1710, 2262, 1710, 734, 174, 21, 1
Offset: 0

Views

Author

F. Chapoton, Mar 25 2025

Keywords

Comments

Every row is symmetric.

Examples

			Triangle begins:
  [0] [1]
  [1] [1,  1]
  [2] [1,  4,   1]
  [3] [1,  7,   7,   1]
  [4] [1,  9,  18,   9,   1]
  [5] [1, 11,  34,  34,  11,   1]
  [6] [1, 13,  54,  86,  54,  13,   1]
  [7] [1, 15,  78, 174, 174,  78,  15,   1]
  [8] [1, 17, 106, 306, 434, 306, 106,  17,  1]
  [9] [1, 19, 138, 490, 914, 914, 490, 138, 19, 1]
  ...
		

Crossrefs

Similar to A008288, A103450 and A382436. Row sums are A265107.

Programs

  • Sage
    y = polygen(QQ, 'y')
    x = y.parent()[['x']].gen()
    gf = (1 + y*x^2 + (y^2 + y)*x^3)/(1 - (1 + y)*x - y*x^2)
    [list(u) for u in list(gf.O(11))]

A350512 Triangle read by rows with T(n,0) = 1 for n >= 0 and T(n,k) = binomial(n-1,k-1)*(2*k*(n-k) + n)/k for 0 < k <= n.

Original entry on oeis.org

1, 1, 1, 1, 4, 1, 1, 7, 7, 1, 1, 10, 18, 10, 1, 1, 13, 34, 34, 13, 1, 1, 16, 55, 80, 55, 16, 1, 1, 19, 81, 155, 155, 81, 19, 1, 1, 22, 112, 266, 350, 266, 112, 22, 1, 1, 25, 148, 420, 686, 686, 420, 148, 25, 1, 1, 28, 189, 624, 1218, 1512, 1218, 624, 189, 28, 1
Offset: 0

Views

Author

Werner Schulte, Jan 02 2022

Keywords

Comments

Depending on some fixed integer m there is a family of number triangles T(m; n,k) for 0 <= k <= n with entries: T(m; n,0) = 1 for n >= 0 and T(m; n,k) = binomial(n-1,k-1)*(m*k*(n-k) + n)/k for 0 < k <= n.
Special cases: m=0 (A007318), m=1 (A103450), and m=2 (this triangle).
Further properties: T(m; n,n) = 1 for n >= 0; T(m; n,k) = T(m; n,n-k) for 0 <= k <= n; T(m; 2*n,n) = A000108(n)*A086270(m,n+1) for n >= 0 and m > 0.
T(m; n,k) = T(m; n-1,k) + T(m; n-1,k-1) + m*binomial(n-2,k-1) for 0 < k < n.
G.f. of column k: (1 + m*k*x) * x^k / (1 - x)^(k+1).
G.f.: A(x, t) = (1 - (1+x)*t + m*x*t^2) / (1 - (1+x)*t)^2.
T(m; n,k) = [x^k] (1 + (m*n - m + 2)*x + x^2) * (1 + x)^(n-2) for n > 0.

Examples

			Triangle T(n, k) for 0 <= k <= n starts:
n\k :  0   1    2    3    4    5    6    7   8  9
=================================================
  0 :  1
  1 :  1   1
  2 :  1   4    1
  3 :  1   7    7    1
  4 :  1  10   18   10    1
  5 :  1  13   34   34   13    1
  6 :  1  16   55   80   55   16    1
  7 :  1  19   81  155  155   81   19    1
  8 :  1  22  112  266  350  266  112   22   1
  9 :  1  25  148  420  686  686  420  148  25  1
  etc.
		

Crossrefs

Row sums are A057711(n+1).

Programs

  • Mathematica
    Flatten[Table[Join[{1},Table[Binomial[n-1,k-1](2*k*(n-k) + n)/k,{k,n}]],{n,0,10}]] (* Stefano Spezia, Jan 06 2022 *)

Formula

T(n, n) = 1; T(n, k) = T(n, n-k).
T(2*n, n) = (n+1)^2 * A000108(n).
T(n, k) = T(n-1, k) + T(n-1, k-1) + 2 * binomial(n-2,k-1) for 0 < k < n.
G.f. of column k: (1 + 2*k*x) * x^k / (1 - x)^(k+1).
G.f.: A(x,t) = (1 - (1 + x) * t + 2 * x * t^2) / (1 - (1 + x) * t)^2.
T(n,k) = [x^k] (1 + 2 * n * x + x^2) * (1 + x)^(n-2) for n > 0.
Showing 1-10 of 11 results. Next