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.

Previous Showing 21-30 of 56 results. Next

A054554 a(n) = 4*n^2 - 10*n + 7.

Original entry on oeis.org

1, 3, 13, 31, 57, 91, 133, 183, 241, 307, 381, 463, 553, 651, 757, 871, 993, 1123, 1261, 1407, 1561, 1723, 1893, 2071, 2257, 2451, 2653, 2863, 3081, 3307, 3541, 3783, 4033, 4291, 4557, 4831, 5113, 5403, 5701, 6007, 6321, 6643, 6973, 7311, 7657, 8011, 8373, 8743
Offset: 1

Views

Author

Keywords

Comments

Move in 1-3 direction in a spiral organized like A068225 etc.
Equals binomial transform of [1, 2, 8, 0, 0, 0, ...]. - Gary W. Adamson, May 03 2008
Ulam's spiral (NE spoke). - Robert G. Wilson v, Oct 31 2011
Number of ternary strings of length 2*(n-1) that have one or no 0's, one or no 1's, and an even number of 2's. For n=2, the 3 strings of length 2 are 01, 10 and 22. For n=3, the 13 strings of length 4 are the 12 permutations of 0122 and 2222. - Enrique Navarrete, Jul 25 2025

Crossrefs

Cf. A014105.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

Formula

a(n) = 8*n + a(n-1) - 14 with n > 1, a(1)=1. - Vincenzo Librandi, Aug 07 2010
G.f.: -x*(7*x^2+1)/(x-1)^3. - Colin Barker, Sep 21 2012
For n > 2, a(n) = A014105(n) + A014105(n-1). - Bruce J. Nicholson, May 07 2017
From Leo Tavares, Feb 21 2022: (Start)
a(n) = A003215(n-2) + 2*A000217(n-1). See Hexagonal Dual Rays illustration in links.
a(n) = A227776(n-1) - 4*A000217(n-1). (End)
a(k+1) = 4k^2 - 2k + 1 in the Numberphile video. - Frank Ellermann, Mar 11 2020
E.g.f.: exp(x)*(7 - 6*x + 4*x^2) - 7. - Stefano Spezia, Apr 24 2024

Extensions

Edited by Frank Ellermann, Feb 24 2002

A054556 a(n) = 4*n^2 - 9*n + 6.

Original entry on oeis.org

1, 4, 15, 34, 61, 96, 139, 190, 249, 316, 391, 474, 565, 664, 771, 886, 1009, 1140, 1279, 1426, 1581, 1744, 1915, 2094, 2281, 2476, 2679, 2890, 3109, 3336, 3571, 3814, 4065, 4324, 4591, 4866, 5149, 5440, 5739, 6046, 6361, 6684, 7015, 7354, 7701, 8056, 8419, 8790
Offset: 1

Views

Author

Keywords

Comments

Move in 1-4 direction in a spiral organized like A068225 etc.
Equals binomial transform of [1, 3, 8, 0, 0, 0, ...]. - Gary W. Adamson, Apr 30 2008
Ulam's spiral (N spoke). - Robert G. Wilson v, Oct 31 2011
Also, numbers of the form m*(4*m+1)+1 for nonpositive m. - Bruno Berselli, Jan 06 2016

Crossrefs

Cf. A266883: m*(4*m+1)+1 for m = 0,-1,1,-2,2,-3,3,...
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

Formula

a(n)^2 = Sum_{i = 0..2*(4*n-5)} (4*n^2-13*n+9+i)^2*(-1)^i = ((n-1)*(4*n-5)+1)^2. - Bruno Berselli, Apr 29 2010
From Harvey P. Dale, Aug 21 2011: (Start)
a(0)=1, a(1)=4, a(2)=15; for n > 2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
G.f.: -x*(6*x^2+x+1)/(x-1)^3. (End)
From Franck Maminirina Ramaharo, Mar 09 2018: (Start)
a(n) = binomial(2*n - 2, 2) + 2*(n - 1)^2 + 1.
a(n) = A000384(n-1) + A058331(n-1).
a(n) = A130883(n-1) + A001105(n-1). (End)
E.g.f.: exp(x)*(6 - 5*x + 4*x^2) - 6. - Stefano Spezia, Apr 24 2024

Extensions

Edited by Frank Ellermann, Feb 24 2002
Incorrect formula deleted by N. J. A. Sloane, Aug 02 2009

A054567 a(n) = 4*n^2 - 7*n + 4.

Original entry on oeis.org

1, 6, 19, 40, 69, 106, 151, 204, 265, 334, 411, 496, 589, 690, 799, 916, 1041, 1174, 1315, 1464, 1621, 1786, 1959, 2140, 2329, 2526, 2731, 2944, 3165, 3394, 3631, 3876, 4129, 4390, 4659, 4936, 5221, 5514, 5815, 6124, 6441, 6766, 7099, 7440, 7789, 8146, 8511, 8884
Offset: 1

Views

Author

Keywords

Comments

The number 1 is placed in the middle of a sheet of squared paper and the numbers 2, 3, 4, 5, 6, etc. are written in a clockwise spiral around 1, as in A068225 etc. This sequence is read off along one of the rays from 1.
Ulam's spiral (W spoke of A054552). - Robert G. Wilson v, Oct 31 2011
Also, numbers of the form m*(4*m+1)+1 for nonnegative m. - Bruno Berselli, Jan 06 2016
The sequence forms the 1x2 diagonal of the square maze arrangement in A081344. - Jarrod G. Sage, Jul 17 2024

Crossrefs

Cf. A266883: m*(4*m+1)+1 for m = 0,-1,1,-2,2,-3,3,...
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

Formula

a(n) = 8*n+a(n-1)-11 for n>1, a(1)=1. - Vincenzo Librandi, Aug 07 2010
a(n) = A204674(n-1) / n. - Reinhard Zumkeller, Jan 18 2012
From Colin Barker, Oct 25 2014: (Start)
a(n) = 3*a(n-1)-3*a(n-2)+a(n-3).
G.f.: -x*(4*x^2+3*x+1) / (x-1)^3. (End)
E.g.f.: exp(x)*(4 - 3*x + 4*x^2) - 4. - Stefano Spezia, Apr 24 2024
a(n) = A016742(n-1) + n. - Jarrod G. Sage, Jul 17 2024

Extensions

Edited by Frank Ellermann, Feb 24 2002
Typo fixed by Charles R Greathouse IV, Oct 28 2009

A056108 Fourth spoke of a hexagonal spiral.

Original entry on oeis.org

1, 5, 15, 31, 53, 81, 115, 155, 201, 253, 311, 375, 445, 521, 603, 691, 785, 885, 991, 1103, 1221, 1345, 1475, 1611, 1753, 1901, 2055, 2215, 2381, 2553, 2731, 2915, 3105, 3301, 3503, 3711, 3925, 4145, 4371, 4603, 4841, 5085, 5335, 5591, 5853, 6121, 6395
Offset: 0

Views

Author

Henry Bottomley, Jun 09 2000

Keywords

Comments

a(n) = sum of (n+1)-th row terms of triangle A134234. - Gary W. Adamson, Oct 14 2007
If Y is a 4-subset of an n-set X then, for n >= 4, a(n-4) is the number of 4-subsets of X having at least two elements in common with Y. - Milan Janjic, Dec 08 2007
Equals binomial transform of [1, 4, 6, 0, 0, 0, ...] - Gary W. Adamson, Apr 30 2008
From A.K. Devaraj, Sep 18 2009: (Start)
Let f(x) be a polynomial in x. Then f(x + n*f(x)) is congruent to 0 (mod f(x)); here n belongs to N.
There is nothing interesting in the quotients f(x + n*f(x))/f(x) when x belongs to Z.
However, when x is irrational these quotients consist of two parts, a) rational integers and b) integer multiples of x.
The present sequence is the integer part when the polynomial is x^2 + x + 1 and x = sqrt(2),
f(x+n*f(x))/f(x) = a(n) + A005563(n)*sqrt(2).
Equals triangle A128229 as an infinite lower triangular matrix * A016777 as a vector, where A016777(n) = (3*n+1). (End)
Numbers of the form ((-h^2+h+1)^2+(h^2-h+1)^2+(h^2+h-1)^2)/(h^2+h+1) for h=n+1. - Bruno Berselli, Mar 13 2013

Crossrefs

Other spirals: A054552.

Programs

Formula

a(n) = 3*n^2 + n + 1.
a(n) = a(n-1) + 6*n - 2 = 2*a(n-1) - a(n-2) + 6
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = A056105(n) + 3*n = A056106(n) + 2*n = A056107(n) + n = A056109(n) - n = A003215(n) - 2*n.
a(n) = A096777(3n+1) . - Reinhard Zumkeller, Dec 29 2007
a(n) = 6*n+a(n-1)-2 with n>0, a(0)=1. - Vincenzo Librandi, Aug 07 2010
G.f.: (1+2*x+3*x^2)/(1-3*x+3*x^2-x^3). - Colin Barker, Jan 04 2012
a(-n) = A056106(n). - Bruno Berselli, Mar 13 2013
E.g.f.: (3*x^2 + 4*x + 1)*exp(x). - G. C. Greubel, Jul 19 2017

A008527 Coordination sequence for body-centered tetragonal lattice.

Original entry on oeis.org

1, 10, 34, 74, 130, 202, 290, 394, 514, 650, 802, 970, 1154, 1354, 1570, 1802, 2050, 2314, 2594, 2890, 3202, 3530, 3874, 4234, 4610, 5002, 5410, 5834, 6274, 6730, 7202, 7690, 8194, 8714, 9250, 9802, 10370, 10954, 11554, 12170, 12802, 13450, 14114, 14794, 15490, 16202, 16930, 17674
Offset: 0

Views

Author

Keywords

Comments

Also sequence found by reading the segment (1, 10) together with the line from 10, in the direction 10, 34, ..., in the square spiral whose vertices are the generalized hexagonal numbers A000217. - Omar E. Pol, Nov 02 2012

Crossrefs

Apart from leading term, same as A108100.
Cf. A206399.
Cf. A016754 (SE), A054554 (NE), A054569 (SW), A053755 (NW), A033951 (S), A054552 (E), A054556 (N), A054567 (W) (Ulam spiral spokes).
A143839 (SSE) + A143855 (ESE) = A143838 (SSW) + A143856 (ENE) = A143854 (WSW) + A143861 (NNE) = A143859 (WNW) + A143860 (NNW) = even bisection = a(2n) = A010021(n).

Programs

  • GAP
    Concatenation([1], List([1..40], n-> 2*(1+4*n^2) )); # G. C. Greubel, Nov 09 2019
  • Magma
    [1] cat [2*(1 + 4*n^2): n in [1..50]]; // G. C. Greubel, Nov 09 2019
    
  • Maple
    1, seq(8*k^2+2, k=1..50);
  • Mathematica
    a[0]:= 1; a[n_]:= 8n^2 +2; Table[a[n], {n,0,50}] (* Alonso del Arte, Sep 06 2011 *)
    LinearRecurrence[{3,-3,1},{1,10,34,74},50] (* Harvey P. Dale, Feb 13 2022 *)
  • PARI
    vector(51, n, if(n==1,1, 2*(1+(2*n-2)^2)) ) \\ G. C. Greubel, Nov 09 2019
    
  • Sage
    [1]+[2*(1+4*n^2) for n in (1..40)] # G. C. Greubel, Nov 09 2019
    

Formula

a(0) = 1; a(n) = 8*n^2+2 for n>0.
From Gary W. Adamson, Dec 27 2007: (Start)
a(n) = (2n+1)^2 + (2n-1)^2 for n>0.
Binomial transform of [1, 9, 15, 1, -1, 1, -1, 1, ...]. (End)
From Colin Barker, Apr 13 2012: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>3.
G.f.: (1+x)*(1+6*x+x^2)/(1-x)^3. (End)
From Bruce J. Nicholson, Jul 31 2019: (Start) Assume n>0.
a(n) = A016754(n) + A016754(n-1).
a(n) = 2 * A053755(n).
a(n) = A054554(n+1) + A054569(n+1).
a(n) = A033951(n) + A054552(n).
a(n) = A054556(n+1) + A054567(n+1). (End)
E.g.f.: -1 + 2*exp(x)*(1 + 2*x)^2. - Stefano Spezia, Aug 02 2019
Sum_{n>=0} 1/a(n) = 3/4+1/8*Pi*coth(Pi/2) = 1.178172.... - R. J. Mathar, May 07 2024

A080335 Diagonal in square spiral or maze arrangement of natural numbers.

Original entry on oeis.org

1, 5, 9, 17, 25, 37, 49, 65, 81, 101, 121, 145, 169, 197, 225, 257, 289, 325, 361, 401, 441, 485, 529, 577, 625, 677, 729, 785, 841, 901, 961, 1025, 1089, 1157, 1225, 1297, 1369, 1445, 1521, 1601, 1681, 1765, 1849, 1937, 2025, 2117, 2209, 2305, 2401, 2501
Offset: 0

Views

Author

Paul Barry, Mar 19 2003

Keywords

Comments

Interleaves the odd squares A016754 with (1+4n^2), A053755.
Squares of positive integers (plus 1 if n is odd). - Wesley Ivan Hurt, Oct 10 2013
a(n) is the maximum total number of queens that can coexist without attacking each other on an [n+3] X [n+3] chessboard, when the lone queen is in the most vulnerable position on the board. Specifically, the lone queen will placed in any center position, facing an opponent's "army" of size a(n)-1 == A137932(n+2). - Bob Selcoe, Feb 12 2015
a(n) is also the edge chromatic number of the complement of the (n+2) X (n+2) rook graph. - Eric W. Weisstein, Jan 31 2024

Crossrefs

Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

Formula

a(n) = (3 + 4*n + 2*n^2 - (-1)^n)/2.
a(2*n) = A016754(n), a(2*n+1) = A053755(n+1).
E.g.f.: exp(x)*(2 + 3*x + x^2) - cosh(x). The sequence 1,1,5,9,... is given by n^2+(1+(-1)^n)/2 with e.g.f. exp(1+x+x^2)*exp(x)-sinh(x). - Paul Barry, Sep 02 2003 and Sep 19 2003
a(0)=1, a(1)=5, a(2)=9, a(3)=17, a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4). - Harvey P. Dale, Jan 29 2012
a(n)+(-1)^n = A137928(n+1). - Philippe Deléham, Feb 17 2012
G.f.: (1 + 3*x - x^2 + x^3)/((1-x)^3*(1+x)). - Colin Barker, Mar 18 2012
a(n) = A000035(n) + A000290(n+1). - Wesley Ivan Hurt, Oct 10 2013
From Bob Selcoe, Feb 12 2015: (Start)
a(n) = A137932(n+2) + 1.
a(n) = (n+1)^2 when n is even; a(n) = (n+1)^2 + 1 when n is odd.
a(n) = A002378(n+2) - A047238(n+3) + 1.
(End)
Sum_{n>=0} 1/a(n) = Pi*coth(Pi/2)/4 + Pi^2/8 - 1/2. - Amiram Eldar, Jul 07 2022

A137932 Terms in an n X n spiral that do not lie on its principal diagonals.

Original entry on oeis.org

0, 0, 0, 4, 8, 16, 24, 36, 48, 64, 80, 100, 120, 144, 168, 196, 224, 256, 288, 324, 360, 400, 440, 484, 528, 576, 624, 676, 728, 784, 840, 900, 960, 1024, 1088, 1156, 1224, 1296, 1368, 1444, 1520, 1600, 1680, 1764, 1848, 1936, 2024, 2116, 2208, 2304, 2400, 2500, 2600, 2704, 2808
Offset: 0

Views

Author

William A. Tedeschi, Feb 29 2008

Keywords

Comments

The count of terms not on the principal diagonals is always even.
The last digit is the repeating pattern 0,0,0,4,8,6,4,6,8,4, which is palindromic if the leading 0's are removed, 4864684.
The sum of the last digits is 40, which is the count of the pattern times 4.
A 4 X 4 spiral is the only spiral, aside from a 0 X 0, whose count of terms that do not lie on its principal diagonals equal the count of terms that do [A137932(4) = A042948(4)] making the 4 X 4 the "perfect spiral".
Yet another property is mod(a(n), A042948(n)) = 0 iff n is even. This is a large family that includes the 4 X 4 spiral.
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+1] X [n+1] chessboard, when the lone queen is in the most vulnerable position on the board, i.e., on a center square. - Bob Selcoe, Feb 12 2015
Also the circumference of the (n-1) X (n-1) grid graph for n > 2. - Eric W. Weisstein, Mar 25 2018
Also the crossing number of the complete bipartite graph K_{5,n}. - Eric W. Weisstein, Sep 11 2018

Examples

			a(0) = 0^2 - (2(0) - mod(0,2)) = 0.
a(3) = 3^2 - (2(3) - mod(3,2)) = 4.
		

Crossrefs

Cf. A042948.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

Formula

a(n) = n^2 - (2*n - mod(n,2)) = n^2 - A042948(n).
a(n) = 2*A007590(n-1). - Enrique Pérez Herrero, Jul 04 2012
G.f.: -4*x^3 / ( (1+x)*(x-1)^3 ). a(n) = 4*A002620(n-1). - R. J. Mathar, Jul 06 2012
From Bob Selcoe, Feb 12 2015: (Start)
a(n) = (n-1)^2 when n is odd; a(n) = (n-1)^2 - 1 when n is even.
a(n) = A002378(n) - A047238(n+1). (End)
From Amiram Eldar, Mar 20 2022: (Start)
Sum_{n>=3} 1/a(n) = Pi^2/24 + 1/4.
Sum_{n>=3} (-1)^(n+1)/a(n) = Pi^2/24 - 1/4. (End)
E.g.f.: x*(x - 1)*cosh(x) + (x^2 - x + 1)*sinh(x). - Stefano Spezia, Oct 17 2022

A156859 The main column of a version of the square spiral.

Original entry on oeis.org

0, 3, 7, 14, 22, 33, 45, 60, 76, 95, 115, 138, 162, 189, 217, 248, 280, 315, 351, 390, 430, 473, 517, 564, 612, 663, 715, 770, 826, 885, 945, 1008, 1072, 1139, 1207, 1278, 1350, 1425, 1501, 1580, 1660, 1743, 1827, 1914, 2002, 2093, 2185, 2280, 2376, 2475, 2575
Offset: 0

Views

Author

Emilio Apricena (emilioapricena(AT)yahoo.it), Feb 17 2009

Keywords

Comments

This spiral is sometimes called an Ulam spiral, but square spiral is a better name. - N. J. A. Sloane, Jul 27 2018
It is easy to see that the only two primes in the sequence are 3, 7. Therefore the primes of the version of Ulam spiral are divided into four parts (see also A035608): northeast (NE), northwest (NW), southwest (SW), and southeast (SE).
Number of pairs (x,y) having x and y of opposite parity with x in {0,...,n} and y in {0,...,2n}. - Clark Kimberling, Jul 02 2012
Partial Sums of A014601(n). - Wesley Ivan Hurt, Oct 11 2013

Crossrefs

Cf. A000290, A000384, A004526, A014601 (first differences), A115258.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

Formula

a(n) = n^2 + n + floor((n+1)/2) = A002378(n) + A004526(n+1) = A002620(n+1) + 3*A002620(n).
From R. J. Mathar, Feb 20 2009: (Start)
G.f.: x*(3+x)/((1+x)*(1-x)^3).
a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4). (End)
a(n-1) = floor(n/(e^(1/n)-1)). - Richard R. Forberg, Jun 19 2013
a(n) = A000290(n+1) + A004526(-n-1). - Wesley Ivan Hurt, Jul 15 2013
a(n) + a(n+1) = A014105(n+1). - R. J. Mathar, Jul 15 2013
a(n) = floor(A000384(n+1)/2). - Bruno Berselli, Nov 11 2013
E.g.f.: (x*(5 + 2*x)*cosh(x) + (1 + 5*x + 2*x^2)*sinh(x))/2. - Stefano Spezia, Apr 24 2024
Sum_{n>=1} 1/a(n) = 4/9 + 2*log(2) - Pi/3. - Amiram Eldar, Apr 26 2024

Extensions

More terms added by Wesley Ivan Hurt, Oct 11 2013

A317186 One of many square spiral sequences: a(n) = n^2 + n - floor((n-1)/2).

Original entry on oeis.org

1, 2, 6, 11, 19, 28, 40, 53, 69, 86, 106, 127, 151, 176, 204, 233, 265, 298, 334, 371, 411, 452, 496, 541, 589, 638, 690, 743, 799, 856, 916, 977, 1041, 1106, 1174, 1243, 1315, 1388, 1464, 1541, 1621, 1702, 1786, 1871, 1959, 2048, 2140, 2233, 2329, 2426
Offset: 0

Views

Author

N. J. A. Sloane, Jul 27 2018

Keywords

Comments

Draw a square spiral on a piece of graph paper, and label the cells starting at the center with the positive (resp. nonnegative) numbers. This produces two versions of the labeled square spiral, shown in the Example section below.
The spiral may proceed clockwise or counterclockwise, and the first arm of the spiral may be along any of the four axes, so there are eight versions of each spiral. However, this has no effect on the resulting sequences, and it is enough to consider just two versions of the square spiral (starting at 1 or starting at 0).
The present sequence is obtained by reading alternate entries on the X-axis (say) of the square spiral started at 1.
The cross-references section lists many sequences that can be read directly off the two spirals. Many other sequences can be obtained from them by using them to extract subsequences from other important sequences. For example, the subsequence of primes indexed by the present sequence gives A317187.
a(n) is also the number of free polyominoes with n + 4 cells whose difference between length and width is n. In this comment the length is the longer of the two dimensions and the width is the shorter of the two dimensions (see the examples of polyominoes). Hence this is also the diagonal 4 of A379625. - Omar E. Pol, Jan 24 2025
From John Mason, Feb 19 2025: (Start)
The sequence enumerates polyominoes of width 2 having precisely 2 horizontal bars. By classifying such polyominoes according to the following templates, it is possible to define a formula that reduces to the one below:
.
OO O O
O OO OO
O O O
O O OO
OO OO O
.
(End)

Examples

			The square spiral when started with 1 begins:
.
  100--99--98--97--96--95--94--93--92--91
                                        |
   65--64--63--62--61--60--59--58--57  90
    |                               |   |
   66  37--36--35--34--33--32--31  56  89
    |   |                       |   |   |
   67  38  17--16--15--14--13  30  55  88
    |   |   |               |   |   |   |
   68  39  18   5---4---3  12  29  54  87
    |   |   |   |       |   |   |   |   |
   69  40  19   6   1---2  11  28  53  86
    |   |   |   |           |   |   |   |
   70  41  20   7---8---9--10  27  52  85
    |   |   |                   |   |   |
   71  42  21--22--23--24--25--26  51  84
    |   |                           |   |
   72  43--44--45--46--47--48--49--50  83
    |                                   |
   73--74--75--76--77--78--79--80--81--82
.
For the square spiral when started with 0, subtract 1 from each entry. In the following diagram this spiral has been reflected and rotated, but of course that makes no difference to the sequences:
.
   99  64--65--66--67--68--69--70--71--72
    |   |                               |
   98  63  36--37--38--39--40--41--42  73
    |   |   |                       |   |
   97  62  35  16--17--18--19--20  43  74
    |   |   |   |               |   |   |
   96  61  34  15   4---5---6  21  44  75
    |   |   |   |   |       |   |   |   |
   95  60  33  14   3   0   7  22  45  76
    |   |   |   |   |   |   |   |   |   |
   94  59  32  13   2---1   8  23  46  77
    |   |   |   |           |   |   |   |
   93  58  31  12--11--10---9  24  47  78
    |   |   |                   |   |   |
   92  57  30--29--28--27--26--25  48  79
    |   |                           |   |
   91  56--55--54--53--52--51--50--49  80
    |                                   |
   90--89--88--87--86--85--84--83--82--81
.
From _Omar E. Pol_, Jan 24 2025: (Start)
For n = 0 there is only one free polyomino with 0 + 4 = 4 cells whose difference between length and width is 0 as shown below, so a(0) = 1.
   _ _
  |_|_|
  |_|_|
.
For n = 1 there are two free polyominoes with 1 + 4 = 5 cells whose difference between length and width is 1 as shown below, so a(1) = 2.
   _ _     _ _
  |_|_|   |_|_|
  |_|_|   |_|_
  |_|     |_|_|
.
(End)
		

Crossrefs

Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.
Filling in these two squares spirals with greedy algorithm: A274640, A274641.
Cf. also A317187.

Programs

  • Mathematica
    a[n_] := n^2 + n - Floor[(n - 1)/2]; Array[a, 50, 0] (* Robert G. Wilson v, Aug 01 2018 *)
    LinearRecurrence[{2, 0, -2 , 1},{1, 2, 6, 11},50] (* or *)
    CoefficientList[Series[(- x^3 - 2 * x^2 - 1) / ((x - 1)^3 * (x + 1)), {x, 0, 50}], x] (* Stefano Spezia, Sep 02 2018 *)

Formula

From Daniel Forgues, Aug 01 2018: (Start)
a(n) = (1/4) * (4 * n^2 + 2 * n + (-1)^n + 3), n >= 0.
a(0) = 1; a(n) = - a(n-1) + 2 * n^2 - n + 2, n >= 1.
a(0) = 1; a(1) = 2; a(2) = 6; a(3) = 11; a(n) = 2 * a(n-1) - 2 * a(n-3) + a(n-4), n >= 4.
G.f.: (- x^3 - 2 * x^2 - 1) / ((x - 1)^3 * (x + 1)). (End)
E.g.f.: ((2 + 3*x + 2*x^2)*cosh(x) + (1 + 3*x + 2*x^2)*sinh(x))/2. - Stefano Spezia, Apr 24 2024
a(n)+a(n+1)=A033816(n). - R. J. Mathar, Mar 21 2025
a(n)-a(n-1) = A042948(n), n>=1. - R. J. Mathar, Mar 21 2025

A056106 Second spoke of a hexagonal spiral.

Original entry on oeis.org

1, 3, 11, 25, 45, 71, 103, 141, 185, 235, 291, 353, 421, 495, 575, 661, 753, 851, 955, 1065, 1181, 1303, 1431, 1565, 1705, 1851, 2003, 2161, 2325, 2495, 2671, 2853, 3041, 3235, 3435, 3641, 3853, 4071, 4295, 4525, 4761, 5003, 5251, 5505, 5765, 6031, 6303
Offset: 0

Views

Author

Henry Bottomley, Jun 09 2000

Keywords

Comments

First differences of A027444. - J. M. Bergot, Jun 04 2012
Numbers of the form ((h^2+h+1)^2+(-h^2+h+1)^2+(h^2+h-1)^2)/(h^2-h+1) for h=n-1. - Bruno Berselli, Mar 13 2013
For n > 0: 2*a(n) = A058331(n) + A001105(n) + A001844(n-1) = A251599(3*n-2) + A251599(3*n-1) + A251599(3*n). - Reinhard Zumkeller, Dec 13 2014
For all n >= 6, a(n+1) expressed in base n is "353". - Mathew Englander, Jan 06 2021

Crossrefs

First differences of A053698, A027444, and A188947.
Cf. A113524 (semiprime terms), A002061.
Other spirals: A054552.

Programs

  • Haskell
    a056106 n = n * (3 * n - 1) + 1  -- Reinhard Zumkeller, Dec 13 2014
  • Magma
    I:=[1,3]; [n le 2 select I[n] else 2*Self(n-1)-Self(n-2)+6: n in [1..50]]; // Vincenzo Librandi, Nov 14 2011
    
  • Mathematica
    Table[3*n^2 - n + 1, {n,0,50}] (* G. C. Greubel, Jul 19 2017 *)
  • PARI
    a(n) = 3*n^2-n+1;
    

Formula

a(n) = 3*n^2 - n + 1.
a(n) = a(n-1) + 6*n - 4 = 2*a(n-1) - a(n-2) + 6.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
E.g.f.: (1+2*x+3*x^2)*exp(x). - Paul Barry, Mar 13 2003
a(n) = A096777(3*n) for n>0. - Reinhard Zumkeller, Dec 29 2007
G.f.: (1+5*x^2)/(1-3*x+3*x^2-x^3). - Colin Barker, Jan 04 2012
a(n) = n*A002061(n+1) - (n-1)*A002061(n). - Bruno Berselli, Jan 15 2013
a(-n) = A056108(n). - Bruno Berselli, Mar 13 2013
Previous Showing 21-30 of 56 results. Next