cp's OEIS Frontend

4*a(n) + 1 are the odd squares A016754(n).

The word "pronic" (used by Dickson) is incorrect. - Michael Somos

According to the 2nd edition of Webster, the correct word is "promic". - R. K. Guy

a(n) is the number of minimal vectors in the root lattice A_n (see Conway and Sloane, p. 109).

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 - a(n)*x - A002415(n)). - Benoit Cloitre, Nov 09 2002

The greatest LCM of all pairs (j, k) for j < k <= n for n > 1. - Robert G. Wilson v, Jun 19 2004

First differences are a(n+1) - a(n) = 2*n + 2 = 2, 4, 6, ... (while first differences of the squares are (n+1)^2 - n^2 = 2*n + 1 = 1, 3, 5, ...). - Alexandre Wajnberg, Dec 29 2005

25 appended to these numbers corresponds to squares of numbers ending in 5 (i.e., to squares of A017329). - Lekraj Beedassy, Mar 24 2006

A rapid (mental) multiplication/factorization technique -- a generalization of Lekraj Beedassy's comment: For all bases b >= 2 and positive integers n, c, d, k with c + d = b^k, we have (n*b^k + c)*(n*b^k + d) = a(n)*b^(2*k) + c*d. Thus the last 2*k base-b digits of the product are exactly those of c*d -- including leading 0(s) as necessary -- with the preceding base-b digit(s) the same as a(n)'s. Examples: In decimal, 113*117 = 13221 (as n = 11, b = 10 = 3 + 7, k = 1, 3*7 = 21, and a(11) = 132); in octal, 61*67 = 5207 (52 is a(6) in octal). In particular, for even b = 2*m (m > 0) and c = d = m, such a product is a square of this type. Decimal factoring: 5609 is immediately seen to be 71*79. Likewise, 120099 = 301*399 (k = 2 here) and 99990000001996 = 9999002*9999998 (k = 3). - Rick L. Shepherd, Jul 24 2021

Number of circular binary words of length n + 1 having exactly one occurrence of 01. Example: a(2) = 6 because we have 001, 010, 011, 100, 101 and 110. Column 1 of A119462. - Emeric Deutsch, May 21 2006

The sequence of iterated square roots sqrt(N + sqrt(N + ...)) has for N = 1, 2, ... the limit (1 + sqrt(1 + 4*N))/2. For N = a(n) this limit is n + 1, n = 1, 2, .... For all other numbers N, N >= 1, this limit is not a natural number. Examples: n = 1, a(1) = 2: sqrt(2 + sqrt(2 + ...)) = 1 + 1 = 2; n = 2, a(2) = 6: sqrt(6 + sqrt(6 + ...)) = 1 + 2 = 3. - Wolfdieter Lang, May 05 2006

Nonsquare integers m divisible by ceiling(sqrt(m)), except for m = 0. - Max Alekseyev, Nov 27 2006

The number of off-diagonal elements of an (n + 1) X (n + 1) matrix. - Artur Jasinski, Jan 11 2007

a(n) is equal to the number of functions f:{1, 2} -> {1, 2, ..., n + 1} such that for a fixed x in {1, 2} and a fixed y in {1, 2, ..., n + 1} we have f(x) <> y. - Aleksandar M. Janjic and Milan Janjic, Mar 13 2007

Numbers m >= 0 such that round(sqrt(m+1)) - round(sqrt(m)) = 1. - Hieronymus Fischer, Aug 06 2007

Numbers m >= 0 such that ceiling(2*sqrt(m+1)) - 1 = 1 + floor(2*sqrt(m)). - Hieronymus Fischer, Aug 06 2007

Numbers m >= 0 such that fract(sqrt(m+1)) > 1/2 and fract(sqrt(m)) < 1/2 where fract(x) is the fractional part (fract(x) = x - floor(x), x >= 0). - Hieronymus Fischer, Aug 06 2007

X values of solutions to the equation 4*X^3 + X^2 = Y^2. To find Y values: b(n) = n(n+1)(2n+1). - Mohamed Bouhamida, Nov 06 2007

Nonvanishing diagonal of A132792, the infinitesimal Lah matrix, so "generalized factorials" composed of a(n) are given by the elements of the Lah matrix, unsigned A111596, e.g., a(1)*a(2)*a(3) / 3! = -A111596(4,1) = 24. - Tom Copeland, Nov 20 2007

If Y is a 2-subset of an n-set X then, for n >= 2, a(n-2) is the number of 2-subsets and 3-subsets of X having exactly one element in common with Y. - Milan Janjic, Dec 28 2007

a(n) coincides with the vertex of a parabola of even width in the Redheffer matrix, directed toward zero. An integer p is prime if and only if for all integer k, the parabola y = kx - x^2 has no integer solution with 1 < x < k when y = p; a(n) corresponds to odd k. - Reikku Kulon, Nov 30 2008

The third differences of certain values of the hypergeometric function 3F2 lead to the squares of the oblong numbers i.e., 3F2([1, n + 1, n + 1], [n + 2, n + 2], z = 1) - 3*3F2([1, n + 2, n + 2], [n + 3, n + 3], z = 1) + 3*3F2([1, n + 3, n + 3], [n + 4, n + 4], z = 1) - 3F2([1, n + 4, n + 4], [n + 5, n + 5], z = 1) = (1/((n+2)*(n+3)))^2 for n = -1, 0, 1, 2, ... . See also A162990. - Johannes W. Meijer, Jul 21 2009

Generalized factorials, [a.(n!)] = a(n)*a(n-1)*...*a(0) = A010790(n), with a(0) = 1 are related to A001263. - Tom Copeland, Sep 21 2011

For n > 1, a(n) is the number of functions f:{1, 2} -> {1, ..., n + 2} where f(1) > 1 and f(2) > 2. Note that there are n + 1 possible values for f(1) and n possible values for f(2). For example, a(3) = 12 since there are 12 functions f from {1, 2} to {1, 2, 3, 4, 5} with f(1) > 1 and f(2) > 2. - Dennis P. Walsh, Dec 24 2011

a(n) gives the number of (n + 1) X (n + 1) symmetric (0, 1)-matrices containing two ones (see [Cameron]). - L. Edson Jeffery, Feb 18 2012

a(n) is the number of positions of a domino in a rectangled triangular board with both legs equal to n + 1. - César Eliud Lozada, Sep 26 2012

a(n) is the number of ordered pairs (x, y) in [n+2] X [n+2] with |x-y| > 1. - Dennis P. Walsh, Nov 27 2012

a(n) is the number of injective functions from {1, 2} into {1, 2, ..., n + 1}. - Dennis P. Walsh, Nov 27 2012

a(n) is the sum of the positive differences of the partition parts of 2n + 2 into exactly two parts (see example). - Wesley Ivan Hurt, Jun 02 2013

a(n)/a(n-1) is asymptotic to e^(2/n). - Richard R. Forberg, Jun 22 2013

Number of positive roots in the root system of type D_{n + 1} (for n > 2). - Tom Edgar, Nov 05 2013

Number of roots in the root system of type A_n (for n > 0). - Tom Edgar, Nov 05 2013

From Felix P. Muga II, Mar 18 2014: (Start)

a(m), for m >= 1, are the only positive integer values t for which the Binet-de Moivre formula for the recurrence b(n) = b(n-1) + t*b(n-2) with b(0) = 0 and b(1) = 1 has a root of a square. PROOF (as suggested by Wolfdieter Lang, Mar 26 2014): The sqrt(1 + 4t) appearing in the zeros r1 and r2 of the characteristic equation is (a positive) integer for positive integer t precisely if 4t + 1 = (2m + 1)^2, that is t = a(m), m >= 1. Thus, the characteristic roots are integers: r1 = m + 1 and r2 = -m.

Let m > 1 be an integer. If b(n) = b(n-1) + a(m)*b(n-2), n >= 2, b(0) = 0, b(1) = 1, then lim_{n->oo} b(n+1)/b(n) = m + 1. (End)

Cf. A130534 for relations to colored forests, disposition of flags on flagpoles, and colorings of the vertices (chromatic polynomial) of the complete graphs (here simply K_2). - Tom Copeland, Apr 05 2014

The set of integers k for which k + sqrt(k + sqrt(k + sqrt(k + sqrt(k + ...) ... is an integer. - Leslie Koller, Apr 11 2014

a(n-1) is the largest number k such that (n*k)/(n+k) is an integer. - Derek Orr, May 22 2014

Number of ways to place a domino and a singleton on a strip of length n - 2. - Ralf Stephan, Jun 09 2014

With offset 1, this appears to give the maximal number of crossings between n nonconcentric circles of equal radius. - Felix Fröhlich, Jul 14 2014

For n > 1, the harmonic mean of the n values a(1) to a(n) is n + 1. The lowest infinite sequence of increasing positive integers whose cumulative harmonic mean is integral. - Ian Duff, Feb 01 2015

a(n) is the maximum number of queens of one color that can coexist without attacking one queen of the opponent's color on an (n+2) X (n+2) chessboard. The lone queen can be placed in any position on the perimeter of the board. - Bob Selcoe, Feb 07 2015

With a(0) = 1, a(n-1) is the smallest positive number not in the sequence such that Sum_{i = 1..n} 1/a(i-1) has a denominator equal to n. - Derek Orr, Jun 17 2015

The positive members of this sequence are a proper subsequence of the so-called 1-happy couple products A007969. See the W. Lang link there, eq. (4), with Y_0 = 1, with a table at the end. - Wolfdieter Lang, Sep 19 2015

For n > 0, a(n) is the reciprocal of the area bounded above by y = x^(n-1) and below by y = x^n for x in the interval [0, 1]. Summing all such areas visually demonstrates the formula below giving Sum_{n >= 1} 1/a(n) = 1. - Rick L. Shepherd, Oct 26 2015

It appears that, except for a(0) = 0, this is the set of positive integers n such that x*floor(x) = n has no solution. (For example, to get 3, take x = -3/2.) - Melvin Peralta, Apr 14 2016

If two independent real random variables, x and y, are distributed according to the same exponential distribution: pdf(x) = lambda * exp(-lambda * x), lambda > 0, then the probability that n - 1 <= x/y < n is given by 1/a(n). - Andres Cicuttin, Dec 03 2016

a(n) is equal to the sum of all possible differences between n different pairs of consecutive odd numbers (see example). - Miquel Cerda, Dec 04 2016

a(n+1) is the dimension of the space of vector fields in the plane with polynomial coefficients up to order n. - Martin Licht, Dec 04 2016

It appears that a(n) + 3 is the area of the largest possible pond in a square (A268311). - Craig Knecht, May 04 2017

Also the number of 3-cycles in the (n+3)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017

Also the Wiener index of the (n+2)-wheel graph. - Eric W. Weisstein, Sep 08 2017

The left edge of a Floyd's triangle that consists of even numbers: 0; 2, 4; 6, 8, 10; 12, 14, 16, 18; 20, 22, 24, 26, 28; ... giving 0, 2, 6, 12, 20, ... The right edge generates A028552. - Waldemar Puszkarz, Feb 02 2018

a(n+1) is the order of rowmotion on a poset obtained by adjoining a unique minimal (or maximal) element to a disjoint union of at least two chains of n elements. - Nick Mayers, Jun 01 2018

From Juhani Heino, Feb 05 2019: (Start)

For n > 0, 1/a(n) = n/(n+1) - (n-1)/n.

For example, 1/6 = 2/3 - 1/2; 1/12 = 3/4 - 2/3.

Corollary of this:

Take 1/2 pill.

Next day, take 1/6 pill. 1/2 + 1/6 = 2/3, so your daily average is 1/3.

Next day, take 1/12 pill. 2/3 + 1/12 = 3/4, so your daily average is 1/4.

And so on. (End)

From Bernard Schott, May 22 2020: (Start)

For an oblong number m >= 6 there exists a Euclidean division m = d*q + r with q < r < d which are in geometric progression, in this order, with a common integer ratio b. For b >= 2 and q >= 1, the Euclidean division is m = qb*(qb+1) = qb^2 * q + qb where (q, qb, qb^2) are in geometric progression.

Some examples with distinct ratios and quotients:

6 | 4 30 | 25 42 | 18

----- ----- -----

2 | 1 , 5 | 1 , 6 | 2 ,

and also:

42 | 12 420 | 100

----- -----

6 | 3 , 20 | 4 .

Some oblong numbers also satisfy a Euclidean division m = d*q + r with q < r < d that are in geometric progression in this order but with a common noninteger ratio b > 1 (see A335064). (End)

For n >= 1, the continued fraction expansion of sqrt(a(n)) is [n; {2, 2n}]. For n=1, this collapses to [1; {2}]. - Magus K. Chu, Sep 09 2022

a(n-2) is the maximum irregularity over all trees with n vertices. The extremal graphs are stars. (The irregularity of a graph is the sum of the differences between the degrees over all edges of the graph.) - Allan Bickle, May 29 2023

For n > 0, number of diagonals in a regular 2*(n+1)-gon that are not parallel to any edge (cf. A367204). - Paolo Xausa, Mar 30 2024

a(n-1) is the maximum Zagreb index over all trees with n vertices. The extremal graphs are stars. (The Zagreb index of a graph is the sum of the squares of the degrees over all vertices of the graph.) - Allan Bickle, Apr 11 2024

For n >= 1, a(n) is the determinant of the distance matrix of a cycle graph on 2*n + 1 vertices (if the length of the cycle is even such a determinant is zero). - Miquel A. Fiol, Aug 20 2024

For n > 1, the continued fraction expansion of sqrt(16*a(n)) is [2n+1; {1, 2n-1, 1, 8n+2}]. - Magus K. Chu, Nov 20 2024

For n>=2, a(n) is the number of faces on a n+1-zone rhombic zonohedron. Each pair of a collection of great circles on a sphere intersects at two points, so there are 2*binomial(n+1,2) intersections. The dual of the implied polyhedron is a rhombic zonohedron, its faces corresponding to the intersections. - Shel Kaphan, Aug 12 2025

Also, 3 times triangular numbers, a(n) = 3*A000217(n).

In the 24-bit RGB color cube, the number of color-lattice-points in r+g+b = n planes at n < 256 equals the triangular numbers. For n = 256, ..., 765 the number of legitimate color partitions is less than A000217(n) because {r,g,b} components cannot exceed 255. For n = 256, ..., 511, the number of non-color partitions are computable with A045943(n-255), while for n = 512, ..., 765, the number of color points in r+g+b planes equals A000217(765-n). - Labos Elemer, Jun 20 2005

If a 3-set Y and an (n-3)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007

a(n) is also the smallest number that may be written both as the sum of n-1 consecutive positive integers and n consecutive positive integers. - Claudio Meller, Oct 08 2010

For n >= 3, a(n) equals 4^(2+n)*Pi^(1 - n) times the coefficient of zeta(3) in the following integral with upper bound Pi/4 and lower bound 0: int x^(n+1) tan x dx. - John M. Campbell, Jul 17 2011

The difference a(n)-a(n-1) = 3*n, for n >= 1. - Stephen Balaban, Jul 25 2011 [Comment clarified by N. J. A. Sloane, Aug 01 2024]

Sequence found by reading the line from 0, in the direction 0, 3, ..., and the same line from 0, in the direction 0, 9, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. This is one of the orthogonal axes of the spiral; the other is A032528. - Omar E. Pol, Sep 08 2011

A005449(a(n)) = A000332(3n + 3) = C(3n + 3, 4), a second pentagonal number of triangular matchstick number index number. Additionally, a(n) - 2n is a pentagonal number (A000326). - Raphie Frank, Dec 31 2012

Sum of the numbers from n to 2n. - Wesley Ivan Hurt, Nov 24 2015

Number of orbits of Aut(Z^7) as function of the infinity norm (n+1) of the representative integer lattice point of the orbit, when the cardinality of the orbit is equal to 5376 or 17920 or 20160. - Philippe A.J.G. Chevalier, Dec 28 2015

Also the number of 4-cycles in the (n+4)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017

Number of terms less than 10^k, k=0,1,2,3,...: 1, 3, 8, 26, 82, 258, 816, 2582, 8165, 25820, 81650, 258199, 816497, 2581989, 8164966, ... - Muniru A Asiru, Jan 24 2018

Numbers of the form 3*m*(2*m + 1) for m = 0, -1, 1, -2, 2, -3, 3, ... - Bruno Berselli, Feb 26 2018

Partial sums of A008585. - Omar E. Pol, Jun 20 2018

Column 1 of A273464. (Number of ways to select a unit lozenge inside an isosceles triangle of side length n; all vertices on a hexagonal lattice.) - R. J. Mathar, Jul 10 2019

Total number of pips in the n-th suit of a double-n domino set. - Ivan N. Ianakiev, Aug 23 2020

The number of edges of a complete tripartite graph of order 3n, K_n,n,n. - Roberto E. Martinez II, Oct 18 2001

From Floor van Lamoen, Jul 21 2001: (Start)

Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0,3,.... The spiral begins:

.

33--32--31--30

/ \

34 16--15--14 29

/ / \ \

35 17 5---4 13 28

/ / / \ \ \

36 18 6 0---3--12--27--48-->

/ / / / / / / /

37 19 7 1---2 11 26 47

\ \ \ / / /

38 20 8---9--10 25 46

\ \ / /

39 21--22--23--24 45

\ /

40--41--42--43--44

(End)

Number of edges of the complete bipartite graph of order 4n, K_n,3n. - Roberto E. Martinez II, Jan 07 2002

Also the number of partitions of 6n + 3 into at most 3 parts. - R. K. Guy, Oct 23 2003

Also the number of partitions of 6n into exactly 3 parts. - Colin Barker, Mar 23 2015

Numbers n such that the imaginary quadratic field Q[sqrt(-n)] has six units. - Marc LeBrun, Apr 12 2006

The denominators of Hoehn's sequence (recalled by G. L. Honaker, Jr.) and the numerators of that sequence reversed. The sequence is 1/3, (1+3)/(5+7), (1+3+5)/(7+9+11), (1+3+5+7)/(9+11+13+15), ...; reduced to 1/3, 4/12, 9/27, 16/48, ... . For the reversal, the reduction is 3/1, 12/4, 27/9, 48/16, ... . - Enoch Haga, Oct 05 2007

Right edge of tables in A200737 and A200741: A200737(n, A000292(n)) = A200741(n, A100440(n)) = a(n). - Reinhard Zumkeller, Nov 21 2011

The Wiener index of the crown graph G(n) (n>=3). The crown graph G(n) is the graph with vertex set {x(1), x(2), ..., x(n), y(1), y(2), ..., y(n)} and edge set {(x(i), y(j)): 1<=i, j<=n, i/=j} (= the complete bipartite graph K(n,n) with horizontal edges removed). Example: a(3)=27 because G(3) is the cycle C(6) and 6*1 + 6*2 + 3*3 = 27. The Hosoya-Wiener polynomial of G(n) is n(n-1)(t+t^2)+nt^3. - Emeric Deutsch, Aug 29 2013

From Michel Lagneau, May 04 2015: (Start)

Integer area A of equilateral triangles whose side lengths are in the commutative ring Z[3^(1/4)] = {a + b*3^(1/4) + c*3^(1/2) + d*3^(3/4), a,b,c and d in Z}.

The area of an equilateral triangle of side length k is given by A = k^2*sqrt(3)/4. In the ring Z[3^(1/4)], if k = q*3^(1/4), then A = 3*q^2/4 is an integer if q is even. Example: 27 is in the sequence because the area of the triangle (6*3^(1/4), 6*3^(1/4), 6*3^(1/4)) is 27. (End)

a(n) is 2*sqrt(3) times the area of a 30-60-90 triangle with short side n. Also, 3 times the area of an n X n square. - Wesley Ivan Hurt, Apr 06 2016

Consider the hexagonal tiling of the plane. Extract any four hexagons adjacent by edge. This can be done in three ways. Fold the four hexagons so that all opposite faces occupy parallel planes. For all parallel projections of the resulting object, at least two correspond to area a(n) for side length of n of the original hexagons. - Torlach Rush, Aug 17 2022

The sequence terms are the exponents in the expansion of Product_{n >= 1} (1 - q^(3*n))/(1 + q^(3*n)) = ( Sum_{n in Z} q^(n*(3*n+1)/2) ) / ( Product_{n >= 1} 1 + q^n ) = 1 - 2*q^3 + 2*q^12 - 2*q^27 + 2*q^48 - 2*q^75 + - .... - Peter Bala, Dec 30 2024

From Floor van Lamoen, Jul 21 2001: (Start)

Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0, 6, ...

The spiral begins:

85--84--83--82--81--80

/ \

86 56--55--54--53--52 79

/ / \ \

87 57 33--32--31--30 51 78

/ / / \ \ \

88 58 34 16--15--14 29 50 77

/ / / / \ \ \ \

89 59 35 17 5---4 13 28 49 76

/ / / / / \ \ \ \ \

<==90==60==36==18===6===0 3 12 27 48 75

/ / / / / / / / / /

61 37 19 7 1---2 11 26 47 74

\ \ \ \ / / / /

62 38 20 8---9--10 25 46 73

\ \ \ / / /

63 39 21--22--23--24 45 72

\ \ / /

64 40--41--42--43--44 71

\ /

65--66--67--68--69--70

(End)

If Y is a 4-subset of an n-set X then, for n >= 5, a(n-5) is the number of (n-4)-subsets of X having exactly two elements in common with Y. - Milan Janjic, Dec 28 2007

a(n) is the maximal number of points of intersection of n+1 distinct triangles drawn in the plane. For example, two triangles can intersect in at most a(1) = 6 points (as illustrated in the Star of David configuration). - Terry Stickels (Terrystickels(AT)aol.com), Jul 12 2008

Also sequence found by reading the line from 0, in the direction 0, 6, ... and the same line from 0, in the direction 0, 18, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. Axis perpendicular to A195143 in the same spiral. - Omar E. Pol, Sep 18 2011

Partial sums of A008588. - R. J. Mathar, Aug 28 2014

Also the number of 5-cycles in the (n+5)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017

a(n-4) is the maximum irregularity over all maximal 3-degenerate graphs with n vertices. The extremal graphs are 3-stars (K_3 joined to n-3 independent vertices). (The irregularity of a graph is the sum of the differences between the degrees over all edges of the graph.) - Allan Bickle, May 29 2023

Write 0,1,2,3,4,... in a triangular spiral; then a(n) is the sequence found by reading from 0 in the vertical upward direction.

Number of edges in the join of two complete graphs of order 2n and n, K_2n * K_n - Roberto E. Martinez II, Jan 07 2002

Column 3 of A048790.

Sequence found by reading the line from 0, in the direction 0, 3, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Sep 08 2011

a(n) is the sum of all perimeters of triangles having two sides of length n. For n=4 one has seven triangles with two sides of length 4 and the other of lengths 1..7. - J. M. Bergot, Mar 26 2014

a(n) is the Wiener index of the complete tripartite graph K_{n,n,n}. - Eric W. Weisstein, Sep 07 2017

Sequence found by reading the line from 0, in the direction 0, 3, ..., in a spiral on an equilateral triangular lattice. - Hans G. Oberlack, Dec 08 2018

a(n) also can be represented as n concentric triangles (see example). - Omar E. Pol, Aug 21 2011

Sequences of numbers of the form n*(n*k - k + 6)/2:

. k from 0 to 10, respectively: A008585, A055998, A005563, A045943, A014105, A005475, A033428, A022264, A033991, A062741, A147874;

. k=11: a(n);

. k=12: A094159;

. k=13: 0, 3, 19, 48, 90, 145, 213, 294, 388, 495, 615, 748, 894, ...;

. k=14: 0, 3, 20, 51, 96, 155, 228, 315, 416, 531, 660, 803, 960, ...;

. k=15: A152773;

. k=16: A139272;

. k=17: 0, 3, 23, 60, 114, 185, 273, 378, 500, 639, 795, 968, ...;

. k=18: A152751;

. k=19: 0, 3, 25, 66, 126, 205, 303, 420, 556, 711, 885, 1078, ...;

. k=20: 0, 3, 26, 69, 132, 215, 318, 441, 584, 747, 930, 1133, ...;

. k=21: A152759;

. k=22: 0, 3, 28, 75, 144, 235, 348, 483, 640, 819, 1020, 1243, ...;

. k=23: 0, 3, 29, 78, 150, 245, 363, 504, 668, 855, 1065, 1298, ...;

. k=24: A152767;

. k=25: 0, 3, 31, 84, 162, 265, 393, 546, 724, 927, 1155, 1408, ...;

. k=26: 0, 3, 32, 87, 168, 275, 408, 567, 752, 963, 1200, 1463, ...;

. k=27: A153783;

. k=28: A195021;

. k=29: 0, 3, 35, 96, 186, 305, 453, 630, 836, 1071, 1335, 1628, ...;

. k=30: A153448;

. k=31: 0, 3, 37, 102, 198, 325, 483, 672, 892, 1143, 1425, 1738, ...;

. k=32: 0, 3, 38, 105, 204, 335, 498, 693, 920, 1179, 1470, 1793, ...;

. k=33: A153875.

Also:

a(n) - n = A180223(n);

a(n) + n = n*(11*n-3)/2 = 0, 4, 19, 45, 82, 130, 189, 259, ...;

a(n) - 2*n = A051865(n);

a(n) + 2*n = A022268(n);

a(n) - 3*n = A152740(n-1);

a(n) + 3*n = A022269(n);

a(n) - 4*n = n*(11*n-13)/2 = 0, -1, 9, 30, 62, 105, 159, 224, ...;

a(n) + 4*n = A254963(n);

a(n) - n*(n-1)/2 = A147874(n+1);

a(n) + n*(n-1)/2 = A094159(n) (case k=12);

a(n) - n*(n-1) = A062741(n) (see above, this is the case k=9);

a(n) + n*(n-1) = n*(13*n-7)/2 (case k=13);

a(n) - n*(n+1)/2 = A135706(n);

a(n) + n*(n+1)/2 = A033579(n);

a(n) - n*(n+1) = A051682(n);

a(n) + n*(n+1) = A186030(n);

a(n) - n^2 = A062708(n);

a(n) + n^2 = n*(13*n-5)/2 = 0, 4, 21, 51, 94, 150, 219, ..., etc.

Sum of reciprocals of a(n), for n > 0: 0.47118857003113149692081665034891...

3*A172078(n) = n*a(n) - Sum_{k=0..n-1} a(k). - Bruno Berselli, Dec 12 2010