A002414 -id:A002414 - cp's OEIS frontend

A000567 Octagonal numbers: n(3n-2). Also called star numbers.

0, 1, 8, 21, 40, 65, 96, 133, 176, 225, 280, 341, 408, 481, 560, 645, 736, 833, 936, 1045, 1160, 1281, 1408, 1541, 1680, 1825, 1976, 2133, 2296, 2465, 2640, 2821, 3008, 3201, 3400, 3605, 3816, 4033, 4256, 4485, 4720, 4961, 5208, 5461

Offset: 0

N. J. A. Sloane

From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0,1,....
The spiral begins:
.
                 85--84--83--82--81--80
                 /                     \
               86  56--55--54--53--52  79
               /   /                 \   \
             87  57  33--32--31--30  51  78
             /   /   /             \   \   \
           88  58  34  16--15--14  29  50  77
           /   /   /   /         \   \   \   \
         89  59  35  17   5---4  13  28  49  76
         /   /   /   /   /     \   \   \   \   \
       90  60  36  18   6   0   3  12  27  48  75
       /   /   /   /   /   /   /   /   /   /   /
     91  61  37  19   7   1---2  11  26  47  74
       \   \   \   \   \ .       /   /   /   /
       92  62  38  20   8---9--10  25  46  73
         \   \   \   \ .           /   /   /
         93  63  39  21--22--23--24  45  72
           \   \   \ .               /   /
           94  64  40--41--42--43--44  71
             \   \ .                   /
             95  65--66--67--68--69--70
               \ .
               96
.
(End)
From Lekraj Beedassy, Oct 02 2003: (Start)
Also the number of distinct three-cell blocks that may be removed out of A000217(n+1) square cells arranged in a stepping triangular array of side (n+1). A 5-layer triangular array of square cells, for instance, has vertices outlined thus:
  x x
  x x x
  x x x x
  x x x x x
  x x x x x x
  x x x x x x  (End)
First derivative at n of A045991. - Ross La Haye, Oct 23 2004
Starting from n=1, the sequence corresponds to the Wiener index of K_{n,n} (the complete bipartite graph wherein each independent set has n vertices). - Kailasam Viswanathan Iyer, Mar 11 2009
Number of divisors of 24^(n-1) for n > 0 (cf A009968). - J. Lowell, Aug 30 2008
a(n) = A001399(6n-5), number of partitions of 6*n - 5 into parts < 4. For example a(2)=8 and partitions of 6*2 - 5 = 7 into parts < 4 are: [1,1,1,1,1,1,1], [1,1,1,1,1,2],[1,1,1,1,3], [1,1,1,2,2], [1,1,2,3], [1,2,2,2], [1,3,3], [2,2,3]. - Adi Dani, Jun 07 2011
Also, sequence found by reading the line from 0 in the direction 0, 8, ..., and the parallel line from 1 in the direction 1, 21, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. - Omar E. Pol, Sep 10 2011
Partial sums give A002414. - Omar E. Pol, Jan 12 2013
Generate a Pythagorean triple using Euclid's formula with (n, n-1) to give A,B,C. a(n) = B + (A + C)/2. - J. M. Bergot, Jul 13 2013
The number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 773", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
For n >= 1, the continued fraction expansion of sqrt(27*a(n)) is [9n-4; {1, 2n-2, 3, 2n-2, 1, 18n-8}]. For n=1, this collapses to [5; {5, 10}]. - Magus K. Chu, Oct 10 2022
a(n)*a(n+1) + 1 = (3n^2 + n - 1)^2. In general, a(n)*a(n+k) + k^2 = (3n^2 + (3k-2)n - k)^2. - Charlie Marion, May 23 2023

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 p. 38.
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.
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 19-20.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 123.

T. D. Noe, Table of n, a(n) for n = 0..1000
Raghavendra N. Bhat, Cristian Cobeli, and Alexandru Zaharescu, A lozenge triangulation of the plane with integers, arXiv:2403.10500 [math.NT], 2024.
Francesco Brenti and Paolo Sentinelli, Wachs permutations, Bruhat order and weak order, arXiv:2212.04932 [math.CO], 2022.
Cesar Ceballos and Viviane Pons, The s-weak order and s-permutahedra II: The combinatorial complex of pure intervals, arXiv:2309.14261 [math.CO], 2023. See p. 42.
C. K. Cook and M. R. Bacon, Some polygonal number summation formulas, Fib. Q., 52 (2014), 336-343.
John Elias, Illustration: Hexagonal spiral grids based on generalized octagonal numbers; Illustration: Generalized Square-Octagonal Grids.
John Elias, Illustration: Nesting Cube-frames; Nesting Cube Animation; Nesting-frames Decomposition; Factorial Compartmentalization.
Ghislain R. Franssens, On a Number Pyramid Related to the Binomial, Deleham, Eulerian, MacMahon and Stirling number triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.4.1.
Lancelot Hogben, Choice and Chance by Cardpack and Chessboard, Vol. 1, Max Parrish and Co, London, 1950, p. 36.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 342.
Milan Janjic and B. Petkovic, A Counting Function, arXiv 1301.4550 [math.CO], 2013.
R. Kemp, On the number of words in the language {w in Sigma* | w = w^R }^2, Discrete Math., 40 (1982), 225-234. See Table 1.
Hyun Kwang Kim, On Regular Polytope Numbers, Proc. Amer. Math. Soc., 131 (2002), 65-75.
Kaie Kubjas, Luca Sodomaco, and Elias Tsigaridas, Exact solutions in low-rank approximation with zeros, arXiv:2010.15636 [math.AG], 2020.
Viktor Levandovskyy, Christoph Koutschan, and Oleksandr Motsak, On Two-generated Non-commutative Algebras Subject to the Affine Relation, arXiv:1108.1108 [cs.SC], 2011.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Omar E. Pol, Illustration of initial terms of A000217, A000290, A000326, A000384, A000566, A000567.
Leo Tavares, Illustration: Square Rays.
Leo Tavares, Illustration: Twin Rectangular Rays.
Leo Tavares, Illustration: Star Rows.
Leo Tavares, Illustration: Split Stars.
Eric Weisstein's World of Mathematics, Complete Bipartite Graph.
Eric Weisstein's World of Mathematics, Octagonal Number.
Eric Weisstein's World of Mathematics, Wiener Index.
Index to sequences related to polygonal numbers.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A014641, A014642, A014793, A014794, A001835, A016777, A045944, A093563 ((6, 1) Pascal, column m=2). A016921 (differences).
Cf. A005408 (the odd numbers).
Cf. A000290, A000217, A046092, A000384, A002378, A003154.

List([0..50], n -> n*(3*n-2)); # G. C. Greubel, Nov 15 2018

a000567 n = n * (3 * n - 2)  -- Reinhard Zumkeller, Dec 20 2012

[n*(3*n-2) : n in [0..50]]; // Wesley Ivan Hurt, Oct 10 2021

A000567 := proc(n)
    n*(3*n-2) ;
end proc:
seq(A000567(n), n=1..50) ;

Table[n (3 n - 2), {n, 0, 50}] (* Harvey P. Dale, May 06 2012 *)
Table[PolygonalNumber[RegularPolygon[8], n], {n, 0, 43}] (* Arkadiusz Wesolowski, Aug 27 2016 *)
PolygonalNumber[8, Range[0, 20]] (* Eric W. Weisstein, Sep 07 2017 *)
LinearRecurrence[{3, -3, 1}, {1, 8, 21}, {0, 20}] (* Eric W. Weisstein, Sep 07 2017 *)

a(n)=n*(3*n-2) \\ Charles R Greathouse IV, Jun 10 2011

vector(50, n, n--; n*(3*n-2)) \\ G. C. Greubel, Nov 15 2018

# 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 + 6, y + 6
A000567 = aList()
print([next(A000567) for i in range(49)]) # Peter Luschny, Aug 04 2019

[n*(3*n-2) for n in range(50)] # Gennady Eremin, Mar 10 2022

[n*(3*n-2) for n in range(50)] # G. C. Greubel, Nov 15 2018

a(n) = n*(3*n-2).
a(n) = (3n-2)*(3n-1)*(3n)/((3n-1) + (3n-2) + (3n)), i.e., (the product of three consecutive numbers)/(their sum). a(1) = 1*2*3/(1+2+3), a(2) = 4*5*6/(4+5+6), etc. - Amarnath Murthy, Aug 29 2002
E.g.f.: exp(x)*(x+3*x^2). - Paul Barry, Jul 23 2003
G.f.: x*(1+5*x)/(1-x)^3. Simon Plouffe in his 1992 dissertation
a(n) = Sum_{k=1..n} (5*n - 4*k). - Paul Barry, Sep 06 2005
a(n) = n + 6*A000217(n-1). - Floor van Lamoen, Oct 14 2005
a(n) = C(n+1,2) + 5*C(n,2).
Starting (1, 8, 21, 40, 65, ...) = binomial transform of [1, 7, 6, 0, 0, 0, ...]. - Gary W. Adamson, Apr 30 2008
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=0, a(1)=1, a(2)=8. - Jaume Oliver Lafont, Dec 02 2008
a(n) = A000578(n) - A007531(n). - Reinhard Zumkeller, Sep 18 2009
a(n) = a(n-1) + 6*n - 5 (with a(0)=0). - Vincenzo Librandi, Nov 20 2010
a(n) = 2*a(n-1) - a(n-2) + 6. - Ant King, Sep 01 2011
a(n) = A000217(n) + 5*A000217(n-1). - Vincenzo Librandi, Nov 20 2010
a(n) = (A185212(n) - 1) / 4. - Reinhard Zumkeller, Dec 20 2012
a(n) = A174709(6n). - Philippe Deléham, Mar 26 2013
a(n) = (2*n-1)^2 - (n-1)^2. - Ivan N. Ianakiev, Apr 10 2013
a(6*a(n) + 16*n + 1) = a(6*a(n) + 16*n) + a(6*n + 1). - Vladimir Shevelev, Jan 24 2014
a(0) = 0, a(n) = Sum_{k=0..n-1} A005408(A051162(n-1,k)), n >= 1. - L. Edson Jeffery, Jul 28 2014
Sum_{n>=1} 1/a(n) = (sqrt(3)*Pi + 9*log(3))/12 = 1.2774090575596367311949534921... . - Vaclav Kotesovec, Apr 27 2016
From Ilya Gutkovskiy, Jul 29 2016: (Start)
Inverse binomial transform of A084857.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(2*sqrt(3)) = A093766. (End)
a(n) = n * A016777(n-1) = A053755(n) - A000290(n+1). - Bruce J. Nicholson, Aug 10 2017
Product_{n>=2} (1 - 1/a(n)) = 3/4. - Amiram Eldar, Jan 21 2021
P(4k+4,n) = ((k+1)*n - k)^2 - (k*n - k)^2 where P(m,n) is the n-th m-gonal number (a generalization of the Apr 10 2013 formula, a(n) = (2*n-1)^2 - (n-1)^2). - Charlie Marion, Oct 07 2021
From Leo Tavares, Oct 31 2021: (Start)
a(n) = A000290(n) + 4*A000217(n-1). See Square Rays illustration.
a(n) = A000290(n) + A046092(n-1)
a(n) = A000384(n) + 2*A000217(n-1). See Twin Rectangular Rays illustration.
a(n) = A000384(n) + A002378(n-1)
a(n) = A003154(n) - A045944(n-1). See Star Rows illustration. (End)

Incorrect example removed by Joerg Arndt, Mar 11 2010

A080851 Square array of pyramidal numbers, read by antidiagonals.

Original entry on oeis.org

1, 1, 3, 1, 4, 6, 1, 5, 10, 10, 1, 6, 14, 20, 15, 1, 7, 18, 30, 35, 21, 1, 8, 22, 40, 55, 56, 28, 1, 9, 26, 50, 75, 91, 84, 36, 1, 10, 30, 60, 95, 126, 140, 120, 45, 1, 11, 34, 70, 115, 161, 196, 204, 165, 55, 1, 12, 38, 80, 135, 196, 252, 288, 285, 220, 66, 1, 13, 42, 90, 155, 231, 308, 372, 405, 385, 286, 78

Offset: 0

Paul Barry, Feb 21 2003

The first row contains the triangular numbers, which are really two-dimensional, but can be regarded as degenerate pyramidal numbers. - N. J. A. Sloane, Aug 28 2015

			Array begins (n>=0, k>=0):
1,  3,  6, 10,  15,  21,  28,  36,  45,   55, ... A000217
1,  4, 10, 20,  35,  56,  84, 120, 165,  220, ... A000292
1,  5, 14, 30,  55,  91, 140, 204, 285,  385, ... A000330
1,  6, 18, 40,  75, 126, 196, 288, 405,  550, ... A002411
1,  7, 22, 50,  95, 161, 252, 372, 525,  715, ... A002412
1,  8, 26, 60, 115, 196, 308, 456, 645,  880, ... A002413
1,  9, 30, 70, 135, 231, 364, 540, 765, 1045, ... A002414
1, 10, 34, 80, 155, 266, 420, 624, 885, 1210, ... A007584

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Numerous sequences in the database are to be found in the array. Rows include the pyramidal numbers A000217, A000292, A000330, A002411, A002412, A002413, A002414, A007584, A007585, A007586.
Columns include or are closely related to A017029, A017113, A017017, A017101, A016777, A017305. Diagonals include A006325, A006484, A002417.
Cf. A057145, A027660 (antidiagonal sums).
See A257199 for another version of this array.

vector(vector(poly_coeff(Taylor((1+kx)/(1-x)^4,x,11),x,n),n,0,11),k,-1,10) VECTOR(VECTOR(comb(k+2,2)+comb(k+2,3)n, k, 0, 11), n, 0, 11)

A080851 := proc(n,k)
    binomial(k+3,3)+(n-1)*binomial(k+2,3) ;
end proc:
seq( seq(A080851(d-k,k),k=0..d),d=0..12) ; # R. J. Mathar, Oct 01 2021

pyramidalFigurative[ ngon_, rank_] := (3 rank^2 + rank^3 (ngon - 2) - rank (ngon - 5))/6; Table[ pyramidalFigurative[n-k-1, k], {n, 4, 15}, {k, n-3}] // Flatten (* Robert G. Wilson v, Sep 15 2015 *)

T(n, k) = binomial(k+3, 3) + (n-1)*binomial(k+2, 3), corrected Oct 01 2021.
T(n, k) = T(n-1, k) + C(k+2, 3) = T(n-1, k) + k*(k+1)*(k+2)/6.
G.f. for rows: (1 + n*x)/(1-x)^4, n>=-1.
T(n,k) = sum_{j=1..k+1} A057145(n+2,j). - R. J. Mathar, Jul 28 2016

A237616 a(n) = n(n + 1)(5*n - 4)/2.

Original entry on oeis.org

0, 1, 18, 66, 160, 315, 546, 868, 1296, 1845, 2530, 3366, 4368, 5551, 6930, 8520, 10336, 12393, 14706, 17290, 20160, 23331, 26818, 30636, 34800, 39325, 44226, 49518, 55216, 61335, 67890, 74896, 82368, 90321, 98770, 107730, 117216, 127243, 137826, 148980, 160720

Offset: 0

Bruno Berselli, Feb 10 2014

Also 17-gonal (or heptadecagonal) pyramidal numbers.

This sequence is related to A226489 by 2*a(n) = n*A226489(n) - Sum_{i=0..n-1} A226489(i).

			After 0, the sequence is provided by the row sums of the triangle:
   1;
   2,  16;
   3,  32,  31;
   4,  48,  62,  46;
   5,  64,  93,  92,  61;
   6,  80, 124, 138, 122,  76;
   7,  96, 155, 184, 183, 152,  91;
   8, 112, 186, 230, 244, 228, 182, 106;
   9, 128, 217, 276, 305, 304, 273, 212, 121;
  10, 144, 248, 322, 366, 380, 364, 318, 242, 136; etc.,
where (r = row index, c = column index):
T(r,r) = T(c,c) = 15*r-14 and T(r,c) = T(r-1,c)+T(r,r) = (r-c+1)*T(r,r), with r>=c>0.

E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93 (fifteenth row of the table).

Bruno Berselli, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Pyramidal Number.
Index to sequences related to pyramidal numbers.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A051869, A104728.

Cf. sequences with formula n*(n+1)*(k*n-k+3)/6: A000217 (k=0), A000292 (k=1), A000330 (k=2), A002411 (k=3), A002412 (k=4), A002413 (k=5), A002414 (k=6), A007584 (k=7), A007585 (k=8), A007586 (k=9), A007587 (k=10), A050441 (k=11), A172073 (k=12), A177890 (k=13), A172076 (k=14), this sequence (k=15), A172078(k=16), A237617 (k=17), A172082 (k=18), A237618 (k=19), A172117(k=20), A256718 (k=21), A256716 (k=22), A256645 (k=23), A256646(k=24), A256647 (k=25), A256648 (k=26), A256649 (k=27), A256650(k=28).

List([0..40], n-> n*(n+1)*(5*n-4)/2); # G. C. Greubel, Aug 30 2019

Magma
```
[n*(n+1)*(5*n-4)/2: n in [0..40]];
```

I:=[0,1,18,66]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // Vincenzo Librandi, Feb 12 2014

seq(n*(n+1)*(5*n-4)/2, n=0..40); # G. C. Greubel, Aug 30 2019

Table[n(n+1)(5n-4)/2, {n, 0, 40}]
CoefficientList[Series[x (1+14x)/(1-x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)
LinearRecurrence[{4,-6,4,-1},{0,1,18,66},50] (* Harvey P. Dale, Jan 11 2015 *)

a(n)=n*(n+1)*(5*n-4)/2 \\ Charles R Greathouse IV, Sep 24 2015

[n*(n+1)*(5*n-4)/2 for n in (0..40)] # G. C. Greubel, Aug 30 2019

G.f.: x*(1 + 14*x)/(1 - x)^4.
For n>0, a(n) = Sum_{i=0..n-1} (n-i)*(15*i+1). More generally, the sequence with the closed form n*(n+1)*(k*n-k+3)/6 is also given by Sum_{i=0..n-1} (n-i)*(k*i+1) for n>0.
a(n) = A104728(A001844(n-1)) for n>0.
Sum_{n>=1} 1/a(n) = (2*sqrt(5*(5 + 2*sqrt(5)))*Pi + 10*sqrt(5)*arccoth(sqrt(5)) + 25*log(5) - 16)/72 = 1.086617842136293176... . - Vaclav Kotesovec, Dec 07 2016
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n >= 4. - Wesley Ivan Hurt, Dec 18 2020
E.g.f.: exp(x)*x*(2 + 16*x + 5*x^2)/2. - Elmo R. Oliveira, Aug 04 2025

A002419 4-dimensional figurate numbers: a(n) = (6n-2)binomial(n+2,3)/4.

Original entry on oeis.org

1, 10, 40, 110, 245, 476, 840, 1380, 2145, 3190, 4576, 6370, 8645, 11480, 14960, 19176, 24225, 30210, 37240, 45430, 54901, 65780, 78200, 92300, 108225, 126126, 146160, 168490, 193285, 220720, 250976, 284240, 320705, 360570, 404040, 451326, 502645, 558220

Offset: 1

N. J. A. Sloane

a(n) is the n-th antidiagonal sum of the convolution array A213761. - Clark Kimberling, Jul 04 2012
Convolution of A000027 with A000567 (excluding 0). - Bruno Berselli, Dec 07 2012
a(n) = the sum of all the ways of adding the k-tuples of A016777(0) to A016777(n-1). For n=4, the terms are 1,4,7,10 giving (1)+(4)+(7)+(10)=22; (1+4)+(4+7)+(7+10)=33; (1+4+7)+(4+7+10)=33; (1+4+7+10)=22; adding 22+33+33+22=110. - J. M. Bergot, Jun 26 2017
Also the number of chordless cycles in the (n+2)-crown graph. - Eric W. Weisstein, Jan 02 2018

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).

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992, arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
Eric Weisstein's World of Mathematics, Chordless Cycle.
Eric Weisstein's World of Mathematics, Crown Graph.
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).
Index to sequences related to pyramidal numbers.

Cf. A093563 ((6, 1) Pascal, column m=4).
Cf. A000027, A000567, A002414 (first differences), A016777, A080852, A213761.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.

List([1..40], n-> n*(n+1)*(n+2)*(3*n-1)/12); # G. C. Greubel, Jul 03 2019

/* A000027 convolved with A000567 (excluding 0): */ A000567:=func; [&+[(n-i+1)*A000567(i): i in [1..n]]: n in [1..40]]; // Bruno Berselli, Dec 07 2012

CoefficientList[Series[(1+5*x)/(1-x)^5, {x,0,40}], x] (* Vincenzo Librandi, Jun 20 2013 *)
LinearRecurrence[{5, -10, 10, -5, 1}, {1, 10, 40, 110, 245}, 40] (* Harvey P. Dale, Nov 30 2014 *)
Table[n(n+1)(n+2)(3n-1)/12, {n, 40}] (* Eric W. Weisstein, Jan 02 2018 *)
Table[Sum[2 x + 3 x^2 - 2 y, {x, 0, g}, {y, x, g}], {g, 1, 20}] (* Horst H. Manninger, Jun 20 2025 *)

a(n)=(3*n-1)*binomial(n+2,3)/2 \\ Charles R Greathouse IV, Sep 24 2015

A002419_list, m = [], [6, 1, 1, 1, 1]
for _ in range(10**2):
    A002419_list.append(m[-1])
    for i in range(4):
        m[i+1] += m[i] # Chai Wah Wu, Jan 24 2016

[n*(n+1)*(n+2)*(3*n-1)/12 for n in (1..40)] # G. C. Greubel, Jul 03 2019

a(n) = (3*n-1)*binomial(n+2, 3)/2.
G.f.: x*(1+5*x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation.
Sum_{n>=1} 1/a(n) = (-24+81*log(3) -9*Pi*sqrt(3))/14 = 1.143929... - R. J. Mathar, Mar 29 2011
a(n) = (3*n^4 + 8*n^3 + 3*n^2 - 2*n)/12. - Chai Wah Wu, Jan 24 2016
a(n) = A080852(6,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: x*(12 + 48*x + 26*x^2 + 3*x^3)*exp(x)/12. - G. C. Greubel, Jul 03 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*(3*sqrt(3)*Pi - 32*log(2) + 8)/7. - Amiram Eldar, Feb 11 2022

A007584 9-gonal (or enneagonal) pyramidal numbers: a(n) = n(n+1)(7*n-4)/6.

Original entry on oeis.org

0, 1, 10, 34, 80, 155, 266, 420, 624, 885, 1210, 1606, 2080, 2639, 3290, 4040, 4896, 5865, 6954, 8170, 9520, 11011, 12650, 14444, 16400, 18525, 20826, 23310, 25984, 28855, 31930, 35216, 38720, 42449, 46410, 50610, 55056, 59755, 64714, 69940, 75440, 81221

Offset: 0

N. J. A. Sloane, R. K. Guy

For n > 1, the digital roots of this sequence A010888(A007584(n)) form the purely periodic 27-cycle 1, 1, 7, 8, 2, 5, 6, 3, 3, 4, 4, 1, 2, 5, 8, 9, 6, 6, 7, 7, 4, 5, 8, 2, 3, 9, 9. For n > 1, the units digits of this sequence A010879(A007584(n)) form the purely periodic 20-cycle 1, 0, 4, 0, 5, 6, 0, 4, 5, 0, 6, 0, 9, 0, 0, 6, 5, 4, 0, 0. - Ant King, Oct 30 2012

Partial sums of A001106. - Joerg Arndt, Jun 10 2013

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.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index to sequences related to pyramidal numbers
Index entries for linear recurrences with constant coefficients, signature (4, -6, 4, -1).

Cf. A093564 ((7, 1) Pascal, column m=3).

Cf. similar sequences listed in A237616.

I:=[0, 1, 10, 34, 80]; [n le 5 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // Vincenzo Librandi, Jun 10 2013

a:=n->sum((n+j)^2-(n+j), j=0..n): seq(a(n)/2, n=0..30); # Zerinvary Lajos, May 26 2008

Table[n*(n+1)(7n-4)/6, {n, 0,100}] (* Vladimir Joseph Stephan Orlovsky, Jun 25 2009 *)
LinearRecurrence[{4,-6,4,-1},{1,10,34,80},30] (* Ant King, Oct 27 2012 *)
CoefficientList[Series[x (1 + 6 x) / (1 - x)^4, {x, 0, 50}], x] (* Vincenzo Librandi, Jun 10 2013 *)

A007584[n]:=n*(n+1)*(7*n-4)/6$
makelist(A007584[n],n,0,30); /* Martin Ettl, Oct 29 2012 */

a(n) = n*(n+1)*(7*n-4)/6; \\ Michel Marcus, Mar 04 2014

a(n) = (7*n-4)*binomial(n+1, 2)/3.
G.f.: x*(1+6*x)/(1-x)^4.
From Ant King, Oct 27 2012: (Start)
a(n) = a(n-1) + n*(7*n-5)/2.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 7.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
a(n) = (n+1)*(2*A001106(n)+n)/6.
a(n) = A000292(n) + 6*A000292(n-1).
a(n) = A002414(n) + A000292(n-1).
a(n) = A000217(n) + 7*A000292(n-1).
a(n) = binomial(n+2,3) + 6*binomial(n+1,3). (End)
a(n) = Sum_{i = 0..n-1} (n-i)*(7*i+1) for n>0. - Bruno Berselli, Feb 10 2014
a(n) = A080851(7,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: (x/6)*(6 + 24*x + 7*x^2)*exp(x). - G. C. Greubel, Oct 29 2017

A093563 (6,1)-Pascal triangle.

Original entry on oeis.org

1, 6, 1, 6, 7, 1, 6, 13, 8, 1, 6, 19, 21, 9, 1, 6, 25, 40, 30, 10, 1, 6, 31, 65, 70, 40, 11, 1, 6, 37, 96, 135, 110, 51, 12, 1, 6, 43, 133, 231, 245, 161, 63, 13, 1, 6, 49, 176, 364, 476, 406, 224, 76, 14, 1, 6, 55, 225, 540, 840, 882, 630, 300, 90, 15, 1, 6, 61, 280, 765, 1380

Offset: 0

Wolfdieter Lang, Apr 22 2004

The array F(6;n,m) gives in the columns m >= 1 the figurate numbers based on A016921, including the octagonal numbers A000567, (see the W. Lang link).
This is the sixth member, d=6, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal), A029653, A093560-2, for d=1..5.
This is an example of a Riordan triangle (see A093560 for a comment and A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group). Therefore the o.g.f. for the row polynomials p(n,x):=Sum_{m=0..n} a(n,m)*x^m is G(z,x)=(1+5*z)/(1-(1+x)*z).
The SW-NE diagonals give A022096(n-1) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with n=0 value 5. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 09 2013

			Triangle begins
  1;
  6,  1;
  6,  7,  1;
  6, 13,  8,  1;
  6, 19, 21,  9,  1;
  6, 25, 40, 30, 10,  1;
  ...

Kurt Hawlitschek, Johann Faulhaber 1580-1635, Veroeffentlichung der Stadtbibliothek Ulm, Band 18, Ulm, Germany, 1995, Ch. 2.1.4. Figurierte Zahlen.
Ivo Schneider: Johannes Faulhaber 1580-1635, Birkhäuser, Basel, Boston, Berlin, 1993, ch.5, pp. 109-122.

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
Wolfdieter Lang, First 10 rows and array of figurate numbers

Row sums: A005009(n-1), n>=1, 1 for n=0, alternating row sums are 1 for n=0, 5 for n=2 and 0 else.
Cf. A007318, A093564 (d=7), A228196, A228576.
The column sequences give for m=1..9: A016921, A000567 (octagonal), A002414, A002419, A051843, A027810, A034265, A054487, A055848.

a093563 n k = a093563_tabl !! n !! k
a093563_row n = a093563_tabl !! n
a093563_tabl = [1] : iterate
               (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [6, 1]
-- Reinhard Zumkeller, Aug 31 2014

lim = 11; s = Series[(1 + 5*x)/(1 - x)^(m + 1), {x, 0, lim}]; t = Table[ CoefficientList[s, x], {m, 0, lim}]; Flatten[ Table[t[[j - k + 1, k]], {j, lim + 1}, {k, j, 1, -1}]] (* Jean-François Alcover, Sep 16 2011, after g.f. *)

from math import comb, isqrt
def A093563(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),a:=n-comb(r+1,2))*(r+5*(r-a))//r if n else 1 # Chai Wah Wu, Nov 12 2024

a(n, m)=F(6;n-m, m) for 0<= m <= n, otherwise 0, with F(6;0, 0)=1, F(6;n, 0)=6 if n>=1 and F(6;n, m):= (6*n+m)*binomial(n+m-1, m-1)/m if m>=1.
Recursion: a(n, m)=0 if m>n, a(0, 0)= 1; a(n, 0)=6 if n>=1; a(n, m)= a(n-1, m) + a(n-1, m-1).
G.f. column m (without leading zeros): (1+5*x)/(1-x)^(m+1), m>=0.
T(n, k) = C(n, k) + 5*C(n-1, k). - Philippe Deléham, Aug 28 2005
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(6 + 13*x + 8*x^2/2! + x^3/3!) = 6 + 19*x + 40*x^2/2! + 70*x^3/3! + 110*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014

cp's OEIS Frontend

A344376 Numbers that are both octagonal numbers (A000567) and octagonal pyramidal numbers (A002414).

Views

Author

Keywords

Comments

Crossrefs

Programs

PARI

Extensions

A322653 Numbers that are sums of consecutive octagonal pyramidal numbers (A002414).

Views

Author

Keywords

Links

Crossrefs

A322856 Number of compositions (ordered partitions) of n into octagonal pyramidal numbers (A002414).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A218327 Even octagonal pyramidal numbers (A002414).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A000567 Octagonal numbers: n*(3*n-2). Also called star numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Python

Sage

Formula

Extensions

A080851 Square array of pyramidal numbers, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Derive

Maple

Mathematica

Formula

A237616 a(n) = n*(n + 1)*(5*n - 4)/2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

A000567 Octagonal numbers: n(3n-2). Also called star numbers.

A237616 a(n) = n(n + 1)(5*n - 4)/2.

A002419 4-dimensional figurate numbers: a(n) = (6n-2)binomial(n+2,3)/4.

A007584 9-gonal (or enneagonal) pyramidal numbers: a(n) = n(n+1)(7*n-4)/6.