A048395 -id:A048395 - cp's OEIS frontend

A000578 The cubes: a(n) = n^3.

0, 1, 8, 27, 64, 125, 216, 343, 512, 729, 1000, 1331, 1728, 2197, 2744, 3375, 4096, 4913, 5832, 6859, 8000, 9261, 10648, 12167, 13824, 15625, 17576, 19683, 21952, 24389, 27000, 29791, 32768, 35937, 39304, 42875, 46656, 50653, 54872, 59319, 64000, 68921, 74088, 79507

Offset: 0

N. J. A. Sloane

a(n) is the sum of the next n odd numbers; i.e., group the odd numbers so that the n-th group contains n elements like this: (1), (3, 5), (7, 9, 11), (13, 15, 17, 19), (21, 23, 25, 27, 29), ...; then each group sum = n^3 = a(n). Also the median of each group = n^2 = mean. As the sum of first n odd numbers is n^2 this gives another proof of the fact that the n-th partial sum = (n(n + 1)/2)^2. - Amarnath Murthy, Sep 14 2002
Total number of triangles resulting from criss-crossing cevians within a triangle so that two of its sides are each n-partitioned. - Lekraj Beedassy, Jun 02 2004. See Propp and Propp-Gubin for a proof.
Also structured triakis tetrahedral numbers (vertex structure 7) (cf. A100175 = alternate vertex); structured tetragonal prism numbers (vertex structure 7) (cf. A100177 = structured prisms); structured hexagonal diamond numbers (vertex structure 7) (cf. A100178 = alternate vertex; A000447 = structured diamonds); and structured trigonal anti-diamond numbers (vertex structure 7) (cf. A100188 = structured anti-diamonds). Cf. A100145 for more on structured polyhedral numbers. - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Schlaefli symbol for this polyhedron: {4, 3}.
Least multiple of n such that every partial sum is a square. - Amarnath Murthy, Sep 09 2005
Draw a regular hexagon. Construct points on each side of the hexagon such that these points divide each side into equally sized segments (i.e., a midpoint on each side or two points on each side placed to divide each side into three equally sized segments or so on), do the same construction for every side of the hexagon so that each side is equally divided in the same way. Connect all such points to each other with lines that are parallel to at least one side of the polygon. The result is a triangular tiling of the hexagon and the creation of a number of smaller regular hexagons. The equation gives the total number of regular hexagons found where n = the number of points drawn + 1. For example, if 1 point is drawn on each side then n = 1 + 1 = 2 and a(n) = 2^3 = 8 so there are 8 regular hexagons in total. If 2 points are drawn on each side then n = 2 + 1 = 3 and a(n) = 3^3 = 27 so there are 27 regular hexagons in total. - Noah Priluck (npriluck(AT)gmail.com), May 02 2007
The solutions of the Diophantine equation: (X/Y)^2 - X*Y = 0 are of the form: (n^3, n) with n >= 1. The solutions of the Diophantine equation: (m^2)*(X/Y)^2k - XY = 0 are of the form: (m*n^(2k + 1), m*n^(2k - 1)) with m >= 1, k >= 1 and n >= 1. The solutions of the Diophantine equation: (m^2)*(X/Y)^(2k + 1) - XY = 0 are of the form: (m*n^(k + 1), m*n^k) with m >= 1, k >= 1 and n >= 1. - Mohamed Bouhamida, Oct 04 2007
Except for the first two terms, the sequence corresponds to the Wiener indices of C_{2n} i.e., the cycle on 2n vertices (n > 1). - K.V.Iyer, Mar 16 2009
Totally multiplicative sequence with a(p) = p^3 for prime p. - Jaroslav Krizek, Nov 01 2009
Sums of rows of the triangle in A176271, n > 0. - Reinhard Zumkeller, Apr 13 2010
One of the 5 Platonic polyhedral (tetrahedral, cube, octahedral, dodecahedral and icosahedral) numbers (cf. A053012). - Daniel Forgues, May 14 2010
Numbers n for which order of torsion subgroup t of the elliptic curve y^2 = x^3 - n is t = 2. - Artur Jasinski, Jun 30 2010
The sequence with the lengths of the Pisano periods mod k is 1, 2, 3, 4, 5, 6, 7, 8, 3, 10, 11, 12, 13, 14, 15, 16, 17, 6, 19, 20, ... for k >= 1, apparently multiplicative and derived from A000027 by dividing every ninth term through 3. Cubic variant of A186646. - R. J. Mathar, Mar 10 2011
The number of atoms in a bcc (body-centered cubic) rhombic hexahedron with n atoms along one edge is n^3 (T. P. Martin, Shells of atoms, eq. (8)). - Brigitte Stepanov, Jul 02 2011
The inverse binomial transform yields the (finite) 0, 1, 6, 6 (third row in A019538 and A131689). - R. J. Mathar, Jan 16 2013
Twice the area of a triangle with vertices at (0, 0), (t(n - 1), t(n)), and (t(n), t(n - 1)), where t = A000217 are triangular numbers. - J. M. Bergot, Jun 25 2013
If n > 0 is not congruent to 5 (mod 6) then A010888(a(n)) divides a(n). - Ivan N. Ianakiev, Oct 16 2013
For n > 2, a(n) = twice the area of a triangle with vertices at points (binomial(n,3),binomial(n+2,3)), (binomial(n+1,3),binomial(n+1,3)), and (binomial(n+2,3),binomial(n,3)). - J. M. Bergot, Jun 14 2014
Determinants of the spiral knots S(4,k,(1,1,-1)). a(k) = det(S(4,k,(1,1,-1))). - Ryan Stees, Dec 14 2014
One of the oldest-known examples of this sequence is shown in the Senkereh tablet, BM 92698, which displays the first 32 terms in cuneiform. - Charles R Greathouse IV, Jan 21 2015
From Bui Quang Tuan, Mar 31 2015: (Start)
We construct a number triangle from the integers 1, 2, 3, ... 2*n-1 as follows. The first column contains all the integers 1, 2, 3, ... 2*n-1. Each succeeding column is the same as the previous column but without the first and last items. The last column contains only n. The sum of all the numbers in the triangle is n^3.
Here is the example for n = 4, where 1 + 2*2 + 3*3 + 4*4 + 3*5 + 2*6 + 7 = 64 = a(4):
  1
  2  2
  3  3  3
  4  4  4  4
  5  5  5
  6  6
  7
(End)
For n > 0, a(n) is the number of compositions of n+11 into n parts avoiding parts 2 and 3. - Milan Janjic, Jan 07 2016
Does not satisfy Benford's law [Ross, 2012]. - N. J. A. Sloane, Feb 08 2017
Number of inequivalent face colorings of the cube using at most n colors such that each color appears at least twice. - David Nacin, Feb 22 2017
Consider A = {a,b,c} a set with three distinct members. The number of subsets of A is 8, including {a,b,c} and the empty set. The number of subsets from each of those 8 subsets is 27. If the number of such iterations is n, then the total number of subsets is a(n-1). - Gregory L. Simay, Jul 27 2018
By Fermat's Last Theorem, these are the integers of the form x^k with the least possible value of k such that x^k = y^k + z^k never has a solution in positive integers x, y, z for that k. - Felix Fröhlich, Jul 27 2018

			For k=3, b(3) = 2 b(2) - b(1) = 4-1 = 3, so det(S(4,3,(1,1,-1))) = 3*3^2 = 27.
For n=3, a(3) = 3 + (3*0^2 + 3*0 + 3*1^2 + 3*1 + 3*2^2 + 3*2) = 27. - _Patrick J. McNab_, Mar 28 2016

Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See p. 191.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 43, 64, 81.
R. L. Graham, D. E. Knuth, and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 255; 2nd. ed., p. 269. Worpitzky's identity (6.37).
Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 292.
T. Aaron Gulliver, "Sequences from cubes of integers", International Mathematical Journal, 4 (2003), no. 5, 439 - 445. See http://www.m-hikari.com/z2003.html for information about this journal. [I expanded the reference to make this easier to find. - N. J. A. Sloane, Feb 18 2019]
J. Propp and A. Propp-Gubin, "Counting Triangles in Triangles", Pi Mu Epsilon Journal (to appear).
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 6-7.
D. Wells, You Are A Mathematician, pp. 238-241, Penguin Books 1995.

N. J. A. Sloane, Table of n, a(n) for n = 0..10000
H. Bottomley, Illustration of initial terms
British National Museum, Tablet 92698
N. Brothers, S. Evans, L. Taalman, L. Van Wyk, D. Witczak, and C. Yarnall, Spiral knots, Missouri J. of Math. Sci., 22 (2010).
M. DeLong, M. Russell, and J. Schrock, Colorability and determinants of T(m,n,r,s) twisted torus knots for n equiv. +/-1(mod m), Involve, Vol. 8 (2015), No. 3, 361-384.
Ralph Greenberg, Math For Poets
R. K. Guy, The strong law of small numbers. Amer. Math. Monthly 95 (1988), no. 8, 697-712. [Annotated scanned copy]
Milan Janjic, Enumerative Formulas for Some Functions on Finite Sets [Cached version at the Wayback Machine]
Hyun Kwang Kim, On Regular Polytope Numbers, Proc. Amer. Math. Soc., 131 (2002), 65-75. - fixed by _Felix Fröhlich_, Jun 16 2014
T. P. Martin, Shells of atoms, Phys. Reports, 273 (1996), 199-241, eq. (8).
Ed Pegg, Jr., Sequence Pictures, Math Games column, Dec 08 2003.
Ed Pegg, Jr., Sequence Pictures, Math Games column, Dec 08 2003 [Cached copy, with permission (pdf only)]
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
James Propp and Adam Propp-Gubin, Counting Triangles in Triangles, arXiv:2409.17117 [math.CO], 25 September 2024.
Kenneth A. Ross, First Digits of Squares and Cubes, Math. Mag. 85 (2012) 36-42. doi:10.4169/math.mag.85.1.36.
Eric Weisstein's World of Mathematics, Cubic Number, and Hex Pyramidal Number
Ronald Yannone, Hilbert Matrix Analyses
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
Index entries for "core" sequences
Index entries for sequences related to Benford's law

(1/12)*t*(n^3-n)+n for t = 2, 4, 6, ... gives A004006, A006527, A006003, A005900, A004068, A000578, A004126, A000447, A004188, A004466, A004467, A007588, A062025, A063521, A063522, A063523.
For sums of cubes, cf. A000537 (partial sums), A003072, A003325, A024166, A024670, A101102 (fifth partial sums).
Cf. A001158 (inverse Möbius transform), A007412 (complement), A030078(n) (cubes of primes), A048766, A058645 (binomial transform), A065876, A101094, A101097.
Subsequence of A145784.
Cf. A260260 (comment). - Bruno Berselli, Jul 22 2015
Cf. A000292 (tetrahedral numbers), A005900 (octahedral numbers), A006566 (dodecahedral numbers), A006564 (icosahedral numbers).
Cf. A098737 (main diagonal).

a000578 = (^ 3)
a000578_list = 0 : 1 : 8 : zipWith (+)
   (map (+ 6) a000578_list)
   (map (* 3) $ tail $ zipWith (-) (tail a000578_list) a000578_list)
-- Reinhard Zumkeller, Sep 05 2015, May 24 2012, Oct 22 2011

[ n^3 : n in [0..50] ]; // Wesley Ivan Hurt, Jun 14 2014

I:=[0,1,8,27]; [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..45]]; // Vincenzo Librandi, Jul 05 2014

A000578 := n->n^3;
seq(A000578(n), n=0..50);
isA000578 := proc(r)
    local p;
    if r = 0 or r =1 then
        true;
    else
        for p in ifactors(r)[2] do
            if op(2, p) mod 3 <> 0 then
                return false;
            end if;
        end do:
        true ;
    end if;
end proc: # R. J. Mathar, Oct 08 2013

Table[n^3, {n, 0, 30}] (* Stefan Steinerberger, Apr 01 2006 *)
CoefficientList[Series[x (1 + 4 x + x^2)/(1 - x)^4, {x, 0, 45}], x] (* Vincenzo Librandi, Jul 05 2014 *)
Accumulate[Table[3n^2+3n+1,{n,0,20}]] (* or *) LinearRecurrence[{4,-6,4,-1},{1,8,27,64},20](* Harvey P. Dale, Aug 18 2018 *)

A000578(n):=n^3$
makelist(A000578(n),n,0,30); /* Martin Ettl, Nov 03 2012 */

A000578(n)=n^3 \\ M. F. Hasler, Apr 12 2008

is(n)=ispower(n,3) \\ Charles R Greathouse IV, Feb 20 2012

A000578_list, m = [], [6, -6, 1, 0]
for _ in range(10**2):
    A000578_list.append(m[-1])
    for i in range(3):
        m[i+1] += m[i] # Chai Wah Wu, Dec 15 2015

(define (A000578 n) (* n n n)) ;; Antti Karttunen, Oct 06 2017

a(n) = Sum_{i=0..n-1} A003215(i).
Multiplicative with a(p^e) = p^(3e). - David W. Wilson, Aug 01 2001
G.f.: x*(1+4*x+x^2)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
Dirichlet generating function: zeta(s-3). - Franklin T. Adams-Watters, Sep 11 2005, Amarnath Murthy, Sep 09 2005
E.g.f.: (1+3*x+x^2)*x*exp(x). - Franklin T. Adams-Watters, Sep 11 2005 - Amarnath Murthy, Sep 09 2005
a(n) = Sum_{i=1..n} (Sum_{j=i..n+i-1} A002024(j,i)). - Reinhard Zumkeller, Jun 24 2007
a(n) = lcm(n, (n - 1)^2) - (n - 1)^2. E.g.: lcm(1, (1 - 1)^2) - (1 - 1)^2 = 0, lcm(2, (2 - 1)^2) - (2 - 1)^2 = 1, lcm(3, (3 - 1)^2) - (3 - 1)^2 = 8, ... - Mats Granvik, Sep 24 2007
Starting (1, 8, 27, 64, 125, ...), = binomial transform of [1, 7, 12, 6, 0, 0, 0, ...]. - Gary W. Adamson, Nov 21 2007
a(n) = A007531(n) + A000567(n). - Reinhard Zumkeller, Sep 18 2009
a(n) = binomial(n+2,3) + 4*binomial(n+1,3) + binomial(n,3). [Worpitzky's identity for cubes. See. e.g., Graham et al., eq. (6.37). - Wolfdieter Lang, Jul 17 2019]
a(n) = n + 6*binomial(n+1,3) = binomial(n,1)+6*binomial(n+1,3). - Ron Knott, Jun 10 2019
A010057(a(n)) = 1. - Reinhard Zumkeller, Oct 22 2011
a(n) = A000537(n) - A000537(n-1), difference between 2 squares of consecutive triangular numbers. - Pierre CAMI, Feb 20 2012
a(n) = A048395(n) - 2*A006002(n). - J. M. Bergot, Nov 25 2012
a(n) = 1 + 7*(n-1) + 6*(n-1)*(n-2) + (n-1)*(n-2)*(n-3). - Antonio Alberto Olivares, Apr 03 2013
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 6. - Ant King Apr 29 2013
a(n) = A000330(n) + Sum_{i=1..n-1} A014105(i), n >= 1. - Ivan N. Ianakiev, Sep 20 2013
a(k) = det(S(4,k,(1,1,-1))) = k*b(k)^2, where b(1)=1, b(2)=2, b(k) = 2*b(k-1) - b(k-2) = b(2)*b(k-1) - b(k-2). - Ryan Stees, Dec 14 2014
For n >= 1, a(n) = A152618(n-1) + A033996(n-1). - Bui Quang Tuan, Apr 01 2015
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Jon Tavasanis, Feb 21 2016
a(n) = n + Sum_{j=0..n-1} Sum_{k=1..2} binomial(3,k)*j^(3-k). - Patrick J. McNab, Mar 28 2016
a(n) = A000292(n-1) * 6 + n. - Zhandos Mambetaliyev, Nov 24 2016
a(n) = n*binomial(n+1, 2) + 2*binomial(n+1, 3) + binomial(n,3). - Tony Foster III, Nov 14 2017
From Amiram Eldar, Jul 02 2020: (Start)
Sum_{n>=1} 1/a(n) = zeta(3) (A002117).
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*zeta(3)/4 (A197070). (End)
From Amiram Eldar, Jan 20 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)/Pi.
Product_{n>=2} (1 - 1/a(n)) = cosh(sqrt(3)*Pi/2)/(3*Pi). (End)
a(n) = Sum_{d|n} sigma_3(d)*mu(n/d) = Sum_{d|n} A001158(d)*A008683(n/d). Moebius transform of sigma_3(n). - Ridouane Oudra, Apr 15 2021

A001844 Centered square numbers: a(n) = 2n(n+1)+1. Sums of two consecutive squares. Also, consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z; then sequence gives Z values.

Original entry on oeis.org

1, 5, 13, 25, 41, 61, 85, 113, 145, 181, 221, 265, 313, 365, 421, 481, 545, 613, 685, 761, 841, 925, 1013, 1105, 1201, 1301, 1405, 1513, 1625, 1741, 1861, 1985, 2113, 2245, 2381, 2521, 2665, 2813, 2965, 3121, 3281, 3445, 3613, 3785, 3961, 4141, 4325, 4513

Offset: 0

N. J. A. Sloane

These are Hogben's central polygonal numbers denoted by
...2...
....P..
...4.n.
Numbers of the form (k^2+1)/2 for k odd.
(y(2x+1))^2 + (y(2x^2+2x))^2 = (y(2x^2+2x+1))^2. E.g., let y = 2, x = 1; (2(2+1))^2 + (2(2+2))^2 = (2(2+2+1))^2, (2(3))^2 + (2(4))^2 = (2(5))^2, 6^2 + 8^2 = 10^2, 36 + 64 = 100. - Glenn B. Cox (igloos_r_us(AT)canada.com), Apr 08 2002
a(n) is also the number of 3 X 3 magic squares with sum 3(n+1). - Sharon Sela (sharonsela(AT)hotmail.com), May 11 2002
For n > 0, a(n) is the smallest k such that zeta(2) - Sum_{i=1..k} 1/i^2 <= zeta(3) - Sum_{i=1..n} 1/i^3. - Benoit Cloitre, May 17 2002
Number of convex polyominoes with a 2 X (n+1) minimal bounding rectangle.
The prime terms are given by A027862. - Lekraj Beedassy, Jul 09 2004
First difference of a(n) is 4n = A008586(n). Any entry k of the sequence is followed by k + 2*(1 + sqrt(2k - 1)). - Lekraj Beedassy, Jun 04 2006
Integers of the form 1 + x + x^2/2 (generating polynomial is Schur's polynomial as in A127876). - Artur Jasinski, Feb 04 2007
If X is an n-set and Y and Z disjoint 2-subsets of X then a(n-4) is equal to the number of 4-subsets of X intersecting both Y and Z. - Milan Janjic, Aug 26 2007
Row sums of triangle A132778. - Gary W. Adamson, Sep 02 2007
Binomial transform of [1, 4, 4, 0, 0, 0, ...]; = inverse binomial transform of A001788: (1, 6, 24, 80, 240, ...). - Gary W. Adamson, Sep 02 2007
Narayana transform (A001263) of [1, 4, 0, 0, 0, ...]. Equals A128064 (unsigned) * [1, 2, 3, ...]. - Gary W. Adamson, Dec 29 2007
k such that the Diophantine equation x^3 - y^3 = x*y + k has a solution with y = x-1. If that solution is (x,y) = (m+1,m) then m^2 + (m+1)^2 = k. Note that this Diophantine equation is an elliptic curve and (m+1,m) is an integer point on it. - James R. Buddenhagen, Aug 12 2008
Numbers k such that (k, k, 2*k-2) are the sides of an isosceles triangle with integer area. Also, k such that 2*k-1 is a square. - James R. Buddenhagen, Oct 17 2008
a(n) is also the least weight of self-conjugate partitions having n+1 different odd parts. - Augustine O. Munagi, Dec 18 2008
Prefaced with a "1": (1, 1, 5, 13, 25, 41, ...) = A153869 * (1, 2, 3, ...). - Gary W. Adamson, Jan 03 2009
Prefaced with a "1": (1, 1, 5, 13, 25, 41, ...) where a(n) = 2n*(n-1)+1, all tuples of square numbers (X-Y, X, X+Y) are produced by ((m*(a(n)-2n))^2, (m*a(n))^2, (m*(a(n)+2n-2))^2) where m is a whole number. - Doug Bell, Feb 27 2009
Equals (1, 2, 3, ...) convolved with (1, 3, 4, 4, 4, ...). E.g., a(3) = 25 = (1, 2, 3, 4) dot (4, 4, 3, 1) = (4 + 8 + 9 + 4). - Gary W. Adamson, May 01 2009
The running sum of squares taken two at a time. - Al Hakanson (hawkuu(AT)gmail.com), May 18 2009
Equals the odd integers convolved with (1, 2, 2, 2, ...). - Gary W. Adamson, May 25 2009
Equals the triangular numbers convolved with [1, 2, 1, 0, 0, 0, ...]. - Gary W. Adamson & Alexander R. Povolotsky, May 29 2009
When the positive integers are written in a square array by diagonals as in A038722, a(n) gives the numbers appearing on the main diagonal. - Joshua Zucker, Jul 07 2009
The finite continued fraction [n,1,1,n] = (2n+1)/(2n^2 + 2n + 1) = (2n+1)/a(n); and the squares of the first two denominators of the convergents = a(n). E.g., the convergents and value of [4,1,1,4] = 1/4, 1/5, 2/9, 9/41 where 4^2 + 5^2 = 41. - Gary W. Adamson, Jul 15 2010
From Keith Tyler, Aug 10 2010: (Start)
Running sum of A008574.
Square open pyramidal number; that is, the number of elements in a square pyramid of height (n) with only surface and no bottom nodes. (End)
For k>0, x^4 + x^2 + k factors over the integers iff sqrt(k) is in this sequence. - James R. Buddenhagen, Aug 15 2010
Create the simple continued fraction from Pythagorean triples to get [2n + 1; 2n^2 + 2n, 2n^2 + 2n + 1]; its value equals the rational number 2n + 1 + a(n) / (4n^4 + 8n^3 + 6n^2 + 2n + 1). - J. M. Bergot, Sep 30 2011
a(n), n >= 1, has in its prime number factorization only primes of the form 4*k+1, i.e., congruent to 1 (mod 4) (see A002144). This follows from the fact that a(n) is a primitive sum of two squares and odd. See Theorem 3.20, p. 164, in the given Niven-Zuckerman-Montgomery reference. E.g., a(3) = 25 = 5^2, a(6) = 85 = 5*17. - Wolfdieter Lang, Mar 08 2012
From Ant King, Jun 15 2012: (Start)
a(n) is congruent to 1 (mod 4) for all n.
The digital roots of the a(n) form a purely periodic palindromic 9-cycle 1, 5, 4, 7, 5, 7, 4, 5, 1.
The units' digits of the a(n) form a purely periodic palindromic 5-cycle 1, 5, 3, 5, 1.
(End)
Number of integer solutions (x,y) of |x| + |y| <= n. Geometrically: number of lattice points inside a square with vertices (n,0), (0,-n), (-n,0), (0,n). - César Eliud Lozada, Sep 18 2012
(a(n)-1)/a(n) = 2*x / (1+x^2) where x = n/(n+1). Note that in this form, this is the velocity-addition formula according to the special theory of relativity (two objects traveling at 1/(n+1) slower than c relative to each other appear to travel at 1/a(n) less than c to a stationary observer). - Christian N. K. Anderson, May 20 2013 [Corrected by Rémi Guillaume, May 22 2025]
A geometric curiosity: the envelope of the circles x^2 + (y-a(n)/2)^2 = ((2n+1)/2)^2 is the parabola y = x^2, the n=0 circle being the osculating circle at the parabola vertex. - Jean-François Alcover, Dec 02 2013
Draw n ellipses in the plane (n>0), any 2 meeting in 4 points; a(n-1) gives the number of internal regions into which the plane is divided (cf. A051890, A046092); a(n-1) = A051890(n) - 1 = A046092(n-1) + 1. - Jaroslav Krizek, Dec 27 2013
a(n) is also, of course, the scalar product of the 2-vector (n, n+1) (or (n+1, n)) with itself. The unique inverse of (n, n+1) as vector in the Clifford algebra over the Euclidean 2-space is (1/a(n))(0, n, n+1, 0) (similarly for the other vector). In general the unique inverse of such a nonzero vector v (odd element in Cl_2) is v^(-1) = (1/|v|^2) v. Note that the inverse with respect to the scalar product is not unique for any nonzero vector. See the P. Lounesto reference, sects. 1.7 - 1.12, pp. 7-14. See also the Oct 15 2014 comment in A147973. - Wolfdieter Lang, Nov 06 2014
Subsequence of A004431, for n >= 1. - Bob Selcoe, Mar 23 2016
Numbers k such that 2k - 1 is a perfect square. - Juri-Stepan Gerasimov, Apr 06 2016
The number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 574", based on the 5-celled von Neumann neighborhood. - Robert Price, May 13 2016
a(n) is the first integer in a sum of (2*n + 1)^2 consecutive integers that equals (2*n + 1)^4. - Patrick J. McNab, Dec 24 2016
Central elements of odd-length rows of the triangular array of positive integers. a(n) is the mean of the numbers in the (2*n + 1)-th row of this triangle. - David James Sycamore, Aug 01 2018
Intersection of A000982 and A080827. - David James Sycamore, Aug 07 2018
An off-diagonal of the array of Delannoy numbers, A008288, (or a row/column when the array is shown as a square). As such, this is one of the crystal ball sequences. - Jack W Grahl, Feb 15 2021 and Shel Kaphan, Jan 18 2023
a(n) appears as a solution to a "Riddler Express" puzzle on the FiveThirtyEight website. The Jan 21 2022 issue (problem) and the Jan 28 2022 issue (solution) present the following puzzle and include a proof. - Fold a square piece of paper in half, obtaining a rectangle. Fold again to obtain a square with 1/4 the size of the original square. Then make n cuts through the folded paper. a(n) is the greatest number of pieces of the unfolded paper after the cutting. - Manfred Boergens, Feb 22 2022
a(n) is (1/6) times the number of 2 X 2 triangles in the n-th order hexagram with 12*n^2 cells. - Donghwi Park, Feb 06 2024
If k is a centered square number, its index in this sequence is n = (sqrt(2k-1)-1)/2. - Rémi Guillaume, Mar 30 2025.
Row sums of the symmetric triangle of odd numbers [1]; [1, 3, 1]; [1, 3, 5, 3, 1]; [1, 3, 5, 7, 5, 3, 1]; .... - Marco Zárate, Jun 15 2025

			G.f.: 1 + 5*x + 13*x^2 + 25*x^3 + 41*x^4 + 61*x^5 + 85*x^6 + 113*x^7 + 145*x^8 + ...
The first few triples are (1,0,1), (3,4,5), (5,12,13), (7,24,25), ...
The first four such partitions, corresponding to n = 0,1,2,3, i.e., to a(n) = 1,5,13,25, are 1, 3+1+1, 5+3+3+1+1, 7+5+5+3+3+1+1. - _Augustine O. Munagi_, Dec 18 2008

T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 3.
A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, p. 125, 1964.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 50.
Pertti Lounesto, Clifford Algebras and Spinors, second edition, Cambridge University Press, 2001.
S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 483.
Ivan Niven, Herbert S. Zuckerman and Hugh L. Montgomery, An Introduction to the Theory Of Numbers, Fifth Edition, John Wiley and Sons, Inc., NY 1991.
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).
Travers et al., The Mysterious Lost Proof, Using Advanced Algebra, (1976), pp. 27.

T. D. Noe, Table of n, a(n) for n = 0..1000
M. Ahmed, J. De Loera and R. Hemmecke, Polyhedral Cones of Magic Cubes and Squares, arXiv:math/0201108 [math.CO], 2002.
U. Alfred, n and n+1 consecutive integers with equal sums of squares, Math. Mag., 35 (1962), 155-164.
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Section 2.3.
Paul Barry, Centered polygon numbers, heptagons and nonagons, and the Robbins numbers, arXiv:2104.01644 [math.CO], 2021.
Matthias Beck, Moshe Cohen, Jessica Cuomo, and Paul Gribelyuk, The number of "magic" squares and hypercubes, arXiv:math/0201013 [math.CO], 2002-2005.
Arthur T. Benjamin and Doron Zeilberger, Pythagorean Primes and Palindromic Continued Fractions, Electronic Journal of Combinatorial Number Theory, 5(1) 2005, #A30.
J. A. De Loera, D. C. Haws and M. Koppe, Ehrhart Polynomials of Matroid Polytopes and Polymatroids, arXiv:0710.4346 [math.CO], 2007; Discrete Comput. Geom., 42 (2009), 670-702.
FiveThirtyEight, "Riddler Express" paper cutting problem and solution, Jan 28 2022.
D. C. Haws, Matroids [Broken link, Oct 30 2017]
D. C. Haws, Matroids [Copy on website of Matthias Koeppe]
D. C. Haws, Matroids [Cached copy, pdf file only]
L. Hogben, Choice and Chance by Cardpack and Chessboard, Vol. 1, Max Parrish and Co, London, 1950, pp. 22 and 36.
Milan Janjic, Two Enumerative Functions. [Broken link; WayBackMachine archive.]
Milan Janjic, On a class of polynomials with integer coefficients, JIS 11 (2008) 08.5.
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
Clark Kimberling, Complementary Equations, Journal of Integer Sequences, Vol. 10 (2007), Article 07.1.4.
Clark Kimberling and John E. Brown, Partial Complements and Transposable Dispersions, J. Integer Seqs., Vol. 7, 2004.
Ron Knott, Pythagorean Triples and Online Calculators.
G. Kreweras, Sur les hiérarchies de segments, Cahiers Bureau Universitaire Recherche Opérationnelle, Cahier 20, Inst. Statistiques, Univ. Paris, 1973.
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973). (Annotated scanned copy)
A. O. Munagi, Pairing conjugate partitions by residue classes, Discrete Math., 308 (2008), 2492-2501.
Mitchell Paukner, Lucy Pepin, Manda Riehl, and Jarred Wieser, Pattern Avoidance in Task-Precedence Posets, arXiv:1511.00080 [math.CO], 2015-2016.
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
John A. Jr. Rochowicz, Harmonic Numbers: Insights, Approximations and Applications, Spreadsheets in Education (eJSiE), 2015, Vol. 8: Iss. 2, Article 4.
Amelia Carolina Sparavigna, Groupoid of OEIS A001844 Numbers (centered square numbers), Politecnico di Torino, Italy (2019).
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
David James Sycamore, Triangular array
Leo Tavares, Illustration: Diamond Rows
B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985), 4545-4558.
Eric Weisstein's World of Mathematics, Centered Polygonal Number, Centered Square Number, Diamond, Pythagorean Triple, and von Neumann Neighborhood.
Index entries for sequences related to centered polygonal numbers
Index entries for two-way infinite sequences.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
Nicolay Avilov, Graphical representation of the sequence members

X values are A005408; Y values are A046092.
Cf. A008586 (first differences), A005900 (partial sums), A254373 (digital roots).
Cf. A069894, A000217, A000290, A001263, A001788, A001845, A002061, A002144, A003215, A005448, A005891, A005917, A008844 (square terms), A027862 (prime terms), A048395, A051890, A056106, A101096, A127876, A128064, A132778, A147973, A153869, A240876, A251599, A000982, A080827.
Subsequence of A004431.
Right edge of A055096; main diagonal of A069480, A078475, A129312.
Row n=2 (or column k=2) of A008288.
Cf. A016754.

a001844 n = 2 * n * (n + 1) + 1
a001844_list = zipWith (+) a000290_list $ tail a000290_list
-- Reinhard Zumkeller, Dec 04 2012

[2*n^2 + 2*n + 1: n in [0..50]]; // Vincenzo Librandi, Jan 19 2013

[n: n in [0..4400] | IsSquare(2*n-1)]; // Juri-Stepan Gerasimov, Apr 06 2016

A001844:=-(z+1)**2/(z-1)**3; # Simon Plouffe in his 1992 dissertation

Table[2n(n + 1) + 1, {n, 0, 50}]
FoldList[#1 + #2 &, 1, 4 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
maxn := 47; Flatten[Table[SeriesCoefficient[Series[(n + (n - 1)*x)/(1 - x)^2, {x, 0, maxn}], k], {n, maxn}, {k, n - 1, n - 1}]] (* L. Edson Jeffery, Aug 24 2014 *)
CoefficientList[ Series[-(x^2 + 2x + 1)/(x - 1)^3, {x, 0, 48}], x] (* or *)
LinearRecurrence[{3, -3, 1}, {1, 5, 13}, 48] (* Robert G. Wilson v, Aug 01 2018 *)
Total/@Partition[Range[0,50]^2,2,1] (* Harvey P. Dale, Dec 05 2020 *)
Table[ j! Coefficient[Series[Exp[x]*(1 + 4*x + 2*x^2), {x, 0, 20}], x,
j], {j, 0, 20}] (* Nikolaos Pantelidis, Feb 07 2023 *)

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

x='x+O('x^200); Vec((1+x)^2/(1-x)^3) \\ Altug Alkan, Mar 23 2016

print([2*n*(n+1)+1 for n in range(48)]) # Michael S. Branicky, Jan 05 2021

[i**2 + (i + 1)**2 for i in range(46)] # Zerinvary Lajos, Jun 27 2008

a(n) = 2*n^2 + 2*n + 1 = n^2 + (n+1)^2.
a(n) = 1 + 3 + 5 + ... + 2*n-1 + 2*n+1 + 2*n-1 + ... + 3 + 1. - Amarnath Murthy, May 28 2001
a(n) = 1/real(z(n+1)) where z(1)=i, (i^2=-1), z(k+1) = 1/(z(k)+2i). - Benoit Cloitre, Aug 06 2002
Nearest integer to 1/Sum_{k>n} 1/k^3. - Benoit Cloitre, Jun 12 2003
G.f.: (1+x)^2/(1-x)^3.
E.g.f.: exp(x)*(1+4x+2x^2).
a(n) = a(n-1) + 4n.
a(-n) = a(n-1).
a(n) = A064094(n+3, n) (fourth diagonal).
a(n) = 1 + Sum_{j=0..n} 4*j. - Xavier Acloque, Oct 08 2003
a(n) = A046092(n)+1 = (A016754(n)+1)/2. - Lekraj Beedassy, May 25 2004
a(n) = Sum_{k=0..n+1} (-1)^k*binomial(n, k)*Sum_{j=0..n-k+1} binomial(n-k+1, j)*j^2. - Paul Barry, Dec 22 2004
a(n) = ceiling((2n+1)^2/2). - Paul Barry, Jul 16 2006
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=1, a(1)=5, a(2)=13. - Jaume Oliver Lafont, Dec 02 2008
a(n)*a(n-1) = 4*n^4 + 1 for n > 0. - Reinhard Zumkeller, Feb 12 2009
Prefaced with a "1" (1, 1, 5, 13, 25, 41, ...): a(n) = 2*n*(n-1)+1. - Doug Bell, Feb 27 2009
a(n) = sqrt((A056220(n)^2 + A056220(n+1)^2) / 2). - Doug Bell, Mar 08 2009
a(n) = floor(2*(n+1)^3/(n+2)). - Gary Detlefs, May 20 2010
a(n) = A000330(n) - A000330(n-2). - Keith Tyler, Aug 10 2010
a(n) = A069894(n)/2. - J. M. Bergot, Jun 11 2012
a(n) = 2*a(n-1) - a(n-2) + 4. - Ant King, Jun 12 2012
Sum_{n>=0} 1/a(n) = (Pi/2)*tanh(Pi/2) = 1.4406595199775... = A228048. - Ant King, Jun 15 2012
a(n) = A209297(2*n+1,n+1). - Reinhard Zumkeller, Jan 19 2013
a(n)^3 = A048395(n)^2 + A048395(-n-1)^2. - Vincenzo Librandi, Jan 19 2013
a(n) = A000217(2n+1) - n. - Ivan N. Ianakiev, Nov 08 2013
a(n) = A251599(3*n+1). - Reinhard Zumkeller, Dec 13 2014
a(n) = A101321(4,n). - R. J. Mathar, Jul 28 2016
From Ilya Gutkovskiy, Jul 30 2016: (Start)
a(n) = Sum_{k=0..n} A008574(k).
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = exp(-1) = A068985. (End)
a(n) = 4 * A000217(n) + 1. - Bruce J. Nicholson, Jul 10 2017
a(n) = A002522(n) + A005563(n) = A002522(n+1) + A005563(n-1). - Bruce J. Nicholson, Aug 05 2017
Sum_{n>=0} a(n)/n! = 7*e. Sum_{n>=0} 1/a(n) = A228048. - Amiram Eldar, Jun 20 2020
a(n) = A000326(n+1) + A000217(n-1). - Charlie Marion, Nov 16 2020
a(n) = Integral_{x=0..2n+2} |1-x| dx. - Pedro Caceres, Dec 29 2020
From Amiram Eldar, Feb 17 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)*sech(Pi/2).
Product_{n>=1} (1 - 1/a(n)) = Pi*csch(Pi)*sinh(Pi/2). (End)
a(n) = A001651(n+1) + 1 - A028242(n). - Charlie Marion, Apr 05 2022
a(n) = A016754(n) - A046092(n). - Leo Tavares, Sep 16 2022
For n>0, a(n) = A101096(n+2) / 30. - Andy Nicol, Feb 06 2025
From Rémi Guillaume, Apr 21 2025: (Start)
a(n) = (2*A003215(n)+1)/3.
a(n) = (4*A005448(n+1)-1)/3.
a(n) + a(n-1) = A001845(n) - A001845(n-1), for n >= 1.
a(n) = (A005917(n+1))/(2n+1). (End)

Partially edited by Joerg Arndt, Mar 11 2010

A000037 Numbers that are not squares (or, the nonsquares).

Original entry on oeis.org

2, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99

Offset: 1

N. J. A. Sloane, Simon Plouffe

Note the remarkable formula for the n-th term (see the FORMULA section)!
These are the natural numbers with an even number of divisors. The number of divisors is odd for the complementary sequence, the squares (sequence A000290) and the numbers for which the number of divisors is divisible by 3 is sequence A059269. - Ola Veshta (olaveshta(AT)my-deja.com), Apr 04 2001
a(n) is the largest integer m not equal to n such that n = (floor(n^2/m) + m)/2. - Alexander R. Povolotsky, Feb 10 2008
Union of A007969 and A007970; A007968(a(n)) > 0. - Reinhard Zumkeller, Jun 18 2011
Terms of even numbered rows in the triangle A199332. - Reinhard Zumkeller, Nov 23 2011
If a(n) and a(n+1) are of the same parity then (a(n)+a(n+1))/2 is a square. - Zak Seidov, Aug 13 2012
Theaetetus of Athens proved the irrationality of the square roots of these numbers in the 4th century BC. - Charles R Greathouse IV, Apr 18 2013
4*a(n) are the even members of A079896, the discriminants of indefinite binary quadratic forms. - Wolfdieter Lang, Jun 14 2013

			For example note that the squares 0, 1, 4, 9, 16 are not included.

Titu Andreescu, Dorin Andrica, and Zuming Feng, 104 Number Theory Problems, Birkhäuser, 2006, 58-60.
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).

Ray Chandler, Table of n, a(n) for n = 1..10000 (first 9900 terms from N. J. A. Sloane)
E. R. Berlekamp, A contribution to mathematical psychometrics, Unpublished Bell Labs Memorandum, Feb 08 1968 [Annotated scanned copy]
A. J. dos Reis and D. M. Silberger, Generating nonpowers by formula, Math. Mag., 63 (1990), 53-55.
Bakir Farhi, An explicit formula generating the non-Fibonacci numbers, arXiv:1105.1127 [math.NT], May 05 2011.
S. R. Finch, Class number theory
Steven R. Finch, Class number theory [Cached copy, with permission of the author]
Henry W. Gould, Letters to N. J. A. Sloane, Oct 1973 and Jan 1974.
S. Kaji, T. Maeno, K. Nuida, and Y. Numata, Polynomial Expressions of Carries in p-ary Arithmetics, arXiv preprint arXiv:1506.02742 [math.CO], 2015-2016.
J. Lambek and L. Moser, Inverse and complementary sequences of natural numbers, Amer. Math. Monthly, 61 (1954), 454-458. doi 10.2307/2308078, see example 4 (includes the formula). [Nicolas Normand (Nicolas.Normand(AT)polytech.univ-nantes.fr), Nov 24 2009]
R. P. Loh, A. G. Shannon, and A. F. Horadam, Divisibility Criteria and Sequence Generators Associated with Fermat Coefficients, Preprint, 1980.
Cristinel Mortici, Remarks on Complementary Sequences, Fibonacci Quart. 48 (2010), no. 4, 343-347.
R. D. Nelson, Sequences which omit powers, The Mathematical Gazette, Number 461, 1988, pages 208-211.
M. A. Nyblom, Some curious sequences involving floor and ceiling functions, Am. Math. Monthly 109 (#6, 2002), 559-564.
Rosetta Code, Sequence of non-squares
J. Scholes, 27th Putnam 1966 Prob. A4
Aaron Snook, Augmented Integer Linear Recurrences, 2012. - From _N. J. A. Sloane_, Dec 19 2012
Eric Weisstein's World of Mathematics, Square Number
Eric Weisstein's World of Mathematics, Continued Fraction

Cf. A007412, A000005, A000290, A059269, A134986, A087153, A172151, A000196, A049068 (subsequence).
Cf. A242401 (subsequence).
Cf. A086849 (partial sums), A048395.

a000037 n = n + a000196 (n + a000196 n)
-- Reinhard Zumkeller, Nov 23 2011

[n : n in [1..1000] | not IsSquare(n) ];

at:=0; for n in [1..10000] do if not IsSquare(n) then at:=at+1; print at, n; end if; end for;

Maple
```
A000037 := n->n+floor(1/2+sqrt(n));
```

a[n_] := (n + Floor[Sqrt[n + Floor[Sqrt[n]]]]); Table[a[n], {n, 71}] (* Robert G. Wilson v, Sep 24 2004 *)
With[{upto=100},Complement[Range[upto],Range[Floor[Sqrt[upto]]]^2]] (* Harvey P. Dale, Dec 02 2011 *)
a[ n_] :=  If[ n < 0, 0, n + Round @ Sqrt @ n]; (* Michael Somos, May 28 2014 *)

A000037(n):=n + floor(1/2 + sqrt(n))$ makelist(A000037(n),n,1,50); /* Martin Ettl, Nov 15 2012 */

{a(n) = if( n<0, 0, n + (1 + sqrtint(4*n)) \ 2)};

from math import isqrt
def A000037(n): return n+isqrt(n+isqrt(n)) # Chai Wah Wu, Mar 31 2022

from math import isqrt
def A000037(n): return n+(k:=isqrt(n))+int(n>=k*(k+1)+1) # Chai Wah Wu, Jun 17 2024

a(n) = n + floor(1/2 + sqrt(n)).
a(n) = n + floor(sqrt( n + floor(sqrt n))).
A010052(a(n)) = 0. - Reinhard Zumkeller, Jan 26 2010
A173517(a(n)) = n; a(n)^2 = A030140(n). - Reinhard Zumkeller, Feb 20 2010
a(n) = A000194(n) + n. - Jaroslav Krizek, Jun 14 2009
a(A002061(n)) = a(n^2-n+1) = A002522(n) = n^2 + 1. - Jaroslav Krizek, Jun 21 2009

Edited by Charles R Greathouse IV, Oct 30 2009

cp's OEIS Frontend

A000578 The cubes: a(n) = n^3.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

Maxima

PARI

PARI

Python

Scheme

Formula

A001844 Centered square numbers: a(n) = 2*n*(n+1)+1. Sums of two consecutive squares. Also, consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z; then sequence gives Z values.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

PARI

Python

Sage

Formula

Extensions

A000037 Numbers that are not squares (or, the nonsquares).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

Maxima

PARI

Python

Python

Formula

Extensions

A086849 Sum of first n nonsquares.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

A199771 Row sums of the triangle in A199332.

Views

Author

A001844 Centered square numbers: a(n) = 2n(n+1)+1. Sums of two consecutive squares. Also, consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z; then sequence gives Z values.