A001105 -id:A001105 - cp's OEIS frontend

A046092 4 times triangular numbers: a(n) = 2n(n+1).

0, 4, 12, 24, 40, 60, 84, 112, 144, 180, 220, 264, 312, 364, 420, 480, 544, 612, 684, 760, 840, 924, 1012, 1104, 1200, 1300, 1404, 1512, 1624, 1740, 1860, 1984, 2112, 2244, 2380, 2520, 2664, 2812, 2964, 3120, 3280, 3444, 3612, 3784, 3960, 4140, 4324

Offset: 0

N. J. A. Sloane and Eric W. Weisstein

Consider all Pythagorean triples (X,Y,Z=Y+1) ordered by increasing Z; sequence gives Y values. X values are 1, 3, 5, 7, 9, ... (A005408), Z values are A001844.
In the triple (X, Y, Z) we have X^2=Y+Z. Actually, the triple is given by {x, (x^2 -+ 1)/2}, where x runs over the odd numbers (A005408) and x^2 over the odd squares (A016754). - Lekraj Beedassy, Jun 11 2004
a(n) is the number of edges in n X n square grid with all horizontal and vertical segments filled in. - Asher Auel, Jan 12 2000 [Corrected by Felix Huber, Apr 09 2024]
a(n) is the only number satisfying an inequality related to zeta(2) and zeta(3): Sum_{i>a(n)+1} 1/i^2 < Sum_{i>n} 1/i^3 < Sum_{i>a(n)} 1/i^2. - Benoit Cloitre, Nov 02 2001
Number of right triangles made from vertices of a regular n-gon when n is even. - Sen-Peng Eu, Apr 05 2001
Number of ways to change two non-identical letters in the word aabbccdd..., where there are n type of letters. - Zerinvary Lajos, Feb 15 2005
a(n) is the number of (n-1)-dimensional sides of an (n+1)-dimensional hypercube (e.g., squares have 4 corners, cubes have 12 edges, etc.). - Freek van Walderveen (freek_is(AT)vanwal.nl), Nov 11 2005
From Nikolaos Diamantis (nikos7am(AT)yahoo.com), May 23 2006: (Start)
Consider a triangle, a pentagon, a heptagon, ..., a k-gon where k is odd. We label a triangle with n=1, a pentagon with n=2, ..., a k-gon with n = floor(k/2). Imagine a player standing at each vertex of the k-gon.
Initially there are 2 frisbees, one held by each of two neighboring players. Every time they throw the frisbee to one of their two nearest neighbors with equal probability. Then a(n) gives the average number of steps needed so that the frisbees meet.
I verified this by simulating the processes with a computer program. For example, a(2) = 12 because in a pentagon that's the expected number of trials we need to perform. That is an exercise in Concrete Mathematics and it can be done using generating functions. (End)
A diagonal of A059056. - Zerinvary Lajos, Jun 18 2007
If X_1,...,X_n is a partition of a 2n-set X into 2-blocks then a(n-1) is equal to the number of 2-subsets of X containing none of X_i, (i=1,...,n). - Milan Janjic, Jul 16 2007
X values of solutions to the equation 2*X^3 + X^2 = Y^2. To find Y values: b(n) = 2n(n+1)(2n+1). - Mohamed Bouhamida, Nov 06 2007
Number of (n+1)-permutations of 3 objects u,v,w, with repetition allowed, containing n-1 u's. Example: a(1)=4 because we have vv, vw, wv and ww; a(2)=12 because we can place u in each of the previous four 2-permutations either in front, or in the middle, or at the end. - Zerinvary Lajos, Dec 27 2007
Sequence found by reading the line from 0, in the direction 0, 4, ... and the same line from 0, in the direction 0, 12, ..., in the square spiral whose vertices are the triangular numbers A000217. - Omar E. Pol, May 03 2008
a(n) is also the least weight of self-conjugate partitions having n different even parts. - Augustine O. Munagi, Dec 18 2008
From Peter Luschny, Jul 12 2009: (Start)
The general formula for alternating sums of powers of even integers is in terms of the Swiss-Knife polynomials P(n,x) A153641 (P(n,1)-(-1)^k P(n,2k+1))/2. Here n=2, thus
a(k) = |(P(2,1) - (-1)^k*P(2,2k+1))/2|. (End)
The sum of squares of n+1 consecutive numbers between a(n)-n and a(n) inclusive equals the sum of squares of n consecutive numbers following a(n). For example, for n = 2, a(2) = 12, and the corresponding equation is 10^2 + 11^2 + 12^2 = 13^2 + 14^2. - Tanya Khovanova, Jul 20 2009
Number of roots in the root system of type D_{n+1} (for n>2). - Tom Edgar, Nov 05 2013
Draw n ellipses in the plane (n>0), any 2 meeting in 4 points; sequence gives number of intersections of these ellipses (cf. A051890, A001844); a(n) = A051890(n+1) - 2 = A001844(n) - 1. - Jaroslav Krizek, Dec 27 2013
a(n) appears also as the second member of the quartet [p0(n), a(n), p2(n), p3(n)] of the square of [n, n+1, n+2, n+3] in the Clifford algebra Cl_2 for n >= 0. p0(n) = -A147973(n+3), p2(n) = A054000(n+1) and p3(n) = A139570(n). See a comment on A147973, also with a reference. - Wolfdieter Lang, Oct 15 2014
a(n) appears also as the third and fourth member of the quartet [p0(n), p0(n), a(n), a(n)] of the square of [n, n, n+1, n+1] in the Clifford algebra Cl_2 for n >= 0. p0(n) = A001105(n). - Wolfdieter Lang, Oct 16 2014
Consider two equal rectangles composed of unit squares. Then surround the 1st rectangle with 1-unit-wide layers to build larger rectangles, and surround the 2nd rectangle just to hide the previous layers. If r(n) and h(n) are the number of unit squares needed for n layers in the 1st case and the 2nd case, then for all rectangles, we have a(n) = r(n) - h(n) for n>=1. - Michel Marcus, Sep 28 2015
When greater than 4, a(n) is the perimeter of a Pythagorean triangle with an even short leg 2*n. - Agola Kisira Odero, Apr 26 2016
Also the number of minimum connected dominating sets in the (n+1)-cocktail party graph. - Eric W. Weisstein, Jun 29 2017
a(n+1) is the harmonic mean of A000384(n+2) and A014105(n+1). - Bob Andriesse, Apr 27 2019
Consider a circular cake from which wedges of equal center angle c are cut out in clockwise succession and turned around so that the bottom comes to the top. This goes on until the cake shows its initial surface again. An interesting case occurs if 360°/c is not an integer. Then, with n = floor(360°/c), the number of wedges which have to be cut out and turned equals a(n). (For the number of cutting line segments see A005408.) - According to Peter Winkler's book "Mathematical Mind-Benders", which presents the problem and its solution (see Winkler, pp. 111, 115) the problem seems to be of French origin but little is known about its history. - Manfred Boergens, Apr 05 2022
a(n-3) is the maximum irregularity over all maximal 2-degenerate graphs with n vertices.  The extremal graphs are 2-stars (K_2 joined to n-2 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
Number of ways of placing a domino on a (n+1)X(n+1) board of squares. - R. J. Mathar, Apr 24 2024
The sequence terms are the exponents in the expansion of (1/(1 + x)) * Sum_{n >= 0} x^n * Product_{k = 1..n} (1 - x^(2*k-1))/(1 + x^(2*k+1)) = 1 - x^4 + x^12 - x^24 + x^40 - x^60 + - ... (Andrews and Berndt, Entry 9.3.3, p. 229). Cf. A153140. - Peter Bala, Feb 15 2025
Number of edges in an (n+1)-dimensional orthoplex. 2D orthoplexes (diamonds) have 4 edges, 3D orthoplexes (octahedrons) have 12 edges, 4D orthoplexes (16-cell) have 24 edges, and so on. - Aaron Franke, Mar 23 2025

			a(7)=112 because 112 = 2*7*(7+1).
The first few triples are (1,0,1), (3,4,5), (5,12,13), (7,24,25), ...
The first such partitions, corresponding to a(n)=1,2,3,4, are 2+2, 4+4+2+2, 6+6+4+4+2+2, 8+8+6+6+4+4+2+2. - _Augustine O. Munagi_, Dec 18 2008

George E. Andrews and Bruce C. Berndt, Ramanujan's Lost Notebook, Part I, Springer, 2005.
Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 3.
Albert H. Beiler, Recreations in the Theory of Numbers. New York: Dover, p. 125, 1964.
Ronald L. Graham, D. E. Knuth and Oren Patashnik, Concrete Mathematics, Reading, Massachusetts: Addison-Wesley, 1994.
Peter Winkler, Mathematical Mind-Benders, Wellesley, Massachusetts: A K Peters, 2007.

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Allan Bickle and Zhongyuan Che, Irregularities of Maximal k-degenerate Graphs, Discrete Applied Math. 331 (2023) 70-87.
Allan Bickle, A Survey of Maximal k-degenerate Graphs and k-Trees, Theory and Applications of Graphs 0 1 (2024) Article 5.
H. J. Brothers, Pascal's Prism: Supplementary Material.
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Exact and asymptotic solutions of the recurrence f(n) = f(floor(n/2)) + f(ceiling(n/2)) + g(n): theory and applications, preprint, 2016.
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Exact and Asymptotic Solutions of a Divide-and-Conquer Recurrence Dividing at Half: Theory and Applications, ACM Transactions on Algorithms, Vol. 13, No. 4 (2017), Article #47.
Z. Janelidze, F. van Niekerk, and J. Viljoen, What is the maximal connected partial symmetry index of a connected graph of a given size?, arXiv:2502.00704 [math.CO], 2025. See p. 3.
Milan Janjic, Two Enumerative Functions
Ron Knott, Pythagorean Triples and Online Calculators
Tanya Khovanova, A Miracle Equation.
Augustine O. Munagi, Pairing conjugate partitions by residue classes, Discrete Math., 308 (2008), 2492-2501. [From _Augustine O. Munagi_, Dec 18 2008]
Enrique Navarrete and Daniel Orellana, Finding Prime Numbers as Fixed Points of Sequences, arXiv:1907.10023 [math.NT], 2019.
Omar E. Pol, Determinacion geometrica de los numeros primos y perfectos.
Amelia Carolina Sparavigna, The groupoids of Mersenne, Fermat, Cullen, Woodall and other Numbers and their representations by means of integer sequences, Politecnico di Torino, Italy (2019), [math.NT].
Amelia Carolina Sparavigna, Some Groupoids and their Representations by Means of Integer Sequences, International Journal of Sciences (2019) Vol. 8, No. 10.
Rusliansyah D. Suprijanto, Observation on Sums of Powers of Integers Divisible by Four, Applied Mathematical Sciences, Vol. 8, 2014, no. 45, 2219 - 2226.
Leo Tavares, Illustration: Diamond Rows
Herman Tulleken, Polyominoes 2.2: How they fit together, (2019).
Eric Weisstein's World of Mathematics, Aztec Diamond.
Eric Weisstein's World of Mathematics, Cocktail Party Graph.
Eric Weisstein's World of Mathematics, Connected Dominating Set.
Eric Weisstein's World of Mathematics, Gear Graph.
Eric Weisstein's World of Mathematics, Hamiltonian Path.
Eric Weisstein's World of Mathematics, Pythagorean Triple.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A045943, A028895, A002943, A054000, A000330, A007290, A002378, A033996, A124080, A028896, A049598, A005563, A000217, A033586, A085250.
Main diagonal of array in A001477.
Equals A033996/2. Cf. A001844. - Augustine O. Munagi, Dec 18 2008
Cf. A078371, A141530 (see Librandi's comment in A078371).
Cf. A097080, A001845.
Cf. similar sequences listed in A299645.
Cf. A005408.
Cf. A016754.
Cf. A002378, A046092, A028896 (irregularities of maximal k-degenerate graphs).

a046092 = (* 2) . a002378  -- Reinhard Zumkeller, Dec 15 2013

[2*n*(n+1): n in [0..50]]; // Vincenzo Librandi, Oct 04 2011

Table[2 n (n + 1), {n, 0, 50}] (* Stefan Steinerberger, Apr 03 2006 *)
LinearRecurrence[{3, -3, 1}, {0, 4, 12}, 50] (* Harvey P. Dale, Jul 25 2011 *)
4*Binomial[Range[50], 2] (* Harvey P. Dale, Jul 25 2011 *)

A046092(n):=2*n*(n+1)$
makelist(A046092(n),n,0,30); /* Martin Ettl, Nov 08 2012 */

a(n)=binomial(n+1,2)<<2 \\ Charles R Greathouse IV, Jun 10 2011

def A046092(n): return n*(n+1)<<1 # Chai Wah Wu, Mar 11 2025

a(n) = A100345(n+1, n-1) for n>0.
a(n) = 2*A002378(n) = 4*A000217(n). - Lekraj Beedassy, May 25 2004
a(n) = C(2n, 2) - n = 4*C(n, 2). - Zerinvary Lajos, Feb 15 2005
From Lekraj Beedassy, Jun 04 2006: (Start)
a(n) - a(n-1)=4*n.
Let k=a(n). Then a(n+1) = k + 2*(1 + sqrt(2k + 1)). (End)
Array read by rows: row n gives A033586(n), A085250(n+1). - Omar E. Pol, May 03 2008
O.g.f.:4*x/(1-x)^3; e.g.f.: exp(x)*(2*x^2+4*x). - Geoffrey Critzer, May 17 2009
From Stephen Crowley, Jul 26 2009: (Start)
a(n) = 1/int(-(x*n+x-1)*(step((-1+x*n)/n)-1)*n*step((x*n+x-1)/(n+1)),x=0..1) where step(x)=piecewise(x<0,0,0<=x,1) is the Heaviside step function.
Sum_{n>=1} 1/a(n) = 1/2. (End)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3); a(0)=0, a(1)=4, a(2)=12. - Harvey P. Dale, Jul 25 2011
For n > 0, a(n) = 1/(Integral_{x=0..Pi/2} (sin(x))^(2*n-1)*(cos(x))^3). - Francesco Daddi, Aug 02 2011
a(n) = A001844(n) - 1. - Omar E. Pol, Oct 03 2011
(a(n) - A000217(k))^2  = A000217(2n-k)*A000217(2n+1+k) - (A002378(n) - A000217(k)), for all k. See also A001105. - Charlie Marion, May 09 2013
From Ivan N. Ianakiev, Aug 30 2013: (Start)
a(n)*(2m+1)^2 + a(m) = a(n*(2m+1)+m), for any nonnegative integers n and m.
t(k)*a(n) + t(k-1)*a(n+1) = a((n+1)*(t(k)-t(k-1)-1)), where k>=2, n>=1, t(k)=A000217(k). (End)
a(n) = A245300(n,n). - Reinhard Zumkeller, Jul 17 2014
2*a(n)+1 = A016754(n) = A005408(n)^2, the odd squares. - M. F. Hasler, Oct 02 2014
Sum_{n>=1} (-1)^(n+1)/a(n) = log(2) - 1/2 = A187832. - Ilya Gutkovskiy, Mar 16 2017
a(n) = lcm(2*n,2*n+2). - Enrique Navarrete, Aug 30 2017
a(n)*a(n+k) + k^2 = m^2 (a perfect square), n >= 1, k >= 0. - Ezhilarasu Velayutham, May 13 2019
From Amiram Eldar, Jan 29 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = cosh(Pi/2)/(Pi/2).
Product_{n>=1} (1 - 1/a(n)) = -2*cos(sqrt(3)*Pi/2)/Pi. (End)
a(n) = A016754(n) - A001844(n). - Leo Tavares, Sep 20 2022

A028387 a(n) = n + (n+1)^2.

Original entry on oeis.org

1, 5, 11, 19, 29, 41, 55, 71, 89, 109, 131, 155, 181, 209, 239, 271, 305, 341, 379, 419, 461, 505, 551, 599, 649, 701, 755, 811, 869, 929, 991, 1055, 1121, 1189, 1259, 1331, 1405, 1481, 1559, 1639, 1721, 1805, 1891, 1979, 2069, 2161, 2255, 2351, 2449, 2549, 2651

Offset: 0

Patrick De Geest

a(n+1) is the least k > a(n) + 1 such that A000217(a(n)) + A000217(k) is a square. - David Wasserman, Jun 30 2005
Values of Fibonacci polynomial n^2 - n - 1 for n = 2, 3, 4, 5, ... - Artur Jasinski, Nov 19 2006
A127701 * [1, 2, 3, ...]. - Gary W. Adamson, Jan 24 2007
Row sums of triangle A135223. - Gary W. Adamson, Nov 23 2007
Equals row sums of triangle A143596. - Gary W. Adamson, Aug 26 2008
a(n-1) gives the number of n X k rectangles on an n X n chessboard (for k = 1, 2, 3, ..., n). - Aaron Dunigan AtLee, Feb 13 2009
sqrt(a(0) + sqrt(a(1) + sqrt(a(2) + sqrt(a(3) + ...)))) = sqrt(1 + sqrt(5 + sqrt(11 + sqrt(19 + ...)))) = 2. - Miklos Kristof, Dec 24 2009
When n + 1 is prime, a(n) gives the number of irreducible representations of any nonabelian group of order (n+1)^3. - Andrew Rupinski, Mar 17 2010
a(n) = A176271(n+1, n+1). - Reinhard Zumkeller, Apr 13 2010
The product of any 4 consecutive integers plus 1 is a square (see A062938); the terms of this sequence are the square roots. - Harvey P. Dale, Oct 19 2011
Or numbers not expressed in the form m + floor(sqrt(m)) with integer m. - Vladimir Shevelev, Apr 09 2012
Left edge of the triangle in A214604: a(n) = A214604(n+1,1). - Reinhard Zumkeller, Jul 25 2012
Another expression involving phi = (1 + sqrt(5))/2 is a(n) = (n + phi)(n + 1 - phi). Therefore the numbers in this sequence, even if they are prime in Z, are not prime in Z[phi]. - Alonso del Arte, Aug 03 2013
a(n-1) = n*(n+1) - 1, n>=0, with a(-1) = -1, gives the values for a*c of indefinite binary quadratic forms [a, b, c] of discriminant D = 5 for b = 2*n+1. In general D = b^2 - 4ac > 0 and the form [a, b, c] is a*x^2 + b*x*y + c*y^2. - Wolfdieter Lang, Aug 15 2013
a(n) has prime factors given by A038872. - Richard R. Forberg, Dec 10 2014
A253607(a(n)) = -1. - Reinhard Zumkeller, Jan 05 2015
An example of a quadratic sequence for which the continued square root map (see A257574) produces the number 2. There are infinitely many sequences with this property - another example is A028387. See Popular Computing link. - N. J. A. Sloane, May 03 2015
Left edge of the triangle in A260910: a(n) = A260910(n+2,1). - Reinhard Zumkeller, Aug 04 2015
Numbers m such that 4m+5 is a square. - Bruce J. Nicholson, Jul 19 2017
The numbers represented as 131 in base n: 131_4 = 29, 131_5 = 41, ... . If 'digits' larger than the base are allowed then 131_2 = 11 and 131_1 = 5 also. - Ron Knott, Nov 14 2017
From Klaus Purath, Mar 18 2019: (Start)
Let m be a(n) or a prime factor of a(n). Then, except for 1 and 5, there are, if m is a prime, exactly two squares y^2 such that the difference y^2 - m contains exactly one pair of factors {x,z} such that the following applies: x*z = y^2 - m, x + y = z with
x < y, where {x,y,z} are relatively prime numbers. {x,y,z} are the initial values of a sequence of the Fibonacci type. Thus each a(n) > 5, if it is a prime, and each prime factor p > 5 of an a(n) can be assigned to exactly two sequences of the Fibonacci type. a(0) = 1 belongs to the original Fibonacci sequence and a(1) = 5 to the Lucas sequence.
But also the reverse assignment applies. From any sequence (f(i)) of the Fibonacci type we get from its 3 initial values by f(i)^2 - f(i-1)*f(i+1) with f(i-1) < f(i) a term a(n) or a prime factor p of a(n). This relation is also valid for any i. In this case we get the absolute value |a(n)| or |p|. (End)
a(n-1) = 2*T(n) - 1, for n>=1, with T = A000217, is a proper subsequence of A089270, and the terms are 0,-1,+1 (mod 5). - Wolfdieter Lang, Jul 05 2019
a(n+1) is the number of wedged n-dimensional spheres in the homotopy of the neighborhood complex of Kneser graph KG_{2,n}. Here, KG_{2,n} is a graph whose vertex set is the collection of subsets of cardinality 2 of set {1,2,...,n+3,n+4} and two vertices are adjacent if and only if they are disjoint. - Anurag Singh, Mar 22 2021
Also the number of squares between (n+2)^2 and (n+2)^4. - Karl-Heinz Hofmann, Dec 07 2021
(x, y, z) = (A001105(n+1), -a(n-1), -a(n)) are solutions of the Diophantine equation x^3 + 4*y^3 + 4*z^3 = 8. - XU Pingya, Apr 25 2022
The least significant digit of terms of this sequence cycles through 1, 5, 1, 9, 9. - Torlach Rush, Jun 05 2024

			From _Ilya Gutkovskiy_, Apr 13 2016: (Start)
Illustration of initial terms:
                                        o               o
                        o           o   o o           o o
            o       o   o o       o o   o o o       o o o
    o   o   o o   o o   o o o   o o o   o o o o   o o o o
o   o o o   o o o o o   o o o o o o o   o o o o o o o o o
n=0  n=1       n=2           n=3               n=4
(End)
From _Klaus Purath_, Mar 18 2019: (Start)
Examples:
a(0) = 1: 1^1-0*1 = 1, 0+1 = 1 (Fibonacci A000045).
a(1) = 5: 3^2-1*4 = 5, 1+3 = 4 (Lucas A000032).
a(2) = 11: 4^2-1*5 = 11, 1+4 = 5 (A000285); 5^2-2*7 = 11, 2+5 = 7 (A001060).
a(3) = 19: 5^2-1*6 = 19, 1+5 = 6 (A022095); 7^2-3*10 = 19, 3+7 = 10 (A022120).
a(4) = 29: 6^2-1*7 = 29, 1+6 = 7 (A022096); 9^2-4*13 = 29, 4+9 = 13 (A022130).
a(11)/5 = 31: 7^2-2*9 = 31, 2+7 = 9 (A022113); 8^2-3*11 = 31, 3+8 = 11 (A022121).
a(24)/11 = 59: 9^2-2*11 = 59, 2+9 = 11 (A022114); 12^2-5*17 = 59, 5+12 = 17 (A022137).
(End)

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Patrick De Geest, World!Of Numbers, Palindromes of the form n+(n+1)^2.
Adalbert Kerber, A matrix of combinatorial numbers related to the symmetric groups, Discrete Math., 21 (1978), 319-321. [Annotated scanned copy]
Clark Kimberling, Complementary Equations, Journal of Integer Sequences, Vol. 10 (2007), Article 07.1.4.
Nandini Nilakantan and Anurag Singh, Homotopy type of neighborhood complexes of Kneser graphs, KG_{2,k}, Proceeding-Mathematical Sciences, 128, Article number: 53(2018).
Yanni Pei and Jiang Zeng, Counting signed derangements with right-to-left minima and excedances, arXiv:2206.11236 [math.CO], 2022.
Popular Computing (Calabasas, CA), The CSR Function, Vol. 4 (No. 34, Jan 1976), pages PC34-10 to PC34-11. Annotated and scanned copy.
Zdzislaw Skupień and Andrzej Żak, Pair-sums packing and rainbow cliques, in Topics In Graph Theory, A tribute to A. A. and T. E. Zykovs on the occasion of A. A. Zykov's 90th birthday, ed. R. Tyshkevich, Univ. Illinois, 2013, pages 131-144, (in English and Russian).
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Complement of A028392. Third column of array A094954.
Cf. A000217, A002522, A062392, A062786, A127701, A135223, A143596, A052905, A162997, A062938 (squares of this sequence).
A110331 and A165900 are signed versions.
Cf. A002327 (primes), A094210.
Frobenius number for k successive numbers: this sequence (k=2), A079326 (k=3), A138984 (k=4), A138985 (k=5), A138986 (k=6), A138987 (k=7), A138988 (k=8).
Cf. also A135223, A176271, A214604, A105728, A005408, A002378, A084990, A253607, A260910, A089270.

a028387 n = n + (n + 1) ^ 2  -- Reinhard Zumkeller, Jul 17 2014

[n + (n+1)^2: n in [0..60]]; // Vincenzo Librandi, Apr 26 2011

FoldList[## + 2 &, 1, 2 Range@ 45] (* Robert G. Wilson v, Feb 02 2011 *)
Table[n + (n + 1)^2, {n, 0, 100}] (* Vincenzo Librandi, Oct 17 2012 *)
Table[ FrobeniusNumber[{n, n + 1}], {n, 2, 30}] (* Zak Seidov, Jan 14 2015 *)

a(n)=n^2+3*n+1 \\ Charles R Greathouse IV, Jun 10 2011

def a(n): return (n**2+3*n+1) # Torlach Rush, May 07 2024

[n+(n+1)^2 for n in range(0,48)] # Zerinvary Lajos, Jul 03 2008

a(n) = sqrt(A062938(n)). - Floor van Lamoen, Oct 08 2001
a(0) = 1, a(1) = 5, a(n) = (n+1)*a(n-1) - (n+2)*a(n-2) for n > 1. - Gerald McGarvey, Sep 24 2004
a(n) = A105728(n+2, n+1). - Reinhard Zumkeller, Apr 18 2005
a(n) = A109128(n+2, 2). - Reinhard Zumkeller, Jun 20 2005
a(n) = 2*T(n+1) - 1, where T(n) = A000217(n). - Gary W. Adamson, Aug 15 2007
a(n) = A005408(n) + A002378(n); A084990(n+1) = Sum_{k=0..n} a(k). - Reinhard Zumkeller, Aug 20 2007
Binomial transform of [1, 4, 2, 0, 0, 0, ...] = (1, 5, 11, 19, ...). - Gary W. Adamson, Sep 20 2007
G.f.: (1+2*x-x^2)/(1-x)^3. a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - R. J. Mathar, Jul 11 2009
a(n) = (n + 2 + 1/phi) * (n + 2 - phi); where phi = 1.618033989... Example: a(3) = 19 = (5 + .6180339...) * (3.381966...). Cf. next to leftmost column in A162997 array. - Gary W. Adamson, Jul 23 2009
a(n) = a(n-1) + 2*(n+1), with n > 0, a(0) = 1. - Vincenzo Librandi, Nov 18 2010
For k < n, a(n) = (k+1)*a(n-k) - k*a(n-k-1) + k*(k+1); e.g., a(5) = 41 = 4*11 - 3*5 + 3*4. - Charlie Marion, Jan 13 2011
a(n) = lower right term in M^2, M = the 2 X 2 matrix [1, n; 1, (n+1)]. - Gary W. Adamson, Jun 29 2011
G.f.: (x^2-2*x-1)/(x-1)^3 = G(0) where G(k) = 1 + x*(k+1)*(k+4)/(1 - 1/(1 + (k+1)*(k+4)/G(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 16 2012
Sum_{n>0} 1/a(n) = 1 + Pi*tan(sqrt(5)*Pi/2)/sqrt(5). - Enrique Pérez Herrero, Oct 11 2013
E.g.f.: exp(x) (1+4*x+x^2). - Tom Copeland, Dec 02 2013
a(n) = A005408(A000217(n)). - Tony Foster III, May 31 2016
From Amiram Eldar, Jan 29 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2).
Product_{n>=1} (1 - 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2)/6. (End)
a(5*n+1)/5 = A062786(n+1). - Torlach Rush, Jun 05 2024

Minor edits by N. J. A. Sloane, Jul 04 2010, following suggestions from the Sequence Fans Mailing List

A005574 Numbers k such that k^2 + 1 is prime.

Original entry on oeis.org

1, 2, 4, 6, 10, 14, 16, 20, 24, 26, 36, 40, 54, 56, 66, 74, 84, 90, 94, 110, 116, 120, 124, 126, 130, 134, 146, 150, 156, 160, 170, 176, 180, 184, 204, 206, 210, 224, 230, 236, 240, 250, 256, 260, 264, 270, 280, 284, 300, 306, 314, 326, 340, 350, 384, 386, 396

Offset: 1

N. J. A. Sloane

Hardy and Littlewood conjectured that the asymptotic number of elements in this sequence not exceeding n is approximately c*sqrt(n)/log(n) for some constant c. - Stefan Steinerberger, Apr 06 2006
Also, nonnegative integers such that a(n)+i is a Gaussian prime. - Maciej Ireneusz Wilczynski, May 30 2011
Apparently Goldbach conjectured that any a > 1 from this sequence can be written as a=b+c where b and c are in this sequence (Lemmermeyer link below). - Jeppe Stig Nielsen, Oct 14 2015
No term > 2 can be both in this sequence and in A001105 because of the Aurifeuillean factorization (2*k^2)^2 + 1 = (2*k^2 - 2*k + 1) * (2*k^2 + 2*k + 1). - Jeppe Stig Nielsen, Aug 04 2019

Harvey Dubner, "Generalized Fermat primes", J. Recreational Math., 18 (1985): 279-280.
R. K. Guy, "Unsolved Problems in Number Theory", 3rd edition, A2.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, p. 15, Thm. 17.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 1..10000
H. Dubner, Generalized Fermat primes, J. Recreational Math. 18.4 (1985-1986), 279. (Annotated scanned copy)
F. Ellermann, Primes of the form (m^2)+1 up to 10^6.
L. Euler, Lettre CXLIX (to Goldbach), 1752.
L. Euler, De numeris primis valde magnis, Novi Commentarii academiae scientiarum Petropolitanae 9 (1764), pp. 99-153. See pp. 123-125.
R. K. Guy, Letter to N. J. A. Sloane, 1988-04-12 (annotated scanned copy).
F. Lemmermeyer, Primes of the form a^2+1, Math Overflow question (2010).
Eric Weisstein's World of Mathematics, Landau's Problems.
Eric Weisstein's World of Mathematics, Power.
Eric Weisstein's World of Mathematics, Near-Square Prime.
Marek Wolf, Search for primes of the form m^2+1, arXiv:0803.1456 [math.NT], 2008-2010.

Cf. A002522, A001912, A002496, A062325, A090693, A000068, A006314, A006313, A006315, A006316, A056994, A056995, A057465, A057002, A088361, A088362, A226528, A226529, A226530, A251597, A253854, A244150, A243959, A321323.
Other sequences of the type "Numbers k such that k^2 + i is prime": this sequence (i=1), A067201 (i=2), A049422 (i=3), A007591 (i=4), A078402 (i=5), A114269 (i=6), A114270 (i=7), A114271 (i=8), A114272 (i=9), A114273 (i=10), A114274 (i=11), A114275 (i=12).
Cf. A010051, A259645, A295405 (characteristic function).

a005574 n = a005574_list !! (n-1)
a005574_list = filter ((== 1) . a010051' . (+ 1) . (^ 2)) [0..]
-- Reinhard Zumkeller, Jul 03 2015

[n: n in [0..400] | IsPrime(n^2+1)]; // Vincenzo Librandi, Nov 18 2010

Select[Range[350], PrimeQ[ #^2 + 1] &] (* Stefan Steinerberger, Apr 06 2006 *)
Join[{1},2Flatten[Position[PrimeQ[Table[x^2+1,{x,2,1000,2}]],True]]]  (* Fred Patrick Doty, Aug 18 2017 *)

isA005574(n) = isprime(n^2+1) \\ Michael B. Porter, Mar 20 2010

for(n=1, 1e3, if(isprime(n^2 + 1), print1(n, ", "))) \\ Altug Alkan, Oct 14 2015

from sympy import isprime; [print(n, end = ', ') for n in range(1, 400) if isprime(n*n+1)] # Ya-Ping Lu, Apr 23 2025

a(n) = A090693(n) - 1.

a(n) = 2*A001912(n-1) for n > 1. - Jeppe Stig Nielsen, Aug 04 2019

A028982 Squares and twice squares.

Original entry on oeis.org

1, 2, 4, 8, 9, 16, 18, 25, 32, 36, 49, 50, 64, 72, 81, 98, 100, 121, 128, 144, 162, 169, 196, 200, 225, 242, 256, 288, 289, 324, 338, 361, 392, 400, 441, 450, 484, 512, 529, 576, 578, 625, 648, 676, 722, 729, 784, 800, 841, 882, 900, 961, 968, 1024

Offset: 1

Patrick De Geest

Numbers n such that sum of divisors of n (A000203) is odd.
Also the numbers with an odd number of run sums (trapezoidal arrangements, number of ways of being written as the difference of two triangular numbers). - Ron Knott, Jan 27 2003
Pell(n)*Sum_{k|n} 1/Pell(k) is odd, where Pell(n) is A000129(n). - Paul Barry, Oct 12 2005
Number of odd divisors of n (A001227) is odd. - Vladeta Jovovic, Aug 28 2007
A071324(a(n)) is odd. - Reinhard Zumkeller, Jul 03 2008
Sigma(a(n)) = A000203(a(n)) = A152677(n). - Jaroslav Krizek, Oct 06 2009
Numbers n such that sum of odd divisors of n (A000593) is odd. - Omar E. Pol, Jul 05 2016
A187793(a(n)) is odd. - Timothy L. Tiffin, Jul 18 2016
If k is odd (k = 2m+1 for m >= 0), then 2^k = 2^(2m+1) = 2*(2^m)^2.  If k is even (k = 2m for m >= 0), then 2^k = 2^(2m) = (2^m)^2.  So, the powers of 2 sequence (A000079) is a subsequence of this one. - Timothy L. Tiffin, Jul 18 2016
Numbers n such that A175317(n) = Sum_{d|n} pod(d) is odd, where pod(m) = the product of divisors of m (A007955). - Jaroslav Krizek, Dec 28 2016
Positions of zeros in A292377 and A292383, positions of ones in A286357 and A292583. (See A292583 for why.) - Antti Karttunen, Sep 25 2017
Numbers of the form A000079(i)*A016754(j), i,j>=0. - R. J. Mathar, May 30 2020
Equivalently, numbers whose odd part is square. Cf. A042968. - Peter Munn, Jul 14 2020
These are the Heinz numbers of the partitions counted by A119620. - Gus Wiseman, Oct 29 2021
Numbers m whose abundance, A033880(m), is odd. - Peter Munn, May 23 2022
Numbers with an odd number of middle divisors (cf. A067742). - Omar E. Pol, Aug 02 2022

T. D. Noe, Table of n, a(n) for n = 1..1000
Tewodros Amdeberhan, Victor H. Moll, Vaishavi Sharma, and Diego Villamizar, Arithmetic properties of the sum of divisors, arXiv:2007.03088 [math.NT], 2020. See p. 5.
J. N. Cooper and A. W. N. Riasanovsky, On the Reciprocal of the Binary Generating Function for the Sum of Divisors, Journal of Integer Sequences, Vol. 16 (2013), #13.1.8.
Patrick De Geest, World!Of Numbers
John S. Rutherford, Sublattice enumeration. IV. Equivalence classes of plane sublattices by parent Patterson symmetry and colour lattice group type, Acta Cryst. (2009). A65, 156-163.
Eric Weisstein's World of Mathematics, Abundance

Complement of A028983.
Characteristic function is A053866, A093709.
Odd terms in A178910.
Cf. A000203, A000290, A000593, A001105, A042968, A187793.
Supersequence of A000079.
Cf. A028260, A033880, A046951, A067742.
Cf. A119620, A119899, A347437, A347438, A348550.

import Data.List.Ordered (union)
a028982 n = a028982_list !! (n-1)
a028982_list = tail $ union a000290_list a001105_list
-- Reinhard Zumkeller, Jun 27 2015

Take[ Sort[ Flatten[ Table[{n^2, 2n^2}, {n, 35}] ]], 57] (* Robert G. Wilson v, Aug 27 2004 *)

list(lim)=vecsort(concat(vector(sqrtint(lim\1),i,i^2), vector(sqrtint(lim\2),i,2*i^2))) \\ Charles R Greathouse IV, Jun 16 2011

from itertools import count, islice
from sympy.ntheory.primetest import is_square
def A028982_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:int(is_square(n) or is_square(n<<1)),count(max(startvalue,1)))
A028982_list = list(islice(A028982_gen(),30)) # Chai Wah Wu, Jan 09 2023

from math import isqrt
def A028982(n):
    def f(x): return n-1+x-isqrt(x)-isqrt(x>>1)
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 22 2024

A001105 UNION A000290.
a(n) is asymptotic to c*n^2 with c = 2/(1+sqrt(2))^2 = 0.3431457.... - Benoit Cloitre, Sep 17 2002
In particular, a(n) = c*n^2 + O(n). - Charles R Greathouse IV, Jan 11 2013
a(A003152(n)) = n^2; a(A003151(n)) = 2*n^2. - Enrique Pérez Herrero, Oct 09 2013
Sum_{n>=1} 1/a(n) = Pi^2/4. - Amiram Eldar, Jun 28 2020

A016742 Even squares: a(n) = (2*n)^2.

Original entry on oeis.org

0, 4, 16, 36, 64, 100, 144, 196, 256, 324, 400, 484, 576, 676, 784, 900, 1024, 1156, 1296, 1444, 1600, 1764, 1936, 2116, 2304, 2500, 2704, 2916, 3136, 3364, 3600, 3844, 4096, 4356, 4624, 4900, 5184, 5476, 5776, 6084, 6400, 6724, 7056, 7396, 7744, 8100, 8464

Offset: 0

N. J. A. Sloane

4 times the squares.
Number of edges in the complete bipartite graph of order 5n, K_{n,4n} - Roberto E. Martinez II, Jan 07 2002
It is conjectured (I think) that a regular Hadamard matrix of order n exists iff n is an even square (cf. Seberry and Yamada, Th. 10.11). A Hadamard matrix is regular if the sum of the entries in each row is the same. - N. J. A. Sloane, Nov 13 2008
Sequence arises from reading the line from 0, in the direction 0, 16, ... and the line from 4, in the direction 4, 36, ... in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008
The entries from a(1) on can be interpreted as pair sums of (2, 2), (8, 8), (18, 18), (32, 32) etc. that arise from a re-arrangement of the subshell orbitals in the periodic table of elements. 8 becomes the maximum number of electrons in the (2s,2p) or (3s,3p) orbitals, 18 the maximum number of electrons in (4s,3d,4p) or (5s,3d,5p) shells, for example. - Julio Antonio Gutiérrez Samanez, Jul 20 2008
The first two terms of the sequence (n=1, 2) give the numbers of chemical elements using only n types of atomic orbitals, i.e., there are a(1)=4 elements (H,He,Li,Be) where electrons reside only on s-orbitals, there are a(2)=16 elements (B,C,N,O,F,Ne,Na,Mg,Al,Si,P,S,Cl,Ar,K,Ca) where electrons reside only on s- and p-orbitals. However, after that, there is 37 (which is one more than a(3)=36) elements (from Sc, Scandium, atomic number 21 to La, Lanthanum, atomic number 57) where electrons reside only on s-, p- and d-orbitals. This is because Lanthanum (with the electron configuration [Xe]5d^1 6s^2) is an exception to the Aufbau principle, which would predict that its electron configuration is [Xe]4f^1 6s^2. - Antti Karttunen, Aug 14 2008.
Number of cycles of length 3 in the king's graph associated with an (n+1) X (n+1) chessboard. - Anton Voropaev (anton.n.voropaev(AT)gmail.com), Feb 01 2009
a(n+1) is the molecular topological index of the n-star graph S_n. - Eric W. Weisstein, Jul 11 2011
a(n) is the sum of two consecutives odd numbers 2*n^2-1 and 2*n^2+1 and the difference of two squares (n^2+1)^2 - (n^2-1)^2. - Pierre CAMI, Jan 02 2012
For n > 3, a(n) is the area of the irregular quadrilateral created by the points ((n-4)*(n-3)/2,(n-3)*(n-2)/2), ((n-2)*(n-1)/2,(n-1)*n/2), ((n+1)*(n+2)/2,n*(n+1)/2), and ((n+3)*(n+4)/2,(n+2)*(n+3)/2). - J. M. Bergot, May 27 2014
Number of terms less than 10^k: 1, 2, 5, 16, 50, 159, 500, 1582, 5000, 15812, 50000, 158114, 500000, ... - Muniru A Asiru, Jan 28 2018
Right-hand side of the binomial coefficient identity Sum_{k = 0..2*n} (-1)^(k+1)* binomial(2*n,k)*binomial(2*n + k,k)*(2*n - k) = a(n). - Peter Bala, Jan 12 2022

R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd ed., 1994, p. 99.
Seberry, Jennifer and Yamada, Mieko; Hadamard matrices, sequences and block designs, in Dinitz and Stinson, eds., Contemporary design theory, pp. 431-560, Wiley-Intersci. Ser. Discrete Math. Optim., Wiley, New York, 1992.
W. D. Wallis, Anne Penfold Street and Jennifer Seberry Wallis, Combinatorics: Room squares, sum-free sets, Hadamard matrices, Lecture Notes in Mathematics, Vol. 292, Springer-Verlag, Berlin-New York, 1972. iv+508 pp.

Vincenzo Librandi, Table of n, a(n) for n = 0..900
R. P. Boas and N. J. A. Sloane, Correspondence, 1974.
Leo Tavares, Illustration: X Squares
Various, Electron Configuration (Discussion in Physics Forums).
Eric Weisstein's World of Mathematics, Graph Cycle.
Eric Weisstein's World of Mathematics, King Graph.
Eric Weisstein's World of Mathematics, Molecular Topological Index.
Wikipedia, Aufbau principle.
Index entries for sequences related to Hadamard matrices
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000290, A001105, A001539, A016754, A016802, A016814, A016826, A016838, A007742, A033991, A245058.
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.
Cf. sequences listed in A254963.
Other n X n king graph cycle counts: A288918 (4-cycles), A288919 (5-cycles), A288920 (6-cycles).
Cf. A000384, A014105.
Cf. A016813.

List([0..100], n -> (2*n)^2); # Muniru A Asiru, Jan 28 2018

a016742 = (* 4) . (^ 2)
a016742_list = 0 : map (subtract 4) (zipWith (+) a016742_list [8, 16 ..])
-- Reinhard Zumkeller, Jun 28 2015, Apr 20 2015

[(2*n)^2: n in [0..50]]; // Vincenzo Librandi, Apr 26 2011

seq((2*n)^2, n=0..100); # Muniru A Asiru, Jan 28 2018

Table[(2n)^2, {n, 0, 46}] (* Alonso del Arte, Apr 26 2011 *)

makelist((2*n)^2,n,0,20); /* Martin Ettl, Jan 22 2013 */

a(n)=4*n^2 \\ Charles R Greathouse IV, Jul 28 2015

O.g.f.: 4*x*(1+x)/(1-x)^3. - R. J. Mathar, Jul 28 2008
a(n) = A000290(n)*4 = A001105(n)*2. - Omar E. Pol, May 21 2008
a(n) = A155955(n,2) for n > 1. - Reinhard Zumkeller, Jan 31 2009
Sum_{n>=1} 1/a(n) = (1/4)*Pi^2/6 = Pi^2/24. - Ant King, Nov 04 2009
a(n) = a(n-1) + 8*n - 4 (with a(0)=0). - Vincenzo Librandi, Nov 19 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) with a(0) = 0, a(1) = 4, a(2) = 16. - Philippe Deléham, Mar 26 2013
a(n) = A118729(8n+3). - Philippe Deléham, Mar 26 2013
Pi = 2*Product_{n>=1} (1 + 1/(a(n)-1)). - Adriano Caroli, Aug 04 2013
Pi = Sum_{n>=0} 8/(a(2n+1)-1). - Adriano Caroli, Aug 06 2013
E.g.f.: exp(x)*(4x^2 + 4x). - Geoffrey Critzer, Oct 07 2013
a(n) = A000384(n) + A014105(n). - Bruce J. Nicholson, Nov 11 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/48 (A245058). - Amiram Eldar, Oct 10 2020
From Amiram Eldar, Jan 25 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = sinh(Pi/2)/(Pi/2) (A308716).
Product_{n>=1} (1 - 1/a(n)) = sin(Pi/2)/(Pi/2) = 2/Pi (A060294). (End)
a(n) = A016754(n) - A016813(n). - Leo Tavares, Feb 24 2022

More terms from Sabir Abdus-Samee (sabdulsamee(AT)prepaidlegal.com), Mar 13 2006

A001951 A Beatty sequence: a(n) = floor(n*sqrt(2)).

Original entry on oeis.org

0, 1, 2, 4, 5, 7, 8, 9, 11, 12, 14, 15, 16, 18, 19, 21, 22, 24, 25, 26, 28, 29, 31, 32, 33, 35, 36, 38, 39, 41, 42, 43, 45, 46, 48, 49, 50, 52, 53, 55, 56, 57, 59, 60, 62, 63, 65, 66, 67, 69, 70, 72, 73, 74, 76, 77, 79, 80, 82, 83, 84, 86, 87, 89, 90, 91, 93, 94, 96, 97, 98, 100

Offset: 0

N. J. A. Sloane

Earliest monotonic sequence greater than 0 satisfying the condition: "a(n) + 2n is not in the sequence". - Benoit Cloitre, Mar 25 2004
Also the integer part of the hypotenuse of isosceles right triangles. The real part of these numbers is irrational. For proof see Jones and Jones.
First differences are 1, 1, 2, 1, 2, 1, 1, 2, 1, 2, 1, 1, 2, ... (A006337 with a 1 in front). - Philippe Deléham, May 29 2006
It appears that the distance between the a(n)-th triangular number and the nearest square is not greater than floor(a(n)/2). - Ralf Stephan, Sep 14 2013
These are the nonnegative integers m satisfying sin(m*Pi/r)*sin((m+1)*Pi/r) <= 0, where r = sqrt(2). In general, the Beatty sequence of an irrational number r > 1 consists of the numbers m satisfying sin(m*x)*sin((m+1)*x) <= 0, where x = Pi/r. Thus the numbers m satisfying sin(m*x)*sin((m+1)*x) > 0 form the Beatty sequence of r/(1-r). - Clark Kimberling, Aug 21 2014
For n > 0: A080764(a(n)) = 1. - Reinhard Zumkeller, Jul 03 2015
From Clark Kimberling, Oct 17 2016: (Start)
We can generate A001951 and A001952 without using sqrt(2).
First write the even positive integers in a row:
  2   4   6   8   10   12   14 . . .
Then put 1 under 2 and add:
  2   4   6   8   10   12   14 . . .
  1
  3
Next, under 4, put the least positive integer that is not yet in rows 2 and 3;
it is 2; and add:
  2   4   6   8   10   12   14 . . .
  1   2
  3   6
Next, under the 6 in row 1, put the least positive integer not yet in rows 2 and 3;
it is 4, and add:
  2   4   6   8   10   12   14 . . .
  1   2   4
  3   6   10
Continue in this manner. (End)
This sequence contains an infinite number of powers of 2 (proof in Crux Mathematicorum link). See A103341. - Bernard Schott, Mar 08 2019
The terms of this sequence generate the multiplicative group of positive rational numbers (observation by Stephen M. Gagola, Jr.; see References). - Allen Stenger, Aug 05 2023
a(n) is also the number of distinct straight cylinders with integer radius and height having the same surface as a sphere with radius n. - Felix Huber, Sep 20 2024
Let P(x,y) be the condition x^2 + y^2 <= n^2, then 4*a(n) is the number of integer points (x,y) such that P(x,y) is true and at least one of P(x+1,y), P(x-1,y), P(x,y+1), P(x,y-1) is false. See LINKS for examples. - Bob de Boisvilliers, May 14 2025

Eric Duchêne, Aviezri S. Fraenkel, Vladimir Gurvich, Nhan Bao Ho, Clark Kimberling, Urban Larsson, Wythoff Visions, Games of No Chance, Vol. 5; MSRI Publications, Vol. 70 (2017), pages 101-153.
Stephen M. Gagola Jr., Solution of Problem 12282, Am. Math. Monthly, 130 (2023), pp. 682-683.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 77.
Gareth A. Jones and J. Mary Jones, Elementary Number Theory, Springer, 1998; pp. 221-222.
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).
Roland Sprague, Recreations in Mathematics, Blackie and Son, (1963).
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, Revised edition (1997), Entry sqrt(2), p. 18.

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
L. Carlitz, Richard Scoville and Verner E. Hoggatt, Jr., Pellian representatives, Fib. Quart., 10 (1972), 449-488.
Ian G. Connell, A generalization of Wythoff's game, Canad. Math. Bull. 2 (1959) 181-190
Bob de Boisvilliers, ASCII circles' boundary pixels count examples from n=0 to n=10
Ed Doolittle, Problem 19, 26th I.M.O. Finland proposed by Romania, Crux Mathematicorum, p. 70, Vol. 14, Mar. 88.
Aviezri S. Fraenkel, How to beat your Wythoff games' opponent on three fronts, Amer. Math. Monthly, 89 (1982), 353-361 (the case a=2).
Aviezri S. Fraenkel, On the recurrence f(m+1)= b(m)*f(m)-f(m-1) and applications, Discrete Mathematics 224 (2000), no. 1-3, pp. 273-279.
Wen An Liu and Xiao Zhao, Adjoining to (s,t)-Wythoff's game its P-positions as moves, Discrete Applied Mathematics, Aug 27 2014; See Table 3.
Luke Schaeffer, Jeffrey Shallit, and Stefan Zorcic, Beatty Sequences for a Quadratic Irrational: Decidability and Applications, arXiv:2402.08331 [math.NT], 2024. See pp. 17-18.
N. J. A. Sloane, Families of Essentially Identical Sequences, Mar 24 2021 (Includes this sequence)
Eric Weisstein's World of Mathematics, Beatty Sequence.
Index entries for sequences related to Beatty sequences

Complement of A001952. Equals A001952(n) - 2*n for n>0.
Equals A003151(n) - n; a bisection of A094077.
Bisections: A022842, A342281.
Cf. A022342, A026250, A080764, A103341.
The following sequences are all essentially the same, in the sense that they are simple transformations of each other, with A003151 as the parent: A003151, A001951, A001952, A003152, A006337, A080763, A082844 (conjectured), A097509, A159684, A188037, A245219 (conjectured), A276862. - N. J. A. Sloane, Mar 09 2021
Partial sums: A194102.

a001951 = floor . (* sqrt 2) . fromIntegral
-- Reinhard Zumkeller, Sep 14 2014

[Floor(n*Sqrt(2)): n in [0..60]]; // Vincenzo Librandi, Oct 22 2011

[Isqrt(2*n^2):n in[0..60]]; // Jason Kimberley, Oct 28 2016

a:=n->floor(n*sqrt(2)): seq(a(n),n=0..80); # Muniru A Asiru, Mar 09 2019

Floor[Range[0, 72] Sqrt[2]] (* Robert G. Wilson v, Oct 17 2012 *)

makelist(floor(n*sqrt(2)), n, 0, 100); /* Martin Ettl, Oct 17 2012 */

f(n) = for(j=1,n,print1(floor(sqrt(2*j^2))","))

a(n)=sqrtint(2*n^2) \\ Charles R Greathouse IV, Oct 19 2016

from sympy import integer_nthroot
def A001951(n): return integer_nthroot(2*n**2,2)[0] # Chai Wah Wu, Mar 16 2021

a(n) = A000196(A001105(n)). - Jason Kimberley, Oct 26 2016
a(n) = floor(csc(1/(sqrt(2)*n))) for n > 0, since sqrt(2)*n < csc(1/(sqrt(2)*n)) < sqrt(2)*n + 1/(3*sqrt(2)*n) < floor(sqrt(2)*n) + 1 for n > 0. - Jianing Song, Sep 07 2021
a(n) = A194102(n) - A194102(n-1) for n > 0. - M. F. Hasler, Apr 23 2022

More terms from David W. Wilson, Sep 20 2000

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

Original entry on oeis.org

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

Offset: 0

N. J. A. Sloane

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

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

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

Vincenzo Librandi, Table of n, a(n) for n = 0..3000
Nikolai Beluhov, Snake paths in king and knight graphs, arXiv:2301.01152 [math.CO], 2023.
M. Benoumhani and M. Kolli, Finite topologies and partitions, JIS 13 (2010) # 10.3.5, t_{N0}(n,4) in theorem 5.
Andrea C. Burgess, Caleb W. Jones, and David A. Pike, Extending Graph Burning to Hypergraphs, arXiv:2403.01001 [math.CO], 2024. See p. 9.
Geoffrey B. Campbell, Vector Partition Identities for 2D, 3D and nD Lattices, arXiv:2302.01091 [math.CO], 2023.
Ronald Cools, Ian H. Sloan, Minimial cubature formulae of trigonometric degree, Math. Comp. 65 (216) (1996) 1583-1600. Table 1 dimension 2.
John Elias, Illustration of Initial Terms: Intersection of a double spaced square grid and centrally aligned triangle.
J. G. Kalbfleisch and R. G. Stanton, A combinatorial problem in matching, J. London Math. Soc. Vol. 1, No. 1 (1969), 60-64. [Corrected by _N. J. A. Sloane_, Feb 08 2019]
J. M. Kantor, Mathématiques venues d'ailleurs: divertissements mathématiques en U.R.S.S., Le cube transpercé, pp. 56-62, Belin, Paris, 1982.
S. Lafortune, A. Ramani, B. Grammaticos, Y. Ohta, and K.M. Tamizhmani, Blending two discrete integrability criteria: ..., arXiv:nlin/0104020 [nlin.SI], 2001.
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
Lara Pudwell, From permutation patterns to the periodic table, Notices of the American Mathematical Society, 67.7 (2020), 994-1001.
Eric Weisstein's World of Mathematics, Black Bishop Graph
Eric Weisstein's World of Mathematics, Edge Cover Number
Eric Weisstein's World of Mathematics, King Graph
Eric Weisstein's World of Mathematics, Topology
Eric Weisstein's World of Mathematics, Vertex Count
Index entries for linear recurrences with constant coefficients, signature (2,0,-2,1).

Cf. A000096, A000217, A001105, A001477, A001844, A002061, A004526, A005843, A007590, A008794, A037270, A081352, A109613, A110654, A116940, A134444, A158946, A168380, A357501.
Column 2 of A195040.
Cf. also A132188, A362931.

a000982 = (`div` 2) . (+ 1) . (^ 2)  -- Reinhard Zumkeller, Jun 27 2013

[(2*n^2 + 1 - (-1)^n) / 4: n in [0..60]]; // Vincenzo Librandi, Jun 16 2011

seq( ceil(n^2/2),n=0..30) ; # R. J. Mathar, Jun 05 2011

Table[Ceiling[n^2/2], {n, 0, 120}] (* Vladimir Joseph Stephan Orlovsky, Apr 02 2011 *)
Accumulate[Join[{0}, (# - Boole[EvenQ[#]] &) /@ Range[80]]] (* Alonso del Arte, Sep 11 2019 *)

a(n)=(n^2+1)\2 \\ Charles R Greathouse IV, Sep 13 2013

x='x+O('x^100); concat([0], Vec(x*(1+x^2)/((1+x)*(1-x)^3))) \\ Altug Alkan, Oct 12 2015

apply( A000982(n)=n^2\/2, [0..55]) \\ M. F. Hasler, Feb 29 2020

def A000982(n): return n**2+1>>1 # Chai Wah Wu, Aug 28 2023

(((1 to 49) by 2) flatMap { List.fill(2)() }).scanLeft(0)( + ) // _Alonso del Arte, Sep 11 2019

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

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

Original entry on oeis.org

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

Offset: 0

Erich Friedman, Dec 12 2000

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

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

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

G. C. Greubel, Table of n, a(n) for n = 0..5000
Steven Edwards and William Griffiths, On Generalized Delannoy Numbers, J. Int. Seq., Vol. 23 (2020), Article 20.3.6.
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.
Leo Tavares, Illustration: Triangular Outlines
Reinhard Zumkeller, Enumerations of Divisors.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000124.
Second row of array A099597.
See A120062 for sequences related to integer-sided triangles with integer inradius n.
Cf. A112295.
Cf. A087113, A002552.
Cf. A005408, A016813, A086514, A000125, A002522, A161701, A161702, A161703, A000127, A161704, A161706, A161707, A161708, A161710, A080856, A161711, A161712, A161713, A161715, A006261.
Cf. A001079, A037270, A071253, A108741, A132592, A146311, A146312, A146313, A173115, A173116, A173121.
Column 2 of array A188645.
Cf. A001105 and A247375. - Bruno Berselli, Sep 16 2014
Cf. A056106, A251599.
Cf. A000384, A000217, A166926.

a058331 = (+ 1) . a001105  -- Reinhard Zumkeller, Dec 13 2014

[2*n^2 + 1 : n in [0..100]]; // Wesley Ivan Hurt, Feb 02 2017

b[g_] := Length[Union[Map[Det, Flatten[ Table[{{i, j}, {k, l}}, {i, 0, g}, {j, 0, g}, {k, 0, g}, {l, 0, g}], 3]]]] Table[b[g], {g, 0, 20}]
2*Range[0, 49]^2 + 1 (* Alonso del Arte, Dec 05 2012 *)

a(n)=2*n^2+1 \\ Charles R Greathouse IV, Jun 16 2011

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

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

A225546 Tek's flip: Write n as the product of distinct factors of the form prime(i)^(2^(j-1)) with i and j integers, and replace each such factor with prime(j)^(2^(i-1)).

Original entry on oeis.org

1, 2, 4, 3, 16, 8, 256, 6, 9, 32, 65536, 12, 4294967296, 512, 64, 5, 18446744073709551616, 18, 340282366920938463463374607431768211456, 48, 1024, 131072, 115792089237316195423570985008687907853269984665640564039457584007913129639936, 24, 81, 8589934592, 36, 768

Offset: 1

Paul Tek, May 10 2013

This is a multiplicative self-inverse permutation of the integers.
A225547 gives the fixed points.
From Antti Karttunen and Peter Munn, Feb 02 2020: (Start)
This sequence operates on the Fermi-Dirac factors of a number. As arranged in array form, in A329050, this sequence reflects these factors about the main diagonal of the array, substituting A329050[j,i] for A329050[i,j], and this results in many relationships including significant homomorphisms.
This sequence provides a relationship between the operations of squaring and prime shift (A003961) because each successive column of the A329050 array is the square of the previous column, and each successive row is the prime shift of the previous row.
A329050 gives examples of how significant sets of numbers can be formed by choosing their factors in relation to rows and/or columns. This sequence therefore maps equivalent derived sets by exchanging rows and columns. Thus odd numbers are exchanged for squares, squarefree numbers for powers of 2 etc.
Alternative construction: For n > 1, form a vector v of length A299090(n), where each element v[i] for i=1..A299090(n) is a product of those distinct prime factors p(i) of n whose exponent e(i) has the bit (i-1) "on", or 1 (as an empty product) if no such exponents are present. a(n) is then Product_{i=1..A299090(n)} A000040(i)^A048675(v[i]). Note that because each element of vector v is squarefree, it means that each exponent A048675(v[i]) present in the product is a "submask" (not all necessarily proper) of the binary string A087207(n).
This permutation effects the following mappings:
A000035(a(n)) = A010052(n), A010052(a(n)) = A000035(n). [Odd numbers <-> Squares]
A008966(a(n)) = A209229(n), A209229(a(n)) = A008966(n). [Squarefree numbers <-> Powers of 2]
(End)
From Antti Karttunen, Jul 08 2020: (Start)
Moreover, we see also that this sequence maps between A016825 (Numbers of the form 4k+2) and A001105 (2*squares) as well as between A008586 (Multiples of 4) and A028983 (Numbers with even sum of the divisors).
(End)

			  7744  = prime(1)^2^(2-1)*prime(1)^2^(3-1)*prime(5)^2^(2-1).
a(7744) = prime(2)^2^(1-1)*prime(3)^2^(1-1)*prime(2)^2^(5-1) = 645700815.

Cf. A225547 (fixed points) and the subsequences listed there.
Transposes A329050, A329332.
An automorphism of positive integers under the binary operations A059895, A059896, A059897, A306697, A329329.
An automorphism of A059897 subgroups: A000379, A003159, A016754, A122132.
Permutes lists where membership is determined by number of Fermi-Dirac factors: A000028, A050376, A176525, A268388.
Also permutes: A036554, A059404, A066427, A066428, A072774, A331593.
Sequences f that satisfy f(a(n)) = f(n): A048675, A064179, A064547, A097248, A302777, A331592.
Pairs of sequences (f,g) that satisfy a(f(n)) = g(a(n)): (A000265,A008833), (A000290,A003961), (A005843,A334747), (A006519,A007913), (A008586,A334748).
Pairs of sequences (f,g) that satisfy a(f(n)) = g(n), possibly with offset change: (A000040,A001146), (A000079,A019565).
Pairs of sequences (f,g) that satisfy f(a(n)) = g(n), possibly with offset change: (A000035, A010052), (A008966, A209229), (A007814, A248663), (A061395, A299090), (A087207, A267116), (A225569, A227291).
Cf. A331287 [= gcd(a(n),n)].
Cf. A331288 [= min(a(n),n)], see also A331301.
Cf. A331309 [= A000005(a(n)), number of divisors].
Cf. A331590 [= a(a(n)*a(n))].
Cf. A331591 [= A001221(a(n)), number of distinct prime factors], see also A331593.
Cf. A331740 [= A001222(a(n)), number of prime factors with multiplicity].
Cf. A331733 [= A000203(a(n)), sum of divisors].
Cf. A331734 [= A033879(a(n)), deficiency].
Cf. A331735 [= A009194(a(n))].
Cf. A331736 [= A000265(a(n)) = a(A008833(n)), largest odd divisor].
Cf. A335914 [= A038040(a(n))].
A self-inverse isomorphism between pairs of A059897 subgroups: (A000079,A005117), (A000244,A062503), (A000290\{0},A005408), (A000302,A056911), (A000351,A113849 U {1}), (A000400,A062838), (A001651,A252895), (A003586,A046100), (A007310,A000583), (A011557,A113850 U {1}), (A028982,A042968), (A053165,A065331), (A262675,A268390).
A bijection between pairs of sets: (A001248,A011764), (A007283,A133466), (A016825, A001105), (A008586, A028983).
Cf. also A336321, A336322 (compositions with another involution, A122111).

Array[If[# == 1, 1, Times @@ Flatten@ Map[Function[{p, e}, Map[Prime[Log2@ # + 1]^(2^(PrimePi@ p - 1)) &, DeleteCases[NumberExpand[e, 2], 0]]] @@ # &, FactorInteger[#]]] &, 28] (* Michael De Vlieger, Jan 21 2020 *)

A019565(n) = factorback(vecextract(primes(logint(n+!n, 2)+1), n));
a(n) = {my(f=factor(n)); for (i=1, #f~, my(p=f[i,1]); f[i,1] = A019565(f[i,2]); f[i,2] = 2^(primepi(p)-1);); factorback(f);} \\ Michel Marcus, Nov 29 2019

A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
A225546(n) = if(1==n,1,my(f=factor(n),u=#binary(vecmax(f[, 2])),prods=vector(u,x,1),m=1,e); for(i=1,u,for(k=1,#f~, if(bitand(f[k,2],m),prods[i] *= f[k,1])); m<<=1); prod(i=1,u,prime(i)^A048675(prods[i]))); \\ Antti Karttunen, Feb 02 2020

from math import prod
from sympy import prime, primepi, factorint
def A225546(n): return prod(prod(prime(i) for i, v in enumerate(bin(e)[:1:-1],1) if v == '1')**(1<Chai Wah Wu, Mar 17 2023

Multiplicative, with a(prime(i)^j) = A019565(j)^A000079(i-1).
a(prime(i)) = 2^(2^(i-1)).
From Antti Karttunen and Peter Munn, Feb 06 2020: (Start)
a(A329050(n,k)) = A329050(k,n).
a(A329332(n,k)) = A329332(k,n).
Equivalently, a(A019565(n)^k) = A019565(k)^n. If n = 1, this gives a(2^k) = A019565(k).
a(A059897(n,k)) = A059897(a(n), a(k)).
The previous formula implies a(n*k) = a(n) * a(k) if A059895(n,k) = 1.
a(A000040(n)) = A001146(n-1); a(A001146(n)) = A000040(n+1).
a(A000290(a(n))) = A003961(n); a(A003961(a(n))) = A000290(n) = n^2.
a(A000265(a(n))) = A008833(n); a(A008833(a(n))) = A000265(n).
a(A006519(a(n))) = A007913(n); a(A007913(a(n))) = A006519(n).
A007814(a(n)) = A248663(n); A248663(a(n)) = A007814(n).
A048675(a(n)) = A048675(n) and A048675(a(2^k * n)) = A048675(2^k * a(n)) = k + A048675(a(n)).
(End)
From Antti Karttunen and Peter Munn, Jul 08 2020: (Start)
For all n >= 1, a(2n) = A334747(a(n)).
In particular, for n = A003159(m), m >= 1, a(2n) = 2*a(n). [Note that A003159 includes all odd numbers]
(End)

Name edited by Peter Munn, Feb 14 2020

"Tek's flip" prepended to the name by Antti Karttunen, Jul 08 2020

A045944 Rhombic matchstick numbers: a(n) = n(3n+2).

Original entry on oeis.org

0, 5, 16, 33, 56, 85, 120, 161, 208, 261, 320, 385, 456, 533, 616, 705, 800, 901, 1008, 1121, 1240, 1365, 1496, 1633, 1776, 1925, 2080, 2241, 2408, 2581, 2760, 2945, 3136, 3333, 3536, 3745, 3960, 4181, 4408, 4641, 4880, 5125, 5376, 5633, 5896, 6165, 6440

Offset: 0

R. K. Guy

From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0, then a(n) is the n-th term of the sequence found by reading the line from 0 in the direction 0,5,.... The spiral begins:
                 .
                 85--84--83--82--81--80
                   .                   \
                   56--55--54--53--52  79
                   / .               \   \
                 57  33--32--31--30  51  78
                 /   / .           \   \   \
               58  34  16--15--14  29  50  77
               /   /   / .       \   \   \   \
             59  35  17   5---4  13  28  49  76
             /   /   /   / .   \   \   \   \   \
           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)
Connection to triangular numbers: a(n) = 4*T_n + S_n where T_n is the n-th triangular number and S_n is the n-th square. - William A. Tedeschi, Sep 12 2010
Also, second octagonal numbers. - Bruno Berselli, Jan 13 2011
Sequence found by reading the line from 0, in the direction 0, 16, ... and the line from 5, in the direction 5, 33, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. - Omar E. Pol, Jul 18 2012
Let P denote the points from the n X n grid. A(n-1) also coincides with the minimum number of points Q needed to "block" P, that is, every line segment spanned by two points from P must contain one point from Q. - Manfred Scheucher, Aug 30 2018
Also the number of internal edges of an (n+1)*(n+1) "square" of hexagons; i.e., n+1 rows, each of n+1 edge-adjacent hexagons, stacked with minimal overhang. - Jon Hart, Sep 29 2019
For n >= 1, the continued fraction expansion of sqrt(27*a(n)) is [9n+2; {1, 2n-1, 1, 1, 1, 2n-1, 1, 18n+4}]. - Magus K. Chu, Oct 13 2022

Ivan Panchenko, Table of n, a(n) for n = 0..1000
Ghislain R. Franssens, On a Number Pyramid Related to the Binomial, Deleham, Eulerian, MacMahon and Stirling number triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.4.1.
M. Janjic and B. Petkovic, A Counting Function, arXiv:1301.4550 [math.CO], 2013.
Leo Tavares, Illustration: Square Stars
Leo Tavares, Illustration: Split Stars
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Bisection of A001859. See Comments of A135713.
Cf. A000217, A000567, A001082, A002378, A016969, A049450, A174709.
Cf. second n-gonal numbers: A005449, A014105, A147875, A179986, A033954, A062728, A135705.
Cf. numbers of the form n*(d*n+10-d)/2: A008587, A056000, A028347, A140090, A014106, A028895, A186029, A007742, A022267, A033429, A022268, A049452, A186030, A135703, A152734, A139273.
Cf. A000290, A022144, A005563, A001105.
Cf. A056109.
Cf. A003154.

[n*(3*n+2) : n in [0..100]]; // Wesley Ivan Hurt, Sep 24 2017

Table[n*(3n+2), {n,0,60}] (* Harvey P. Dale, May 05 2011 *)
LinearRecurrence[{3,-3,1},{0,5,16},60] (* Harvey P. Dale, Jan 19 2016 *)
CoefficientList[Series[x*(5 + x)/(1 - x)^3,{x, 0, 60}], x] (* Stefano Spezia, Sep 01 2018 *)

a(n)=n*(3*n+2) \\ Charles R Greathouse IV, Nov 20 2012

O.g.f.: x*(5+x)/(1-x)^3. - R. J. Mathar, Jan 07 2008
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), with a(0)=0, a(1)=5, a(2)=16. - Harvey P. Dale, May 06 2011
a(n) = a(n-1) + 6*n - 1 (with a(0)=0). - Vincenzo Librandi, Nov 18 2010
For n > 0, a(n)^3 + (a(n)+1)^3 + ... + (a(n)+n)^3 + 2*A000217(n)^2 = (a(n) + n + 1)^3 + ... + (a(n) + 2n)^3; see also A033954. - Charlie Marion, Dec 08 2007
a(n) = Sum_{i=0..n-1} A016969(i) for n > 0. - Bruno Berselli, Jan 13 2011
a(n) = A174709(6*n+4). - Philippe Deléham, Mar 26 2013
a(n) = A001082(2*n). - Michael Turniansky, Aug 24 2013
Sum_{n>=1} 1/a(n) = (9 + sqrt(3)*Pi - 9*log(3))/12 = 0.3794906245574721941... . - Vaclav Kotesovec, Apr 27 2016
a(n) = A002378(n) + A014105(n). - J. M. Bergot, Apr 24 2018
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/sqrt(12) - 3/4. - Amiram Eldar, Jul 03 2020
E.g.f.: exp(x)*x*(5 + 3*x). - Stefano Spezia, Jun 08 2021
From Leo Tavares, Oct 14 2021: (Start)
a(n) = A000290(n) + 4*A000217(n). See Square Stars illustration.
a(n) = A000567(n+2) - A022144(n+1)
a(n) = A005563(n) + A001105(n).
a(n) = A056109(n) - 1. (End)
From Leo Tavares, Oct 06 2022: (Start)
a(n) = A003154(n+1) - A000567(n+1). See Split Stars illustration.
a(n) = A014105(n) + 2*A000217(n). (End)

cp's OEIS Frontend

A046092 4 times triangular numbers: a(n) = 2*n*(n+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

Maxima

PARI

Python

Formula

A028387 a(n) = n + (n+1)^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

Python

Sage

Formula

Extensions

A005574 Numbers k such that k^2 + 1 is prime.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

PARI

Python

Formula

A028982 Squares and twice squares.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Python

Formula

A016742 Even squares: a(n) = (2*n)^2.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

Maxima

A046092 4 times triangular numbers: a(n) = 2n(n+1).