A010879 -id:A010879 - cp's OEIS frontend

A000330 Square pyramidal numbers: a(n) = 0^2 + 1^2 + 2^2 + ... + n^2 = n(n+1)(2*n+1)/6.

0, 1, 5, 14, 30, 55, 91, 140, 204, 285, 385, 506, 650, 819, 1015, 1240, 1496, 1785, 2109, 2470, 2870, 3311, 3795, 4324, 4900, 5525, 6201, 6930, 7714, 8555, 9455, 10416, 11440, 12529, 13685, 14910, 16206, 17575, 19019, 20540, 22140, 23821, 25585, 27434, 29370

Offset: 0

N. J. A. Sloane

The sequence contains exactly one square greater than 1, namely 4900 (according to Gardner). - Jud McCranie, Mar 19 2001, Mar 22 2007 [This is a result from Watson. - Charles R Greathouse IV, Jun 21 2013] [See A351830 for further related comments and references.]
Number of rhombi in an n X n rhombus. - Matti De Craene (Matti.DeCraene(AT)rug.ac.be), May 14 2000
Number of acute triangles made from the vertices of a regular n-polygon when n is odd (cf. A007290). - Sen-Peng Eu, Apr 05 2001
Gives number of squares with sides parallel to the axes formed from an n X n square. In a 1 X 1 square, one is formed. In a 2 X 2 square, five squares are formed. In a 3 X 3 square, 14 squares are formed and so on. - Kristie Smith (kristie10spud(AT)hotmail.com), Apr 16 2002; edited by Eric W. Weisstein, Mar 05 2025
a(n-1) = B_3(n)/3, where B_3(x) = x(x-1)(x-1/2) is the third Bernoulli polynomial. - Michael Somos, Mar 13 2004
Number of permutations avoiding 13-2 that contain the pattern 32-1 exactly once.
Since 3*r = (r+1) + r + (r-1) = T(r+1) - T(r-2), where T(r) = r-th triangular number r*(r+1)/2, we have 3*r^2 = r*(T(r+1) - T(r-2)) = f(r+1) - f(r-1) ... (i), where f(r) = (r-1)*T(r) = (r+1)*T(r-1). Summing over n, the right hand side of relation (i) telescopes to f(n+1) + f(n) = T(n)*((n+2) + (n-1)), whence the result Sum_{r=1..n} r^2 = n*(n+1)*(2*n+1)/6 immediately follows. - Lekraj Beedassy, Aug 06 2004
Also as a(n) = (1/6)*(2*n^3 + 3*n^2 + n), n > 0: structured trigonal diamond numbers (vertex structure 5) (cf. A006003 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Number of triples of integers from {1, 2, ..., n} whose last component is greater than or equal to the others.
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Jun 12 2005
Sum of the first n positive squares. - Cino Hilliard, Jun 18 2007
Maximal number of cubes of side 1 in a right pyramid with a square base of side n and height n. - Pasquale CUTOLO (p.cutolo(AT)inwind.it), Jul 09 2007
If a 2-set Y and an (n-2)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 4-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
We also have the identity 1 + (1+4) + (1+4+9) + ... + (1+4+9+16+ ... + n^2) = n(n+1)(n+2)(n+(n+1)+(n+2))/36; ... and in general the k-fold nested sum of squares can be expressed as n(n+1)...(n+k)(n+(n+1)+...+(n+k))/((k+2)!(k+1)/2). - Alexander R. Povolotsky, Nov 21 2007
The terms of this sequence are coefficients of the Engel expansion of the following converging sum: 1/(1^2) + (1/1^2)*(1/(1^2+2^2)) + (1/1^2)*(1/(1^2+2^2))*(1/(1^2+2^2+3^2)) + ... - Alexander R. Povolotsky, Dec 10 2007
Convolution of A000290 with A000012. - Sergio Falcon, Feb 05 2008
Hankel transform of binomial(2*n-3, n-1) is -a(n). - Paul Barry, Feb 12 2008
Starting (1, 5, 14, 30, ...) = binomial transform of [1, 4, 5, 2, 0, 0, 0, ...]. - Gary W. Adamson, Jun 13 2008
Starting (1,5,14,30,...) = second partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = Sum_{i=0..n} binomial(n+2,i+2)*b(i), where b(i)=1,2,0,0,0,... - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
Convolution of A001477 with A005408: a(n) = Sum_{k=0..n} (2*k+1)*(n-k). - Reinhard Zumkeller, Mar 07 2009
Sequence of the absolute values of the z^1 coefficients of the polynomials in the GF1 denominators of A156921. See A157702 for background information. - Johannes W. Meijer, Mar 07 2009
The sequence is related to A000217 by a(n) = n*A000217(n) - Sum_{i=0..n-1} A000217(i) and this is the case d = 1 in the identity n^2*(d*n-d+2)/2 - Sum_{i=0..n-1} i*(d*i-d+2)/2 = n*(n+1)(2*d*n-2*d+3)/6, or also the case d = 0 in n^2*(n+2*d+1)/2 - Sum_{i=0..n-1} i*(i+2*d+1)/2 = n*(n+1)*(2*n+3*d+1)/6. - Bruno Berselli, Apr 21 2010, Apr 03 2012
a(n)/n = k^2 (k = integer) for n = 337; a(337) = 12814425, a(n)/n = 38025, k = 195, i.e., the number k = 195 is the quadratic mean (root mean square) of the first 337 positive integers. There are other such numbers -- see A084231 and A084232. - Jaroslav Krizek, May 23 2010
Also the number of moves to solve the "alternate coins game": given 2n+1 coins (n+1 Black, n White) set alternately in a row (BWBW...BWB) translate (not rotate) a pair of adjacent coins at a time (1 B and 1 W) so that at the end the arrangement shall be BBBBB..BW...WWWWW (Blacks separated by Whites). Isolated coins cannot be moved. - Carmine Suriano, Sep 10 2010
From J. M. Bergot, Aug 23 2011: (Start)
Using four consecutive numbers n, n+1, n+2, and n+3 take all possible pairs (n, n+1), (n, n+2), (n, n+3), (n+1, n+2), (n+1, n+3), (n+2, n+3) to create unreduced Pythagorean triangles. The sum of all six areas is 60*a(n+1).
Using three consecutive odd numbers j, k, m, (j+k+m)^3 - (j^3 + k^3 + m^3) equals 576*a(n) = 24^2*a(n) where n = (j+1)/2. (End)
From Ant King, Oct 17 2012: (Start)
For n > 0, the digital roots of this sequence A010888(a(n)) form the purely periodic 27-cycle {1, 5, 5, 3, 1, 1, 5, 6, 6, 7, 2, 2, 9, 7, 7, 2, 3, 3, 4, 8, 8, 6, 4, 4, 8, 9, 9}.
For n > 0, the units' digits of this sequence A010879(a(n)) form the purely periodic 20-cycle {1, 5, 4, 0, 5, 1, 0, 4, 5, 5, 6, 0, 9, 5, 0, 6, 5, 9, 0, 0}. (End)
Length of the Pisano period of this sequence mod n, n>=1: 1, 4, 9, 8, 5, 36, 7, 16, 27, 20, 11, 72, 13, 28, 45, 32, 17, 108, 19, 40, ... . - R. J. Mathar, Oct 17 2012
Sum of entries of n X n square matrix with elements min(i,j). - Enrique Pérez Herrero, Jan 16 2013
The number of intersections of diagonals in the interior of regular n-gon for odd n > 1 divided by n is a square pyramidal number; that is, A006561(2*n+1)/(2*n+1) = A000330(n-1) = (1/6)*n*(n-1)*(2*n-1). - Martin Renner, Mar 06 2013
For n > 1, a(n)/(2n+1) = A024702(m), for n such that 2n+1 = prime, which results in 2n+1 = A000040(m). For example, for n = 8, 2n+1 = 17 = A000040(7), a(8) = 204, 204/17 = 12 = A024702(7). - Richard R. Forberg, Aug 20 2013
A formula for the r-th successive summation of k^2, for k = 1 to n, is (2*n+r)*(n+r)!/((r+2)!*(n-1)!) (H. W. Gould). - Gary Detlefs, Jan 02 2014
The n-th square pyramidal number = the n-th triangular dipyramidal number (Johnson 12), which is the sum of the n-th + (n-1)-st tetrahedral numbers. E.g., the 3rd tetrahedral number is 10 = 1+3+6, the 2nd is 4 = 1+3. In triangular "dipyramidal form" these numbers can be written as 1+3+6+3+1 = 14. For "square pyramidal form", rebracket as 1+(1+3)+(3+6) = 14. - John F. Richardson, Mar 27 2014
Beukers and Top prove that no square pyramidal number > 1 equals a tetrahedral number A000292. - Jonathan Sondow, Jun 21 2014
Odd numbered entries are related to dissections of polygons through A100157. - Tom Copeland, Oct 05 2014
From Bui Quang Tuan, Apr 03 2015: (Start)
We construct a number triangle from the integers 1, 2, 3, ..., n as follows. The first column contains 2*n-1 integers 1. The second column contains 2*n-3 integers 2, ... The last column contains only one integer n. The sum of all the numbers in the triangle is a(n).
Here is an example with n = 5:
  1
  1  2
  1  2  3
  1  2  3  4
  1  2  3  4  5
  1  2  3  4
  1  2  3
  1  2
  1
(End)
The Catalan number series A000108(n+3), offset 0, gives Hankel transform revealing the square pyramidal numbers starting at 5, A000330(n+2), offset 0 (empirical observation). - Tony Foster III, Sep 05 2016; see Dougherty et al. link p. 2. - Andrey Zabolotskiy, Oct 13 2016
Number of floating point additions in the factorization of an (n+1) X (n+1) real matrix by Gaussian elimination as e.g. implemented in LINPACK subroutines sgefa.f or dgefa.f. The number of multiplications is given by A007290. - Hugo Pfoertner, Mar 28 2018
The Jacobi polynomial P(n-1,-n+2,2,3) or equivalently the sum of dot products of vectors from the first n rows of Pascal's triangle (A007318) with the up-diagonal Chebyshev T coefficient vector (1,3,2,0,...) (A053120) or down-diagonal vector (1,-7,32,-120,400,...) (A001794). a(5) = 1 + (1,1).(1,3) + (1,2,1).(1,3,2) + (1,3,3,1).(1,3,2,0) + (1,4,6,4,1).(1,3,2,0,0) = (1 + (1,1).(1,-7) + (1,2,1).(1,-7,32) + (1,3,3,1).(1,-7,32,-120) + (1,4,6,4,1).(1,-7,32,-120,400))*(-1)^(n-1) = 55. - Richard Turk, Jul 03 2018
Coefficients in the terminating series identity 1 - 5*n/(n + 4) + 14*n*(n - 1)/((n + 4)*(n + 5)) - 30*n*(n - 1)*(n - 2)/((n + 4)*(n + 5)*(n + 6)) + ... = 0 for n = 1,2,3,.... Cf. A002415 and A108674. - Peter Bala, Feb 12 2019
n divides a(n) iff n == +- 1 (mod 6) (see A007310). (See De Koninck reference.) Examples: a(11) = 506 = 11 * 46, and a(13) = 819 = 13 * 63. - Bernard Schott, Jan 10 2020
For n > 0, a(n) is the number of ternary words of length n+2 having 3 letters equal to 2 and 0 only occurring as the last letter. For example, for n=2, the length 4 words are 2221,2212,2122,1222,2220. - Milan Janjic, Jan 28 2020
Conjecture: Every integer can be represented as a sum of three generalized square pyramidal numbers. A related conjecture is given in A336205 corresponding to pentagonal case. A stronger version of these conjectures is that every integer can be expressed as a sum of three generalized r-gonal pyramidal numbers for all r >= 3. In here "generalized" means negative indices are included. - Altug Alkan, Jul 30 2020
The natural number y is a term if and only if y = a(floor((3 * y)^(1/3))). - Robert Israel, Dec 04 2024
Also the number of directed bishop moves on an n X n chessboard, where two moves are considered the same if one can be obtained from the other by a rotation of the board. Reflections are ignored. Equivalently, number of directed bishop moves on an n X n chessboard, where two moves are considered the same if one can be obtained from the other by an axial reflection of the board (horizontal or vertical). Rotations and diagonal reflections are ignored. - Hilko Koning, Aug 22 2025

			G.f. = x + 5*x^2 + 14*x^3 + 30*x^4 + 55*x^5 + 91*x^6 + 140*x^7 + 204*x^8 + ...

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 813.
A. H. Beiler, Recreations in the Theory of Numbers, Dover Publications, NY, 1964, p. 194.
A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 215,223.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 122, see #19 (3(1)), I(n); p. 155.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 47-49.
H. S. M. Coxeter, Polyhedral numbers, pp. 25-35 of R. S. Cohen, J. J. Stachel and M. W. Wartofsky, eds., For Dirk Struik: Scientific, historical and political essays in honor of Dirk J. Struik, Reidel, Dordrecht, 1974.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p.165).
J. M. De Koninck and A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 310, pp. 46-196, Ellipses, Paris, 2004.
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.
M. Gardner, Fractal Music, Hypercards and More, Freeman, NY, 1991, p. 293.
Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 293.
M. Holt, Math puzzles and games, Walker Publishing Company, 1977, p. 2 and p. 89.
Simon Singh, The Simpsons and Their Mathematical Secrets. London: Bloomsbury Publishing PLC (2013): 188.
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).
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 126.

Felix Fröhlich, Table of n, a(n) for n = 0..10000 (first 1001 terms from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
L. Ancora, Quadrature of the Parabola with the Square Pyramidal Number, Mondadori Education, Archimede 66, No. 3, 139-144 (2014).
Jack Anderson, Amy Woodall, and Alexandru Zaharescu, Arithmetic Polygons and Sums of Consecutive Squares, arXiv:2411.08398 [math.NT], 2024.
Ben Babcock and Adam Van Tuyl, Revisiting the spreading and covering numbers, arXiv preprint arXiv:1109.5847 [math.AC], 2011.
Joshua L. Bailey, Jr., A table to facilitate the fitting of certain logistic curves, Annals Math. Stat., 2 (1931), 355-359.
Joshua L. Bailey, A table to facilitate the fitting of certain logistic curves, Annals Math. Stat., 2 (1931), 355-359. [Annotated scanned copy]
Michael A. Bennett, Lucas' square pyramid problem revisited, Acta Arithmetica 105 (2002), 341-347.
Bruno Berselli, A description of the recursive method in Comments lines: website Matem@ticamente (in Italian).
Fritz Beukers and Jaap Top, On oranges and integral points on certain plane cubic curves, Nieuw Arch. Wiskd., IV (1988), Ser. 6, No. 3, 203-210.
Henry Bottomley, Illustration of initial terms.
Steve Butler and Pavel Karasik, A note on nested sums, J. Int. Seq. 13 (2010), 10.4.4, p=1 in first displayed equation page 4.
Bikash Chakraborty, Proof Without Words: Sums of Powers of Natural numbers, arXiv:2012.11539 [math.HO], 2020.
Robert Dawson, On Some Sequences Related to Sums of Powers, J. Int. Seq., Vol. 21 (2018), Article 18.7.6.
Alexander Diaz-Lopez, Pamela E. Harris, Erik Insko, and Darleen Perez-Lavin, Peaks Sets of Classical Coxeter Groups, arXiv preprint arXiv:1505.04479 [math.GR], 2015.
Anji Dong, Katerina Saettone, Kendra Song, and Alexandru Zaharescu, Cannonball Polygons with Multiplicities, arXiv:2507.18057 [math.NT], 2025. See p. 1.
Michael Dougherty, Christopher French, Benjamin Saderholm, and Wenyang Qian, Hankel Transforms of Linear Combinations of Catalan Numbers, Journal of Integer Sequences, Vol. 14 (2011), Article 11.5.1.
David Galvin and Courtney Sharpe, Independent set sequence of linear hyperpaths, arXiv:2409.15555 [math.CO], 2024. See p. 7.
Manfred Goebel, Rewriting Techniques and Degree Bounds for Higher Order Symmetric Polynomials, Applicable Algebra in Engineering, Communication and Computing (AAECC), Volume 9, Issue 6 (1999), 559-573.
T. Aaron Gulliver, Sequences from hexagonal pyramid of integers, International Mathematical Forum, Vol. 6, 2011, no. 17, p. 821-827.
Tian-Xiao He, Peter J.-S. Shiue, Zihan Nie, and Minghao Chen, Recursive sequences and Girard-Waring identities with applications in sequence transformation, Electronic Research Archive (2020) Vol. 28, No. 2, 1049-1062.
Milan Janjic, Two Enumerative Functions.
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
Milan Janjic and B. Petkovic, A Counting Function, arXiv preprint arXiv:1301.4550 [math.CO], 2013.
R. Jovanovic, First 2500 Pyramidal numbers.
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See pp. 9, 13-15, 24.
R. P. Loh, A. G. Shannon, and A. F. Horadam, Divisibility Criteria and Sequence Generators Associated with Fermat Coefficients, Preprint, 1980.
Toufik Mansour, Restricted permutations by patterns of type 2-1, arXiv:math/0202219 [math.CO], 2002.
Mircea Merca, A Special Case of the Generalized Girard-Waring Formula, J. Integer Sequences, Vol. 15 (2012), Article 12.5.7.
Cleve Moler, LINPACK subroutine sgefa.f, University of New Mexico, Argonne National Lab, 1978.
Michael Penn, Counting on a chessboard., YouTube video, 2021.
Claudio de J. Pita Ruiz V., Some Number Arrays Related to Pascal and Lucas Triangles, J. Int. Seq. 16 (2013) #13.5.7.
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.
Torsten Sillke, Square Counting.
Think Twice, Sum of n squares | explained visually |, video (2017).
Herman Tulleken, Polyominoes 2.2: How they fit together, (2019).
G. N. Watson, The problem of the square pyramid, Messenger of Mathematics 48 (1918), pp. 1-22.
Eric Weisstein's World of Mathematics, Faulhaber's Formula.
Eric Weisstein's World of Mathematics, Square Pyramidal Number.
Eric Weisstein's World of Mathematics, Square Tiling.
Eric Weisstein's World of Mathematics, Power Sum.
Wikipedia, Faulhaber's formula.
G. Xiao, Sigma Server, Operate on"n^2".
Index entries for "core" sequences.
Index to sequences related to pyramidal numbers.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
Index entries for two-way infinite sequences.

Cf. A000217, A000292, A000537, A005408, A006003, A006331, A033994, A033999, A046092, A050409, A050446, A050447, A060493, A100157, A132124, A132112, A156921, A157702, A258708, A351105, A351830.
Sums of 2 consecutive terms give A005900.
Column 0 of triangle A094414.
Column 1 of triangle A008955.
Right side of triangle A082652.
Row 2 of array A103438.
Partial sums of A000290.
Cf. similar sequences listed in A237616 and A254142.
Cf. |A084930(n, 1)|.
Cf. A253903 (characteristic function).
Cf. A034705 (differences of any two terms).

List([0..30], n-> n*(n+1)*(2*n+1)/6); # G. C. Greubel, Dec 31 2019

a000330 n = n * (n + 1) * (2 * n + 1) `div` 6
a000330_list = scanl1 (+) a000290_list
-- Reinhard Zumkeller, Nov 11 2012, Feb 03 2012

[n*(n+1)*(2*n+1)/6: n in [0..50]]; // Wesley Ivan Hurt, Jun 28 2014

[0] cat [((2*n+3)*Binomial(n+2,2))/3: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014

A000330 := n -> n*(n+1)*(2*n+1)/6;
a := n->(1/6)*n*(n+1)*(2*n+1): seq(a(n),n=0..53); # Emeric Deutsch
with(combstruct): ZL:=[st, {st=Prod(left, right), left=Set(U, card=r), right=Set(U, card=r), U=Sequence(Z, card>=1)}, unlabeled]: subs(r=1, stack): seq(count(subs(r=2, ZL), size=m*2), m=1..45) ; # Zerinvary Lajos, Jan 02 2008
nmax := 44; for n from 0 to nmax do fz(n) := product( (1-(2*m-1)*z)^(n+1-m) , m=1..n); c(n) := abs(coeff(fz(n),z,1)); end do: a := n-> c(n): seq(a(n), n=0..nmax); # Johannes W. Meijer, Mar 07 2009

Table[Binomial[w+2, 3] + Binomial[w+1, 3], {w, 0, 30}]
CoefficientList[Series[x(1+x)/(1-x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 30 2014 *)
Accumulate[Range[0,50]^2] (* Harvey P. Dale, Sep 25 2014 *)

A000330(n):=binomial(n+2,3)+binomial(n+1,3)$
makelist(A000330(n),n,0,20); /* Martin Ettl, Nov 12 2012 */

PARI
```
{a(n) = n * (n+1) * (2*n+1) / 6};
```

upto(n) = [x*(x+1)*(2*x+1)/6 | x<-[0..n]] \\ Cino Hilliard, Jun 18 2007, edited by M. F. Hasler, Jan 02 2024

a=lambda n: (n*(n+1)*(2*n+1))//6 # Indranil Ghosh, Jan 04 2017

[n*(n+1)*(2*n+1)/6 for n in (0..30)] # G. C. Greubel, Dec 31 2019

G.f.: x*(1+x)/(1-x)^4. - Simon Plouffe (in his 1992 dissertation: generating function for sequence starting at a(1))
E.g.f.: (x + 3*x^2/2 + x^3/3)*exp(x).
a(n) = n*(n+1)*(2*n+1)/6 = binomial(n+2, 3) + binomial(n+1, 3).
2*a(n) = A006331(n). - N. J. A. Sloane, Dec 11 1999
Can be extended to Z with a(n) = -a(-1-n) for all n in Z.
a(n) = A002492(n)/4. - Paul Barry, Jul 19 2003
a(n) = (((n+1)^4 - n^4) - ((n+1)^2 - n^2))/12. - Xavier Acloque, Oct 16 2003
From Alexander Adamchuk, Oct 26 2004: (Start)
a(n) = sqrt(A271535(n)).
a(n) = (Sum_{k=1..n} Sum_{j=1..n} Sum_{i=1..n} (i*j*k)^2)^(1/3). (End)
a(n) = Sum_{i=1..n} i*(2*n-2*i+1); sum of squares gives 1 + (1+3) + (1+3+5) + ... - Jon Perry, Dec 08 2004
a(n+1) = A000217(n+1) + 2*A000292(n). - Creighton Dement, Mar 10 2005
Sum_{n>=1} 1/a(n) = 6*(3-4*log(2)); Sum_{n>=1} (-1)^(n+1)*1/a(n) = 6*(Pi-3). - Philippe Deléham, May 31 2005
Sum of two consecutive tetrahedral (or pyramidal) numbers a(n) = A000292(n-1) + A000292(n). - Alexander Adamchuk, May 17 2006
Euler transform of length-2 sequence [ 5, -1 ]. - Michael Somos, Sep 04 2006
a(n) = a(n-1) + n^2. - Rolf Pleisch, Jul 22 2007
a(n) = A132121(n,0). - Reinhard Zumkeller, Aug 12 2007
a(n) = binomial(n, 2) + 2*binomial(n, 3). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009, corrected by M. F. Hasler, Jan 02 2024
a(n) = A168559(n) + 1 for n > 0. - Reinhard Zumkeller, Feb 03 2012
a(n) = Sum_{i=1..n} J_2(i)*floor(n/i), where J_2 is A007434. - Enrique Pérez Herrero, Feb 26 2012
a(n) = s(n+1, n)^2 - 2*s(n+1, n-1), where s(n, k) are Stirling numbers of the first kind, A048994. - Mircea Merca, Apr 03 2012
a(n) = A001477(n) + A000217(n) + A007290(n+2) + 1. - J. M. Bergot, May 31 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 2. - Ant King, Oct 17 2012
a(n) = Sum_{i = 1..n} Sum_{j = 1..n} min(i,j). - Enrique Pérez Herrero, Jan 15 2013
a(n) = A000217(n) + A007290(n+1). - Ivan N. Ianakiev, May 10 2013
a(n) = (A047486(n+2)^3 - A047486(n+2))/24. - Richard R. Forberg, Dec 25 2013
a(n) = Sum_{i=0..n-1} (n-i)*(2*i+1), with a(0) = 0. After 0, row sums of the triangle in A101447. - Bruno Berselli, Feb 10 2014
a(n) = n + 1 + Sum_{i=1..n+1} (i^2 - 2i). - Wesley Ivan Hurt, Feb 25 2014
a(n) = A000578(n+1) - A002412(n+1). - Wesley Ivan Hurt, Jun 28 2014
a(n) = Sum_{i = 1..n} Sum_{j = i..n} max(i,j). - Enrique Pérez Herrero, Dec 03 2014
a(n) = A055112(n)/6, see Singh (2013). - Alonso del Arte, Feb 20 2015
For n >= 2, a(n) = A028347(n+1) + A101986(n-2). - Bui Quang Tuan, Apr 03 2015
For n > 0: a(n) = A258708(n+3,n-1). - Reinhard Zumkeller, Jun 23 2015
a(n) = A175254(n) + A072481(n), n >= 1. - Omar E. Pol, Aug 12 2015
a(n) = A000332(n+3) - A000332(n+1). - Antal Pinter, Dec 27 2015
Dirichlet g.f.: zeta(s-3)/3 + zeta(s-2)/2 + zeta(s-1)/6. - Ilya Gutkovskiy, Jun 26 2016
a(n) = A080851(2,n-1). - R. J. Mathar, Jul 28 2016
a(n) = (A005408(n) * A046092(n))/12 = (2*n+1)*(2*n*(n+1))/12. - Bruce J. Nicholson, May 18 2017
12*a(n) = (n+1)*A001105(n) + n*A001105(n+1). - Bruno Berselli, Jul 03 2017
a(n) = binomial(n-1, 1) + binomial(n-1, 2) + binomial(n, 3) + binomial(n+1, 2) + binomial(n+1, 3). - Tony Foster III, Aug 24 2018
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Nathan Fox, Dec 04 2019
Let T(n) = A000217(n), the n-th triangular number.  Then a(n) = (T(n)+1)^2 + (T(n)+2)^2 + ... + (T(n)+n)^2 - (n+2)*T(n)^2. - Charlie Marion, Dec 31 2019
a(n) = 2*n - 1 - a(n-2) + 2*a(n-1). - Boštjan Gec, Nov 09 2023
a(n) = 2/(2*n)! * Sum_{j = 1..n} (-1)^(n+j) * j^(2*n+2) * binomial(2*n, n-j). Cf. A060493. - Peter Bala, Mar 31 2025

Partially edited by Joerg Arndt, Mar 11 2010

A010888 Digital root of n (repeatedly add the digits of n until a single digit is reached).

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 4, 5

Offset: 0

N. J. A. Sloane

This is sometimes also called the additive digital root of n.
n mod 9 (A010878) is a very similar sequence.
Partial sums are given by A130487(n-1) + n (for n > 0). - Hieronymus Fischer, Jun 08 2007
Decimal expansion of 13717421/111111111 is 0.123456789123456789123456789... with period 9. - Eric Desbiaux, May 19 2008
Decimal expansion of 13717421 / 1111111110 = 0.0[123456789] (periodic) - Daniel Forgues, Feb 27 2017
a(A005117(n)) < 9. - Reinhard Zumkeller, Mar 30 2010
My friend Jahangeer Kholdi has found that 19 is the smallest prime p such that for each number n, a(p*n) = a(n). In fact we have: a(m*n) = a(a(m)*a(n)) so all numbers with digital root 1 (numbers of the form 9k + 1) have this property. See comment lines of A017173. Also we have a(m+n) = a(a(m) + a(n)). - Farideh Firoozbakht, Jul 23 2010

			The digits of 37 are 3 and 7, and 3 + 7 = 10. And the digits of 10 are 1 and 0, and 1 + 0 = 1, so a(37) = 1.

Martin Gardner, Mathematics, Magic and Mystery, 1956.

N. J. A. Sloane, Table of n, a(n) for n = 0..10000
Eric Weisstein's World of Mathematics, Digitaddition
Eric Weisstein's World of Mathematics, Digital Root
Wikipedia, Vedic square
Index entries for Colombian or self numbers and related sequences
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,1).

Cf. A007953, A007954, A031347, A113217, A113218, A010878 (n mod 9), A010872, A010873, A010874, A010875, A010876, A010877, A010879, A004526, A002264, A002265, A002266, A017173, A031286 (additive persistence of n), (multiplicative digital root of n), A031346 (multiplicative persistence of n).

a010888 = until (< 10) a007953
-- Reinhard Zumkeller, Oct 17 2011, May 12 2011

[n eq 0 select 0 else 1+(n-1) mod 9: n in [0..110]]; // Bruno Berselli, Mar 18 2016

A010888 := n->if n=0 then 0 else ((n-1) mod 9) + 1; fi; # N. J. A. Sloane, Feb 20 2013

Join[{0}, Array[Mod[ # - 1, 9] + 1 &, 104]] (* Robert G. Wilson v, Jan 04 2006 *)
Join[Range[0, 1], Table[n - 9 Floor[(n - 1) / 9], {n, 2, 100}]] (* José de Jesús Camacho Medina, Nov 10 2014 *) (* Corrected by Vincenzo Librandi, Nov 11 2014 *)
Join[{0},LinearRecurrence[{0, 0, 0, 0, 0, 0, 0, 0, 1},{1, 2, 3, 4, 5, 6, 7, 8, 9},104]] (* Ray Chandler, Aug 26 2015 *)
Table[FixedPoint[Total[IntegerDigits[#, 10]] &, n], {n, 0, 104}] (* IWABUCHI Yu(u)ki, Jun 03 2016 *)

A010888(n)=if(n,(n-1)%9+1) \\ M. F. Hasler, Jan 04 2011

def A010888(n):
    return 1 + (n - 1) % 9 if n else 0 # Chai Wah Wu, Aug 23 2014, Apr 23 2023

0 :: List.fill(10)(1 to 9).flatten // Alonso del Arte, Feb 01 2020

If n = 0 then a(n) = 0; otherwise a(n) = (n reduced mod 9), but if the answer is 0 change it to 9.
Equivalently, if n = 0 then a(n) = 0, otherwise a(n) = (n - 1 reduced mod 9) + 1.
If the initial 0 term is ignored, the sequence is periodic with period 9.
From Hieronymus Fischer, Jun 08 2007: (Start)
a(n) = A010878(n-1) + 1 (for n > 0).
G.f.: g(x) = x*(Sum_{k = 0..8}(k+1)*x^k)/(1 - x^9). Also: g(x) = x(9x^10 - 10x^9 + 1)/((1 - x^9)(1 - x)^2). (End)
a(n) = n - 9*floor((n-1)/9), for n > 0. - José de Jesús Camacho Medina, Nov 10 2014

A000030 Initial digit of n.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1

Offset: 0

N. J. A. Sloane, Simon Plouffe

When n - a(n)*10^[log_10 n] >= 10^[(log_10 n) - 1], where [] denotes floor, or when n < 100 and 10|n, n is the concatenation of a(n) and A217657(n). - Reinhard Zumkeller, Oct 10 2012, improved by M. F. Hasler, Nov 17 2018, and corrected by Glen Whitney, Jul 01 2022

Equivalent definition: The initial a(0) = 0 is followed by each digit in S = {1,...,9} once. Thereafter, repeat 10 times each digit in S. Then, repeat 100 times each digit in S, etc.

			23 begins with a 2, so a(23) = 2.

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

David W. Wilson, Table of n, a(n) for n = 0..10000
A. Cobham, Uniform Tag Sequences, Mathematical Systems Theory, 6 (1972), 164-192.

Cf. A010879 (final digit of n).
Cf. A061681, A130571, A109453, A134010, A052038.
Cf. A002993, A089951, A002994, A143464, A098174, A098175, A072543, A072544, A073600, A073601, A037904. - Reinhard Zumkeller, Aug 17 2008

a000030 = until (< 10) (`div` 10) -- Reinhard Zumkeller, Feb 20 2012, Feb 11 2011

[Intseq(n)[#Intseq(n)]: n in [1..100]]; // Vincenzo Librandi, Nov 17 2018

A000030 := proc(n)
    if n = 0 then
        0;
    else
        convert(n,base,10) ;
        %[-1] ;
    end if;
end proc:
seq(A000030(n),n=0..200) ;# N. J. A. Sloane, Feb 10 2017

Join[{0},First[IntegerDigits[#]]&/@Range[90]] (* Harvey P. Dale, Mar 01 2011 *)
Table[Floor[n/10^(Floor[Log10[n]])], {n, 1, 50}] (* G. C. Greubel, May 16 2017 *)
Table[NumberDigit[n,IntegerLength[n]-1],{n,0,100}] (* Harvey P. Dale, Aug 29 2021 *)

a(n)=if(n<10,n,a(n\10)) \\ Mainly for illustration.

A000030(n)=n\10^logint(n+!n,10) \\ Twice as fast as a(n)=digits(n)[1]. Before digits() was added in PARI v.2.6.0 (2013), one could use, e.g., Vecsmall(Str(n))[1]-48. - M. F. Hasler, Nov 17 2018

def a(n): return int(str(n)[0])
print([a(n) for n in range(85)]) # Michael S. Branicky, Jul 01 2022

a(n) = [n / 10^([log_10(n)])] where [] denotes floor and log_10(n) is the logarithm is base 10. - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 07 2001

a(n) = k for k*10^j <= n < (k+1)*10^j for some j. - M. F. Hasler, Mar 23 2015

A010785 Repdigit numbers, or numbers whose digits are all equal.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 22, 33, 44, 55, 66, 77, 88, 99, 111, 222, 333, 444, 555, 666, 777, 888, 999, 1111, 2222, 3333, 4444, 5555, 6666, 7777, 8888, 9999, 11111, 22222, 33333, 44444, 55555, 66666, 77777, 88888, 99999, 111111, 222222, 333333, 444444, 555555, 666666

Offset: 0

N. J. A. Sloane

Complement of A139819. - David Wasserman, May 21 2008
Subsequence of A134336 and of A178403. - Reinhard Zumkeller, May 27 2010
Subsequence of A193460. - Reinhard Zumkeller, Jul 26 2011
Intersection of A009994 and A009996. - David F. Marrs, Sep 29 2018
Beiler (1964) called these numbers "monodigit numbers". The term "repdigit numbers" was used by Trigg (1974). - Amiram Eldar, Jan 21 2022

Albert H. Beiler, Recreations in the Theory of Numbers, Dover, New York, 1964, p. 83.

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Eric F. Bravo, Carlos A. Gómez and Florian Luca, Product of Consecutive Tribonacci Numbers With Only One Distinct Digit, J. Int. Seq., Vol. 22 (2019), Article 19.6.3.
Eric Fernando Bravo, On concatenations of Padovan and Perrin numbers, Math. Commun. (2023) Vol 28, 105-119.
Mahadi Ddamulira, Repdigits as sums of three balancing numbers, Mathematica Slovaca, (2019), hal-02405969.
Mahadi Ddamulira, Padovan numbers that are concatenations of two distinct repdigits, arXiv:2003.10705 [math.NT], 2020.
Mahadi Ddamulira, Tribonacci numbers that are concatenations of two repdigits, hal-02547159, Mathematics [math] / Number Theory [math.NT], 2020.
Mahadi Ddamulira, Padovan numbers that are concatenations of two distinct repdigits, Mathematica Slovaca, Vol. 71, No. 2 (2021), pp. 275-284.
Bart Goddard and Jeremy Rouse, Sum of two repdigits a square, arXiv:1607.06681 [math.NT], 2016. Mentions this sequence.
Bir Kafle, Florian Luca and Alain Togbé, Triangular Repblocks, Fibonacci Quart., Vol. 56, No. 4 (2018), pp. 325-328.
Bir Kafle, Florian Luca and Alain Togbé, Pentagonal and heptagonal repdigits, Annales Mathematicae et Informaticae, Vol. 52 (2020), pp. 137-145.
Benedict Vasco Normenyo, Bir Kafle, and Alain Togbé, Repdigits as Sums of Two Fibonacci Numbers and Two Lucas Numbers, Integers, Vol. 19 (2019), Article A55.
Salah Eddine Rihane and Alain Togbé, Repdigits as products of consecutive Padovan or Perrin numbers, Arab. J. Math., Vol. 10 (2021), pp. 469-480.
Charles W. Trigg, Infinite sequences of palindromic triangular numbers, The Fibonacci Quarterly, Vol. 12, No. 2 (1974), pp. 209-212.
Eric Weisstein's World of Mathematics, Repdigit.
Wikipedia, Repdigit.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,11,0,0,0,0,0,0,0,0,-10).

Cf. A000005, A009994, A009996, A010879, A037904, A047842, A059995, A065444, A134336, A139819, A151949, A178401, A178403, A180410, A193459, A193460, A202022, A227362, A244112.

a010785 n = a010785_list !! n
a010785_list = 0 : r [1..9] where
   r (x:xs) = x : r (xs ++ [10*x + x `mod` 10])
-- Reinhard Zumkeller, Jul 26 2011

[(n-9*Floor((n-1)/9))*(10^Floor((n+8)/9)-1)/9: n in [0..50]]; // Vincenzo Librandi, Nov 10 2014

A010785 := proc(n)
    (n-9*floor(((n-1)/9)))*((10^(floor(((n+8)/9)))-1)/9) ;
end proc:
seq(A010785(n), n = 0 .. 100); # Robert Israel, Nov 09 2014

fQ[n_]:=Module[{id=IntegerDigits[n]}, Length[Union[id]]==1]; Select[Range[0,10000], fQ] (* Vladimir Joseph Stephan Orlovsky, Dec 29 2010 *)
Union[FromDigits/@Flatten[Table[PadRight[{},i,n],{n,0,9},{i,6}],1]] (* or *) LinearRecurrence[{0,0,0,0,0,0,0,0,11,0,0,0,0,0,0,0,0,-10}, {0,1,2,3,4,5,6,7,8,9,11,22,33,44,55,66,77,88},40] (* Harvey P. Dale, Dec 28 2011 *)
Union@ Flatten@ Table[k (10^n - 1)/9, {k, 0, 9}, {n, 6}] (* Robert G. Wilson v, Oct 09 2014 *)
Table[(n - 9 Floor[(n-1)/9]) (10^Floor[(n+8)/9] - 1)/9, {n, 0, 50}] (* José de Jesús Camacho Medina, Nov 06 2014 *)

a(n)=10^((n+8)\9)\9*((n-1)%9+1) \\ Charles R Greathouse IV, Jun 15 2011

nxt(n,t=n%10)=if(t<9,n*(t+1),n*10+9)\t \\ Yields the term a(k+1) following a given term a(k)=n. M. F. Hasler, Jun 24 2016

is(n)={1==#Set(digits(n))}
inv(n) = 9*#Str(n) + n%10 - 9 \\ David A. Corneth, Jun 24 2016

def a(n): return 0 if n == 0 else int(str((n-1)%9+1)*((n-1)//9+1))
print([a(n) for n in range(55)]) # Michael S. Branicky, Dec 29 2021

print([0]+[int(d*r) for r in range(1, 7) for d in "123456789"]) # Michael S. Branicky, Dec 29 2021

# without string operations
def a(n): return 0 if n == 0 else (10**((n-1)//9+1)-1)//9*((n-1)%9+1)
print([a(n) for n in range(55)]) # Michael S. Branicky, Nov 03 2023

A037904(a(n)) = 0. - Reinhard Zumkeller, Dec 14 2007
A178401(a(n)) > 0. - Reinhard Zumkeller, May 27 2010
From Reinhard Zumkeller, Jul 26 2011: (Start)
For n > 0: A193459(a(n)) = A000005(a(n)).
for n > 10: a(n) mod 10 = floor(a(n)/10) mod 10.
A010879(n) = A010879(A059995(n)). (End)
A202022(a(n)) = 1. - Reinhard Zumkeller, Dec 09 2011
a(0)=0, a(1)=1, a(2)=2, a(3)=3, a(4)=4, a(5)=5, a(6)=6, a(7)=7, a(8)=8, a(9)=9, a(10)=11, a(11)=22, a(12)=33, a(13)=44, a(14)=55, a(15)=66, a(16)=77, a(17)=88, a(n) = 11*a(n-9) - 10*a(n-18). - Harvey P. Dale, Dec 28 2011
A151949(a(n)) = 0; A180410(a(n)) = A227362(a(n)). - Reinhard Zumkeller, Jul 09 2013
a(n) = (n - 9*floor((n-1)/9))*(10^floor((n+8)/9) - 1)/9. - José de Jesús Camacho Medina, Nov 06 2014
G.f.: x*(1+2*x+3*x^2+4*x^3+5*x^4+6*x^5+7*x^6+8*x^7+9*x^8)/((1-x^9)*(1-10*x^9)). - Robert Israel, Nov 09 2014
A047842(a(n)) = A244112(a(n)). - Reinhard Zumkeller, Nov 11 2014
Sum_{n>=1} 1/a(n) = (7129/2520) * A065444 = 3.11446261209177581335... - Amiram Eldar, Jan 21 2022

Name clarified by Jon E. Schoenfield, Nov 10 2023

A002411 Pentagonal pyramidal numbers: a(n) = n^2*(n+1)/2.

Original entry on oeis.org

0, 1, 6, 18, 40, 75, 126, 196, 288, 405, 550, 726, 936, 1183, 1470, 1800, 2176, 2601, 3078, 3610, 4200, 4851, 5566, 6348, 7200, 8125, 9126, 10206, 11368, 12615, 13950, 15376, 16896, 18513, 20230, 22050, 23976, 26011, 28158, 30420, 32800, 35301, 37926, 40678

Offset: 0

N. J. A. Sloane

a(n) = n^2(n+1)/2 is half the number of colorings of three points on a line with n+1 colors. - R. H. Hardin, Feb 23 2002
Sum of n smallest multiples of n. - Amarnath Murthy, Sep 20 2002
a(n) = number of (n+6)-bit binary sequences with exactly 7 1's none of which is isolated. A 1 is isolated if its immediate neighbor(s) are 0. - David Callan, Jul 15 2004
Also as a(n) = (1/6)*(3*n^3+3*n^2), n > 0: structured trigonal prism numbers (cf. A100177 - structured prisms; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Nov 18 2005
If Y is a 3-subset of an n-set X then, for n >= 5, a(n-4) is the number of 5-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
a(n-1), n >= 2, is the number of ways to have n identical objects in m=2 of altogether n distinguishable boxes (n-2 boxes stay empty). - Wolfdieter Lang, Nov 13 2007
a(n+1) is the convolution of (n+1) and (3n+1). - Paul Barry, Sep 18 2008
The number of 3-character strings from an alphabet of n symbols, if a string and its reversal are considered to be the same.
Partial sums give A001296. - Jonathan Vos Post, Mar 26 2011
a(n-1):=N_1(n), n >= 1, is the number of edges of n planes in generic position in three-dimensional space. See a comment under A000125 for general arrangement. Comment to Arnold's problem 1990-11, see the Arnold reference, p.506. - Wolfdieter Lang, May 27 2011
Partial sums of pentagonal numbers A000326. - Reinhard Zumkeller, Jul 07 2012
From Ant King, Oct 23 2012: (Start)
For n > 0, the digital roots of this sequence A010888(A002411(n)) form the purely periodic 9-cycle {1,6,9,4,3,9,7,9,9}.
For n > 0, the units' digits of this sequence A010879(A002411(n)) form the purely periodic 20-cycle {1,6,8,0,5,6,6,8,5,0,6,6,3,0,0,6,1,8,0,0}.
(End)
a(n) is the number of inequivalent ways to color a path graph having 3 nodes using at most n colors.  Note, here there is no restriction on the color of adjacent nodes as in the above comment by R. H. Hardin (Feb 23 2002).  Also, here the structures are counted up to graph isomorphism, where as in the above comment the "three points on a line" are considered to be embedded in the plane. - Geoffrey Critzer, Mar 20 2013
After 0, row sums of the triangle in A101468. - Bruno Berselli, Feb 10 2014
Latin Square Towers: Take a Latin square of order n, with symbols from 1 to n, and replace each symbol x with a tower of height x. Then the total number of unit cubes used is a(n). - Arun Giridhar, Mar 29 2015
This is the case k = n+4 of b(n,k) = n*((k-2)*n-(k-4))/2, which is the n-th k-gonal number. Therefore, this is the 3rd upper diagonal of the array in A139600. - Luciano Ancora, Apr 11 2015
For n > 0, a(n) is the number of compositions of n+7 into n parts avoiding the part 2. - Milan Janjic, Jan 07 2016
Also the Wiener index of the n-antiprism graph. - Eric W. Weisstein, Sep 07 2017
For n > 0, a(2n+1) is the number of non-isomorphic 5C_m-snakes, where m = 2n+1 or m = 2n (for n >= 2). A kC_n-snake is a connected graph in which the k >= 2 blocks are isomorphic to the cycle C_n and the block-cutpoint graph is a path. - Christian Barrientos, May 15 2019
For n >= 1, a(n-1) is the number of 0°- and 45°-tilted squares that can be drawn by joining points in an n X n lattice. - Paolo Xausa, Apr 13 2021
a(n) is the number of all possible products of n rolls of a six-sided die. This can be easily seen by the recursive formula a(n) = a(n - 1) + 2 * binomial(n, 2) + binomial(n + 1, 2). - Rafal Walczak, Jun 15 2024
a(n) is the number of all triples consisting of nonnegative integers smaller than n such that the sum of the first two integers is less than n. - Ruediger Jehn, Aug 17 2025

			a(3)=18 because 4 identical balls can be put into m=2 of n=4 distinguishable boxes in binomial(4,2)*(2!/(1!*1!) + 2!/2!) = 6*(2+1) = 18 ways. The m=2 part partitions of 4, namely (1,3) and (2,2), specify the filling of each of the 6 possible two-box choices. - _Wolfdieter Lang_, Nov 13 2007

V. I. Arnold (ed.), Arnold's Problems, Springer, 2004, comments on Problem 1990-11 (p. 75), pp. 503-510. Numbers N_1.
Christian Barrientos, Graceful labelings of cyclic snakes, Ars Combin., Vol. 60 (2001), pp. 85-96.
Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see p. 166, Table 10.4/I/5).
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.
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).

T. D. Noe, Table of n, a(n) for n = 0..1000
Somaya Barati, Beáta Bényi, Abbas Jafarzadeh and Daniel Yaqubi, Mixed restricted Stirling numbers, arXiv:1812.02955 [math.CO], 2018.
Phyllis Chinn and Silvia Heubach, Integer Sequences Related to Compositions without 2's, J. Integer Seq., Vol. 6 (2003), Article 03.2.3.
Robert Coquereaux and Jean-Bernard Zuber, Counting partitions by genus. II. A compendium of results, arXiv:2305.01100 [math.CO], 2023. See p. 17.
Lancelot Hogben, Choice and Chance by Cardpack and Chessboard, Vol. 1, Max Parrish and Co, London, 1950, p. 36.
Milan Janjic, Binomial Coefficients and Enumeration of Restricted Words, Journal of Integer Sequences, Vol 19 (2016), Article 16.7.3.
C. Krishnamachaki, The operator (xD)^n, J. Indian Math. Soc., Vol. 15 (1923), pp. 3-4. [Annotated scanned copy]
S. M. Losanitsch, Die Isomerie-Arten bei den Homologen der Paraffin-Reihe, Chem. Ber., Vol. 30 (1897), pp. 1917-1926.
S. M. Losanitsch, Die Isomerie-Arten bei den Homologen der Paraffin-Reihe, Chem. Ber., Vol. 30 (1897), pp. 1917-1926. (Annotated scanned copy)
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, Pentagonal Pyramidal Number.
Eric Weisstein's World of Mathematics, Wiener Index.
Index entries for two-way infinite sequences.
Index to sequences related to polygonal numbers.
Index to sequences related to pyramidal numbers.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A001296, A011379, A015223, A015224, A014799, A014800, A127739, A132118, A139600.
A006002(n) = -a(-1-n).
a(n) = A093560(n+2, 3), (3, 1)-Pascal column.
A row or column of A132191.
Second column of triangle A103371.
Cf. similar sequences listed in A237616.

List([0..45], n->n^2*(n+1)/2); # Muniru A Asiru, Feb 19 2018

a002411 n = n * a000217 n  -- Reinhard Zumkeller, Jul 07 2012

[n^2*(n+1)/2: n in [0..40]]; // Wesley Ivan Hurt, May 25 2014

Maple
```
seq(n^2*(n+1)/2, n=0..40);
```

Table[n^2 (n + 1)/2, {n, 0, 40}]
LinearRecurrence[{4, -6, 4, -1}, {0, 1, 6, 18}, 50] (* Harvey P. Dale, Oct 20 2011 *)
Nest[Accumulate, Range[1, 140, 3], 2] (* Vladimir Joseph Stephan Orlovsky, Jan 21 2012 *)
CoefficientList[Series[x (1 + 2 x) / (1 - x)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Jan 08 2016 *)

PARI
```
a(n)=n^2*(n+1)/2
```

concat(0, Vec(x*(1+2*x)/(1-x)^4 + O(x^100))) \\ Altug Alkan, Jan 07 2016

Average of n^2 and n^3.
G.f.: x*(1+2*x)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
a(n) = n*Sum_{k=0..n} (n-k) = n*Sum_{k=0..n} k. - Paul Barry, Jul 21 2003
a(n) = n*A000217(n). - Xavier Acloque, Oct 27 2003
a(n) = (1/2)*Sum_{j=1..n} Sum_{i=1..n} (i+j) = (1/2)*(n^2+n^3) = (1/2)*A011379(n). - Alexander Adamchuk, Apr 13 2006
Row sums of triangle A127739, triangle A132118; and binomial transform of [1, 5, 7, 3, 0, 0, 0, ...] = (1, 6, 18, 40, 75, ...). - Gary W. Adamson, Aug 10 2007
G.f.: x*F(2,3;1;x). - Paul Barry, Sep 18 2008
Sum_{j>=1} 1/a(j) = hypergeom([1, 1, 1], [2, 3], 1) = -2 + 2*zeta(2) = A195055 - 2. - Stephen Crowley, Jun 28 2009
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4); a(0)=0, a(1)=1, a(2)=6, a(3)=18. - Harvey P. Dale, Oct 20 2011
From Ant King, Oct 23 2012: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 3.
a(n) = (n+1)*(2*A000326(n)+n)/6 = A000292(n) + 2*A000292(n-1).
a(n) = A000330(n)+A000292(n-1) = A000217(n) + 3*A000292(n-1).
a(n) = binomial(n+2,3) + 2*binomial(n+1,3).
(End)
a(n) = (A000330(n) + A002412(n))/2 = (A000292(n) + A002413(n))/2. - Omar E. Pol, Jan 11 2013
a(n) = (24/(n+3)!)*Sum_{j=0..n} (-1)^(n-j)*binomial(n,j)*j^(n+3). - Vladimir Kruchinin, Jun 04 2013
Sum_{n>=1} a(n)/n! = (7/2)*exp(1). - Richard R. Forberg, Jul 15 2013
E.g.f.: x*(2 + 4*x + x^2)*exp(x)/2. - Ilya Gutkovskiy, May 31 2016
From R. J. Mathar, Jul 28 2016: (Start)
a(n) = A057145(n+4,n).
a(n) = A080851(3,n-1). (End)
For n >= 1, a(n) = (Sum_{i=1..n} i^2) + Sum_{i=0..n-1} i^2*((i+n) mod 2). - Paolo Xausa, Apr 13 2021
a(n) = Sum_{k=1..n} GCD(k,n) * LCM(k,n). - Vaclav Kotesovec, May 22 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = 2 + Pi^2/6 - 4*log(2). - Amiram Eldar, Jan 03 2022

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

N. J. A. Sloane

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

			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

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

T. D. Noe and William A. Tedeschi, Table of n, a(n) for n = 0..10000 (first 1000 terms computed by T. D. Noe)
Abdullah Atmaca and A. Yavuz Oruç, On the size of two families of unlabeled bipartite graphs, AKCE International Journal of Graphs and Combinatorics, Vo. 16, No. 2 (2019), pp. 222-229.
Bruno Berselli, A description of the transform in Comments lines: website Matem@ticamente (in Italian).
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.
Luis Verde-Star, A Matrix Approach to Generalized Delannoy and Schröder Arrays, J. Int. Seq., Vol. 24 (2021), Article 21.4.1.
Eric Weisstein's World of Mathematics, Hexagonal Pyramidal Number.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Bisection of A002623. Equals A000578(n) - A000330(n-1).
Cf. A000292, A016061.
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.

List([0..40],n->n*(n+1)*(4*n-1)/6); # Muniru A Asiru, Mar 18 2019

[n*(n+1)*(4*n-1)/6: n in [0..40]]; // Vincenzo Librandi, Nov 28 2015

seq(sum(i*(2*k-i), i=1..k), k=0..100); # Wesley Ivan Hurt, Sep 25 2013

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

A002412(n):=n*(n+1)*(4*n-1)/6$ makelist(A002412(n),n,0,20); /* Martin Ettl, Dec 12 2012 */

v=vector(40,i,(i*(i+1))\2); s=0; print1(s","); forstep(i=1,40,2,s+=v[i]; print1(s","))

print([n*(n+1)*(4*n-1)//6 for n in range(40)]) # Michael S. Branicky, Mar 28 2022

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

cp's OEIS Frontend

A130488 a(n) = Sum_{k=0..n} (k mod 10) (Partial sums of A010879).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Python

Sage

Formula

A000330 Square pyramidal numbers: a(n) = 0^2 + 1^2 + 2^2 + ... + n^2 = n*(n+1)*(2*n+1)/6.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Magma

Maple

Mathematica

Maxima

PARI

PARI

Python

Sage

Formula

Extensions

A010888 Digital root of n (repeatedly add the digits of n until a single digit is reached).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Scala

Formula

A000030 Initial digit of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Formula

A010785 Repdigit numbers, or numbers whose digits are all equal.

Views

Author

Keywords

A000330 Square pyramidal numbers: a(n) = 0^2 + 1^2 + 2^2 + ... + n^2 = n(n+1)(2*n+1)/6.