A016802 -id:A016802 - cp's OEIS frontend

A016754 Odd squares: a(n) = (2n+1)^2. Also centered octagonal numbers.

1, 9, 25, 49, 81, 121, 169, 225, 289, 361, 441, 529, 625, 729, 841, 961, 1089, 1225, 1369, 1521, 1681, 1849, 2025, 2209, 2401, 2601, 2809, 3025, 3249, 3481, 3721, 3969, 4225, 4489, 4761, 5041, 5329, 5625, 5929, 6241, 6561, 6889, 7225, 7569, 7921, 8281, 8649, 9025

Offset: 0

N. J. A. Sloane

The brown rat (rattus norwegicus) breeds very quickly. It can give birth to other rats 7 times a year, starting at the age of three months. The average number of pups is 8. The present sequence gives the total number of rats, when the intervals are 12/7 of a year and a young rat starts having offspring at 24/7 of a year. - Hans Isdahl, Jan 26 2008
Numbers n such that tau(n) is odd where tau(x) denotes the Ramanujan tau function (A000594). - Benoit Cloitre, May 01 2003
If Y is a fixed 2-subset of a (2n+1)-set X then a(n-1) is the number of 3-subsets of X intersecting Y. - Milan Janjic, Oct 21 2007
Binomial transform of [1, 8, 8, 0, 0, 0, ...]; Narayana transform (A001263) of [1, 8, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
All terms of this sequence are of the form 8k+1. For numbers 8k+1 which aren't squares see A138393. Numbers 8k+1 are squares iff k is a triangular number from A000217. And squares have form 4n(n+1)+1. - Artur Jasinski, Mar 27 2008
Sequence arises from reading the line from 1, in the direction 1, 25, ... and the line from 9, in the direction 9, 49, ..., in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008
Equals the triangular numbers convolved with [1, 6, 1, 0, 0, 0, ...]. - Gary W. Adamson & Alexander R. Povolotsky, May 29 2009
First differences: A008590(n) = a(n) - a(n-1) for n>0. - Reinhard Zumkeller, Nov 08 2009
Central terms of the triangle in A176271; cf. A000466, A053755. - Reinhard Zumkeller, Apr 13 2010
Odd numbers with odd abundance. Odd numbers with even abundance are in A088828. Even numbers with odd abundance are in A088827. Even numbers with even abundance are in A088829. - Jaroslav Krizek, May 07 2011
Appear as numerators in the non-simple continued fraction expansion of Pi-3: Pi-3 = K_{k>=1} (1-2*k)^2/6 = 1/(6+9/(6+25/(6+49/(6+...)))), see also the comment in A007509. - Alexander R. Povolotsky, Oct 12 2011
Ulam's spiral (SE spoke). - Robert G. Wilson v, Oct 31 2011
All terms end in 1, 5 or 9. Modulo 100, all terms are among { 1, 9, 21, 25, 29, 41, 49, 61, 69, 81, 89 }. - M. F. Hasler, Mar 19 2012
Right edge of both triangles A214604 and A214661: a(n) = A214604(n+1,n+1) = A214661(n+1,n+1). - Reinhard Zumkeller, Jul 25 2012
Also: Odd numbers which have an odd sum of divisors (= sigma = A000203). - M. F. Hasler, Feb 23 2013
Consider primitive Pythagorean triangles (a^2 + b^2 = c^2, gcd(a, b) = 1) with hypotenuse c (A020882) and respective even leg b (A231100); sequence gives values c-b, sorted with duplicates removed. - K. G. Stier, Nov 04 2013
For n>1 a(n) is twice the area of the irregular quadrilateral created by the points ((n-2)*(n-1),(n-1)*n/2), ((n-1)*n/2,n*(n+1)/2), ((n+1)*(n+2)/2,n*(n+1)/2), and ((n+2)*(n+3)/2,(n+1)*(n+2)/2). - J. M. Bergot, May 27 2014
Number of pairs (x, y) of Z^2, such that max(abs(x), abs(y)) <= n. - Michel Marcus, Nov 28 2014
Except for a(1)=4, the number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 737", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
a(n) is the sum of 2n+1 consecutive numbers, the first of which is n+1. - Ivan N. Ianakiev, Dec 21 2016
a(n) is the number of 2 X 2 matrices with all elements in {0..n} with determinant = 2*permanent. - Indranil Ghosh, Dec 25 2016
Engel expansion of Pi*StruveL_0(1)/2 where StruveL_0(1) is A197037. - Benedict W. J. Irwin, Jun 21 2018
Consider all Pythagorean triples (X,Y,Z=Y+1) ordered by increasing Z; the segments on the hypotenuse {p = a(n)/A001844(n), q = A060300(n)/A001844(n) = A001844(n) - p} and their ratio p/q = a(n)/A060300(n) are irreducible fractions in Q\Z. X values are A005408, Y values are A046092, Z values are A001844. - Ralf Steiner, Feb 25 2020
a(n) is the number of large or small squares that are used to tile primitive squares of type 2 (A344332). - Bernard Schott, Jun 03 2021
Also, positive odd integers with an odd number of odd divisors (for similar sequence with 'even', see A348005). - Bernard Schott, Nov 21 2021
a(n) is the least odd number k = x + y, with 0 < x < y, such that there are n distinct pairs (x,y) for which x*y/k is an integer; for example, a(2) = 25 and the two corresponding pairs are (5,20) and (10,15). The similar sequence with 'even' is A016742 (see Comment of Jan 26 2018). - Bernard Schott, Feb 24 2023
From Peter Bala, Jan 03 2024: (Start)
The sequence terms are the exponents of q in the series expansions of the following infinite products:
1) q*Product_{n >= 1} (1 - q^(16*n))*(1 + q^(8*n)) = q + q^9 + q^25 + q^49 + q^81 + q^121 + q^169 + ....
2) q*Product_{n >= 1} (1 + q^(16*n))*(1 - q^(8*n)) = q - q^9 - q^25 + q^49 + q^81 - q^121 - q^169 + + - - ....
3) q*Product_{n >= 1} (1 - q^(8*n))^3 = q - 3*q^9 + 5*q^25 - 7*q^49 + 9*q^81 - 11*q^121 + 13*q^169 - + ....
4) q*Product_{n >= 1} ( (1 + q^(8*n))*(1 - q^(16*n))/(1 + q^(16*n)) )^3 = q + 3*q^9 - 5*q^25 - 7*q^49 + 9*q^81 + 11*q^121 - 13*q^169 - 15*q^225 + + - - .... (End)

L. Lorentzen and H. Waadeland, Continued Fractions with Applications, North-Holland 1992, p. 586.

Paolo Xausa, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe).
Jeremiah Bartz, Bruce Dearden, and Joel Iiams, Classes of Gap Balancing Numbers, arXiv:1810.07895 [math.NT], 2018.
Bruce C. Berndt and Ken Ono, Ramanujan's unpublished manuscript on the partition and tau functions with proofs and commentary, Séminaire Lotharingien de Combinatoire, B42c (1999), 63 pp.
John Elias, Illustration: 8-fold Square Progression of Ulam's Spiral.
Milan Janjic, Two Enumerative Functions.
Scientific American, Cover of the March 1964 issue.
Amelia Carolina Sparavigna, Groupoids of OEIS A002378 and A016754 Numbers (oblong and odd square numbers), Politecnico di Torino (Italy, 2019).
Leo Tavares, Illustration: Diamond Triangles.
Leo Tavares, Illustration: Diamond Stars.
Eric Weisstein's World of Mathematics, Moore Neighborhood.
R. Yin, J. Mu, and T. Komatsu, The p-Frobenius Number for the Triple of the Generalized Star Numbers, Preprints 2024, 2024072280. See p. 2.
Index entries for sequences related to centered polygonal numbers.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000290, A000384, A001263, A001539, A001844, A003881, A005408, A006752, A014105, A016742, A016802, A016814, A016826, A016838, A033996, A046092, A060300, A138393, A167661, A167700.
Cf. A000447 (partial sums).
Cf. A005917, A344330, A344332.
Cf. A348005, A379481 [= a(A048673(n)-1)].
Partial sums of A022144.
Positions of odd terms in A341528.
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. A014634, A003154.

a016754 n = a016754_list !! n
a016754_list = scanl (+) 1 $ tail a008590_list
-- Reinhard Zumkeller, Apr 02 2012

[n^2: n in [1..100 by 2]]; // Vincenzo Librandi, Jan 03 2017

Range[1, 100, 2]^2 (* Paolo Xausa, Feb 13 2025 *)

A016754(n):=(n+n+1)^2$
makelist(A016754(n),n,0,20); /* Martin Ettl, Nov 12 2012 */

A016754(n)=(n<<1+1)^2 \\ Charles R Greathouse IV, Jun 16 2011, corrected and edited by M. F. Hasler, Apr 11 2023

def A016754(n): return ((n<<1)|1)**2 # Chai Wah Wu, Jul 06 2023

a(n) = 1 + Sum_{i=1..n} 8*i = 1 + 8*A000217(n). - Xavier Acloque, Jan 21 2003; Zak Seidov, May 07 2006; Robert G. Wilson v, Dec 29 2010
O.g.f.: (1+6*x+x^2)/(1-x)^3. - R. J. Mathar, Jan 11 2008
a(n) = 4*n*(n + 1) + 1 = 4*n^2 + 4*n + 1. - Artur Jasinski, Mar 27 2008
a(n) = A061038(2+4n). - Paul Curtz, Oct 26 2008
Sum_{n>=0} 1/a(n) = Pi^2/8 = A111003. - Jaume Oliver Lafont, Mar 07 2009
a(n) = A000290(A005408(n)). - Reinhard Zumkeller, Nov 08 2009
a(n) = a(n-1) + 8*n with n>0, a(0)=1. - Vincenzo Librandi, Aug 01 2010
a(n) = A033951(n) + n. - Reinhard Zumkeller, May 17 2009
a(n) = A033996(n) + 1. - Omar E. Pol, Oct 03 2011
a(n) = (A005408(n))^2. - Zak Seidov, Nov 29 2011
From George F. Johnson, Sep 05 2012: (Start)
a(n+1) = a(n) + 4 + 4*sqrt(a(n)).
a(n-1) = a(n) + 4 - 4*sqrt(a(n)).
a(n+1) = 2*a(n) - a(n-1) + 8.
a(n+1) = 3*a(n) - 3*a(n-1) + a(n-2).
(a(n+1) - a(n-1))/8 = sqrt(a(n)).
a(n+1)*a(n-1) = (a(n)-4)^2.
a(n) = 2*A046092(n) + 1 = 2*A001844(n) - 1 = A046092(n) + A001844(n).
Limit_{n -> oo} a(n)/a(n-1) = 1. (End)
a(n) = binomial(2*n+2,2) + binomial(2*n+1,2). - John Molokach, Jul 12 2013
E.g.f.: (1 + 8*x + 4*x^2)*exp(x). - Ilya Gutkovskiy, May 23 2016
a(n) = A101321(8,n). - R. J. Mathar, Jul 28 2016
Product_{n>=1} A033996(n)/a(n) = Pi/4. - Daniel Suteu, Dec 25 2016
a(n) = A014105(n) + A000384(n+1). - Bruce J. Nicholson, Nov 11 2017
a(n) = A003215(n) + A002378(n). - Klaus Purath, Jun 09 2020
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=0} a(n)/n! = 13*e.
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = 3/e. (End)
Sum_{n>=0} (-1)^n/a(n) = A006752. - Amiram Eldar, Oct 10 2020
From Amiram Eldar, Jan 28 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = cosh(Pi/2).
Product_{n>=1} (1 - 1/a(n)) = Pi/4 (A003881). (End)
From Leo Tavares, Nov 24 2021: (Start)
a(n) = A014634(n) - A002943(n). See Diamond Triangles illustration.
a(n) = A003154(n+1) - A046092(n). See Diamond Stars illustration. (End)
From Peter Bala, Mar 11 2024: (Start)
Sum_{k = 1..n+1} 1/(k*a(k)*a(k-1)) = 1/(9 - 3/(17 - 60/(33 - 315/(57 - ... - n^2*(4*n^2 - 1)/((2*n + 1)^2 + 2*2^2 ))))).
3/2 - 2*log(2) = Sum_{k >= 1} 1/(k*a(k)*a(k-1)) = 1/(9 - 3/(17 - 60/(33 - 315/(57 - ... - n^2*(4*n^2 - 1)/((2*n + 1)^2 + 2*2^2 - ... ))))).
Row 2 of A142992. (End)
From Peter Bala, Mar 26 2024: (Start)
8*a(n) = (2*n + 1)*(a(n+1) - a(n-1)).
Sum_{n >= 0} (-1)^n/(a(n)*a(n+1)) = 1/2 - Pi/8 = 1/(9 + (1*3)/(8 + (3*5)/(8 + ... + (4*n^2 - 1)/(8 + ... )))). For the continued fraction use Lorentzen and Waadeland, p. 586, equation 4.7.9 with n = 1. Cf. A057813. (End)

Additional description from Terrel Trotter, Jr., Apr 06 2002

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

A016814 a(n) = (4*n + 1)^2.

Original entry on oeis.org

1, 25, 81, 169, 289, 441, 625, 841, 1089, 1369, 1681, 2025, 2401, 2809, 3249, 3721, 4225, 4761, 5329, 5929, 6561, 7225, 7921, 8649, 9409, 10201, 11025, 11881, 12769, 13689, 14641, 15625, 16641, 17689, 18769, 19881, 21025, 22201, 23409, 24649, 25921, 27225, 28561, 29929

Offset: 0

N. J. A. Sloane

A bisection of A016754. Sequence arises from reading the line from 1, in the direction 1, 25, ..., in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008

G. C. Greubel, Table of n, a(n) for n = 0..10000 (terms 0..200 from Ivan Panchenko).
Leo Tavares, Illustration: Square trapeziums
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Sequences of the form (m*n+1)^2: A000012 (m=0), A000290 (m=1), A016754 (m=2), A016778 (m-3), this sequence (m=4), A016862 (m=5), A016922 (m=6), A016994 (m=7), A017078 (m=8), A017174 (m=9), A017282 (m=10), A017402 (m=11), A017534 (m=12), A134934 (m=14).
Cf. A001539, A006752, A016742, A016802, A016813.
Cf. A016826, A016838, A017126, A068467, A087197.
Cf. A272399, A014105.

[(4*n+1)^2: n in [0..40]]; // G. C. Greubel, Dec 28 2022

A016814:=n->(4*n+1)^2; seq(A016814(k), k=0..100); # Wesley Ivan Hurt, Nov 02 2013

(4*Range[0,40] +1)^2 (* or *) LinearRecurrence[{3,-3,1}, {1,25,81}, 40] (* Harvey P. Dale, Nov 20 2012 *)
Accumulate[32Range[0, 47] - 8] + 9 (* Alonso del Arte, Aug 19 2017 *)

a(n)=(4*n+1)^2 \\ Charles R Greathouse IV, Oct 07 2015

[(4*n+1)^2 for n in range(41)] # G. C. Greubel, Dec 28 2022

a(n) = a(n-1) + 32*n - 8, n > 0. - Vincenzo Librandi, Dec 15 2010
From George F. Johnson, Sep 28 2012: (Start)
G.f.: (1 + 22*x + 9*x^2)/(1 - x)^3.
a(n+1) = a(n) + 16 + 8*sqrt(a(n)).
a(n+1) = 2*a(n) - a(n-1) + 32 = 3*a(n) - 3*a(n-1) + a(n-2).
a(n-1)*a(n+1) = (a(n) - 16)^2 ; a(n+1) - a(n-1) = 16*sqrt(a(n)).
a(n) = A016754(2*n) = (A016813(n))^2. (End)
Sum_{n>=0} 1/a(n) = G/2 + Pi^2/16, where G is the Catalan constant (A006752). - Amiram Eldar, Jun 28 2020
Product_{n>=1} (1 - 1/a(n)) = 2*Gamma(5/4)^2/sqrt(Pi) = 2 * A068467^2 * A087197. - Amiram Eldar, Feb 01 2021
From G. C. Greubel, Dec 28 2022: (Start)
a(2*n) = A017078(n).
a(2*n+1) = A017126(n).
E.g.f.: (1 + 24*x + 16*x^2)*exp(x). (End)
a(n) = A272399(n+1) - A014105(n). - Leo Tavares, Dec 24 2023

A244630 a(n) = 17*n^2.

Original entry on oeis.org

0, 17, 68, 153, 272, 425, 612, 833, 1088, 1377, 1700, 2057, 2448, 2873, 3332, 3825, 4352, 4913, 5508, 6137, 6800, 7497, 8228, 8993, 9792, 10625, 11492, 12393, 13328, 14297, 15300, 16337, 17408, 18513, 19652, 20825, 22032, 23273, 24548, 25857, 27200, 28577, 29988

Offset: 0

Vincenzo Librandi, Jul 03 2014

First bisection of A195047. - Bruno Berselli, Jul 03 2014

Norms of purely imaginary numbers in Z[sqrt(-17)] (for example, 3*sqrt(-17) has norm 153). - Alonso del Arte, Jun 23 2018

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A008599, A195047.

Cf. similar sequences of the type k*n^2: A000290 (k = 1), A001105 (k = 2), A033428 (k = 3), A016742 (k = 4), A033429 (k = 5), A033581 (k = 6), A033582 (k = 7), A139098 (k = 8), A016766 (k = 9), A033583 (k = 10), A033584 (k = 11), A135453 (k = 12), A152742 (k = 13), A144555 (k = 14), A064761 (k = 15), A016802 (k = 16), this sequence (k = 17), A195321 (k = 18), A244631 (k = 19), A195322 (k = 20), A064762 (k = 21), A195323 (k = 22), A244632 (k = 23), A195824 (k = 24), A016850 (k = 25), A244633 (k = 26), A244634 (k = 27), A064763 (k = 28), A244635 (k = 29), A244636 (k = 30).

List([0..45],n->17*n^2); # Muniru A Asiru, Jun 29 2018

Magma
```
[17*n^2: n in [0..40]];
```

seq(17*n^2,n=0..45); # Muniru A Asiru, Jun 29 2018

Mathematica
```
Table[17 n^2, {n, 0, 40}]
```

a(n)=17*n^2 \\ Charles R Greathouse IV, Oct 07 2015

for (i <- 0 to 50) yield 17 * (i * i) // Alonso del Arte, Jun 29 2018

G.f.: 17*x*(1 + x)/(1 - x)^3. [corrected by Bruno Berselli, Jul 03 2014]
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2.
a(n) = 17*A000290(n). - Omar E. Pol, Jul 03 2014
a(n) = a(-n). - Muniru A Asiru, Jun 29 2018
From Elmo R. Oliveira, Dec 02 2024: (Start)
E.g.f.: 17*x*(1 + x)*exp(x).
a(n) = n*A008599(n) = A195047(2*n). (End)

A195146 Concentric 16-gonal numbers.

Original entry on oeis.org

0, 1, 16, 33, 64, 97, 144, 193, 256, 321, 400, 481, 576, 673, 784, 897, 1024, 1153, 1296, 1441, 1600, 1761, 1936, 2113, 2304, 2497, 2704, 2913, 3136, 3361, 3600, 3841, 4096, 4353, 4624, 4897, 5184, 5473, 5776, 6081, 6400, 6721, 7056, 7393, 7744, 8097, 8464

Offset: 0

Omar E. Pol, Sep 17 2011

Concentric hexadecagonal numbers or concentric hexakaidecagonal numbers.

Sequence found by reading the line from 0, in the direction 0, 16, ..., and the same line from 1, in the direction 1, 33, ..., in the square spiral whose vertices are the generalized decagonal numbers A074377. Main axis, perpendicular to A033996 in the same spiral.

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Index entries for linear recurrences with constant coefficients, signature (2,0,-2,1).

A016802 and A195315 interleaved.
Cf. A032528, A033996, A074377, A077221, A195142, A195143, A195145, A195147, A195148, A195149.
Column 16 of A195040.

[(8*n^2+3*(-1)^n-3)/2: n in [0..50]]; // Vincenzo Librandi, Sep 27 2011

LinearRecurrence[{2, 0, -2, 1}, {0, 1, 16, 33}, 50] (* Amiram Eldar, Jan 16 2023 *)

a(n)=(8*n^2+3*(-1)^n-3)/2 \\ Charles R Greathouse IV, Oct 07 2015

From Vincenzo Librandi, Sep 27 2011: (Start)
a(n) = (8*n^2 + 3*(-1)^n - 3)/2;
a(n) = -a(n-1) + 8*n^2 - 8*n + 1. (End)
G.f. -x*(1+14*x+x^2) / ( (1+x)*(x-1)^3 ). - R. J. Mathar, Sep 18 2011
Sum_{n>=1} 1/a(n) = Pi^2/96 + tan(sqrt(3)*Pi/4)*Pi/(8*sqrt(3)). - Amiram Eldar, Jan 16 2023

A016826 a(n) = (4n + 2)^2.

Original entry on oeis.org

4, 36, 100, 196, 324, 484, 676, 900, 1156, 1444, 1764, 2116, 2500, 2916, 3364, 3844, 4356, 4900, 5476, 6084, 6724, 7396, 8100, 8836, 9604, 10404, 11236, 12100, 12996, 13924, 14884, 15876, 16900, 17956

Offset: 0

N. J. A. Sloane

A bisection of A016742. Sequence arises from reading 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

Harvey P. Dale, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Equals A001539 + 1.

Cf. A000290, A010503, A016742, A016754, A016802, A016814, A016838.

A016826:=n->(4*n + 2)^2; seq(A016826(n), n=0..40); # Wesley Ivan Hurt, Feb 24 2014

(4*Range[0,40]+2)^2 (* or *) LinearRecurrence[{3,-3,1},{4,36,100},40] (* Harvey P. Dale, Nov 24 2011 *)

a(n)=(4*n+2)^2 \\ Charles R Greathouse IV, Oct 07 2015

a(n) = a(n-1) + 32*n (with a(0)=4). - Vincenzo Librandi, Dec 15 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), with a(0)=4, a(1)=36, a(2)=100. - Harvey P. Dale, Nov 24 2011
G.f.: -((4*(x^2+6*x+1))/(x-1)^3). - Harvey P. Dale, Nov 24 2011
a(n) = A000290(A016825(n)). - Wesley Ivan Hurt, Feb 24 2014
From Amiram Eldar, Jun 28 2020: (Start)
Sum_{n>=0} 1/a(n) = Pi^2/32.
Sum_{n>=0} (-1)^n/a(n) = G/4, where G is the Catalan constant (A006752). (End)
From Amiram Eldar, Jan 29 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = cosh(Pi/4).
Product_{n>=0} (1 - 1/a(n)) = 1/sqrt(2) (A010503). (End)

A016838 a(n) = (4n + 3)^2.

Original entry on oeis.org

9, 49, 121, 225, 361, 529, 729, 961, 1225, 1521, 1849, 2209, 2601, 3025, 3481, 3969, 4489, 5041, 5625, 6241, 6889, 7569, 8281, 9025, 9801, 10609, 11449, 12321, 13225, 14161, 15129, 16129, 17161, 18225

Offset: 0

N. J. A. Sloane

If Y is a fixed 2-subset of a (4n+1)-set X then a(n-1) is the number of 3-subsets of X intersecting Y. - Milan Janjic, Oct 21 2007
A bisection of A016754. Sequence arises from reading the line from 9, in the direction 9, 49, ... in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008
Using (n,n+1) to generate a Pythagorean triangle with sides of lengths xJ. M. Bergot, Jul 17 2013

Vincenzo Librandi, Table of n, a(n) for n = 0..860
Milan Janjic, Two Enumerative Functions
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000290, A001539, A016742, A016754, A016802, A016814, A016826, A068465, A087197.

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

A016838:=n->(4*n + 3)^2; seq(A016838(n), n=0..50); # Wesley Ivan Hurt, Feb 24 2014

Table[(4*n + 3)^2, {n, 0, 40}] (* Vaclav Kotesovec, Nov 14 2017 *)

a(n)=(4*n+3)^2 \\ Charles R Greathouse IV, Jun 17 2017

Denominators of first differences Zeta[2,(4n-1)/4]-Zeta[2,(4(n+1)-1)/4]. - Artur Jasinski, Mar 03 2010
From George F. Johnson, Oct 03 2012: (Start)
G.f.: (9+22*x+x^2)/(1-x)^3.
a(n+1) = a(n) + 16 + 8*sqrt(a(n)).
a(n+1) = 2*a(n) - a(n-1) + 32 = 3*a(n) - 3*a(n-1) + a(n-2).
a(n-1)*a(n+1) = (a(n)-16)^2; a(n+1) - a(n-1) = 16*sqrt(a(n)).
a(n) = A016754(2*n+1) = (A004767(n))^2.
(End)
Sum_{n>=0} 1/a(n) = Pi^2/16 - G/2, where G is the Catalan constant (A006752). - Amiram Eldar, Jun 28 2020
Product_{n>=0} (1 - 1/a(n)) = Gamma(3/4)^2/sqrt(Pi) = A068465^2 * A087197. - Amiram Eldar, Feb 01 2021

A188824 T(n,k) = number of n X k binary arrays without the pattern 0 1 diagonally or antidiagonally.

Original entry on oeis.org

2, 4, 4, 8, 9, 8, 16, 25, 16, 16, 32, 64, 48, 25, 32, 64, 169, 144, 81, 36, 64, 128, 441, 432, 256, 120, 49, 128, 256, 1156, 1296, 841, 400, 169, 64, 256, 512, 3025, 3888, 2704, 1360, 576, 224, 81, 512, 1024, 7921, 11664, 8836, 4624, 2025, 784, 289, 100, 1024, 2048

Offset: 1

R. H. Hardin, Apr 11 2011

Table starts
....2...4...8...16...32....64...128....256.....512....1024.....2048.....4096
....4...9..25...64..169...441..1156...3025....7921...20736....54289...142129
....8..16..48..144..432..1296..3888..11664...34992..104976...314928...944784
...16..25..81..256..841..2704..8836..28561...93025..301401...980100..3179089
...32..36.120..400.1360..4624.15776..53824..183744..627264..2141568..7311616
...64..49.169..576.2025..7056.24964..87616..310249.1092025..3865156.13623481
..128..64.224..784.2800.10000.36000.129600..468000.1690000..6110000.22090000
..256..81.289.1024.3721.13456.49284.179776..660969.2421136..8916196.32729841
..512.100.360.1296.4752.17424.64416.238144..884256.3283344.12220128.45481536
.1024.121.441.1600.5929.21904.81796.304704.1142761.4276624.16080100.60341824

			Some solutions for 5X3
..1..1..1....1..1..1....1..1..1....1..1..1....1..1..1....0..1..0....1..1..1
..1..1..1....0..1..1....1..0..1....1..1..1....1..1..1....1..0..1....1..1..1
..1..1..1....0..0..1....0..1..0....1..1..1....0..1..0....0..1..0....0..1..1
..1..1..1....0..0..0....0..0..0....0..0..0....0..0..0....1..0..1....1..0..1
..1..0..0....0..0..0....0..0..0....0..0..0....0..0..0....0..0..0....0..1..0

R. H. Hardin, Table of n, a(n) for n = 1..1627

Column 2 is A000290(n+1).
Column 4 is A016802.
Row 2 is A007598(n+2).

Empirical: a(n,1) = 2^n
Empirical: a(n,2) = n^2 + 2*n + 1
Empirical: a(n,3) = 2*a(n-1,3) - 2*a(n-3,3) + a(n-4,3)
Empirical: a(n,4) = 16*n^2
Empirical: a(n,5) = 2*a(n-1,5) - 2*a(n-3,5) + a(n-4,5) for n>5
Empirical: a(n,6) = 256*n^2 - 384*n + 144 for n>2
Empirical: a(n,7) = 2*a(n-1,7) - 2*a(n-3,7) + a(n-4,7) for n>7
Empirical: a(n,8) = 4096*n^2 - 11264*n + 7744 for n>4

A195321 a(n) = 18*n^2.

Original entry on oeis.org

0, 18, 72, 162, 288, 450, 648, 882, 1152, 1458, 1800, 2178, 2592, 3042, 3528, 4050, 4608, 5202, 5832, 6498, 7200, 7938, 8712, 9522, 10368, 11250, 12168, 13122, 14112, 15138, 16200, 17298, 18432, 19602, 20808, 22050, 23328, 24642, 25992, 27378, 28800, 30258, 31752

Offset: 0

Omar E. Pol, Sep 16 2011

Sequence found by reading the line from 0, in the direction 0, 18, ..., in the square spiral whose vertices are the generalized hendecagonal numbers A195160. Semi-axis opposite to A195316 in the same spiral.
Area of a square with diagonal 6n. - Wesley Ivan Hurt, Jun 19 2014
Number of identical tessellation tiles that are composed of 48 equilateral edge joined triangles that can be formed into a order n hexagon. The example tiles shown in the link below are tessellated with eight sphinx tiles. See A291582. - Craig Knecht, Sep 02 2017

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Craig Knecht, Hexagon tessellation.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Bisection of A195147.
Cf. A016802, A033581, A033583, A135453, A139098, A144555, A195322, A195323.
Cf. A000290, A001105, A008600, A016766, A033428, A195160, A195316, A291582.

[18*n^2:n in [0..40]]; // Vincenzo Librandi, Sep 20 2011

A195321:=n->18*n^2; seq(A195321(n), n=0..50); # Wesley Ivan Hurt, Jun 19 2014

18 Range[0, 50]^2 (* or *) CoefficientList[Series[18 x*(1 + x)/(1 - x)^3, {x, 0, 30}], x] (* Wesley Ivan Hurt, Jun 20 2014 *)
LinearRecurrence[{3,-3,1},{0,18,72},50] (* Harvey P. Dale, Mar 26 2023 *)

a(n)=18*n^2 \\ Charles R Greathouse IV, Jun 17 2017

a(n) = 18*A000290(n) = 9*A001105(n) = 6*A033428(n) = 3*A033581(n) = 2*A016766(n).
G.f.: 18*x*(1+x)/(1-x)^3. - Wesley Ivan Hurt, Jun 20 2014
From Elmo R. Oliveira, Dec 01 2024: (Start)
E.g.f.: 18*x*(1 + x)*exp(x).
a(n) = n*A008600(n) = A195147(2*n).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2. (End)

A195322 a(n) = 20*n^2.

Original entry on oeis.org

0, 20, 80, 180, 320, 500, 720, 980, 1280, 1620, 2000, 2420, 2880, 3380, 3920, 4500, 5120, 5780, 6480, 7220, 8000, 8820, 9680, 10580, 11520, 12500, 13520, 14580, 15680, 16820, 18000, 19220, 20480, 21780, 23120, 24500, 25920, 27380, 28880, 30420, 32000, 33620, 35280

Offset: 0

Omar E. Pol, Sep 16 2011

Sequence found by reading the line from 0, in the direction 0, 20, ..., in the square spiral whose vertices are the generalized dodecagonal numbers A195162. Semiaxis opposite to A195317 in the same spiral.
a(n) is the sum of all the integers less than 10*n which are not multiple of 2 or 5. a(2) = (1 + 3 + 7 + 9) + (11 + 13 + 17 + 19) = 20 + 60 = 80 = 20 * 2^2. (Link Crux Mathematicorum). - Bernard Schott, May 15 2017
Number of terms less than 10^k (k=0, 1, 2, ...): 1, 1, 3, 8, 23, 71, 224, 708, 2237, 7072, 22361, 70711, ... - Muniru A Asiru, Feb 01 2018

			From _Muniru A Asiru_, Feb 01 2018: (Start)
n=0, a(0) = 20*0^2 = 0.
n=1, a(1) = 20*1^2 = 20.
n=1, a(2) = 20*2^2 = 80.
n=1, a(3) = 20*3^2 = 180.
n=1, a(4) = 20*4^2 = 320.
...
(End)

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Léo Sauvé, Problem 53, Crux Mathematicorum, Vol. 1, Nov. 1975, page 88.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Bisection of A195148.
Cf. A016802, A033581, A033583, A135453, A139098, A144555, A195321, A195323.
Cf. A000290, A001105, A005899, A010014, A016742, A033429, A195162, A195317.
Cf. A008602.

List([0..10^3],n->20*n^2); # Muniru A Asiru, Feb 01 2018

[20*n^2: n in [0..40]]; // Vincenzo Librandi, Sep 20 2011

a := n -> 20*n^2; seq(a(n), n=0..10^3); # Muniru A Asiru, Feb 01 2018

20 Range[0, 40]^2 (* or *) LinearRecurrence[{3, -3, 1}, {0, 20, 80}, 50] (* Harvey P. Dale, Jan 18 2013 *)

a(n) = 20*n^2 \\ Charles R Greathouse IV, Oct 07 2015

a(n) = 20*A000290(n) = 10*A001105(n) = 5*A016742(n) = 4*A033429(n) = 2*A033583(n).
a(0)=0, a(1)=20, a(2)=80; for n > 2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Jan 18 2013
a(n) = A010014(n) - A005899(n) for n > 0. - R. J. Cano, Sep 29 2015
From Elmo R. Oliveira, Nov 30 2024: (Start)
G.f.: 20*x*(1 + x)/(1-x)^3.
E.g.f.: 20*x*(1 + x)*exp(x).
a(n) = n*A008602(n) = A195148(2*n). (End)

cp's OEIS Frontend

A016754 Odd squares: a(n) = (2n+1)^2. Also centered octagonal numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

Maxima

PARI

Python

Formula

Extensions

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

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

Maxima

PARI

Formula

Extensions

A016814 a(n) = (4*n + 1)^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A244630 a(n) = 17*n^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Scala

Formula

A195146 Concentric 16-gonal numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A016826 a(n) = (4n + 2)^2.

Views

Author