cp's OEIS Frontend

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

Showing 1-10 of 51 results. Next

A317714 Chessboard rectangles sequence (see Comments), also A037270 interleaved with A163102.

Original entry on oeis.org

0, 0, 1, 2, 10, 18, 45, 72, 136, 200, 325, 450, 666, 882, 1225, 1568, 2080, 2592, 3321, 4050, 5050, 6050, 7381, 8712, 10440, 12168, 14365, 16562, 19306, 22050, 25425, 28800, 32896, 36992, 41905, 46818, 52650, 58482, 65341, 72200, 80200, 88200, 97461, 106722, 117370
Offset: 1

Views

Author

Ivan N. Ianakiev, Aug 05 2018

Keywords

Comments

Take a chessboard of n X n unit squares in which the a1 square is black. a(n) is the number of composite rectangles of p X q unit squares whose vertices are covered by black unit squares (1 < p <= n, 1 < q <= n).

Examples

			In a 4 X 4 chessboard there are two such rectangles (for both p = q = 3) and the coordinates of their lower left vertices are a1 and b2. Therefore, a(4) = 2.
		

Crossrefs

Programs

  • Magma
    [(5-5*(-1)^n-12*n+12*(-1)^n*n+14*n^2-6*(-1)^n*n^2-8*n^3+2*n^4)/64: n in [1..50]]; // Vincenzo Librandi, Aug 05 2018
    
  • Mathematica
    CoefficientList[Series[-((x^2 (1+4 x^2+x^4))/((-1+x)^5 (1+x)^3)),{x,0,44}],x]
    LinearRecurrence[{2, 2, -6, 0, 6, -2, -2, 1}, {0, 0, 1, 2, 10, 18, 45, 72}, 80] (* Vincenzo Librandi, Aug 06 2018 *)
  • PARI
    a(n) = sum(i = 1, n-1, floor(i/2)^3); \\ Jinyuan Wang, Aug 12 2019
  • Python
    n, a = 0, 0
    while n < 10:
        print(n,a)
    n, a = n+1, a+((n+1)//2)**3 # A.H.M. Smeets, Aug 09 2019
    

Formula

a(n) = 2*a(n-1) + 2*a(n-2) - 6*a(n-3) + 6*a(n-5) - 2*a(n-6) - 2*a(n-7) + a(n-8), with a(1)=0, a(2)=0, a(3)=1, a(4)=2, a(5)=10, a(6)=18, a(7)=45, a(8)=72.
G.f.: -(x^3*(1+ 4*x^2 + x^4))/((-1+x)^5*(1+x)^3).
a(n) = (5 - 5*(-1)^n - 12*n + 12*(-1)^n*n + 14*n^2 - 6*(-1)^n*n^2 - 8*n^3 + 2*n^4)/64.
a(n) = Sum_{i=1..n-1} floor(i/2)^3. - Ridouane Oudra, Jul 24 2019
E.g.f.: (1/64)*exp(-x)*(-5-6*x-6*x^2+exp(2*x)*(5-4*x+4*x^2+4*x^3+2*x^4)). - Stefano Spezia, Aug 14 2019
a(2*n) = A163102(n-1) and a(2*n+1) = A037270(n). - Ridouane Oudra, Mar 24 2024
Sum_{n>=3} 1/a(n) = Pi^2 - Pi*coth(Pi) - 5. - Amiram Eldar, Jul 04 2025

A006003 a(n) = n*(n^2 + 1)/2.

Original entry on oeis.org

0, 1, 5, 15, 34, 65, 111, 175, 260, 369, 505, 671, 870, 1105, 1379, 1695, 2056, 2465, 2925, 3439, 4010, 4641, 5335, 6095, 6924, 7825, 8801, 9855, 10990, 12209, 13515, 14911, 16400, 17985, 19669, 21455, 23346, 25345, 27455, 29679, 32020, 34481, 37065, 39775
Offset: 0

Views

Author

Keywords

Comments

Write the natural numbers in groups: 1; 2,3; 4,5,6; 7,8,9,10; ... and add the groups. In other words, "sum of the next n natural numbers". - Felice Russo
Number of rhombi in an n X n rhombus, if 'crossformed' rhombi are allowed. - Matti De Craene (Matti.DeCraene(AT)rug.ac.be), May 14 2000
Also the sum of the integers between T(n-1)+1 and T(n), the n-th triangular number (A000217). Sum of n-th row of A000027 regarded as a triangular array.
Unlike the cubes which have a similar definition, it is possible for 2 terms of this sequence to sum to a third. E.g., a(36) + a(37) = 23346 + 25345 = 48691 = a(46). Might be called 2nd-order triangular numbers, thus defining 3rd-order triangular numbers (A027441) as n(n^3+1)/2, etc. - Jon Perry, Jan 14 2004
Also as a(n)=(1/6)*(3*n^3+3*n), n > 0: structured trigonal diamond numbers (vertex structure 4) (cf. A000330 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers). - James A. Record (james.record(AT)gmail.com), Nov 07 2004
The sequence M(n) of magic constants for n X n magic squares (numbered 1 through n^2) from n=3 begins M(n) = 15, 34, 65, 111, 175, 260, ... - Lekraj Beedassy, Apr 16 2005 [comment corrected by Colin Hall, Sep 11 2009]
The sequence Q(n) of magic constants for the n-queens problem in chess begins 0, 0, 0, 0, 34, 65, 111, 175, 260, ... - Paul Muljadi, Aug 23 2005
Alternate terms of A057587. - Jeremy Gardiner, Apr 10 2005
Also partial differences of A063488(n) = (2*n-1)*(n^2-n+2)/2. a(n) = A063488(n) - A063488(n-1) for n>1. - Alexander Adamchuk, Jun 03 2006
In an n X n grid of numbers from 1 to n^2, select -- in any manner -- one number from each row and column. Sum the selected numbers. The sum is independent of the choices and is equal to the n-th term of this sequence. - F.-J. Papp (fjpapp(AT)umich.edu), Jun 06 2006
Nonnegative X values of solutions to the equation (X-Y)^3 - (X+Y) = 0. To find Y values: b(n) = (n^3-n)/2. - Mohamed Bouhamida, May 16 2006
For the equation: m*(X-Y)^k - (X+Y) = 0 with X >= Y, k >= 2 and m is an odd number the X values are given by the sequence defined by a(n) = (m*n^k+n)/2. The Y values are given by the sequence defined by b(n) = (m*n^k-n)/2. - Mohamed Bouhamida, May 16 2006
If X is an n-set and Y a fixed 3-subset of X then a(n-3) is equal to the number of 4-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
(m*(2n)^k+n, m*(2n)^k-n) solves the Diophantine equation: 2m*(X-Y)^k - (X+Y) = 0 with X >= Y, k >= 2 where m is a positive integer. - Mohamed Bouhamida, Oct 02 2007
Also c^(1/2) in a^(1/2) + b^(1/2) = c^(1/2) such that a^2 + b = c. - Cino Hilliard, Feb 09 2008
a(n) = n*A000217(n) - Sum_{i=0..n-1} A001477(i). - Bruno Berselli, Apr 25 2010
a(n) is the number of triples (w,x,y) having all terms in {0,...,n} such that at least one of these inequalities fails: x+y < w, y+w < x, w+x < y. - Clark Kimberling, Jun 14 2012
Sum of n-th row of the triangle in A209297. - Reinhard Zumkeller, Jan 19 2013
The sequence starting with "1" is the third partial sum of (1, 2, 3, 3, 3, ...). - Gary W. Adamson, Sep 11 2015
a(n) is the largest eigenvalue of the matrix returned by the MATLAB command magic(n) for n > 0. - Altug Alkan, Nov 10 2015
a(n) is the number of triples (x,y,z) having all terms in {1,...,n} such that all these triangle inequalities are satisfied: x+y > z, y+z > x, z+x > y. - Heinz Dabrock, Jun 03 2016
Shares its digital root with the stella octangula numbers (A007588). See A267017. - Peter M. Chema, Aug 28 2016
Can be proved to be the number of nonnegative solutions of a system of three linear Diophantine equations for n >= 0 even: 2*a_{11} + a_{12} + a_{13} = n, 2*a_{22} + a_{12} + a_{23} = n and 2*a_{33} + a_{13} + a_{23} = n. The number of solutions is f(n) = (1/16)*(n+2)*(n^2 + 4n + 8) and a(n) = n*(n^2 + 1)/2 is obtained by remapping n -> 2*n-2. - Kamil Bradler, Oct 11 2016
For n > 0, a(n) coincides with the trace of the matrix formed by writing the numbers 1...n^2 back and forth along the antidiagonals (proved, see A078475 for the examples of matrix). - Stefano Spezia, Aug 07 2018
The trace of an n X n square matrix where the elements are entered on the ascending antidiagonals. The determinant is A069480. - Robert G. Wilson v, Aug 07 2018
Bisections are A317297 and A005917. - Omar E. Pol, Sep 01 2018
Number of achiral colorings of the vertices (or faces) of a regular tetrahedron with n available colors. An achiral coloring is identical to its reflection. - Robert A. Russell, Jan 22 2020
a(n) is the n-th centered triangular pyramidal number. - Lechoslaw Ratajczak, Nov 02 2021
a(n) is the number of words of length n defined on 4 letters {b,c,d,e} that contain one or no b's, one c or two d's, and any number of e's. For example, a(3) = 15 since the words are (number of permutations in parentheses): bce (6), bdd (3), cee (3), and dde (3). - Enrique Navarrete, Jun 21 2025

Examples

			G.f. = x + 5*x^2 + 15*x^3 + 34*x^4 + 65*x^5 + 111*x^6 + 175*x^7 + 260*x^8 + ...
For a(2)=5, the five tetrahedra have faces AAAA, AAAB, AABB, ABBB, and BBBB with colors A and B. - _Robert A. Russell_, Jan 31 2020
		

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 15, p. 5, Ellipses, Paris 2008.
  • F.-J. Papp, Colloquium Talk, Department of Mathematics, University of Michigan-Dearborn, March 6, 2005.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000330, A000537, A066886, A057587, A027480, A002817 (partial sums).
Cf. A000578 (cubes).
(1/12)*t*(n^3-n)+n for t = 2, 4, 6, ... gives A004006, A006527, this sequence, A005900, A004068, A000578, A004126, A000447, A004188, A004466, A004467, A007588, A062025, A063521, A063522, A063523.
Antidiagonal sums of array in A000027. Row sums of the triangular view of A000027.
Cf. A063488 (sum of two consecutive terms), A005917 (bisection), A317297 (bisection).
Cf. A105374 / 8.
Tetrahedron colorings: A006008 (oriented), A000332(n+3) (unoriented), A000332 (chiral), A037270 (edges).
Other polyhedron colorings: A337898 (cube faces, octahedron vertices), A337897 (octahedron faces, cube vertices), A337962 (dodecahedron faces, icosahedron vertices), A337960 (icosahedron faces, dodecahedron vertices).
Row 3 of A325001 (simplex vertices and facets) and A337886 (simplex faces and peaks).

Programs

  • GAP
    a_n:=List([0..nmax], n->n*(n^2 + 1)/2); # Stefano Spezia, Aug 12 2018
    
  • Haskell
    a006003 n = n * (n ^ 2 + 1) `div` 2
    a006003_list = scanl (+) 0 a005448_list
    -- Reinhard Zumkeller, Jun 20 2013
    
  • MATLAB
    % Also works with FreeMat.
    for(n=0:nmax); tm=n*(n^2 + 1)/2; fprintf('%d\t%0.f\n', n, tm); end
    % Stefano Spezia, Aug 12 2018
    
  • Magma
    [n*(n^2 + 1)/2 : n in [0..50]]; // Wesley Ivan Hurt, Sep 11 2015
    
  • Magma
    [Binomial(n,3)+Binomial(n-1,3)+Binomial(n-2,3): n in [2..60]]; // Vincenzo Librandi, Sep 12 2015
    
  • Mathematica
    Table[ n(n^2 + 1)/2, {n, 0, 45}]
    LinearRecurrence[{4,-6,4,-1}, {0,1,5,15},50] (* Harvey P. Dale, May 16 2012 *)
    CoefficientList[Series[x (1 + x + x^2)/(x - 1)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Sep 12 2015 *)
    With[{n=50},Total/@TakeList[Range[(n(n^2+1))/2],Range[0,n]]] (* Requires Mathematica version 11 or later *) (* Harvey P. Dale, Nov 28 2017 *)
  • Maxima
    a(n):=n*(n^2 + 1)/2$ makelist(a(n), n, 0, nmax); /* Stefano Spezia, Aug 12 2018 */
    
  • PARI
    {a(n) = n * (n^2 + 1) / 2}; /* Michael Somos, Dec 24 2011 */
    
  • PARI
    concat(0, Vec(x*(1+x+x^2)/(x-1)^4 + O(x^20))) \\ Felix Fröhlich, Oct 11 2016
    
  • Python
    def A006003(n): return n*(n**2+1)>>1 # Chai Wah Wu, Mar 25 2024

Formula

a(n) = binomial(n+2, 3) + binomial(n+1, 3) + binomial(n, 3). [corrected by Michel Marcus, Jan 22 2020]
G.f.: x*(1+x+x^2)/(x-1)^4. - Floor van Lamoen, Feb 11 2002
Partial sums of A005448. - Jonathan Vos Post, Mar 16 2006
Binomial transform of [1, 4, 6, 3, 0, 0, 0, ...] = (1, 5, 15, 34, 65, ...). - Gary W. Adamson, Aug 10 2007
a(n) = -a(-n) for all n in Z. - Michael Somos, Dec 24 2011
a(n) = Sum_{k = 1..n} A(k-1, k-1-n) where A(i, j) = i^2 + i*j + j^2 + i + j + 1. - Michael Somos, Jan 02 2012
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)=5, a(3)=15. - Harvey P. Dale, May 16 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 3. - Ant King, Jun 13 2012
a(n) = A000217(n) + n*A000217(n-1). - Bruno Berselli, Jun 07 2013
a(n) = A057145(n+3,n). - Luciano Ancora, Apr 10 2015
E.g.f.: (1/2)*(2*x + 3*x^2 + x^3)*exp(x). - G. C. Greubel, Dec 18 2015; corrected by Ilya Gutkovskiy, Oct 12 2016
a(n) = T(n) + T(n-1) + T(n-2), where T means the tetrahedral numbers, A000292. - Heinz Dabrock, Jun 03 2016
From Ilya Gutkovskiy, Oct 11 2016: (Start)
Convolution of A001477 and A008486.
Convolution of A000217 and A158799.
Sum_{n>=1} 1/a(n) = H(-i) + H(i) = 1.343731971048019675756781..., where H(k) is the harmonic number, i is the imaginary unit. (End)
a(n) = A000578(n) - A135503(n). - Miquel Cerda, Dec 25 2016
Euler transform of length 3 sequence [5, 0, -1]. - Michael Somos, Dec 25 2016
a(n) = A037270(n)/n for n > 0. - Kritsada Moomuang, Dec 15 2018
a(n) = 3*A000292(n-1) + n. - Bruce J. Nicholson, Nov 23 2019
a(n) = A011863(n) - A011863(n-2). - Bruce J. Nicholson, Dec 22 2019
From Robert A. Russell, Jan 22 2020: (Start)
a(n) = C(n,1) + 3*C(n,2) + 3*C(n,3), where the coefficient of C(n,k) is the number of tetrahedron colorings using exactly k colors.
a(n) = C(n+3,4) - C(n,4).
a(n) = 2*A000332(n+3) - A006008(n) = A006008(n) - 2*A000332(n) = A000332(n+3) - A000332(n).
a(n) = A325001(3,n). (End)
From Amiram Eldar, Aug 21 2023: (Start)
Sum_{n>=1} 1/a(n) = 2 * (A248177 + A001620).
Product_{n>=2} (1 - 1/a(n)) = cosh(sqrt(7)*Pi/2)*cosech(Pi)/4.
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(7)*Pi/2)*cosech(Pi). (End)

Extensions

Better description from Albert Rich (Albert_Rich(AT)msn.com), Mar 1997

A002415 4-dimensional pyramidal numbers: a(n) = n^2*(n^2-1)/12.

Original entry on oeis.org

0, 0, 1, 6, 20, 50, 105, 196, 336, 540, 825, 1210, 1716, 2366, 3185, 4200, 5440, 6936, 8721, 10830, 13300, 16170, 19481, 23276, 27600, 32500, 38025, 44226, 51156, 58870, 67425, 76880, 87296, 98736, 111265, 124950, 139860, 156066, 173641, 192660, 213200, 235340
Offset: 0

Views

Author

Keywords

Comments

Also number of ways to legally insert two pairs of parentheses into a string of m := n-1 letters. (There are initially 2C(m+4,4) (A034827) ways to insert the parentheses, but we must subtract 2(m+1) for illegal clumps of 4 parentheses, 2m(m+1) for clumps of 3 parentheses, C(m+1,2) for 2 clumps of 2 parentheses and (m-1)C(m+1,2) for 1 clump of 2 parentheses, giving m(m+1)^2(m+2)/12 = n^2*(n^2-1)/12.) See also A000217.
E.g., for n=2 there are 6 ways: ((a))b, ((a)b), ((ab)), (a)(b), (a(b)), a((b)).
Let M_n denote the n X n matrix M_n(i,j)=(i+j); then the characteristic polynomial of M_n is x^(n-2) * (x^2-A002378(n)*x - a(n)). - Benoit Cloitre, Nov 09 2002
Let M_n denote the n X n matrix M_n(i,j)=(i-j); then the characteristic polynomial of M_n is x^n + a(n)x^(n-2). - Michael Somos, Nov 14 2002 [See A114327 for the infinite matrix M in triangular form. - Wolfdieter Lang, Feb 05 2018]
Number of permutations of [n] which avoid the pattern 132 and have exactly 2 descents. - Mike Zabrocki, Aug 26 2004
Number of tilings of a <2,n,2> hexagon.
a(n) is the number of squares of side length at least 1 having vertices at the points of an n X n unit grid of points (the vertices of an n-1 X n-1 chessboard). [For a proof, see Comments in A051602. - N. J. A. Sloane, Sep 29 2021] For example, on the 3 X 3 grid (the vertices of a 2 X 2 chessboard) there are four 1 X 1 squares, one (skew) sqrt(2) X sqrt(2) square, and one 3 X 3 square, so a(3)=6. On the 4 X 4 grid (the vertices of a 3 X 3 chessboard) there are 9 1 X 1 squares, 4 2 X 2 squares, 1 3 X 3 square, 4 sqrt(2) X sqrt(2) squares, and 2 sqrt(5) X sqrt(5) squares, so a(4) = 20. See also A024206, A108279. [Comment revised by N. J. A. Sloane, Feb 11 2015]
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Jun 12 2005
Number of distinct components of the Riemann curvature tensor. - Gene Ward Smith, Apr 24 2006
a(n) is the number of 4 X 4 matrices (symmetrical about each diagonal) M = [a,b,c,d;b,e,f,c;c,f,e,b;d,c,b,a] with a+b+c+d=b+e+f+c=n+2; (a,b,c,d,e,f natural numbers). - Philippe Deléham, Apr 11 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 5-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
a(n) is the number of Dyck (n+1)-paths with exactly n-1 peaks. - David Callan, Sep 20 2007
Starting (1,6,20,50,...) = third partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = Sum_{i=0..n} C(n+3,i+3)*b(i), where b(i)=[1,2,0,0,0,...]. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
4-dimensional square numbers. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
Equals row sums of triangle A177877; a(n), n > 1 = (n-1) terms in (1,2,3,...) dot (...,3,2,1) with additive carryovers. Example: a(4) = 20 = (1,2,3) dot (3,2,1) with carryovers = (1*3) + (2*2 + 3) + (3*1 + 7) = (3 + 7 + 10).
Convolution of the triangular numbers A000217 with the odd numbers A004273.
a(n+2) is the number of 4-tuples (w,x,y,z) with all terms in {0,...,n} and w-x=max{w,x,y,z}-min{w,x,y,z}. - Clark Kimberling, May 28 2012
The second level of finite differences is a(n+2) - 2*a(n+1) + a(n) = (n+1)^2, the squares. - J. M. Bergot, May 29 2012
Because the differences of this sequence give A000330, this is also the number of squares in an n+1 X n+1 grid whose sides are not parallel to the axes.
a(n+2) gives the number of 2*2 arrays that can be populated with 0..n such that rows and columns are nondecreasing. - Jon Perry, Mar 30 2013
For n consecutive numbers 1,2,3,...,n, the sum of all ways of adding the k-tuples of consecutive numbers for n=a(n+1). As an example, let n=4: (1)+(2)+(3)+(4)=10; (1+2)+(2+3)+(3+4)=15; (1+2+3)+(2+3+4)=15; (1+2+3+4)=10 and the sum of these is 50=a(4+1)=a(5). - J. M. Bergot, Apr 19 2013
If P(n,k) = n*(n+1)*(k*n-k+3)/6 is the n-th (k+2)-gonal pyramidal number, then a(n) = P(n,k)*P(n-1,k-1) - P(n-1,k)*P(n,k-1). - Bruno Berselli, Feb 18 2014
For n > 1, a(n) = 1/6 of the area of the trapezoid created by the points (n,n+1), (n+1,n), (1,n^2+n), (n^2+n,1). - J. M. Bergot, May 14 2014
For n > 3, a(n) is twice the area of a triangle with vertices at points (C(n,4),C(n+1,4)), (C(n+1,4),C(n+2,4)), and (C(n+2,4),C(n+3,4)). - J. M. Bergot, Jun 03 2014
a(n) is the dimension of the space of metric curvature tensors (those having the symmetries of the Riemann curvature tensor of a metric) on an n-dimensional real vector space. - Daniel J. F. Fox, Dec 15 2018
Coefficients in the terminating series identity 1 - 6*n/(n + 5) + 20*n*(n - 1)/((n + 5)*(n + 6)) - 50*n*(n - 1)*(n - 2)/((n + 5)*(n + 6)*(n + 7)) + ... = 0 for n = 1,2,3,.... Cf. A000330 and A005585. - Peter Bala, Feb 18 2019

Examples

			a(7) = 6*21 - (6*0 + 4*1 + 2*3 + 0*6 - 2*10 - 4*15) = 196. - _Bruno Berselli_, Jun 22 2013
G.f. = x^2 + 6*x^3 + 20*x^4 + 50*x^5 + 105*x^6 + 196*x^7 + 336*x^8 + ...
		

References

  • O. D. Anderson, Find the next sequence, J. Rec. Math., 8 (No. 4, 1975-1976), 241.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p.165).
  • R. Euler and J. Sadek, "The Number of Squares on a Geoboard", Journal of Recreational Mathematics, 251-5 30(4) 1999-2000 Baywood Pub. NY
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 238.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = ((-1)^n)*A053120(2*n, 4)/8 (one-eighth of fifth unsigned column of Chebyshev T-triangle, zeros omitted). Cf. A001296.
Second row of array A103905.
Third column of Narayana numbers A001263.
Partial sums of A000330.
The expression binomial(m+n-1,n)^2-binomial(m+n,n+1)*binomial(m+n-2,n-1) for the values m = 2 through 14 produces sequences A000012, A000217, A002415, A006542, A006857, A108679, A134288, A134289, A134290, A134291, A140925, A140935, A169937.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers (A000027) with the k-gonal numbers.

Programs

  • GAP
    List([0..45],n->Binomial(n^2,2)/6); # Muniru A Asiru, Dec 15 2018
  • Magma
    [n^2*(n^2-1)/12: n in [0..50]]; // Wesley Ivan Hurt, May 14 2014
    
  • Maple
    A002415 := proc(n) binomial(n^2,2)/6 ; end proc: # Zerinvary Lajos, Jan 07 2008
  • Mathematica
    Table[(n^4 - n^2)/12, {n, 0, 40}] (* Zerinvary Lajos, Mar 21 2007 *)
    LinearRecurrence[{5,-10,10,-5,1},{0,0,1,6,20},40] (* Harvey P. Dale, Nov 29 2011 *)
  • PARI
    a(n) = n^2 * (n^2 - 1) / 12;
    
  • PARI
    x='x+O('x^200); concat([0, 0], Vec(x^2*(1+x)/(1-x)^5)) \\ Altug Alkan, Mar 23 2016
    

Formula

G.f.: x^2*(1+x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = Sum_{i=0..n} (n-i)*i^2 = a(n-1) + A000330(n-1) = A000217(n)*A000292(n-2)/n = A000217(n)*A000217(n-1)/3 = A006011(n-1)/3, convolution of the natural numbers with the squares. - Henry Bottomley, Oct 19 2000
a(n)+1 = A079034(n). - Mario Catalani (mario.catalani(AT)unito.it), Feb 12 2003
a(n) = 2*C(n+2, 4) - C(n+1, 3). - Paul Barry, Mar 04 2003
a(n) = C(n+2, 4) + C(n+1, 4). - Paul Barry, Mar 13 2003
a(n) = Sum_{k=1..n} A000330(n-1). - Benoit Cloitre, Jun 15 2003
a(n) = n*C(n+1,3)/2 = C(n+1,3)*C(n+1,2)/(n+1). - Mitch Harris, Jul 06 2006
a(n) = A006011(n)/3 = A008911(n)/2 = A047928(n-1)/12 = A083374(n)/6. - Zerinvary Lajos, May 09 2007
a(n) = (1/2)*Sum_{1 <= x_1, x_2 <= n} (det V(x_1,x_2))^2 = (1/2)*Sum_{1 <= i,j <= n} (i-j)^2, where V(x_1,x_2) is the Vandermonde matrix of order 2. - Peter Bala, Sep 21 2007
a(n) = C(n+1,3) + 2*C(n+1,4). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
a(n) = (1/48)*sinh(2*arccosh(n))^2. - Artur Jasinski, Feb 10 2010
a(n) = n*A000292(n-1)/2. - Tom Copeland, Sep 13 2011
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5), n > 4. - Harvey P. Dale, Nov 29 2011
a(n) = (n-1)*A000217(n-1) - Sum_{i=0..n-2} (n-1-2*i)*A000217(i) for n > 1. - Bruno Berselli, Jun 22 2013
a(n) = C(n,2)*C(n+1,3) - C(n,3)*C(n+1,2). - J. M. Bergot, Sep 17 2013
a(n) = Sum_{k=1..n} ( (2k-n)* k(k+1)/2 ). - Wesley Ivan Hurt, Sep 26 2013
a(n) = floor(n^2/3) + 3*Sum_{k=1..n} k^2*floor((n-k+1)/3). - Mircea Merca, Feb 06 2014
Euler transform of length 2 sequence [6, -1]. - Michael Somos, May 28 2014
G.f. x^2*2F1(3,4;2;x). - R. J. Mathar, Aug 09 2015
Sum_{n>=2} 1/a(n) = 21 - 2*Pi^2 = 1.260791197821282762331... . - Vaclav Kotesovec, Apr 27 2016
a(n) = A080852(2,n-2). - R. J. Mathar, Jul 28 2016
a(n) = A046092(n) * A046092(n-1)/48 = A000217(n) * A000217(n-1)/3. - Bruce J. Nicholson, Jun 06 2017
E.g.f.: (1/12)*exp(x)*x^2*(6 + 6*x + x^2). - Stefano Spezia, Dec 07 2018
Sum_{n>=2} (-1)^n/a(n) = Pi^2 - 9 (See A002388). - Amiram Eldar, Jun 28 2020

Extensions

Typo in link fixed by Matthew Vandermast, Nov 22 2010
Redundant comment deleted and more detail on relationship with A000330 added by Joshua Zucker, Jan 01 2013

A000982 a(n) = ceiling(n^2/2).

Original entry on oeis.org

0, 1, 2, 5, 8, 13, 18, 25, 32, 41, 50, 61, 72, 85, 98, 113, 128, 145, 162, 181, 200, 221, 242, 265, 288, 313, 338, 365, 392, 421, 450, 481, 512, 545, 578, 613, 648, 685, 722, 761, 800, 841, 882, 925, 968, 1013, 1058, 1105, 1152, 1201, 1250, 1301, 1352, 1405
Offset: 0

Views

Author

Keywords

Comments

a(n) = number of pairs (i,j) in [1..n] X [1..n] with integral arithmetic mean. Cf. A132188, A362931. - N. J. A. Sloane, Aug 28 2023
Also, floor( (n^2+1)/2 ). - N. J. A. Sloane, Feb 08 2019
Floor(arithmetic mean of next n numbers). - Amarnath Murthy, Mar 11 2003
Pairwise sums of repeated squares (A008794).
Also, number of topologies on n+1 unlabeled elements with exactly 4 elements in the topology. a(3) gives 4 elements a,b,c,d; the valid topologies are (0,a,ab,abcd), (0,a,abc,abcd), (0,ab,abc,abcd), (0,a,bcd,abcd) and (0,ab,cd,abcd), with a count of 5. - Jon Perry, Mar 05 2004
Partition n into two parts, say, r and s, so that r^2 + s^2 is minimal, then a(n) = r^2 + s^2. Geometrical significance: folding a rod with length n units at right angles in such a way that the end points are at the least distance, which is given by a(n)^(1/2) as the hypotenuse of a right triangle with the sum of the base and height = n units. - Amarnath Murthy, Apr 18 2004
Convolution of A002061(n)-0^n and (-1)^n. Convolution of n (A001477) with {1,0,2,0,2,0,2,...}. Partial sums of repeated odd numbers {0,1,1,3,3,5,5,...}. - Paul Barry, Jul 22 2004
The ratio of the sum of terms over the total number of terms in an n X n spiral. The sum of terms of an n X n spiral is A037270, or Sum_{k=0..n^2} k = (n^4 + n^2)/2 and the total number of terms is n^2. - William A. Tedeschi, Feb 27 2008
Starting with offset 1 = row sums of triangle A158946. - Gary W. Adamson, Mar 31 2009
Partial sums of A109613. - Reinhard Zumkeller, Dec 05 2009
Also the number of compositions of even natural numbers into 2 parts < n. For example a(3)=5 are the compositions (0,0), (0,2), (2,0), (1,1), (2,2) of even natural numbers into 2 parts < 3. a(4)=8 are the compositions (0,0), (0,2), (2,0), (1,1), (2,2), (1,3), (3,1), (3,3) of even natural numbers into 2 parts < 4. - Adi Dani, Jun 05 2011
A001105 and A001844 interleaved. - Omar E. Pol, Sep 18 2011
Number of (w,x,y) having all terms in {0,...,n} and w=average(x,y). - Clark Kimberling, May 15 2012
For n > 0, minimum number of lines necessary to get through all unit cubes of an n X n X n cube (see Kantor link). - Michel Marcus, Apr 13 2013
Sum_{n > 0} 1/a(n) = Sum_{n > 0} 1/(2*n^2) + Sum_{n >= 0} 1/(2*n + 2*n^2 + 1) = (zeta(2) + (Pi* tanh(Pi/2)))/2 = 2.26312655.... - Enrique Pérez Herrero, Jun 17 2013
For n > 1, a(n) is the edge cover number of the n X n king graph. - Eric W. Weisstein, Jun 20 2017
Also the number of vertices in the n X n black bishop graph. - Eric W. Weisstein, Jun 26 2017
The same sequence arises in the triangular array of the integers >= 1, according to a simple "zig-zag" rule for selection of terms. a(n-1) lies in the (n-1)-th row of the array, and the second row of that sub-array (with apex a(n-1)) contains just two numbers, one odd, one even. The one with opposite parity to a(n-1) is a(n). - David James Sycamore, Jul 29 2018
Size of minimal ternary 1-covering code with code length n, i.e., K_n(3,1). See Kalbfleisch and Stanton. - Patrick Wienhöft, Jan 29 2019
For n > 1, a(n-1) is the maximum number of inversions in a permutation consisting of a single n-cycle on n symbols. - M. Ryan Julian Jr., Sep 10 2019
Also the number of classes of convex inscribed polyominoes in a (2,n) rectangular grid; two polyominoes are in the same class if one of them can be obtained by a reflection or 180-degree rotation of the other. - Jean-Luc Manguin, Jan 29 2020
a(n) is the number of pairs (p,q) such that 1 <= p, p+1 < q <= n+2 and q <> 2*p. - César Eliud Lozada, Oct 25 2020
a(n) is the maximum number of copies of a 12 permutation pattern in an alternating (or zig-zag) permutation of length n+1. The maximum number of copies of 123 in an alternating permutation is motivated in the Notices reference, and the argument here is analogous. - Lara Pudwell, Dec 01 2020
It appears that a(n) is the largest number of nodes of an induced path in the n X n king graph. An induced path going in a simple spiraling pattern, starting in a corner, has a(n) nodes. For even n this is optimal, because an induced path can have at most two nodes in any 2 X 2 subsquare. For odd n, I cannot see how to prove that (n^2+1)/2 is best possible. See also A357501. - Pontus von Brömssen, Oct 02 2022 [Proved by Beluhov (2023). - Pontus von Brömssen, Jan 30 2023]
a(n) = n + 2*(n-2) + 2*(n-4) + 2*(n-6) + ... number of black squares on an n X n chessboard. - R. J. Mathar, Dec 03 2022

Examples

			G.f. = x + 2*x^2 + 5*x^3 + 8*x^4 + 13*x^5 + 18*x^6 + 25*x^7 + 32*x^8 + ...
Centrosymmetric 3 X 3 matrix: [[a,b,c],[d,e,d],[c,b,a]], a(3) = 3*(3-1)/2 + (3-1)/2 + 1 = (3^2+1)/2 = 5 from a,b,c,d,e. 4 X 4 case: [[a,b,c,d],[e,f,g,h],[h,g,f,e],[d,c,b,a]], a(4) = 4*4/2 = 8. - _Wolfdieter Lang_, Oct 12 2015
a(3) = 5. The alternating permutation of length 3 + 1 = 4 with the maximum number of copies of 123 is 1324. The five copies are 12, 13, 14, 23, and 24. - _Lara Pudwell_, Dec 01 2020
		

References

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

Crossrefs

Programs

Formula

a(2*n) = 2*n^2, a(2*n+1) = 2*n^2 + 2*n + 1.
G.f.: -x*(1+x^2) / ( (1+x)*(x-1)^3 ). - Simon Plouffe in his 1992 dissertation
From Benoit Cloitre, Nov 06 2002: (Start)
a(n) = (2*n^2 + 1 - (-1)^n) / 4.
a(0)=0, a(1)=1; for n>1, a(n+1) = n + 1 + max(2*floor(a(n)/2), 3*floor(a(n)/3)). (End)
G.f.: (x + x^2 + x^3 + x^4)/((1 - x)*(1 - x^2)^2), not reduced. - Len Smiley
a(n) = a(n-2) + 2n - 2. - Paul Barry, Jul 17 2004
From Paul Barry, Jul 22 2004: (Start)
G.f.: x*(1+x^2)/((1-x^2)*(1-x)^2) = x*(1+x^2)/((1+x)*(1-x)^3);
a(n) = Sum_{k=0..n} (k^2 - k + 1 - 0^k)*(-1)^(n-k);
a(n) = Sum_{k=0..n} (1 + (-1)^(n-k) - 0^(n-k))*k. (End)
From Reinhard Zumkeller, Feb 27 2006: (Start)
a(0) = 0, a(n+1) = a(n) + 2*floor(n/2) + 1.
a(n) = A116940(n) - A005843(n). (End)
Starting with offset 1, = row sums of triangle A134444. Also, with offset 1, = binomial transform of [1, 1, 2, -2, 4, -8, 16, -32, ...]. - Gary W. Adamson, Oct 25 2007
a(n) = floor((n^2+1)/2). - William A. Tedeschi, Feb 27 2008
a(n) = A004526(n+1) + A000217(n-1). - Yosu Yurramendi, Sep 12 2008, corrected by Klaus Purath, Jun 15 2021
From Jaume Oliver Lafont, Dec 05 2008: (Start)
a(n) = a(n-1) + a(n-2) - a(n-3) + 2.
a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4). (End)
a(n) = A004526(n)^2 + A110654(n)^2. - Philippe Deléham, Mar 12 2009
a(n) = n^2 - floor(n^2/2). - Wesley Ivan Hurt, Jun 14 2013
Euler transform is length 4 sequence [2, 2, 0, -1].
a(n) = a(-n) for all n in Z. - Michael Somos, May 05 2015
a(n) is also the number of independent entries in a centrosymmetric n X n matrix: M(i, j) = M(n-i+1, n-j+1). - Wolfdieter Lang, Oct 12 2015
For n > 1, a(n+1)/a(n) = 3 - A081352(n-2)/a(n). - Miko Labalan, Mar 26 2016
E.g.f.: (1/2)*(x*(1 + x)*cosh(x) + (1 + x + x^2)*sinh(x)). - Stefano Spezia, Feb 03 2020
a(n) = binomial(n+1,2) - floor(n/2). - César Eliud Lozada, Oct 25 2020
From Klaus Purath, Jun 15 2021: (Start)
a(n-1) + a(n) = A002061(n).
a(n) = (a(n-1)^2 + 1) / a(n-2), n >= 3 odd.
a(n) = (a(n-1)^2 - (n-1)^2) / a(n-2), n >= 4 even. (End)

A005898 Centered cube numbers: n^3 + (n+1)^3.

Original entry on oeis.org

1, 9, 35, 91, 189, 341, 559, 855, 1241, 1729, 2331, 3059, 3925, 4941, 6119, 7471, 9009, 10745, 12691, 14859, 17261, 19909, 22815, 25991, 29449, 33201, 37259, 41635, 46341, 51389, 56791, 62559, 68705, 75241, 82179, 89531, 97309, 105525, 114191, 123319, 132921
Offset: 0

Views

Author

Keywords

Comments

Write the natural numbers in groups: 1; 2,3,4; 5,6,7,8,9; 10,11,12,13,14,15,16; ..... and add the groups, i.e., a(n) = Sum_{j=n^2-2(n-1)..n^2} j. - Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Sep 05 2001
The numbers 1, 9, 35, 91, etc. are divisible by 1, 3, 5, 7, etc. Therefore there are no prime numbers in this list. 9 is divisible by 3 and every third number after 9 is also divisible by 3. 35 is divisible by 5 and 7 and every fifth number after 35 is also divisible by 5 and every seventh number after 35 is also divisible by 7. This pattern continues indefinitely. - Howard Berman (howard_berman(AT)hotmail.com), Nov 07 2008
n^3 + (n+1)^3 = (2n+1)*(n^2+n+1), hence all terms are composite. - Zak Seidov, Feb 08 2011
This is the order of an n-ball centered at a node in the Kronecker product (or direct product) of three cycles, each of whose lengths is at least 2n+2. - Pranava K. Jha, Oct 10 2011
Positive y values of 4*x^3 - 3*x^2 = y^2. - Bruno Berselli, Apr 28 2018

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 52.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

(1/12)*t*(2*n^3 - 3*n^2 + n) + 2*n - 1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A005917, A063493, A063494, A063495, A063496.
Partial sums of A005897.
The 28 uniform 3D tilings: cab: A299266, A299267; crs: A299268, A299269; fcu: A005901, A005902; fee: A299259, A299265; flu-e: A299272, A299273; fst: A299258, A299264; hal: A299274, A299275; hcp: A007899, A007202; hex: A005897, A005898; kag: A299256, A299262; lta: A008137, A299276; pcu: A005899, A001845; pcu-i: A299277, A299278; reo: A299279, A299280; reo-e: A299281, A299282; rho: A008137, A299276; sod: A005893, A005894; sve: A299255, A299261; svh: A299283, A299284; svj: A299254, A299260; svk: A010001, A063489; tca: A299285, A299286; tcd: A299287, A299288; tfs: A005899, A001845; tsi: A299289, A299290; ttw: A299257, A299263; ubt: A299291, A299292; bnn: A007899, A007202. See the Proserpio link in A299266 for overview.

Programs

Formula

a(n) = Sum_{i=0..n} A005897(i), partial sums. - Jonathan Vos Post, Feb 06 2011
G.f.: (x^2+4*x+1)*(1+x)/(1-x)^3. - Simon Plouffe (see MAPLE section) and Colin Barker, Jan 02 2012; edited by N. J. A. Sloane, Feb 07 2018
a(n) = A037270(n+1) - A037270(n). - Ivan N. Ianakiev, May 13 2012
a(n) = A000217(n+1)^2 - A000217(n-1)^2. - Bob Selcoe, Mar 25 2016
a(n) = A005408(n) * A002061(n+1). - Miquel Cerda, Oct 05 2016
From Ilya Gutkovskiy, Oct 06 2016: (Start)
E.g.f.: (1 + 8*x + 9*x^2 + 2*x^3)*exp(x).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). (End)
a(n) = (A081435(n))^2 - (A081435(n) - 1)^2. - Sergey Pavlov, Mar 01 2017

A058331 a(n) = 2*n^2 + 1.

Original entry on oeis.org

1, 3, 9, 19, 33, 51, 73, 99, 129, 163, 201, 243, 289, 339, 393, 451, 513, 579, 649, 723, 801, 883, 969, 1059, 1153, 1251, 1353, 1459, 1569, 1683, 1801, 1923, 2049, 2179, 2313, 2451, 2593, 2739, 2889, 3043, 3201, 3363, 3529, 3699, 3873, 4051
Offset: 0

Views

Author

Erich Friedman, Dec 12 2000

Keywords

Comments

Maximal number of regions in the plane that can be formed with n hyperbolas.
Also the number of different 2 X 2 determinants with integer entries from 0 to n.
Number of lattice points in an n-dimensional ball of radius sqrt(2). - David W. Wilson, May 03 2001
Equals A112295(unsigned) * [1, 2, 3, ...]. - Gary W. Adamson, Oct 07 2007
Binomial transform of A166926. - Gary W. Adamson, May 03 2008
a(n) = longest side a of all integer-sided triangles with sides a <= b <= c and inradius n >= 1. Triangle has sides (2n^2 + 1, 2n^2 + 2, 4n^2 + 1).
{a(k): 0 <= k < 3} = divisors of 9. - Reinhard Zumkeller, Jun 17 2009
Number of ways to partition a 3*n X 2 grid into 3 connected equal-area regions. - R. H. Hardin, Oct 31 2009
Let A be the Hessenberg matrix of order n defined by: A[1, j] = 1, A[i, i] := 2, (i > 1), A[i, i - 1] = -1, and A[i, j] = 0 otherwise. Then, for n >= 3, a(n - 1) = coeff(charpoly(A, x), x^(n - 2)). - Milan Janjic, Jan 26 2010
Except for the first term of [A002522] and [A058331] if X = [A058331], Y = [A087113], A = [A002522], we have, for all other terms, Pell's equation: [A058331]^2 - [A002522]*[A087113]^2 = 1; (X^2 - A*Y^2 = 1); e.g., 3^2 -2*2^2 = 1; 9^2 - 5*4^2 = 1; 129^2 - 65*16^2 = 1, and so on. - Vincenzo Librandi, Aug 07 2010
Niven (1961) gives this formula as an example of a formula that does not contain all odd integers, in contrast to 2n + 1 and 2n - 1. - Alonso del Arte, Dec 05 2012
Numbers m such that 2*m-2 is a square. - Vincenzo Librandi, Apr 10 2015
Number of n-tuples from the set {1,0,-1} where at most two elements are nonzero. - Michael Somos, Oct 19 2022
a(n) gives the x-value of the integral solution (x,y) of the Pellian equation x^2 - (n^2 + 1)*y^2 = 1. The y-value is given by 2*n (see Tattersall). - Stefano Spezia, Jul 23 2025

Examples

			a(1) = 3 since (0 0 / 0 0), (1 0 / 0 1) and (0 1 / 1 0) have different determinants.
G.f. = 1 + 3*x + 9*x^2 + 19*x^3 + 33*x^4 + 51*x^5 + 73*x^6 + ... - _Michael Somos_, Oct 19 2022
		

References

  • Ivan Niven, Numbers: Rational and Irrational, New York: Random House for Yale University (1961): 17.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 256.

Crossrefs

Cf. A000124.
Second row of array A099597.
See A120062 for sequences related to integer-sided triangles with integer inradius n.
Cf. A112295.
Column 2 of array A188645.
Cf. A001105 and A247375. - Bruno Berselli, Sep 16 2014

Programs

  • Haskell
    a058331 = (+ 1) . a001105  -- Reinhard Zumkeller, Dec 13 2014
    
  • Magma
    [2*n^2 + 1 : n in [0..100]]; // Wesley Ivan Hurt, Feb 02 2017
  • Mathematica
    b[g_] := Length[Union[Map[Det, Flatten[ Table[{{i, j}, {k, l}}, {i, 0, g}, {j, 0, g}, {k, 0, g}, {l, 0, g}], 3]]]] Table[b[g], {g, 0, 20}]
    2*Range[0, 49]^2 + 1 (* Alonso del Arte, Dec 05 2012 *)
  • PARI
    a(n)=2*n^2+1 \\ Charles R Greathouse IV, Jun 16 2011
    

Formula

G.f.: (1 + 3x^2)/(1 - x)^3. - Paul Barry, Apr 06 2003
a(n) = M^n * [1 1 1], leftmost term, where M = the 3 X 3 matrix [1 1 1 / 0 1 4 / 0 0 1]. a(0) = 1, a(1) = 3; a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). E.g., a(4) = 33 since M^4 *[1 1 1] = [33 17 1]. - Gary W. Adamson, Nov 11 2004
a(n) = cosh(2*arccosh(n)). - Artur Jasinski, Feb 10 2010
a(n) = 4*n + a(n-1) - 2 for n > 0, a(0) = 1. - Vincenzo Librandi, Aug 07 2010
a(n) = (((n-1)^2 + n^2))/2 + (n^2 + (n+1)^2)/2. - J. M. Bergot, May 31 2012
a(n) = A251599(3*n) for n > 0. - Reinhard Zumkeller, Dec 13 2014
a(n) = sqrt(8*(A000217(n-1)^2 + A000217(n)^2) + 1). - J. M. Bergot, Sep 03 2015
E.g.f.: (2*x^2 + 2*x + 1)*exp(x). - G. C. Greubel, Jul 14 2017
a(n) = A002378(n) + A002061(n). - Bruce J. Nicholson, Aug 06 2017
From Amiram Eldar, Jul 15 2020: (Start)
Sum_{n>=0} 1/a(n) = (1 + (Pi/sqrt(2))*coth(Pi/sqrt(2)))/2.
Sum_{n>=0} (-1)^n/a(n) = (1 + (Pi/sqrt(2))*csch(Pi/sqrt(2)))/2. (End)
From Amiram Eldar, Feb 05 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = sqrt(2)*csch(Pi/sqrt(2))*sinh(Pi).
Product_{n>=1} (1 - 1/a(n)) = (Pi/sqrt(2))*csch(Pi/sqrt(2)). (End)
From Leo Tavares, May 23 2022: (Start)
a(n) = A000384(n+1) - 3*n.
a(n) = 3*A000217(n) + A000217(n-2). (End)
a(n) = a(-n) for all n in Z and A037235(n) = Sum_{k=0..n-1} a(k). - Michael Somos, Oct 19 2022

Extensions

Revised description from Noam Katz (noamkj(AT)hotmail.com), Jan 28 2001

A027441 a(n) = (n^4 + n)/2 (Row sums of an n X n X n magic cube, when it exists).

Original entry on oeis.org

0, 1, 9, 42, 130, 315, 651, 1204, 2052, 3285, 5005, 7326, 10374, 14287, 19215, 25320, 32776, 41769, 52497, 65170, 80010, 97251, 117139, 139932, 165900, 195325, 228501, 265734, 307342, 353655, 405015, 461776, 524304, 592977, 668185, 750330, 839826, 937099
Offset: 0

Views

Author

Keywords

Comments

Starting with offset 1 = binomial transform of (1, 8, 25, 30, 12, 0, 0, 0, ...). - Gary W. Adamson, May 20 2009

Crossrefs

Subsequence of A057590.

Programs

Formula

O.g.f.: x*(1+4*x+7*x^2)/(1-x)^5. - R. J. Mathar, Feb 13 2008
a(n) = Sum_{k=n..n^2} k; for n>0: a(n) = A037270(n) - A000217(n-1). - Reinhard Zumkeller, Jul 06 2010
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Wesley Ivan Hurt, Aug 13 2014
a(n) = A071232(n) / n^2, for n > 0. - Wesley Ivan Hurt, Aug 13 2014
a(n) = A002061(n)*A000217(n). - Anton Zakharov, Dec 16 2016
a(n) = (n+1)*(a(n-1)/(n-1) + n*(n-1)), a(0)=0, a(1)=1. - Vladimir Kruchinin, Oct 10 2018

Extensions

More terms from Wesley Ivan Hurt, Aug 13 2014

A083374 a(n) = n^2 * (n^2 - 1)/2.

Original entry on oeis.org

0, 6, 36, 120, 300, 630, 1176, 2016, 3240, 4950, 7260, 10296, 14196, 19110, 25200, 32640, 41616, 52326, 64980, 79800, 97020, 116886, 139656, 165600, 195000, 228150, 265356, 306936, 353220, 404550, 461280, 523776, 592416, 667590, 749700, 839160, 936396
Offset: 1

Views

Author

Alan Sutcliffe (alansut(AT)ntlworld.com), Jun 05 2003

Keywords

Comments

Triangular numbers t_n as n runs through the squares.
Partial sums of A055112: If one generated Pythagorean primitive triangles from n, n+1, then the collective areas of n of them would be equal to the numbers in this sequence. The sum of the first three triangles is 6+30+84 = 120 which is the third nonzero term of the sequence. - J. M. Bergot, Jul 14 2011
Second leg of Pythagorean triangles with smallest side a cube: A000578(n)^2 + a(n)^2 = A037270(n)^2. - Martin Renner, Nov 12 2011
a(n) is the number of segments on an n X n grid or geoboard. - Martin Renner, Apr 17 2014
Consider the partitions of 2n into two parts (p,q). Then a(n) is the total volume of the family of rectangular prisms with dimensions p, q and |q-p|. - Wesley Ivan Hurt, Apr 15 2018

References

  • Albert H. Beiler, Recreations in the theory of numbers, New York: Dover, (2nd ed.) 1966, p. 106, table 55.

Crossrefs

Programs

Formula

a(n) = (n + 1) * A006002(n).
a(n) = A047928(n)/2 = A002415(n+1)*6 = A006011(n+1)*2 = A008911(n+1)*3. - Zerinvary Lajos, May 09 2007
a(n) = binomial(n^2,2), n>=1. - Zerinvary Lajos, Jan 07 2008
a(n) = 5*a(n-1)-10*a(n-2)+10*a(n-3)-5*a(n-4)+a(n-5) for n>5. - R. J. Mathar, Apr 10 2009
G.f.: -6*x^2*(1+x)/(x-1)^5. - R. J. Mathar, Apr 10 2009
Sum_{n>1} 1/a(n) = (21 - 2*Pi^2)/6. - Enrique Pérez Herrero, Apr 01 2013
a(n) = Sum_{k=0..n-1} k*A000217(2*k+1). - Bruno Berselli, Sep 04 2013
a(n) = 2*A000217(n-1)*A000217(n). - Gionata Neri, May 04 2015
a(n) = Sum_{i=1..n^2-1} i. - Wesley Ivan Hurt, Nov 24 2015
E.g.f.: exp(x)*x^2*(6 + 6*x + x^2)/2. - Stefano Spezia, Jun 06 2021
Sum_{n>=2} (-1)^n/a(n) = Pi^2/6 - 3/2. - Amiram Eldar, Nov 02 2021

Extensions

Corrected and extended by T. D. Noe, Oct 25 2006

A277504 Array read by descending antidiagonals: T(n,k) is the number of unoriented strings with n beads of k or fewer colors.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 3, 1, 0, 1, 4, 6, 6, 1, 0, 1, 5, 10, 18, 10, 1, 0, 1, 6, 15, 40, 45, 20, 1, 0, 1, 7, 21, 75, 136, 135, 36, 1, 0, 1, 8, 28, 126, 325, 544, 378, 72, 1, 0, 1, 9, 36, 196, 666, 1625, 2080, 1134, 136, 1, 0, 1, 10, 45, 288, 1225, 3996, 7875, 8320, 3321, 272, 1, 0
Offset: 0

Author

Jean-François Alcover, Oct 18 2016

Keywords

Comments

From Petros Hadjicostas, Jul 07 2018: (Start)
Column k of this array is the "BIK" (reversible, indistinct, unlabeled) transform of k,0,0,0,....
Consider the input sequence (c_k(n): n >= 1) with g.f. C_k(x) = Sum_{n>=1} c_k(n)*x^n. Let a_k(n) = BIK(c_k(n): n >= 1) be the output sequence under Bower's BIK transform. It can proved that the g.f. of BIK(c_k(n): n >= 1) is A_k(x) = (1/2)*(C_k(x)/(1-C_k(x)) + (C_k(x^2) + C_k(x))/(1-C_k(x^2))). (See the comments for sequence A001224.)
For column k of this two-dimensional array, the input sequence is defined by c_k(1) = k and c_k(n) = 0 for n >= 1. Thus, C_k(x) = k*x, and hence the g.f. of column k is (1/2)*(C_k(x)/(1-C_k(x)) + (C_k(x^2) + C_k(x))/(1-C_k(x^2))) = (1/2)*(k*x/(1-k*x) + (k*x^2 + k*x)/(1-k*x^2)) = (2 + (1-k)*x - 2*k*x^2)*k*x/(2*(1-k*x^2)*(1-k*x)).
Using the first form the g.f. above and the expansion 1/(1-y) = 1 + y + y^2 + ..., we can easily prove J.-F. Alcover's formula T(n,k) = (k^n + k^((n + mod(n,2))/2))/2.
(End)

Examples

			Array begins with T(0,0):
1 1   1     1      1       1        1         1         1          1 ...
0 1   2     3      4       5        6         7         8          9 ...
0 1   3     6     10      15       21        28        36         45 ...
0 1   6    18     40      75      126       196       288        405 ...
0 1  10    45    136     325      666      1225      2080       3321 ...
0 1  20   135    544    1625     3996      8575     16640      29889 ...
0 1  36   378   2080    7875    23436     58996    131328     266085 ...
0 1  72  1134   8320   39375   140616    412972   1050624    2394765 ...
0 1 136  3321  32896  195625   840456   2883601   8390656   21526641 ...
0 1 272  9963 131584  978125  5042736  20185207  67125248  193739769 ...
0 1 528 29646 524800 4884375 30236976 141246028 536887296 1743421725 ...
...
		

References

Crossrefs

Columns 0-6 are A000007, A000012, A005418(n+1), A032120, A032121, A032122, A056308.
Rows 0-20 are A000012, A001477, A000217 (triangular numbers), A002411 (pentagonal pyramidal numbers), A037270, A168178, A071232, A168194, A071231, A168372, A071236, A168627, A071235, A168663, A168664, A170779, A170780, A170790, A170791, A170801, A170802.
Main diagonal is A275549.
Transpose is A284979.
Cf. A003992 (oriented), A293500 (chiral), A321391 (achiral).

Programs

  • Magma
    [[n le 0 select 1 else ((n-k)^k + (n-k)^Ceiling(k/2))/2: k in [0..n]]: n in [0..15]]; // G. C. Greubel, Nov 15 2018
  • Mathematica
    Table[If[n>0, ((n-k)^k + (n-k)^Ceiling[k/2])/2, 1], {n, 0, 15}, {k, 0, n}] // Flatten (* updated Jul 10 2018 *) (* Adapted to T(0,k)=1 by Robert A. Russell, Nov 13 2018 *)
  • PARI
    for(n=0,15, for(k=0,n, print1(if(n==0,1, ((n-k)^k + (n-k)^ceil(k/2))/2), ", "))) \\ G. C. Greubel, Nov 15 2018
    
  • PARI
    T(n,k) = {(k^n + k^ceil(n/2)) / 2} \\ Andrew Howroyd, Sep 13 2019
    

Formula

T(n,k) = [n==0] + [n>0] * (k^n + k^ceiling(n/2)) / 2. [Adapted to T(0,k)=1 by Robert A. Russell, Nov 13 2018]
G.f. for column k: (1 - binomial(k+1,2)*x^2) / ((1-k*x)*(1-k*x^2)). - Petros Hadjicostas, Jul 07 2018 [Adapted to T(0,k)=1 by Robert A. Russell, Nov 13 2018]
From Robert A. Russell, Nov 13 2018: (Start)
T(n,k) = (A003992(k,n) + A321391(n,k)) / 2.
T(n,k) = A003992(k,n) - A293500(n,k) = A293500(n,k) + A321391(n,k).
G.f. for row n: (Sum_{j=0..n} S2(n,j)*j!*x^j/(1-x)^(j+1) + Sum_{j=0..ceiling(n/2)} S2(ceiling(n/2),j)*j!*x^j/(1-x)^(j+1)) / 2, where S2 is the Stirling subset number A008277.
G.f. for row n>0: x*Sum_{k=0..n-1} A145882(n,k) * x^k / (1-x)^(n+1).
E.g.f. for row n: (Sum_{k=0..n} S2(n,k)*x^k + Sum_{k=0..ceiling(n/2)} S2(ceiling(n/2),k)*x^k) * exp(x) / 2, where S2 is the Stirling subset number A008277.
T(0,k) = 1; T(1,k) = k; T(2,k) = binomial(k+1,2); for n>2, T(n,k) = k*(T(n-3,k)+T(n-2,k)-k*T(n-1,k)).
For k>n, T(n,k) = Sum_{j=1..n+1} -binomial(j-n-2,j) * T(n,k-j). (End)

Extensions

Array transposed for greater consistency by Andrew Howroyd, Apr 04 2017
Origin changed to T(0,0) by Robert A. Russell, Nov 13 2018

A132592 X-values of solutions to the equation X*(X + 1) - 8*Y^2 = 0.

Original entry on oeis.org

0, 8, 288, 9800, 332928, 11309768, 384199200, 13051463048, 443365544448, 15061377048200, 511643454094368, 17380816062160328, 590436102659356800, 20057446674355970888, 681362750825443653408, 23146276081390728245000, 786292024016459316676608, 26710782540478226038759688
Offset: 0

Author

Mohamed Bouhamida, Nov 14 2007

Keywords

Comments

Equivalently, numbers k such that both k/2 and k+1 are squares. - Karl-Heinz Hofmann, Sep 20 2022
Equivalently, numbers k such that the k-dimensional volume and total (k-1)-dimensional volume are equal, with side length being a positive integer, for all regular polyhedra constructible in k dimensions. - Matt Moir, Jul 09 2024

Programs

Formula

a(0)=0, a(1)=8 and a(n) = 34*a(n-1) - a(n-2) + 16.
a(n) = (A056771(n) - 1)/2. - Max Alekseyev, Nov 13 2009
a(n) = sinh(2*n*arccosh(sqrt(2))^2) (n=0,1,2,3,...). - Artur Jasinski, Feb 10 2010
G.f.: -8*x*(x+1)/((x-1)*(x^2-34*x+1)). - Colin Barker, Oct 24 2012
a(n) = A055792(n+1)-1 = A001541(n)^2 - 1. - Antti Karttunen, Oct 03 2016

Extensions

More terms from Max Alekseyev, Nov 13 2009
Showing 1-10 of 51 results. Next