A042948 -id:A042948 - cp's OEIS frontend

A035608 Expansion of g.f. x(1 + 3x)/((1 + x)*(1 - x)^3).

0, 1, 5, 10, 18, 27, 39, 52, 68, 85, 105, 126, 150, 175, 203, 232, 264, 297, 333, 370, 410, 451, 495, 540, 588, 637, 689, 742, 798, 855, 915, 976, 1040, 1105, 1173, 1242, 1314, 1387, 1463, 1540, 1620, 1701, 1785, 1870, 1958, 2047, 2139, 2232, 2328, 2425, 2525, 2626

Offset: 0

N. J. A. Sloane

Maximum value of Voronoi's principal quadratic form of the first type when variables restricted to {-1,0,1}. - Michael Somos, Mar 10 2004
This is the main row of a version of the "square spiral" when read alternatively from left to right (see link). See also A001107, A007742, A033954, A033991. It is easy to see that the only prime in the sequence is 5. - Emilio Apricena (emilioapricena(AT)yahoo.it), Feb 08 2009
From Mitch Phillipson, Manda Riehl, Tristan Williams, Mar 06 2009: (Start)
a(n) gives the number of elements of S_2 \wr C_k that avoid the pattern 12, using the following ordering:
In S_j, a permutation p avoids a pattern q if it has no subsequence that is order-isomorphic to q. For example, p avoids the pattern 132 if it has no subsequence abc with a < c < b. We extend this notion to S_j \wr C_n as follows. Element psi =[ alpha_1^beta_1, ... alpha_j^beta_j ] avoids tau = [ a_1 ... a_m ] (tau in S_m) if psi' = [ alpha_1*beta_1 ... alpha_j*beta_j ] avoids tau in the usual sense. For n=2, there are 5 elements of S_2 \wr C_2 that avoid the pattern 12. They are: [ 2^1,1^1 ], [ 2^2,1^1 ], [ 2^2,1^2 ], [ 2^1,1^2 ], [ 1^2,2^1 ].
For example, if psi = [2^1,1^2], then psi'=[2,2] which avoids tau=[1,2] because no subsequence ab of psi' has a < b. (End)

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 115.

William A. Tedeschi, Table of n, a(n) for n = 0..10000
Emilio Apricena, A version of the Ulam spiral (This is called the square spiral. - _N. J. A. Sloane_, Jul 27 2018)
Emanuele Munarini, Topological indices for the antiregular graphs, Le Mathematiche, Vol. 76, No. 1 (2021), see p. 283.
Index entries for linear recurrences with constant coefficients, signature (2,0,-2,1).

Partial sums of A042948.
Cf. A002265, A004526, A006578, A011848, A133983, A139592, A173562, A253146.
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.

[n^2 + n - 1 - Floor((n-1)/2): n in [0..25]]; // G. C. Greubel, Oct 29 2017

A035608:=n->floor((n + 1/4)^2): seq(A035608(n), n=0..100); # Wesley Ivan Hurt, Oct 29 2017

Table[n^2 + Floor[n/2], {n, 0, 100}] (* Vladimir Joseph Stephan Orlovsky, Apr 12 2011 *)
CoefficientList[Series[x (1 + 3 x)/((1 + x) (1 - x)^3), {x, 0, 60}], x] (* or *) LinearRecurrence[{2, 0, -2, 1}, {0, 1, 5, 10}, 60] (* Harvey P. Dale, Feb 21 2013 *)

PARI
```
a(n)=n^2+n-1-(n-1)\2
```

a(n) = n^2 + n - 1 - floor((n-1)/2).
a(n) = A011848(2*n+1).
a(n) = A002378(n) - A004526(n+1). - Reinhard Zumkeller, Jan 27 2010
a(n) = 2*A006578(n) - A002378(n)/2 = A139592(n)/2. - Reinhard Zumkeller, Feb 07 2010
a(n) = A002265(n+2) + A173562(n). - Reinhard Zumkeller, Feb 21 2010
a(n) = floor((n + 1/4)^2). - Reinhard Zumkeller, Jan 27 2010
a(n) = (-1)^n*Sum_{i=0..n} (-1)^i*(2*i^2 + 3*i + 1). Omits the leading 0. - William A. Tedeschi, Aug 25 2010
a(n) = n^2 + floor(n/2), from Mathematica section. - Vladimir Joseph Stephan Orlovsky, Apr 12 2011
a(0)=0, a(1)=1, a(2)=5, a(3)=10; for n > 3, a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4). - Harvey P. Dale, Feb 21 2013
For n > 1: a(n) = a(n-2) + 4*n - 3; see also row sums of triangle A253146. - Reinhard Zumkeller, Dec 27 2014
a(n) = 3*A002620(n) + A002620(n+1). - R. J. Mathar, Jul 18 2015
From Amiram Eldar, Mar 20 2022: (Start)
Sum_{n>=1} 1/a(n) = 4 - 2*log(2) - Pi/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*Pi/3 - 4*(1-log(2)). (End)
E.g.f.: (x*(2*x + 3)*cosh(x) + (2*x^2 + 3*x - 1)*sinh(x))/2. - Stefano Spezia, Apr 24 2024

A079896 Discriminants of indefinite binary quadratic forms.

Original entry on oeis.org

5, 8, 12, 13, 17, 20, 21, 24, 28, 29, 32, 33, 37, 40, 41, 44, 45, 48, 52, 53, 56, 57, 60, 61, 65, 68, 69, 72, 73, 76, 77, 80, 84, 85, 88, 89, 92, 93, 96, 97, 101, 104, 105, 108, 109, 112, 113, 116, 117, 120, 124, 125, 128, 129, 132, 133, 136, 137, 140, 141, 145, 148

Offset: 1

Wolfdieter Lang, Jan 31 2003

Numbers n such that n == 0 (mod 4) or n == 1 (mod 4), but n is not a square.
For an indefinite binary quadratic form over the integers a*x^2 + b*x*y + c*y^2 the discriminant is D = b^2 - 4*a*c > 0; and D not a square is assumed.
Also, a superset of A227453. - Ralf Stephan, Sep 22 2013
For the period length of the continued fraction of sqrt(a(n)) see A267857(n). - Wolfdieter Lang, Feb 18 2016
[I changed the offset to 1, since this is an important list. Many parts of the entry, including the b-file, will need to be changed. - N. J. A. Sloane, Mar 14 2023]

McMullen, Curtis. "Billiards and Teichmüller curves." Bulletin of the American Mathematical Society, 60:2 (2023), 195-250. See Table C.1.
A. Scholz and B. Schoeneberg, Einführung in die Zahlentheorie, 5. Aufl., de Gruyter, Berlin, New York, 1973, p. 112.

Robin Visser, Table of n, a(n) for n = 1..10000 (terms n = 1..2001 from Vincenzo Librandi).
S. R. Finch, Class number theory
Steven R. Finch, Class number theory [Cached copy, with permission of the author]

Cf. A014601, A042948 (with squares), A087048 (class numbers), A267857.

Select[ Range[148], (Mod[ #, 4] == 0 || Mod[ #, 4] == 1) && !IntegerQ[ Sqrt[ # ]] & ]

seq(N) = {
  my(n = 1, v = vector(N), top = 0);
  while (top < N,
    if (n%4 < 2 && !issquare(n), v[top++] = n); n++;);
  return(v);
};
seq(62) \\ Gheorghe Coserea, Nov 07 2016

a(2*k^2 + 2*k + 1) = 4*(k+1)^2 + 1 for k >= 0. - Gheorghe Coserea, Nov 07 2016

a(2*k^2 + 4*k + 2 + (k+1)*(-1)^k) = (2*k + 3)*(2*k + 3 + (-1)^k) for k >= 0. - Bruno Berselli, Nov 10 2016

More terms from Robert G. Wilson v, Mar 26 2003

Offset changed to 1 (since this is a list). - N. J. A. Sloane, Mar 14 2023

A042964 Numbers that are congruent to 2 or 3 mod 4.

Original entry on oeis.org

2, 3, 6, 7, 10, 11, 14, 15, 18, 19, 22, 23, 26, 27, 30, 31, 34, 35, 38, 39, 42, 43, 46, 47, 50, 51, 54, 55, 58, 59, 62, 63, 66, 67, 70, 71, 74, 75, 78, 79, 82, 83, 86, 87, 90, 91, 94, 95, 98, 99, 102, 103, 106, 107, 110, 111, 114, 115, 118, 119, 122, 123, 126, 127

Offset: 1

N. J. A. Sloane

Also numbers m such that binomial(m+2, m) mod 2 = 0. - Hieronymus Fischer, Oct 20 2007
Also numbers m such that floor(1+(m/2)) mod 2 = 0. - Hieronymus Fischer, Oct 20 2007
Partial sums of the sequence 2, 1, 3, 1, 3, 1, 3, 1, 3, 1, ... which has period 2. - Hieronymus Fischer, Oct 20 2007
In groups of four add and divide by two the odd and even numbers. - George E. Antoniou, Dec 12 2001
From Jeremy Gardiner, Jan 22 2006: (Start)
Comments on the "mystery calculator". There are 6 cards.
Card 0: 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, 35, 37, 39, ... (A005408 sequence).
Card 1: 2, 3, 6, 7, 10, 11, 14, 15, 18, 19, 22, 23, 26, 27, 30, 31, 34, 35, 38, 39, ... (this sequence).
Card 2: 4, 5, 6, 7, 12, 13, 14, 15, 20, 21, 22, 23, 28, 29, 30, 31, 36, 37, 38, 39, ... (A047566).
Card 3: 8, 9, 10, 11, 12, 13, 14, 15, 24, 25, 26, 27, 28, 29, 30, 31, 40, 41, 42, ... (A115419).
Card 4: 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 48, 49, 50, ... (A115420).
Card 5: 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, ... (A115421).
The trick: You secretly select a number between 1 and 63 from one of the cards. You indicate to me the cards on which that number appears; I tell you the number you selected!
The solution: I add together the first term from each of the indicated cards. The total equals the selected number. The numbers in each sequence all have a "1" in the same position in their binary expansion. Example: You indicate cards 1, 3 and 5. Your selected number is 2 + 8 + 32 = 42.
Numbers having a 1 in position 1 of their binary expansion. One of the mystery calculator sequences: A005408, A042964, A047566, A115419, A115420, A115421. (End)
Complement of A042948. - Reinhard Zumkeller, Oct 03 2008
Also the 2nd Witt transform of A040000 [Moree]. - R. J. Mathar, Nov 08 2008
In general, sequences of numbers congruent to {a,a+i} mod k will have a closed form of (k-2*i)*(2*n-1+(-1)^n)/4+i*n+a, from offset 0. - Gary Detlefs, Oct 29 2013
Union of A004767 and A016825; Fixed points of A098180. - Wesley Ivan Hurt, Jan 14 2014, Oct 13 2015

David Lovler, Table of n, a(n) for n = 1..10000
Maths Magic, Mystery Calculator.
Pieter Moree, The formal series Witt transform, Discr. Math. no. 295 vol. 1-3 (2005) 143-160.
Index entries for linear recurrences with constant coefficients, signature (1,1,-1).

Cf. A000040, A133620, A133621, A133622, A133630, A133635.
Cf. A133872, A133882, A133890, A133900, A133910.
Card trick: A005408, A047566, A115419, A115420, A115421.
Cf. A004767, A016825 A040000, A042948, A047406, A098180.

[2*n+((-1)^(n-1)-1)/2 : n in [1..100]]; // Wesley Ivan Hurt, Oct 13 2015

[n: n in [1..150] | n mod 4 in [2, 3]]; // Vincenzo Librandi, Oct 13 2015

A042964:=n->2*n+((-1)^(n-1)-1)/2; seq(A042964(n), n=1..100); # Wesley Ivan Hurt, Jan 07 2014

Flatten[Table[4n + {2, 3}, {n, 0, 31}]] (* Alonso del Arte, Feb 07 2013 *)
Select[Range[200],MemberQ[{2,3},Mod[#,4]]&] (* or *) LinearRecurrence[ {1,1,-1},{2,3,6},90] (* Harvey P. Dale, Nov 28 2018 *)

PARI
```
a(n)=2*n+2-n%2
```

Vec((2+x+x^2)/((1-x)*(1-x^2)) + O(x^100)) \\ Altug Alkan, Oct 13 2015

a(n) = A047406(n)/2.
From Michael Somos, Jan 12 2000: (Start)
G.f.: x*(2+x+x^2)/((1-x)*(1-x^2)).
a(n) = a(n-1) + 2 + (-1)^n. (End)
a(n) = 2n if n is odd, otherwise n = 2n - 1. - Amarnath Murthy, Oct 16 2003
a(n) = (3 + (-1)^(n-1))/2 + 2*(n-1) = 2n + 2 - (n mod 2). - Hieronymus Fischer, Oct 20 2007
A133872(a(n)) = 0. - Reinhard Zumkeller, Oct 03 2008
a(n) = 4*n - a(n-1) - 3 (with a(1) = 2). - Vincenzo Librandi, Nov 17 2010
a(n) = 2*n + ((-1)^(n-1) - 1)/2. - Gary Detlefs, Oct 29 2013
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/8 - log(2)/4. - Amiram Eldar, Dec 05 2021
E.g.f.: 1 + ((4*x - 1)*exp(x) - exp(-x))/2. - David Lovler, Aug 08 2022

Edited by N. J. A. Sloane, Jun 30 2008 at the suggestion of R. J. Mathar

Corrected by Jaroslav Krizek, Dec 18 2009

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

William A. Tedeschi, Feb 29 2008

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

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

Enrique Pérez Herrero, Table of n, a(n) for n = 0..5000
Kival Ngaokrajang, Illustration of initial terms.
Eric Weisstein's World of Mathematics, Complete Bipartite Graph.
Eric Weisstein's World of Mathematics, Graph Circumference.
Eric Weisstein's World of Mathematics, Graph Crossing Number.
Eric Weisstein's World of Mathematics, Grid Graph.
Index entries for linear recurrences with constant coefficients, signature (2,0,-2,1).

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.

A137932:=n->2*floor((n-1)^2/2); seq(A137932(n), n=0..50); # Wesley Ivan Hurt, Apr 19 2014

Table[2 Floor[(n - 1)^2/2], {n, 0, 20}] (* Enrique Pérez Herrero, Jul 04 2012 *)
2 Floor[(Range[20] - 1)^2/2] (* Eric W. Weisstein, Sep 11 2018 *)
Table[n^2 - 2 n + (1 - (-1)^n)/2, {n, 20}] (* Eric W. Weisstein, Sep 11 2018 *)
LinearRecurrence[{2, 0, -2, 1}, {0, 0, 4, 8}, 20] (* Eric W. Weisstein, Sep 11 2018 *)
CoefficientList[Series[-((4 x^2)/((-1 + x)^3 (1 + x))), {x, 0, 20}], x] (* Eric W. Weisstein, Sep 11 2018 *)

A137932(n)={ return(n^2 - (2*n-n%2))} ;
print(vector(30,n,A137932(n-1))); /* R. J. Mathar, May 12 2014 */

Python
```
a = lambda n: n**2 - (2*n - (n%2))
```

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

A193356 If n is even then 0, otherwise n.

Original entry on oeis.org

1, 0, 3, 0, 5, 0, 7, 0, 9, 0, 11, 0, 13, 0, 15, 0, 17, 0, 19, 0, 21, 0, 23, 0, 25, 0, 27, 0, 29, 0, 31, 0, 33, 0, 35, 0, 37, 0, 39, 0, 41, 0, 43, 0, 45, 0, 47, 0, 49, 0, 51, 0, 53, 0, 55, 0, 57, 0, 59, 0, 61, 0, 63, 0, 65, 0, 67, 0, 69, 0, 71, 0, 73, 0, 75

Offset: 1

José María Grau Ribas, Jul 24 2011

Multiplicative with a(2^e)=0 if e>0 and a(p^e)=p^e for odd primes p. - R. J. Mathar, Aug 01 2011
A005408 and A000004 interleaved (the usual OEIS policy is not to include sequences like this where alternate terms are zero; this is an exception). - Omar E. Pol, Feb 02 2013
Row sums of A057211. - Omar E. Pol, Mar 05 2014
Column k=2 of triangle A196020. - Omar E. Pol, Aug 07 2015
a(n) is the determinant of the (n+2) X (n+2) circulant matrix with the first row [0,0,1,1,...,1]. This matrix is closely linked with the famous ménage problem (see also comments of Vladimir Shevelev in sequence A000179). Namely it defines the class of permutations p of 1,2,...,n+2 such that p(i)<>i and p(i)<>i+1 for i=1,2,...,n+1, and p(n+2)<>1,n+2. And a(n) is also the difference between the number of even and odd such permutations. - Dmitry Efimov, Feb 02 2016

Franz Lemmermeyer, Reciprocity Laws. From Euler to Eisenstein, Springer, 2000, p. 237, eq. (8.5).

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
C. Kravvaritis, Determinant evaluations for binary circulant matrices, Special Matrices, V2(1) (2014), 187-199.
Index entries for linear recurrences with constant coefficients, signature (0,2,0,-1).

Cf. A000004, A000010, A005408, A027656, A042948, A057211, A196020.

I:=[1,0,3,0]; [n le 4 select I[n] else 2*Self(n-2)-Self(n-4): n in [1..80]]; // Vincenzo Librandi, Feb 24 2014

A193356:=n->(1-(-1)^n)*n/2: seq(A193356(n), n=1..100); # Wesley Ivan Hurt, Aug 07 2015

Mathematica
```
Table[PowerMod[n,n,2*n], {n,200}]
```

a(n)=if(n%2,n) \\ Charles R Greathouse IV, Jul 24 2011

A193356 = [i if i%2 else 0 for i in range(1, 101)] # Jwalin Bhatt, Jun 17 2025

a(n) = n^k mod 2n, for any k>=2, also for k=n. [extended by Wolfdieter Lang, Dec 21 2011]
Dirichlet g.f.: (1-2^(1-s))*zeta(s-1). - R. J. Mathar, Aug 01 2011
G.f.: x*(1+x^2)/(1-x^2)^2. - Philippe Deléham, Feb 13 2012
a(n) = A027656(A042948(n-1)) = (1-(-1)^n)*n/2. - Bruno Berselli, Feb 19 2012
a(n) = n * (n mod 2). - Wesley Ivan Hurt, Jun 29 2013
G.f.: Sum_{n >= 1} A000010(n)*x^n/(1 + x^n). - Mircea Merca, Feb 22 2014
a(n) = 2*a(n-2)-a(n-4), for n>4. - Wesley Ivan Hurt, Aug 07 2015
E.g.f.: x*cosh(x). - Robert Israel, Feb 03 2016
a(n) = Product_{k=1..floor(n/2)}(sin(2*Pi*k/n))^2, for n >= 1 (with the empty product put to 1). Trivial for even n from the factor 0 for k = n/2. For odd n see, e.g., the Lemmermeyer reference, eq. (8.5) on p. 237. - Wolfdieter Lang, Aug 29 2016
a(n) = Sum_{k=1..n} (-1)^((n-k)*k). - Rick L. Shepherd, Sep 18 2020
a(n) = Sum_{k=1..n} (-1)^(1+gcd(k,n)) = Sum_{d | n} (-1)^(d+1)*phi(n/d), where phi(n) = A000010(n). - Peter Bala, Jan 14 2024
Dirichlet g.f.: DirichletLambda(s-1). - Michael Shamos, Jun 13 2025

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

N. J. A. Sloane, Jul 27 2018

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)

			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)

Index entries for linear recurrences with constant coefficients, signature (2,0,-2,1).
Index entries for sequences related to polyominoes.

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.
Cf. A000105, A379625.

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

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

A007417 If k appears, 3k does not.

Original entry on oeis.org

1, 2, 4, 5, 7, 8, 9, 10, 11, 13, 14, 16, 17, 18, 19, 20, 22, 23, 25, 26, 28, 29, 31, 32, 34, 35, 36, 37, 38, 40, 41, 43, 44, 45, 46, 47, 49, 50, 52, 53, 55, 56, 58, 59, 61, 62, 63, 64, 65, 67, 68, 70, 71, 72, 73, 74, 76, 77, 79, 80, 81, 82, 83, 85, 86, 88, 89, 90, 91, 92, 94, 95, 97, 98, 99, 100

Offset: 1

N. J. A. Sloane and Simon Plouffe

The characteristic function of this sequence is given by A014578. - Philippe Deléham, Mar 21 2004
Numbers whose ternary representation ends in even number of zeros. - Philippe Deléham, Mar 25 2004
Numbers for which 3 is not an infinitary divisor. - Vladimir Shevelev, Mar 18 2013
Where odd terms occur in A051064. - Reinhard Zumkeller, May 23 2013

			From _Gary W. Adamson_, Mar 02 2010: (Start)
Given the following multiplication table: top row = "not multiples of 3", left column = powers of 3; we get:
   1   2   4   5   7   8   10   11   13
   3   6  12  15  21  24   30   33   39
   9  18  36  45  63  72   90   99  114
  27  54 108
  81
If rows are labeled (1, 2, 3, ...) then odd-indexed rows are in the set; but evens not. Examples: 9 is in the set since 3 is not, but 27 in row 4 can't be. (End)

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

Iain Fox, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Aviezri S. Fraenkel, The vile, dopey, evil and odious game players, Discrete Math. 312 (2012), no. 1, 42-46.
Simon Plouffe, Email to N. J. A. Sloane, Jun. 1994
David Wakeham and David R. Wood, On multiplicative Sidon sets, INTEGERS, 13 (2013), #A26.
Index entries for 3-automatic sequences.

Complement of A145204. - Reinhard Zumkeller, Oct 04 2008

Cf. A007949, A014578 (characteristic function), A042948, A051064, A052330, A092400, A092401.

import Data.List (delete)
a007417 n = a007417_list !! (n-1)
a007417_list = s [1..] where
   s (x:xs) = x : s (delete (3*x) xs)

Select[ Range[100], (# // IntegerDigits[#, 3]& // Split // Last // Count[#, 0]& // EvenQ)&] (* Jean-François Alcover, Mar 01 2013, after Philippe Deléham *)
Select[Range[100], EvenQ@ IntegerExponent[#, 3] &] (* Michael De Vlieger, Sep 01 2020 *)

is(n) = { my(i = 0); while(n%3==0, n/=3; i++); i%2==0; } \\ Iain Fox, Nov 17 2017

is(n)=valuation(n,3)%2==0; \\ Joerg Arndt, Aug 08 2020

from sympy import integer_log
def A007417(n):
    def f(x): return n+x-sum(((m:=x//9**i)-2)//3+(m-1)//3+2 for i in range(integer_log(x,9)[0]+1))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Feb 15 2025

Limit_{n->infinity} a(n)/n = 4/3. - Philippe Deléham, Mar 21 2004
Partial sums of A092400. Indices of even numbers in A007949. Indices of odd numbers in A051064. a(n) = A092401(2n-1). - Philippe Deléham, Mar 29 2004
{a(n)} = A052330({A042948(n)}), where {a(n)} denotes the set of integers in the sequence. - Peter Munn, Aug 31 2019

More terms from Philippe Deléham, Mar 29 2004

Typo corrected by Philippe Deléham, Apr 15 2010

A367338 Comma-successor to n: second term of commas sequence if initial term is n, or -1 if there is no second term.

Original entry on oeis.org

12, 24, 36, 48, 61, 73, 85, 97, 100, 11, 23, 35, 47, 59, 72, 84, 96, -1, 110, 22, 34, 46, 58, 71, 83, 95, -1, 109, 120, 33, 45, 57, 69, 82, 94, -1, 108, 119, 130, 44, 56, 68, 81, 93, -1, 107, 118, 129, 140, 55, 67, 79, 92, -1, 106, 117, 128, 139, 150, 66, 78, 91, -1, 105, 116

Offset: 1

N. J. A. Sloane, Nov 15 2023

Construct the commas sequence as in A121805, but take the first term to be n. Then a(n), the comma-successor to n, is the second term, or -1 if no second term exists.
More generally, we define a comma-child of n to be any number m with the property that m-n = 10*x+y, where x is the least significant digit of n and y is the most significant digit of m.
A positive number can have 0, 1, or 2 comma-children. In accordance with the Law of Primogeniture, the first-born child (i.e. the smallest), if there is one, is the comma-successor.
Comment from N. J. A. Sloane, Nov 19 2023: (Start)
The following is a proof of a slight modification of a conjecture made by Ivan N. Ianakiev in A367341.
The Comma-Successor Theorem.
Let D(b) denote the set of numbers k which have no comma-successor in base b ("comma-successor" is the base-b generalization of the rule that defines A121805). If a commas sequence reaches a number in D(b) it will end there.
Then D(b) consists precisely of the numbers which when written in base b have the form
  cc...cxy = (b^i-1)*b^2/(b-1) + b*x + y,
   with i >= 0 copies of c = b-1, where x and y are in the range [1..b-2] and satisfy x+y = b-1.  .... (*)
For b = 10 the numbers D(10) are listed in A367341.
For an outline of the proof, see the attached text-file.
Note that in base b = 2, no values of x satisfying (*) exist, and the theorem asserts that D(2) is empty. In fact it is easy to check directly that every commas sequence in base 2 is infinite.  If the initial term is 0 or 1 mod 4 then the sequence will merge with A042948, and if the initial term is 2 or 3 mod 4 then the sequence will merge with A042964.
(End)

			a(1) = A121803(2) = 12,
a(2) = A139284(2) = 24,
a(3) = 36, since the full commas sequence starting with 3 is [3, 36] (which also implies a(36) = -1),
a(4) = A366492(2) = 48, and so on.
60 is the first number that is a comma-child (a member of A367312) but is missing from the present sequence (it is a comma-child but not a comma-successor, since it loses out to 59).

Michael S. Branicky, Table of n, a(n) for n = 1..10000
Eric Angelini, Michael S. Branicky, Giovanni Resta, N. J. A. Sloane, and David W. Wilson, The Comma Sequence: A Simple Sequence With Bizarre Properties, arXiv:2401.14346, Fibonacci Quarterly 62:3 (2024), 215-232.
Eric Angelini, Michael S. Branicky, Giovanni Resta, N. J. A. Sloane, and David W. Wilson, The Comma Sequence: A Simple Sequence With Bizarre Properties, Local copy.
N. J. A. Sloane, Eric Angelini's Comma Sequence, Experimental Math Seminar, Rutgers Univ., January 18, 2024, Youtube video; Slides

A367346 lists those n for which there is more than one choice for the second term.
A367612 lists the numbers that are comma-children of some number k.
Cf. A121805, A139284, A366492, A367337.
See also A042948, A042964, A367339, A367340, A367341.

Ldigit:=proc(n) local v; v:=convert(n, base, 10); v[-1]; end;
A367338 := proc(n) local f,i,d;
f := (n mod 10);
d:=10*f;
for i from 1 to 9 do
d := d+1;
if Ldigit(n+d) = i then return(n+d); fi;
od:
return(-1);
end;
for n from 1 to 50 do lprint(n, A367338(n)); od: # N. J. A. Sloane, Dec 06 2023

a[n_] := a[n] = Module[{l = n, y = 1, d}, While[y < 10, l = l + 10*(Mod[l, 10]); y = 1; While[y < 10, d = IntegerDigits[l + y][[1]]; If[d == y, l = l + y; Break[];]; y++;]; If[y < 10, Return[l]];]; Return[-1];];
Table[a[n], {n, 1, 65}] (* Robert P. P. McKone, Dec 18 2023 *)

from itertools import islice
def a(n):
    an, y = n, 1
    while y < 10:
        an, y = an + 10*(an%10), 1
        while y < 10:
            if str(an+y)[0] == str(y):
                an += y
                break
            y += 1
        if y < 10:
            return an
    return -1
print([a(n) for n in range(1, 66)]) # Michael S. Branicky, Nov 15 2023

A047470 Numbers that are congruent to {0, 3} mod 8.

Original entry on oeis.org

0, 3, 8, 11, 16, 19, 24, 27, 32, 35, 40, 43, 48, 51, 56, 59, 64, 67, 72, 75, 80, 83, 88, 91, 96, 99, 104, 107, 112, 115, 120, 123, 128, 131, 136, 139, 144, 147, 152, 155, 160, 163, 168, 171, 176, 179, 184, 187, 192, 195, 200, 203, 208, 211, 216, 219, 224, 227, 232

Offset: 1

N. J. A. Sloane

Maximum number of squares attacked by a queen on an n X n chessboard. - Stewart Gordon, Mar 23 2001
Equivalently, maximum vertex degree in the n X n queen graph. - Eric W. Weisstein, Jun 20 2017
Number of squares attacked by a queen on a toroidal chessboard. - Diego Torres (torresvillarroel(AT)hotmail.com), May 19 2001
List of squared distances between points of diamond 'lattice' with minimal distance sqrt(3). - Arnold Neumaier (Arnold.Neumaier(AT)univie.ac.at), Aug 01 2003
Draw a figure-eight knot diagram on the plane and assign a list of nonnegative numbers at each crossing as follows. Start with 0 and choose a crossing on the knot. Pick a direction and walk around the knot, appending the following nonnegative number everytime a crossing is visited. Two series of sequences are obtained: this sequence, A047535, A047452, A047617 and A047615, A047461, A047452, A047398 (see example). - Franck Maminirina Ramaharo, Jul 22 2018

			From _Franck Maminirina Ramaharo_, Jul 22 2018: (Start)
Consider the following equivalent figure-eight knot diagrams:
+---------------------+           +-----------------n
|                     |           |                 |
|           +---------B-----+     |           w-----A---e
|           |         |     |     |           |     |   |
|     n-----C---+     |     |     |           |     |   |
|     |     |   |     |     | <=> |   +-------B-----s   |
|     |     +---D-----+     |     |   |       |         |
|     |         |           |     |   |       |         |
w-----A---------e           |     +---C-------D---------+
      |                     |         |       |
      s---------------------+         +-------+
Uppercases A,B,C,D denote crossings, and lowercases n,s,w,e denote directions. Due to symmetry and ambient isotopy, all possible sequences are obtained by starting from crossing A and choose either direction 'n' or 's'.
Direction 'n':
A: 0, 3,  8, 11, 16, 19, 24, 27, 32, 35, 40, ... (this sequence);
B: 4, 7, 12, 15, 20, 23, 28, 31, 36, 39, 44, ... A047535;
C: 1, 6,  9, 14, 17, 22, 25, 30, 33, 38, 41, ... A047452;
D: 2, 5, 10, 13, 18, 21, 26, 29, 34, 37, 42, ... A047617.
Direction 's':
A: 0, 5,  8, 13, 16, 21, 24, 29, 32, 37, 40, ... A047615;
B: 1, 4,  9, 12, 17, 20, 25, 28, 33, 36, 41, ... A047461;
C: 2, 7, 10, 15, 18, 23, 26, 31, 34, 39, 42, ... A047524;
D: 3, 6, 11, 14, 19, 22, 27, 30, 35, 38, 43, ... A047398.
(End)

Muniru A Asiru, Table of n, a(n) for n = 1..5000
Eric Weisstein's World of Mathematics, Maximum Vertex Degree.
Eric Weisstein's World of Mathematics, Queen Graph.
Index entries for linear recurrences with constant coefficients, signature (1,1,-1).

Cf. A042948, A047398, A047461, A047452, A047524, A047535, A047615, A047617.

a:=[0,3,8];; for n in [4..50] do a[n]:=a[n-1]+a[n-2]-a[n-3]; od; a; # Muniru A Asiru, Jul 23 2018

a:=n->add(4+(-1)^j,j=1..n):seq(a(n),n=0..64); # Zerinvary Lajos, Dec 13 2008

With[{c = 8 Range[0, 30]}, Sort[Join[c, c + 3]]] (* Harvey P. Dale, Oct 11 2011 *)
Table[(8 n - 9 - (-1)^n)/2, {n, 20}] (* Eric W. Weisstein, Jun 20 2017 *)
LinearRecurrence[{1, 1, -1}, {0, 3, 8}, 20] (* Eric W. Weisstein, Jun 20 2017 *)
CoefficientList[Series[(x (3 + 5 x))/((-1 + x)^2 (1 + x)), {x, 0, 20}], x]  (* Eric W. Weisstein, Jun 20 2017 *)

forstep(n=0,200,[3,5],print1(n", ")) \\ Charles R Greathouse IV, Oct 17 2011

def A047470(n): return (n-1<<2)-(n&1^1) # Chai Wah Wu, Mar 30 2024

a(n) = a(n-1) + 4 + (-1)^n.
a(n) = a(n-1) + a(n-2) - a(n-3).
a(n) = A042948(n) + A005843(n).
G.f.: (3x+5*x^2)/((1-x)*(1-x^2)).
a(n) = 8*n - a(n-1) - 13 (with a(1)=0). - Vincenzo Librandi, Aug 06 2010
a(n+1) = Sum_{k>=0} A030308(n,k)*A171497(k). - Philippe Deléham, Oct 17 2011
a(n) = 4*n -(9 + (-1)^n)/2. - Arkadiusz Wesolowski, Sep 18 2012
E.g.f: (10 - exp(-x) + (8*x - 9)*exp(x))/2. - Franck Maminirina Ramaharo, Jul 22 2018
Sum_{n>=2} (-1)^n/a(n) = (sqrt(2)-1)*Pi/16 + log(2)/2 - sqrt(2)*log(sqrt(2)+1)/8. - Amiram Eldar, Dec 11 2021

More terms from Vincenzo Librandi, Aug 06 2010

A047467 Numbers that are congruent to {0, 2} mod 8.

Original entry on oeis.org

0, 2, 8, 10, 16, 18, 24, 26, 32, 34, 40, 42, 48, 50, 56, 58, 64, 66, 72, 74, 80, 82, 88, 90, 96, 98, 104, 106, 112, 114, 120, 122, 128, 130, 136, 138, 144, 146, 152, 154, 160, 162, 168, 170, 176, 178, 184, 186, 192, 194, 200, 202, 208, 210, 216, 218, 224, 226, 232

Offset: 1

N. J. A. Sloane

David Lovler, Table of n, a(n) for n = 1..10000
Index entries for linear recurrences with constant coefficients, signature (1,1,-1).

Union of A008590 and A017089.

Cf. A030308, A042948.

{#,#+2}&/@(8*Range[0,30])//Flatten (* or *) LinearRecurrence[{1,1,-1},{0,2,8},60] (* Harvey P. Dale, Nov 30 2019 *)

forstep(n=0,200,[2,6],print1(n", ")) \\ Charles R Greathouse IV, Oct 17 2011

a(n) = 4*n - 5 - (-1)^n; \\ David Lovler, Jul 25 2022

From R. J. Mathar, Sep 19 2008: (Start)
a(n) = 4*n - 5 - (-1)^n = 2*A042948(n-1).
G.f.: 2*x^2*(1+3x)/((1-x)^2*(1+x)). (End)
a(n) = 8*n - a(n-1) - 14 with a(1)=0. - Vincenzo Librandi, Aug 06 2010
a(n+1) = Sum_{k>=0} A030308(n,k)*b(k) with b(0)=2 and b(k)=2^(k+2)for k > 0. - Philippe Deléham, Oct 17 2011
a(n) = floor((8/3)*floor(3*n/2)). - Clark Kimberling, Jul 04 2012
Sum_{n>=2} (-1)^n/a(n) = Pi/16 + 3*log(2)/8. - Amiram Eldar, Dec 18 2021
E.g.f.: 6 + (4*x - 5)*exp(x) - exp(-x). - David Lovler, Jul 22 2022

More terms from Vincenzo Librandi, Aug 06 2010

cp's OEIS Frontend

A035608 Expansion of g.f. x*(1 + 3*x)/((1 + x)*(1 - x)^3).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A079896 Discriminants of indefinite binary quadratic forms.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A042964 Numbers that are congruent to 2 or 3 mod 4.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Magma

Maple

Mathematica

PARI

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A193356 If n is even then 0, otherwise n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

Formula

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

Views

Author

Keywords

Comments

Examples

Links

A035608 Expansion of g.f. x(1 + 3x)/((1 + x)*(1 - x)^3).