cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-7 of 7 results.

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

Examples

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

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 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.

Crossrefs

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

Programs

  • GAP
    List([0..30], n-> n*(n+1)*(2*n+1)/6); # G. C. Greubel, Dec 31 2019
  • Haskell
    a000330 n = n * (n + 1) * (2 * n + 1) `div` 6
    a000330_list = scanl1 (+) a000290_list
    -- Reinhard Zumkeller, Nov 11 2012, Feb 03 2012
    
  • Magma
    [n*(n+1)*(2*n+1)/6: n in [0..50]]; // Wesley Ivan Hurt, Jun 28 2014
    
  • Magma
    [0] cat [((2*n+3)*Binomial(n+2,2))/3: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014
    
  • Maple
    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
  • Mathematica
    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 *)
  • Maxima
    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};
    
  • PARI
    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
    
  • Python
    a=lambda n: (n*(n+1)*(2*n+1))//6 # Indranil Ghosh, Jan 04 2017
    
  • Sage
    [n*(n+1)*(2*n+1)/6 for n in (0..30)] # G. C. Greubel, Dec 31 2019
    

Formula

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

Extensions

Partially edited by Joerg Arndt, Mar 11 2010

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

Original entry on oeis.org

2, 11, 32, 70, 130, 217, 336, 492, 690, 935, 1232, 1586, 2002, 2485, 3040, 3672, 4386, 5187, 6080, 7070, 8162, 9361, 10672, 12100, 13650, 15327, 17136, 19082, 21170, 23405, 25792, 28336, 31042, 33915, 36960, 40182, 43586, 47177, 50960, 54940
Offset: 1

Views

Author

Barry E. Williams, Dec 16 1999

Keywords

Comments

Partial sums of A005476.
a(n) is the dot product of the vectors of the first n positive integers and the next n integers. - Michel Marcus, Sep 02 2020

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 194-196.

Crossrefs

Programs

  • GAP
    a:=List([1..40],n->n*(n+1)*(5*n+1)/6);; Print(a); # Muniru A Asiru, Jan 01 2019
  • Magma
    [n*(n+1)*(5*n+1)/6 : n in [1..40]]; // Vincenzo Librandi, Jan 01 2019
    
  • Maple
    [n*(n+1)*(5*n+1)/6$n=1..40]; # Muniru A Asiru, Jan 01 2019
  • Mathematica
    Table[Range[x].Range[x+1,2x],{x,40}] (* or *) LinearRecurrence[{4,-6,4,-1},{2,11,32,70},40] (* Harvey P. Dale, Jun 01 2018 *)
  • PARI
    a(n) = n*(n+1)*(5*n+1)/6;
    

Formula

G.f.: x*(2+3*x)/(1-x)^4.
a(n) = A132121(n,1). - Reinhard Zumkeller, Aug 12 2007
a(n) = A000292(n) + A002412(n) = A000330(n) + A002411(n). - Omar E. Pol, Jan 11 2013
a(n) = Sum_{i=1..n} Sum_{j=1..n} i+min(i,j). - Enrique Pérez Herrero, Jan 15 2013
a(n) = Sum_{i=1..n} i*(n+i). - Charlie Marion, Apr 10 2013
Sum_{n>=1} 1/a(n) = 36 - 3*Pi*5^(3/4)*phi^(3/2)/4 - 15*sqrt(5)*log(phi)/4 - 75*log(5)/8 = 0.66131826232008423794478..., where phi = A001622 = (1 + sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Jun 01 2018
E.g.f.: exp(x)*x*(12 + 21*x + 5*x^2)/6. - Stefano Spezia, Feb 21 2024

Extensions

More terms from James Sellers, Jan 19 2000

A050409 Truncated square pyramid numbers: a(n) = Sum_{k = n..2*n} k^2.

Original entry on oeis.org

0, 5, 29, 86, 190, 355, 595, 924, 1356, 1905, 2585, 3410, 4394, 5551, 6895, 8440, 10200, 12189, 14421, 16910, 19670, 22715, 26059, 29716, 33700, 38025, 42705, 47754, 53186, 59015, 65255, 71920, 79024, 86581, 94605, 103110, 112110, 121619
Offset: 0

Views

Author

Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Dec 22 1999

Keywords

Crossrefs

Cf. A225144. - Bruno Berselli, Jun 06 2013
Cf. A045943: Sum_{k = n..2*n} k.
Cf. A304993: Sum_{k = n..2*n} k*(k+1)/2.

Programs

  • GAP
    List([0..40], n-> n*(n+1)*(14*n+1)/6); # G. C. Greubel, Oct 30 2019
  • Magma
    [&+[k^2: k in [n..2*n]]: n in [0..40]]; // Bruno Berselli, Feb 11 2011
    
  • Magma
    I:=[0, 5, 29, 86]; [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..40]]; // Vincenzo Librandi, Jun 22 2012
    
  • Maple
    seq(add((n+k)^2,k=0..n),n=0..40); # Zerinvary Lajos, Dec 01 2006
  • Mathematica
    LinearRecurrence[{4,-6,4,-1},{0,5,29,86},40] (* Vincenzo Librandi, Jun 22 2012 *)
    Table[(n(n+1)(14n+1))/6,{n,0,40}] (* Harvey P. Dale, Mar 08 2020 *)
  • PARI
    a(n)=sum(k=n,n+n,k^2)
    
  • PARI
    vector(40, n, n*(n-1)*(14*n-13)/6) \\ G. C. Greubel, Oct 30 2019
    
  • Sage
    [n*(n+1)*(14*n+1)/6 for n in (0..40)] # G. C. Greubel, Oct 30 2019
    

Formula

a(n) = n*(n+1)*(14*n+1)/6.
a(n) = A132121(n,4) for n>3. - Reinhard Zumkeller, Aug 12 2007
From Bruno Berselli, Feb 11 2011: (Start)
G.f.: x*(5+9*x)/(1-x)^4.
a(n) = A129371(2*n). (End)
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Vincenzo Librandi, Jun 22 2012
E.g.f.: x*(30 + 57*x + 14*x^2)*exp(x)/6. - G. C. Greubel, Oct 30 2019

A132111 Triangle read by rows: T(n,k) = n^2 + k*n + k^2, 0 <= k <= n.

Original entry on oeis.org

0, 1, 3, 4, 7, 12, 9, 13, 19, 27, 16, 21, 28, 37, 48, 25, 31, 39, 49, 61, 75, 36, 43, 52, 63, 76, 91, 108, 49, 57, 67, 79, 93, 109, 127, 147, 64, 73, 84, 97, 112, 129, 148, 169, 192, 81, 91, 103, 117, 133, 151, 171, 193, 217, 243, 100, 111, 124, 139, 156, 175, 196, 219
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 10 2007

Keywords

Comments

Permutation of A003136, the Loeschian numbers. [This is false - some terms are repeated, the first being 49. - Joerg Arndt, Dec 18 2015]
Row sums give A132112.
Central terms give A033582.
T(n,k+1) = T(n,k) + n + 2*k + 1;
T(n+1,k) = T(n,k) + 2*n + k + 1;
T(n+1,k+1) = T(n,k) + 3*(n+k+1);
T(n,0) = A000290(n);
T(n,1) = A002061(n+1) for n>0;
T(n,2) = A117950(n+1) for n>1;
T(n,n-2) = A056107(n-1) for n>1;
T(n,n-1) = A003215(n-1) for n>0;
T(n,n) = A033428(n).
T(n,k) is the norm N(alpha) of the integer alpha = n*1 - k*omega, where omega = exp(2*Pi*i/3) = (-1 + i*sqrt(3))/2 in the imaginary quadratic number field Q(sqrt(-3)): N = |alpha|^2 = (n + k/2)^2 + (3/4)*k^2 = n^2 + n*k + k^2 = T(n,k), with n >= 0, and k <= n. See also triangle A073254 for T(n,-k). - Wolfdieter Lang, Jun 13 2021

Examples

			From _Philippe Deléham_, Apr 16 2014: (Start)
Triangle begins:
   0;
   1,  3;
   4,  7,  12;
   9, 13,  19,  27;
  16, 21,  28,  37,  48;
  25, 31,  39,  49,  61,  75;
  36, 43,  52,  63,  76,  91, 108;
  49, 57,  67,  79,  93, 109, 127, 147;
  64, 73,  84,  97, 112, 129, 148, 169, 192;
  81, 91, 103, 117, 133, 151, 171, 193, 217, 243;
  ...
(End)
		

Crossrefs

Cf. A073254.

Programs

  • Mathematica
    Flatten[Table[n^2+k*n+k^2,{n,0,10},{k,0,n}]] (* Harvey P. Dale, Jun 10 2013 *)

A132124 a(n) = n*(n+1)*(8*n + 1)/6.

Original entry on oeis.org

0, 3, 17, 50, 110, 205, 343, 532, 780, 1095, 1485, 1958, 2522, 3185, 3955, 4840, 5848, 6987, 8265, 9690, 11270, 13013, 14927, 17020, 19300, 21775, 24453, 27342, 30450, 33785, 37355, 41168, 45232, 49555, 54145, 59010, 64158, 69597, 75335, 81380, 87740, 94423
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 12 2007

Keywords

Comments

Convolution of the sequences (0,3,5,0,0,0,...) and (binomial(n+3, 3)), n >= 0. - Emeric Deutsch, Aug 30 2007

Crossrefs

Programs

  • Maple
    seq((1/6)*n*(n+1)*(8*n+1),n=0..40); # Emeric Deutsch, Aug 30 2007
  • Mathematica
    a[n_] := n*(n + 1)*(8*n + 1)/6; Array[a, 42, 0] (* Amiram Eldar, May 20 2023 *)

Formula

a(n) = A132121(n,2) for n > 1.
G.f.: x*(3+5*x)/(1-x)^4. - Emeric Deutsch, Aug 30 2007
From Bruno Berselli, Nov 25 2010: (Start)
a(n) = n*A014105(n) - A016061(n-1), since A016061(n-1) = Sum_{k=0..n-1} A014105(k) (n > 0).
Also a(n) = A002412(n) + A006331(n) = A007585(n) + A002378(n). (End)
Sum_{n>=1} 1/a(n) = 54 - 24*(sqrt(2)+1)*Pi/7 - 24*(sqrt(2)+8)*log(2)/7 + 48*sqrt(2)*log(2-sqrt(2))/7. - Amiram Eldar, May 20 2023
E.g.f.: exp(x)*x*(18 + 33*x + 8*x^2)/6. - Stefano Spezia, Feb 21 2024

A132121 Triangle read by rows: T(n,k)=n*(n+1)*((3*k+2)*n+1)/6, 0<=k<=n.

Original entry on oeis.org

0, 1, 2, 5, 11, 17, 14, 32, 50, 68, 30, 70, 110, 150, 190, 55, 130, 205, 280, 355, 430, 91, 217, 343, 469, 595, 721, 847, 140, 336, 532, 728, 924, 1120, 1316, 1512, 204, 492, 780, 1068, 1356, 1644, 1932, 2220, 2508, 285, 690, 1095, 1500, 1905, 2310, 2715, 3120, 3525, 3930
Offset: 0

Views

Author

Reinhard Zumkeller, Aug 12 2007

Keywords

Comments

Row sums give A132122; central terms give A132123.
T(n,0) = A000330(n);
T(n,1) = A033994(n) for n>0;
T(n,2) = A132124(n) for n>1;
T(n,3) = A132112(n) for n>2;
T(n,4) = A050409(n) for n>3.

Examples

			0;
1, 2;
5, 11, 17;
14, 32, 50, 68;
30, 70, 110, 150, 190;
55, 130, 205, 280, 355, 430;
91, 217, 343, 469, 595, 721, 847;
140, 336, 532, 728, 924, 1120, 1316, 1512;
204, 492, 780, 1068, 1356, 1644, 1932, 2220, 2508;
		

Crossrefs

Programs

  • Maple
    A132121 := proc(n,k)
        n*(n+1)*((3*k+2)*n+1)/6 ;
    end proc:
    seq(seq(A132121(n,k),k=0..n),n=0..13) ; # R. J. Mathar, Feb 19 2020
  • Mathematica
    Table[n(n+1)((3k+2)n+1)/6,{n,0,9},{k,0,n}]//Flatten (* James C. McMahon, Mar 04 2025 *)

Formula

G.f.: Sum_{n>=0} Sum_{k>=0} T(n,k)*x^n*y^k = x*(x*y+1+x)/((1-x)^4*(1-y)^2). - R. J. Mathar, Jul 28 2016. Note that this generates a full array, not just the triangular subspace.

Extensions

a(53)-a(54) from James C. McMahon, Mar 04 2025

A254407 a(n) = n*(n+1)*(11*n +10)/6.

Original entry on oeis.org

0, 7, 32, 86, 180, 325, 532, 812, 1176, 1635, 2200, 2882, 3692, 4641, 5740, 7000, 8432, 10047, 11856, 13870, 16100, 18557, 21252, 24196, 27400, 30875, 34632, 38682, 43036, 47705, 52700, 58032, 63712, 69751, 76160, 82950, 90132, 97717, 105716, 114140, 123000
Offset: 0

Views

Author

Bruno Berselli, Jan 30 2015

Keywords

Comments

Similar sequences of the type m*P(s,m) - Sum_{i=1..m} P(s-1,i), where P(s,m) is the m-th s-gonal number:
s=3: A027480(n) = (n+1)*A000217(n+1) - Sum_{i=1..n+1} i;
s=4: A162148(n) = (n+1)*A000290(n+1) - Sum_{i=1..n+1} A000217(i);
s=5: A245301(n) = (n+1)*A000326(n+1) - Sum_{i=1..n+1} A000290(i);
s=6: A085788(n) = (n+1)*A000384(n+1) - Sum_{i=1..n+1} A000326(i);
s=7: a(n) = (n+1)*A000566(n+1) - Sum_{i=1..n+1} A000384(i).

Examples

			532 is the 7th term because A000566(7)=112 and Sum_{i=1..7} A000384(i)=252, therefore 7*112-252 = 532.
		

Crossrefs

Programs

  • Magma
    [n*(n+1)*(11*n+10)/6: n in [0..40]];
    
  • Maple
    A254407:= n-> n*(n+1)*(11*n+10)/6; seq(A254407(n), n=0..50); # G. C. Greubel, Mar 31 2021
  • Mathematica
    Table[n (n + 1) (11 n + 10)/6, {n, 0, 40}]
    Column[CoefficientList[Series[x (7 + 4 x) / (1 - x)^4, {x, 0, 60}], x]] (* Vincenzo Librandi, Jan 31 2015 *)
  • Maxima
    makelist(n*(n+1)*(11*n+10)/6, n, 0, 40);
  • PARI
    vector(40, n, n--; n*(n+1)*(11*n+10)/6)
    
  • Sage
    [n*(n+1)*(11*n+10)/6 for n in (0..40)]
    

Formula

G.f.: x*(7 + 4*x)/(1 - x)^4.
a(-n) = -A132112(n-1).
a(n) = Sum_{k=0..n} A011875(11*k+2).
Equivalently, partial sums of A254963.
E.g.f.: x*(42 + 54*x + 11*x^2)*exp(x)/6. - G. C. Greubel, Mar 31 2021
Showing 1-7 of 7 results.