A061037 -id:A061037 - cp's OEIS frontend

A247004 Denominator of (n+4)/gcd(n, 4)^2, a 16-periodic sequence that associates A061037 with A106617.

4, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1, 4, 1, 2, 1, 2

Offset: 0

Jean-François Alcover and Paul Curtz, Sep 09 2014

This sequence may also be defined as the denominators of A061037(n+3)/(n+1), or also as A060819 / A109008.

One can notice that the analog numerators [numerators of (n+4)/gcd(n, 4)^2] are A106617 left-shifted 4 places.

			Fractions begin:
1/4,  5,  3/2,  7, 1/2,  9,  5/2, 11, 3/4, 13,  7/2, 15, 1, 17,  9/2, 19,
5/4, 21, 11/2, 23, 3/2, 25, 13/2, 27, 7/4, 29, 15/2, 31, 2, 33, 17/2, 35,
...
Numerators begin:
1,  5,  3,  7, 1,  9,  5, 11, 3, 13,  7, 15, 1, 17,  9, 19,
5, 21, 11, 23, 3, 25, 13, 27, 7, 29, 15, 31, 2, 33, 17, 35,
...
Periodic part = [4, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1];

Antti Karttunen, Table of n, a(n) for n = 0..65537
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1).

Cf. A060819, A061037, A106617, A109008, A188134, A213268.

[Denominator((n+4)/Gcd(n,4)^2): n in [0..100]]; // G. C. Greubel, Aug 05 2018

a[n_] := (n+4)/GCD[n, 4]^2 // Denominator;  Table[a[n], {n, 0, 100}]
(* or: *)
Table[{1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 1, 1, 1, 2, 1, 4}[[Mod[n, 16, 1]]], {n, 0, 100}]

for(n=0,100, print1(denominator((n+4)/gcd(n,4)^2), ", ")) \\ G. C. Greubel, Aug 05 2018

A247004 = A060819 / A109008.
(n+4) / gcd(n, 4)^2 = A188134(n+4) / 4. - Michael Somos, Sep 12 2014
a(n) = a(n+16) = a(-n), a(2*n + 1) = 1 for all n in Z. - Michael Somos, Sep 13 2014

A144392 Inverse binomial transform of A061037 (read as with offset 0).

Original entry on oeis.org

0, 5, -7, 27, -84, 240, -630, 1554, -3702, 8694, -20310, 47190, -108684, 247572, -557928, 1246320, -2765760, 6106944, -13427424, 29404704, -64138560, 139366080, -301749888, 651236352, -1401444864, 3007987200, -6440578560, 13759041024

Offset: 0

Paul Curtz, Oct 02 2008

Empirical g.f.: x*(5 +53*x +273*x^2 +878*x^3 +1934*x^4 +3042*x^5 +3462*x^6 +2814*x^7 +1554*x^8 +518*x^9 +74*x^10) / ((1 +2*x)^3*(1 +2*x +2*x^2)^3). - Colin Barker, Mar 17 2016

Extended by R. J. Mathar, Dec 02 2008

A002378 Oblong (or promic, pronic, or heteromecic) numbers: a(n) = n*(n+1).

Original entry on oeis.org

0, 2, 6, 12, 20, 30, 42, 56, 72, 90, 110, 132, 156, 182, 210, 240, 272, 306, 342, 380, 420, 462, 506, 552, 600, 650, 702, 756, 812, 870, 930, 992, 1056, 1122, 1190, 1260, 1332, 1406, 1482, 1560, 1640, 1722, 1806, 1892, 1980, 2070, 2162, 2256, 2352, 2450, 2550

Offset: 0

N. J. A. Sloane

4*a(n) + 1 are the odd squares A016754(n).
The word "pronic" (used by Dickson) is incorrect. - Michael Somos
According to the 2nd edition of Webster, the correct word is "promic". - R. K. Guy
a(n) is the number of minimal vectors in the root lattice A_n (see Conway and Sloane, p. 109).
Let M_n denote the n X n matrix M_n(i, j) = (i + j); then the characteristic polynomial of M_n is x^(n-2) * (x^2 - a(n)*x - A002415(n)). - Benoit Cloitre, Nov 09 2002
The greatest LCM of all pairs (j, k) for j < k <= n for n > 1. - Robert G. Wilson v, Jun 19 2004
First differences are a(n+1) - a(n) = 2*n + 2 = 2, 4, 6, ... (while first differences of the squares are (n+1)^2 - n^2 = 2*n + 1 = 1, 3, 5, ...). - Alexandre Wajnberg, Dec 29 2005
25 appended to these numbers corresponds to squares of numbers ending in 5 (i.e., to squares of A017329). - Lekraj Beedassy, Mar 24 2006
A rapid (mental) multiplication/factorization technique -- a generalization of Lekraj Beedassy's comment: For all bases b >= 2 and positive integers n, c, d, k with c + d = b^k, we have (n*b^k + c)*(n*b^k + d) = a(n)*b^(2*k) + c*d. Thus the last 2*k base-b digits of the product are exactly those of c*d -- including leading 0(s) as necessary -- with the preceding base-b digit(s) the same as a(n)'s. Examples: In decimal, 113*117 = 13221 (as n = 11, b = 10 = 3 + 7, k = 1, 3*7 = 21, and a(11) = 132); in octal, 61*67 = 5207 (52 is a(6) in octal). In particular, for even b = 2*m (m > 0) and c = d = m, such a product is a square of this type. Decimal factoring: 5609 is immediately seen to be 71*79. Likewise, 120099 = 301*399 (k = 2 here) and 99990000001996 = 9999002*9999998 (k = 3). - Rick L. Shepherd, Jul 24 2021
Number of circular binary words of length n + 1 having exactly one occurrence of 01. Example: a(2) = 6 because we have 001, 010, 011, 100, 101 and 110. Column 1 of A119462. - Emeric Deutsch, May 21 2006
The sequence of iterated square roots sqrt(N + sqrt(N + ...)) has for N = 1, 2, ... the limit (1 + sqrt(1 + 4*N))/2. For N = a(n) this limit is n + 1, n = 1, 2, .... For all other numbers N, N >= 1, this limit is not a natural number. Examples: n = 1, a(1) = 2: sqrt(2 + sqrt(2 + ...)) = 1 + 1 = 2; n = 2, a(2) = 6: sqrt(6 + sqrt(6 + ...)) = 1 + 2 = 3. - Wolfdieter Lang, May 05 2006
Nonsquare integers m divisible by ceiling(sqrt(m)), except for m = 0. - Max Alekseyev, Nov 27 2006
The number of off-diagonal elements of an (n + 1) X (n + 1) matrix. - Artur Jasinski, Jan 11 2007
a(n) is equal to the number of functions f:{1, 2} -> {1, 2, ..., n + 1} such that for a fixed x in {1, 2} and a fixed y in {1, 2, ..., n + 1} we have f(x) <> y. - Aleksandar M. Janjic and Milan Janjic, Mar 13 2007
Numbers m >= 0 such that round(sqrt(m+1)) - round(sqrt(m)) = 1. - Hieronymus Fischer, Aug 06 2007
Numbers m >= 0 such that ceiling(2*sqrt(m+1)) - 1 = 1 + floor(2*sqrt(m)). - Hieronymus Fischer, Aug 06 2007
Numbers m >= 0 such that fract(sqrt(m+1)) > 1/2 and fract(sqrt(m)) < 1/2 where fract(x) is the fractional part (fract(x) = x - floor(x), x >= 0). - Hieronymus Fischer, Aug 06 2007
X values of solutions to the equation 4*X^3 + X^2 = Y^2. To find Y values: b(n) = n(n+1)(2n+1). - Mohamed Bouhamida, Nov 06 2007
Nonvanishing diagonal of A132792, the infinitesimal Lah matrix, so "generalized factorials" composed of a(n) are given by the elements of the Lah matrix, unsigned A111596, e.g., a(1)*a(2)*a(3) / 3! = -A111596(4,1) = 24. - Tom Copeland, Nov 20 2007
If Y is a 2-subset of an n-set X then, for n >= 2, a(n-2) is the number of 2-subsets and 3-subsets of X having exactly one element in common with Y. - Milan Janjic, Dec 28 2007
a(n) coincides with the vertex of a parabola of even width in the Redheffer matrix, directed toward zero. An integer p is prime if and only if for all integer k, the parabola y = kx - x^2 has no integer solution with 1 < x < k when y = p; a(n) corresponds to odd k. - Reikku Kulon, Nov 30 2008
The third differences of certain values of the hypergeometric function 3F2 lead to the squares of the oblong numbers i.e., 3F2([1, n + 1, n + 1], [n + 2, n + 2], z = 1) - 3*3F2([1, n + 2, n + 2], [n + 3, n + 3], z = 1) + 3*3F2([1, n + 3, n + 3], [n + 4, n + 4], z = 1) - 3F2([1, n + 4, n + 4], [n + 5, n + 5], z = 1) = (1/((n+2)*(n+3)))^2 for n = -1, 0, 1, 2, ... . See also A162990. - Johannes W. Meijer, Jul 21 2009
Generalized factorials, [a.(n!)] = a(n)*a(n-1)*...*a(0) = A010790(n), with a(0) = 1 are related to A001263. - Tom Copeland, Sep 21 2011
For n > 1, a(n) is the number of functions f:{1, 2} -> {1, ..., n + 2} where f(1) > 1 and f(2) > 2. Note that there are n + 1 possible values for f(1) and n possible values for f(2). For example, a(3) = 12 since there are 12 functions f from {1, 2} to {1, 2, 3, 4, 5} with f(1) > 1 and f(2) > 2. - Dennis P. Walsh, Dec 24 2011
a(n) gives the number of (n + 1) X (n + 1) symmetric (0, 1)-matrices containing two ones (see [Cameron]). - L. Edson Jeffery, Feb 18 2012
a(n) is the number of positions of a domino in a rectangled triangular board with both legs equal to n + 1. - César Eliud Lozada, Sep 26 2012
a(n) is the number of ordered pairs (x, y) in [n+2] X [n+2] with |x-y| > 1. - Dennis P. Walsh, Nov 27 2012
a(n) is the number of injective functions from {1, 2} into {1, 2, ..., n + 1}. - Dennis P. Walsh, Nov 27 2012
a(n) is the sum of the positive differences of the partition parts of 2n + 2 into exactly two parts (see example). - Wesley Ivan Hurt, Jun 02 2013
a(n)/a(n-1) is asymptotic to e^(2/n). - Richard R. Forberg, Jun 22 2013
Number of positive roots in the root system of type D_{n + 1} (for n > 2). - Tom Edgar, Nov 05 2013
Number of roots in the root system of type A_n (for n > 0). - Tom Edgar, Nov 05 2013
From Felix P. Muga II, Mar 18 2014: (Start)
a(m), for m >= 1, are the only positive integer values t for which the Binet-de Moivre formula for the recurrence b(n) = b(n-1) + t*b(n-2) with b(0) = 0 and b(1) = 1 has a root of a square. PROOF (as suggested by Wolfdieter Lang, Mar 26 2014): The sqrt(1 + 4t) appearing in the zeros r1 and r2 of the characteristic equation is (a positive) integer for positive integer t precisely if 4t + 1 = (2m + 1)^2, that is t = a(m), m >= 1. Thus, the characteristic roots are integers: r1 = m + 1 and r2 = -m.
Let m > 1 be an integer. If b(n) = b(n-1) + a(m)*b(n-2), n >= 2, b(0) = 0, b(1) = 1, then lim_{n->oo} b(n+1)/b(n) = m + 1. (End)
Cf. A130534 for relations to colored forests, disposition of flags on flagpoles, and colorings of the vertices (chromatic polynomial) of the complete graphs (here simply K_2). - Tom Copeland, Apr 05 2014
The set of integers k for which k + sqrt(k + sqrt(k + sqrt(k + sqrt(k + ...) ... is an integer. - Leslie Koller, Apr 11 2014
a(n-1) is the largest number k such that (n*k)/(n+k) is an integer. - Derek Orr, May 22 2014
Number of ways to place a domino and a singleton on a strip of length n - 2. - Ralf Stephan, Jun 09 2014
With offset 1, this appears to give the maximal number of crossings between n nonconcentric circles of equal radius. - Felix Fröhlich, Jul 14 2014
For n > 1, the harmonic mean of the n values a(1) to a(n) is n + 1. The lowest infinite sequence of increasing positive integers whose cumulative harmonic mean is integral. - Ian Duff, Feb 01 2015
a(n) is the maximum number of queens of one color that can coexist without attacking one queen of the opponent's color on an (n+2) X (n+2) chessboard. The lone queen can be placed in any position on the perimeter of the board. - Bob Selcoe, Feb 07 2015
With a(0) = 1, a(n-1) is the smallest positive number not in the sequence such that Sum_{i = 1..n} 1/a(i-1) has a denominator equal to n. - Derek Orr, Jun 17 2015
The positive members of this sequence are a proper subsequence of the so-called 1-happy couple products A007969. See the W. Lang link there, eq. (4), with Y_0 = 1, with a table at the end. - Wolfdieter Lang, Sep 19 2015
For n > 0, a(n) is the reciprocal of the area bounded above by y = x^(n-1) and below by y = x^n for x in the interval [0, 1]. Summing all such areas visually demonstrates the formula below giving Sum_{n >= 1} 1/a(n) = 1. - Rick L. Shepherd, Oct 26 2015
It appears that, except for a(0) = 0, this is the set of positive integers n such that x*floor(x) = n has no solution. (For example, to get 3, take x = -3/2.) - Melvin Peralta, Apr 14 2016
If two independent real random variables, x and y, are distributed according to the same exponential distribution: pdf(x) = lambda * exp(-lambda * x), lambda > 0, then the probability that n - 1 <= x/y < n is given by 1/a(n). - Andres Cicuttin, Dec 03 2016
a(n) is equal to the sum of all possible differences between n different pairs of consecutive odd numbers (see example). - Miquel Cerda, Dec 04 2016
a(n+1) is the dimension of the space of vector fields in the plane with polynomial coefficients up to order n. - Martin Licht, Dec 04 2016
It appears that a(n) + 3 is the area of the largest possible pond in a square (A268311). - Craig Knecht, May 04 2017
Also the number of 3-cycles in the (n+3)-triangular honeycomb acute knight graph. - Eric W. Weisstein, Jul 27 2017
Also the Wiener index of the (n+2)-wheel graph. - Eric W. Weisstein, Sep 08 2017
The left edge of a Floyd's triangle that consists of even numbers: 0; 2, 4; 6, 8, 10; 12, 14, 16, 18; 20, 22, 24, 26, 28; ... giving 0, 2, 6, 12, 20, ... The right edge generates A028552. - Waldemar Puszkarz, Feb 02 2018
a(n+1) is the order of rowmotion on a poset obtained by adjoining a unique minimal (or maximal) element to a disjoint union of at least two chains of n elements. - Nick Mayers, Jun 01 2018
From Juhani Heino, Feb 05 2019: (Start)
For n > 0, 1/a(n) = n/(n+1) - (n-1)/n.
For example, 1/6 = 2/3 - 1/2; 1/12 = 3/4 - 2/3.
Corollary of this:
Take 1/2 pill.
Next day, take 1/6 pill. 1/2 + 1/6 = 2/3, so your daily average is 1/3.
Next day, take 1/12 pill. 2/3 + 1/12 = 3/4, so your daily average is 1/4.
And so on. (End)
From Bernard Schott, May 22 2020: (Start)
For an oblong number m >= 6 there exists a Euclidean division m = d*q + r with q < r < d which are in geometric progression, in this order, with a common integer ratio b. For b >= 2 and q >= 1, the Euclidean division is m = qb*(qb+1) = qb^2 * q + qb where (q, qb, qb^2) are in geometric progression.
Some examples with distinct ratios and quotients:
   6 | 4            30 | 25          42 | 18
     -----             -----            -----
   2 | 1      ,      5 |  1   ,       6 |  2    ,
and also:
  42 | 12          420 | 100
     -----             -----
   6 |  3     ,     20 |   4  .
Some oblong numbers also satisfy a Euclidean division m = d*q + r with q < r < d that are in geometric progression in this order but with a common noninteger ratio b > 1 (see A335064). (End)
For n >= 1, the continued fraction expansion of sqrt(a(n)) is [n; {2, 2n}]. For n=1, this collapses to [1; {2}]. - Magus K. Chu, Sep 09 2022
a(n-2) is the maximum irregularity over all trees with n vertices. The extremal graphs are stars. (The irregularity of a graph is the sum of the differences between the degrees over all edges of the graph.) - Allan Bickle, May 29 2023
For n > 0, number of diagonals in a regular 2*(n+1)-gon that are not parallel to any edge (cf. A367204). - Paolo Xausa, Mar 30 2024
a(n-1) is the maximum Zagreb index over all trees with n vertices. The extremal graphs are stars. (The Zagreb index of a graph is the sum of the squares of the degrees over all vertices of the graph.) - Allan Bickle, Apr 11 2024
For n >= 1, a(n) is the determinant of the distance matrix of a cycle graph on 2*n + 1 vertices (if the length of the cycle is even such a determinant is zero). - Miquel A. Fiol, Aug 20 2024
For n > 1, the continued fraction expansion of sqrt(16*a(n)) is [2n+1; {1, 2n-1, 1, 8n+2}]. - Magus K. Chu, Nov 20 2024
For n>=2, a(n) is the number of faces on a n+1-zone rhombic zonohedron. Each pair of a collection of great circles on a sphere intersects at two points, so there are 2*binomial(n+1,2) intersections. The dual of the implied polyhedron is a rhombic zonohedron, its faces corresponding to the intersections. - Shel Kaphan, Aug 12 2025

			a(3) = 12, since 2(3)+2 = 8 has 4 partitions with exactly two parts: (7,1), (6,2), (5,3), (4,4). Taking the positive differences of the parts in each partition and adding, we get: 6 + 4 + 2 + 0 = 12. - _Wesley Ivan Hurt_, Jun 02 2013
G.f. = 2*x + 6*x^2 + 12*x^3 + 20*x^4 + 30*x^5 + 42*x^6 + 56*x^7 + ... - _Michael Somos_, May 22 2014
From _Miquel Cerda_, Dec 04 2016: (Start)
a(1) = 2, since 45-43 = 2;
a(2) = 6, since 47-45 = 2 and 47-43 = 4, then 2+4 = 6;
a(3) = 12, since 49-47 = 2, 49-45 = 4, and 49-43 = 6, then 2+4+6 = 12. (End)

W. W. Berman and D. E. Smith, A Brief History of Mathematics, 1910, Open Court, page 67.
J. H. Conway and R. K. Guy, The Book of Numbers, 1996, p. 34.
J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag.
L. E. Dickson, History of the Theory of Numbers, Vol. 1: Divisibility and Primality. New York: Chelsea, p. 357, 1952.
L. E. Dickson, History of the Theory of Numbers, Vol. 2: Diophantine Analysis. New York: Chelsea, pp. 6, 232-233, 350 and 407, 1952.
H. Eves, An Introduction to the History of Mathematics, revised, Holt, Rinehart and Winston, 1964, page 72.
Nicomachus of Gerasa, Introduction to Arithmetic, translation by Martin Luther D'Ooge, Ann Arbor, University of Michigan Press, 1938, p. 254.
Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.6 Figurate Numbers, p. 291.
Granino A. Korn and Theresa M. Korn, Mathematical Handbook for Scientists and Engineers, McGraw-Hill Book Company, New York (1968), pp. 980-981.
C. S. Ogilvy and J. T. Anderson, Excursions in Number Theory, Oxford University Press, 1966, pp. 61-62.
Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, pages 54-55.
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).
F. J. Swetz, From Five Fingers to Infinity, Open Court, 1994, p. 219.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 2-6.

T. D. Noe, Table of n, a(n) for n = 0..1000
D. Applegate, M. LeBrun and N. J. A. Sloane, Dismal Arithmetic, J. Int. Seq. 14 (2011) # 11.9.8.
R. Bapat, S. J. Kirkland, and M. Neumann, On distance matrices and Laplacians, Linear Algebra Appl. 401 (2005), 193-209.
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.
Allan Bickle, Zagreb Indices of Maximal k-degenerate Graphs, Australas. J. Combin. 89 1 (2024) 167-178.
Alin Bostan, Frédéric Chyzak, and Vincent Pilaud, Refined product formulas for Tamari intervals, arXiv:2303.10986 [math.CO], 2023.
H. Bottomley, Illustration of initial terms of A000217, A002378.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 410
P. Cameron, T. Prellberg and D. Stark, Asymptotics for incidence matrix classes, Electron. J. Combin. 13 (2006), #R85, p. 11.
S. Crowley, Two new zeta constants: fractal string and hypergeometric aspects of the Riemann zeta function, viXra:1202.0066, 2012.
J. Estes and B. Wei, Sharp bounds of the Zagreb indices of k-trees, J Comb Optim 27 (2014), 271-291.
Ricardo Gómez Aíza, Trees with flowers: A catalog of integer partition and integer composition trees with their asymptotic analysis, arXiv:2402.16111 [math.CO], 2024. See p. 23.
I. Gutman and K. Das, The first Zagreb index 30 years after, MATCH Commun. Math. Comput. Chem. 50 (2004), 83-92.
Milan Janjic, Enumerative Formulas for Some Functions on Finite Sets.
L. B. W. Jolley, Summation of Series, Dover, 1961.
Refik Keskin and Olcay Karaatli, Some New Properties of Balancing Numbers and Square Triangular Numbers, Journal of Integer Sequences, Vol. 15 (2012), Article #12.1.4.
Craig Knecht, Largest pond by area in a square.
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See p. 27.
Enrique Navarrete and Daniel Orellana, Finding Prime Numbers as Fixed Points of Sequences, arXiv:1907.10023 [math.NT], 2019.
László Németh, The trinomial transform triangle, J. Int. Seqs., Vol. 21 (2018), Article #18.7.3. Also arXiv:1807.07109 [math.NT], 2018.
Lee Melvin Peralta, Solutions to the Equation [x]x = n, The Mathematics Teacher, Vol. 111, No. 2 (October 2017), pp. 150-154.
Aleksandar Petojević, A Note about the Pochhammer Symbol, Mathematica Moravica, Vol. 12-1 (2008), 37-42.
A. Petojevic and N. Dapic, The vAm(a,b,c;z) function, Preprint 2013.
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.
John D. Roth, David A. Garren, and R. Clark Robertson, Integer Carrier Frequency Offset Estimation in OFDM with Zadoff-Chu Sequences, IEEE Int'l Conference on Acoustics, Speech and Signal Processing (ICASSP 2021) 4850-4854.
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014-2015.
Michelle Rudolph-Lilith, On the Product Representation of Number Sequences, with Application to the Fibonacci Family, arXiv preprint arXiv:1508.07894 [math.NT], 2015.
Amelia Carolina Sparavigna, Groupoids of OEIS A002378 and A016754 Numbers (oblong and odd square numbers), Politecnico di Torino (Italy, 2019).
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].
J. Striker and N. Williams, Promotion and Rowmotion , arXiv preprint arXiv:1108.1172 [math.CO], 2011-2012.
D. Suprijanto and Rusliansyah, Observation on Sums of Powers of Integers Divisible by Four, Applied Mathematical Sciences, Vol. 8, 2014, no. 45, 2219 - 2226.
Leo Tavares, Square illustration
R. Tijdeman, Some applications of Diophantine approximation, pp. 261-284 of Surveys in Number Theory (Urbana, May 21, 2000), ed. M. A. Bennett et al., Peters, 2003.
G. Villemin's Almanach of Numbers, Nombres Proniques.
Eric Weisstein's World of Mathematics, Crown Graph.
Eric Weisstein's World of Mathematics, Graph Cycle.
Eric Weisstein's World of Mathematics, Pronic Number.
Eric Weisstein's World of Mathematics, Leibniz Harmonic Triangle.
Eric Weisstein's World of Mathematics, Wheel Graph.
Eric Weisstein's World of Mathematics, Wiener Index.
Wikipedia, Pronic number.
Wolfram Research, Hypergeometric Function 3F2, The Wolfram Functions site.
Index entries for "core" sequences
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Partial sums of A005843 (even numbers). Twice triangular numbers (A000217).
1/beta(n, 2) in A061928.
Cf. A035106, A087811, A119462, A127235, A049598, A124080, A033996, A028896, A046092, A000217, A005563, A046092, A001082, A059300, A059297, A059298, A166373, A002943 (bisection), A002939 (bisection), A078358 (complement).
A036689 and A036690 are subsequences. Cf. numbers of the form n*(n*k-k+4)/2 listed in A226488. - Bruno Berselli, Jun 10 2013
Row n=2 of A185651.
Cf. A007745, A169810, A213541, A005369 (characteristic function).
Cf. A281026. - Bruno Berselli, Jan 16 2017
Cf. A045943 (4-cycles in triangular honeycomb acute knight graph), A028896 (5-cycles), A152773 (6-cycles).
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.
A335064 is a subsequence.
Second column of A003506.
Cf. A002378, A046092, A028896 (irregularities of maximal k-degenerate graphs).
Cf. A347213 (Dgf at s=4).
Cf. A002378, A152811, A371912 (Zagreb indices of maximal k-degenerate graphs).
Cf. A016742, A049450.

a002378 n = n * (n + 1)
a002378_list = zipWith (*) [0..] [1..]
-- Reinhard Zumkeller, Aug 27 2012, Oct 12 2011

[n*(n+1) : n in [0..100]]; // Wesley Ivan Hurt, Oct 26 2015

A002378 := proc(n)
        n*(n+1) ;
end proc:
seq(A002378(n),n=0..100) ;

Table[n(n + 1), {n, 0, 50}] (* Robert G. Wilson v, Jun 19 2004 *)
oblongQ[n_] := IntegerQ @ Sqrt[4 n + 1]; Select[Range[0, 2600], oblongQ] (* Robert G. Wilson v, Sep 29 2011 *)
2 Accumulate[Range[0, 50]] (* Harvey P. Dale, Nov 11 2011 *)
LinearRecurrence[{3, -3, 1}, {2, 6, 12}, {0, 20}] (* Eric W. Weisstein, Jul 27 2017 *)

PARI
```
{a(n) = n*(n+1)};
```

concat(0, Vec(2*x/(1-x)^3 + O(x^100))) \\ Altug Alkan, Oct 26 2015

is(n)=my(m=sqrtint(n)); m*(m+1)==n \\ Charles R Greathouse IV, Nov 01 2018

is(n)=issquare(4*n+1) \\ Charles R Greathouse IV, Mar 16 2022

def a(n): return n*(n+1)
print([a(n) for n in range(51)]) # Michael S. Branicky, Jan 13 2022

(2 to 100 by 2).scanLeft(0)( + ) // Alonso del Arte, Sep 12 2019

G.f.: 2*x/(1-x)^3. - Simon Plouffe in his 1992 dissertation.
a(n) = a(n-1) + 2*n, a(0) = 0.
Sum_{n >= 1} a(n) = n*(n+1)*(n+2)/3 (cf. A007290, partial sums).
Sum_{n >= 1} 1/a(n) = 1. (Cf. Tijdeman)
Sum_{n >= 1} (-1)^(n+1)/a(n) = log(4) - 1 = A016627 - 1 [Jolley eq (235)].
1 = 1/2 + Sum_{n >= 1} 1/(2*a(n)) = 1/2 + 1/4 + 1/12 + 1/24 + 1/40 + 1/60 + ... with partial sums: 1/2, 3/4, 5/6, 7/8, 9/10, 11/12, 13/14, ... - Gary W. Adamson, Jun 16 2003
a(n)*a(n+1) = a(n*(n+2)); e.g., a(3)*a(4) = 12*20 = 240 = a(3*5). - Charlie Marion, Dec 29 2003
Sum_{k = 1..n} 1/a(k) = n/(n+1). - Robert G. Wilson v, Feb 04 2005
a(n) = A046092(n)/2. - Zerinvary Lajos, Jan 08 2006
Log 2 = Sum_{n >= 0} 1/a(2n+1) = 1/2 + 1/12 + 1/30 + 1/56 + 1/90 + ... = (1 - 1/2) + (1/3 - 1/4) + (1/5 - 1/6) + (1/7 - 1/8) + ... = Sum_{n >= 0} (-1)^n/(n+1) = A002162. - Gary W. Adamson, Jun 22 2003
a(n) = A110660(2*n). - N. J. A. Sloane, Sep 21 2005
a(n-1) = n^2 - n = A000290(n) - A000027(n) for n >= 1. a(n) is the inverse (frequency distribution) sequence of A000194(n). - Mohammad K. Azarian, Jul 26 2007
(2, 6, 12, 20, 30, ...) = binomial transform of (2, 4, 2). - Gary W. Adamson, Nov 28 2007
a(n) = 2*Sum_{i=0..n} i = 2*A000217(n). - Artur Jasinski, Jan 09 2007, and Omar E. Pol, May 14 2008
a(n) = A006503(n) - A000292(n). - Reinhard Zumkeller, Sep 24 2008
a(n) = A061037(4*n) = (n+1/2)^2 - 1/4 = ((2n+1)^2 - 1)/4 = (A005408(n)^2 - 1)/4. - Paul Curtz, Oct 03 2008 and Klaus Purath, Jan 13 2022
a(0) = 0, a(n) = a(n-1) + 1 + floor(x), where x is the minimal positive solution to fract(sqrt(a(n-1) + 1 + x)) = 1/2. - Hieronymus Fischer, Dec 31 2008
E.g.f.: (x+2)*x*exp(x). - Geoffrey Critzer, Feb 06 2009
Product_{i >= 2} (1-1/a(i)) = -2*sin(Pi*A001622)/Pi = -2*sin(A094886)/A000796 = 2*A146481. - R. J. Mathar, Mar 12 2009, Mar 15 2009
E.g.f.: ((-x+1)*log(-x+1)+x)/x^2 also Integral_{x = 0..1} ((-x+1)*log(-x+1) + x)/x^2 = zeta(2) - 1. - Stephen Crowley, Jul 11 2009
a(A007018(n)) = A007018(n+1), i.e., A007018(n+1) = A007018(n)-th oblong numbers. - Jaroslav Krizek, Sep 13 2009
a(n) = floor((n + 1/2)^2). a(n) = A035608(n) + A004526(n+1). - Reinhard Zumkeller, Jan 27 2010
a(n) = 2*(2*A006578(n) - A035608(n)). - Reinhard Zumkeller, Feb 07 2010
a(n-1) = floor(n^5/(n^3 + n^2 + 1)). - Gary Detlefs, Feb 11 2010
For n > 1: a(n) = A173333(n+1, n-1). - Reinhard Zumkeller, Feb 19 2010
a(n) = A004202(A000217(n)). - Reinhard Zumkeller, Feb 12 2011
a(n) = A188652(2*n+1) + 1. - Reinhard Zumkeller, Apr 13 2011
For n > 0 a(n) = 1/(Integral_{x=0..Pi/2} 2*(sin(x))^(2*n-1)*(cos(x))^3). - Francesco Daddi, Aug 02 2011
a(n) = A002061(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(0) = 0, a(n) = A005408(A034856(n)) - A005408(n-1). - Ivan N. Ianakiev, Dec 06 2012
a(n) = A005408(A000096(n)) - A005408(n). - Ivan N. Ianakiev, Dec 07 2012
a(n) = A001318(n) + A085787(n). - Omar E. Pol, Jan 11 2013
Sum_{n >= 1} 1/(a(n))^(2s) = Sum_{t = 1..2*s} binomial(4*s - t - 1, 2*s - 1) * ( (1 + (-1)^t)*zeta(t) - 1). See Arxiv:1301.6293. - R. J. Mathar, Feb 03 2013
a(n)^2 + a(n+1)^2 = 2 * a((n+1)^2), for n > 0. - Ivan N. Ianakiev, Apr 08 2013
a(n) = floor(n^2 * e^(1/n)) and a(n-1) = floor(n^2 / e^(1/n)). - Richard R. Forberg, Jun 22 2013
a(n) = 2*C(n+1, 2), for n >= 0. - Felix P. Muga II, Mar 11 2014
A005369(a(n)) = 1. - Reinhard Zumkeller, Jul 05 2014
Binomial transform of [0, 2, 2, 0, 0, 0, ...]. - Alois P. Heinz, Mar 10 2015
a(2n) = A002943(n) for n >= 0, a(2n-1) = A002939(n) for n >= 1. - M. F. Hasler, Oct 11 2015
For n > 0, a(n) = 1/(Integral_{x=0..1} (x^(n-1) - x^n) dx). - Rick L. Shepherd, Oct 26 2015
a(n) = A005902(n) - A007588(n). - Peter M. Chema, Jan 09 2016
For n > 0, a(n) = lim_{m -> oo} (1/m)*1/(Sum_{i=m*n..m*(n+1)} 1/i^2), with error of ~1/m. - Richard R. Forberg, Jul 27 2016
From Ilya Gutkovskiy, Jul 28 2016: (Start)
Dirichlet g.f.: zeta(s-2) + zeta(s-1).
Convolution of nonnegative integers (A001477) and constant sequence (A007395).
Sum_{n >= 0} a(n)/n! = 3*exp(1). (End)
From Charlie Marion, Mar 06 2020: (Start)
a(n)*a(n+2k-1) + (n+k)^2 = ((2n+1)*k + n^2)^2.
a(n)*a(n+2k) + k^2 = ((2n+1)*k + a(n))^2. (End)
Product_{n>=1} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)/Pi. - Amiram Eldar, Jan 20 2021
A generalization of the Dec 29 2003 formula, a(n)*a(n+1) = a(n*(n+2)), follows. a(n)*a(n+k) = a(n*(n+k+1)) + (k-1)*n*(n+k+1). - Charlie Marion, Jan 02 2023
a(n) = A016742(n) - A049450(n). - Leo Tavares, Mar 15 2025

Additional comments from Michael Somos

Comment and cross-reference added by Christopher Hunt Gribble, Oct 13 2009

A060819 a(n) = n / gcd(n,4).

Original entry on oeis.org

1, 1, 3, 1, 5, 3, 7, 2, 9, 5, 11, 3, 13, 7, 15, 4, 17, 9, 19, 5, 21, 11, 23, 6, 25, 13, 27, 7, 29, 15, 31, 8, 33, 17, 35, 9, 37, 19, 39, 10, 41, 21, 43, 11, 45, 23, 47, 12, 49, 25, 51, 13, 53, 27, 55, 14, 57, 29, 59, 15, 61, 31, 63, 16, 65, 33, 67, 17, 69, 35, 71, 18, 73, 37, 75, 19

Offset: 1

Len Smiley, Apr 30 2001

From Peter Bala, Feb 19 2019: (Start)
We make some general remarks about the sequence a(n) = numerator(n/(n + k)) = n/gcd(n,k) for k a fixed positive integer. The present sequence is the case k = 4. Several other cases are listed in the Crossrefs. In addition to being multiplicative these sequences are also strong divisibility sequences, that is, gcd(a(n),a(m)) = a(gcd(n, m)) for n, m >= 1. In particular, it follows that a(n) is a divisibility sequence: if n divides m then a(n) divides a(m).
By the multiplicativeness and strong divisibility property of the sequence a(n) it follows that if gcd(n, m) = 1 then a(a(n)*a(m) ) = a(a(n)) * a(a(m)), a(a(a(n))*a(a(m)) ) = a(a(a(n))) * a(a(a(m))) and so on.
The sequence a(n) has the rational generating function Sum_{d divides k} f(d)*x^d/(1 - x^d)^2, where f(n) is the Dirichlet inverse of the Euler totient function A000010. f(n) is a multiplicative function defined on prime powers p^k by f(p^k) = 1 - p. See A023900. Cf. A181318. (End)

			From _Peter Bala_, Feb 21 2019: (Start)
Sum_{n >= 1} n*a(n)*x^n = G(x) - 2*G(x^2) - 4*G(x^4), where G(x) = x*(1 + x)/(1 - x)^3.
Sum_{n >= 1} (1/n)*a(n)*x^n = H(x) - (1/2)*H(x^2) - (1/4)*H(x^4), where H(x) = x/(1 - x).
Sum_{n >= 1} (1/n^2)*a(n)*x^n = L(x) - (1/2^2)*L(x^2) - (1/4^2)*L(x^4), where L(x) = Log(1/(1 - x)).
Sum_{n >= 1} (1/a(n))*x^n = L(x) + (1/2)*L(x^2) + (1/2)*L(x^4). (End)

Harry J. Smith, Table of n, a(n) for n = 1..1000
Peter Bala, A note on the sequence of numerators of a rational function, 2019.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,2,0,0,0,-1).

Cf. A026741, A051176, A060791, A060789. Cf. Other sequences given by the formula numerator(n/(n + k)): A106608 thru A106612 (k = 7 thru 11), A051724 (k = 12), A106614 thru A106621 (k = 13 thru 20).

Cf. A000010, A023900, A061037, A061038, A109008, A145979, A167192, A176895, A181318, A220466.

List([1..80],n->n/Gcd(n,4)); # Muniru A Asiru, Feb 20 2019

a060819 n = n `div` a109008 n  -- Reinhard Zumkeller, Nov 25 2013

[n/GCD(n,4): n in [1..80]]; // G. C. Greubel, Sep 19 2018

A060819 := n -> numer(1/2-1/(n+2)): seq(A060819(n),n=1..75); # Gary Detlefs, Sep 16 2011

Mathematica
```
f[n_]:= n/GCD[n, 4]; Array[f, 80]
```

a(n) = { n / gcd(n, 4) } \\ Harry J. Smith, Jul 12 2009

[lcm(n,4)/4 for n in (1..80)] # Zerinvary Lajos, Jun 07 2009

G.f.: x*(1 +x +3*x^2 +x^3 +3*x^4 +x^5 +x^6)/(1 - x^4)^2.
a(n) = 2*a(n-4) - a(n-8).
a(n) = (n/16)*(11 - 5*(-1)^n - i^n - (-i)^n). - Ralf Stephan, Mar 15 2003
a(2*n+1) = a(4*n+2) = 2*n+1, a(4*n+4) = n+1. - Ralf Stephan, Jun 10 2005
Multiplicative with a(2^e) = 2^max(0, e-2), a(p^e) = p^e, p >= 3. - Mitch Harris, Jun 29 2005
a(n) = A167192(n+4,4). - Reinhard Zumkeller, Oct 30 2009
From R. J. Mathar, Apr 18 2011: (Start)
a(n) = A109045(n)/4.
Dirichlet g.f.: zeta(s-1)*(1-1/2^s-1/2^(2s)). (End)
a(n+4) - a(n) = A176895(n).  - Paul Curtz, Apr 05 2011
a(n) = numerator(Sum_{k=1..n} 1/((k+1)*(k+2))). This summation has a closed form of 1/2 - 1/(n+2) and denominator of A145979(n). - Gary Detlefs, Sep 16 2011
a((2*n-1)*2^p) = ceiling(2^(p-2))*(2*n-1), p >= 0 and n >= 1. - Johannes W. Meijer, Feb 06 2013
a(n) = n / A109008(n). - Reinhard Zumkeller, Nov 25 2013
a(n) = denominator((2n-4)/n). - Wesley Ivan Hurt, Dec 22 2016
From Peter Bala, Feb 21 2019: (Start)
O.g.f.: Sum_{n >= 0} a(n)*x^n = F(x) - F(x^2) - F(x^4), where F(x) = x/(1 - x)^2.
More generally, Sum_{n >= 0} (a(n)^m)*x^n = F(m,x) + (1 - 2^m)*( F(m,x^2) + F(m,x^4) ), where F(m,x) = A(m,x)/(1 - x)^(m+1) with A(m,x) the m_th Eulerian polynomial: A(1,x) = x, A(2,x) = x*(1 + x), A(3,x) = x*(1 + 4*x + x^2) - see A008292.
Repeatedly applying the Euler operator x*d/dx or its inverse operator to the o.g.f. for the sequence a(n) produces generating functions for the sequences ((n^m)*a(n))n>=1 for m in Z. Some examples are given below.
(End)
Sum_{k=1..n} a(k) ~ (11/32) * n^2. - Amiram Eldar, Nov 25 2022
E.g.f.: x*(8*cosh(x) + sin(x) + 3*sinh(x))/8. - Stefano Spezia, Dec 02 2023

A061038 Denominator of 1/4 - 1/n^2.

Original entry on oeis.org

1, 36, 16, 100, 9, 196, 64, 324, 25, 484, 144, 676, 49, 900, 256, 1156, 81, 1444, 400, 1764, 121, 2116, 576, 2500, 169, 2916, 784, 3364, 225, 3844, 1024, 4356, 289, 4900, 1296, 5476, 361, 6084, 1600, 6724, 441, 7396, 1936, 8100, 529, 8836

Offset: 2

N. J. A. Sloane, May 26 2001

Harry J. Smith, Table of n, a(n) for n = 2..1000
Index entries for linear recurrences with constant coefficients, signature (0,0,0,3,0,0,0,-3,0,0,0,1).

See A061037 for comments, references, links.

Cf. A145979. - Bruno Berselli, Mar 21 2011

import Data.Ratio ((%), denominator)
a061038 n = denominator (1%4 - 1%n^2)  -- Reinhard Zumkeller, Jan 22 2012

[ Denominator(1/4-1/n^2): n in [2..50] ]; // Vincenzo Librandi, Feb 10 2011

Table[Denominator[1/4 - 1/n^2], {n, 2, 60}] (* Stefan Steinerberger, Apr 08 2006 *)
LinearRecurrence[{0,0,0,3,0,0,0,-3,0,0,0,1},{1,36,16,100,9,196,64,324, 25,484,144,676},50] (* Harvey P. Dale, Aug 05 2018 *)

a(n) = { denominator(1/4 - 1/n^2) } \\ Harry J. Smith, Jul 17 2009

def A061038(n): return denominator(1/4 - 1/n^2)
[A061038(n) for n in range(2,51)] # G. C. Greubel, Apr 18 2023

a(4n+2) = (2n+1)^2, a(2n+3) = (4n+6)^2, a(4n+4) = (4n+4)^2. - Ralf Stephan, Jun 10 2005
a(n) = 3*a(n-4) - 3*a(n-8) + a(n-12). - Paul Curtz, Feb 25 2011
From Bruno Berselli, Mar 21 2011: (Start)
G.f.: x^2*(1 +36*x +16*x^2 +100*x^3 +6*x^4 +88*x^5 +16*x^6 +24*x^7 +x^8 +4*x^9 +4*x^11)/(1-x^4)^3.
a(n) = (1/64)*( n*(16 - (1+(-1)^n)*(5-i^n)) )^2 with i=sqrt(-1).
a(n) = (n/(n-4))^2 * a(n-4) for n>5. (End)
a(n) = 4*n^2 / gcd(4*n^2, (n^2-4)). - Colin Barker, Jan 13 2014
Sum_{n>=2} 1/a(n) = Pi^2/6 - 1/4. - Amiram Eldar, Aug 12 2022

More terms from Stefan Steinerberger, Apr 08 2006

A078371 a(n) = (2n+5)(2*n+1).

Original entry on oeis.org

5, 21, 45, 77, 117, 165, 221, 285, 357, 437, 525, 621, 725, 837, 957, 1085, 1221, 1365, 1517, 1677, 1845, 2021, 2205, 2397, 2597, 2805, 3021, 3245, 3477, 3717, 3965, 4221, 4485, 4757, 5037, 5325, 5621, 5925, 6237, 6557, 6885, 7221, 7565, 7917, 8277, 8645

Offset: 0

Wolfdieter Lang, Nov 29 2002

This is the generic form of D in the (nontrivially) solvable Pell equation x^2 - D*y^2 = +4. See A077428 and A078355.
Consider all primitive Pythagorean triples (a,b,c) with c-a=8, sequence gives values of a. (Corresponding values for b are A017113(n), while c follows A078370(n).) - Lambert Klasen (Lambert.Klasen(AT)gmx.net), Nov 19 2004
From Vincenzo Librandi, Aug 08 2010: (Start)
The identity (4*n^3 + 18*n^2 + 24*n + 9)^2 - (4*n^2 + 12*n + 5)*(2*n^2 + 6*n + 4)^2 = 1 (see Ramasamy's paper in link) can be written as A141530(n+2)^2 - a(n)*A046092(n+1)^2 = 1.
a(n)^3 + 6*a(n)^2 + 9*a(n) + 4 is a square: in fact, a(n)^3 + 6*a(n)^2 + 9*a(n) + 4 = (a(n) + 1)^2*(a(n) + 4), where a(n) + 4 = (2*n+3)^2. (End)
Products of two positive odd integers with difference 4 (i.e., 1*5, 3*7, 5*9, 7*11, 9*13, ...). - Wesley Ivan Hurt, Nov 19 2013
Starting with stage 1, the number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 675", based on the 5-celled von Neumann neighborhood. - Robert Price, May 21 2016
The continued fraction expansion of (sqrt(a(n))-1)/2 is [n; {1,2*n+1}] with periodic part of length 2: repeat {1,2*n+1}. - Ron Knott, May 11 2017
a(n) is the sum of 2*n+5 consecutive integers starting from n-1. - Bruno Berselli, Jan 16 2018
The continued fraction expansion of sqrt(a(n)) is [2n+2; {1, n, 2, n, 1, 4n+4}]. For n=0, this collapses to [2; {4}]. - Magus K. Chu, Aug 26 2022

Bruno Berselli, Table of n, a(n) for n = 0..1000
Soren Laing Aletheia-Zomlefer, Lenny Fukshansky, and Stephan Ramon Garcia, The Bateman-Horn Conjecture: Heuristics, History, and Applications, arXiv:1807.08899 [math.NT], 2018-2019. See Example 6.6.5 p. 34.
A. M. S. Ramasamy, Polynomial solutions for the Pell's equation, Indian Journal of Pure and Applied Mathematics 25 (1994), p. 578.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Subsequence of A077425 (D values (not a square) for which Pell x^2 - D*y^2 = +4 is solvable in positive integers).
Cf. A017113, A046092, A061037, A077428, A078355, A078370.
Supersequence of A143206.

[(2*n+5)*(2*n+1): n in [0..100]]; // G. C. Greubel, Sep 19 2018

seq((2*n+5)*(2*n+1), n=0..48); # Emeric Deutsch, Feb 24 2005

Table[(2 n + 5) (2 n + 1), {n, 0, 100}] (* Wesley Ivan Hurt, Nov 19 2013 *)
LinearRecurrence[{3,-3,1},{5,21,45},50] (* Harvey P. Dale, Oct 18 2020 *)

lista(nn) = {for (n=0, nn, print1((2*n+1)*(2*n+5), ", "));} \\ Michel Marcus, Nov 21 2013

a(n) = 8*(binomial(n+2, 2)-1)+5, hence subsequence of A004770 (5 (mod 8) numbers).
G.f.: (5 + 6*x - 3*x^2)/(1-x)^3.
a(n) = A061037(2*n+1) = (2*n+3)^2 - 4. For A061037: a(2*n+1) = (2*n+1)*(2*n+5) = (2*n+3)^2-4. - Paul Curtz, Sep 24 2008
a(n) = 8*(n+1) + a(n-1) for n > 0, a(0)=5. - Vincenzo Librandi, Aug 08 2010
From Ilya Gutkovskiy, May 22 2016: (Start)
E.g.f.: (5 + 4*x*(4 + x))*exp(x).
Sum_{n>=0} 1/a(n) = 1/3. (End)
Sum_{n>=0} (-1)^n/a(n) = 1/6. - Amiram Eldar, Oct 08 2023

More terms from Emeric Deutsch, Feb 24 2005

A061035 Triangle T(m,n) = numerator of 1/m^2 - 1/n^2, n >= 1, m=n,n-1,n-2,...,1.

Original entry on oeis.org

0, 0, 3, 0, 5, 8, 0, 7, 3, 15, 0, 9, 16, 21, 24, 0, 11, 5, 1, 2, 35, 0, 13, 24, 33, 40, 45, 48, 0, 15, 7, 39, 3, 55, 15, 63, 0, 17, 32, 5, 56, 65, 8, 77, 80, 0, 19, 9, 51, 4, 3, 21, 91, 6, 99, 0, 21, 40, 57, 72, 85, 96, 105, 112, 117, 120, 0, 23, 11, 7, 5, 95, 1, 119, 1, 5, 35, 143, 0

Offset: 1

N. J. A. Sloane, May 26 2001

Wavelengths in hydrogen spectrum are given by Rydberg's formula 1/wavelength = constant*(1/m^2 - 1/n^2).

			Triangle 1/m^2-1/n^2, m >= 1, 1<=n<=m, (i.e. with rows reversed) begins
0
3/4, 0
8/9, 5/36, 0
15/16, 3/16, 7/144, 0
24/25, 21/100, 16/225, 9/400, 0
35/36, 2/9, 1/12, 5/144, 11/900, 0

J. E. Brady and G. E. Humiston, General Chemistry, 3rd. ed., Wiley; p. 77.

Reinhard Zumkeller, Rows n=1..100 of triangle, flattened
J. J. O'Connor and E. F. Robertson, Johannes Robert Rydberg
Eric Weisstein's World of Physics, Balmer Formula

Cf. A061036. Rows give A061037-A061050.

Cf. A126252.

import Data.Ratio ((%), numerator)
a061035 n k = a061035_tabl !! (n-1) !! (k-1)
a061035_row = map numerator . balmer where
   balmer n = map (subtract (1 % n ^ 2) . (1 %) . (^ 2)) [n, n-1 .. 1]
a061035_tabl = map a061035_row [1..]
-- Reinhard Zumkeller, Apr 12 2012

t[m_, n_] := Numerator[1/m^2 - 1/n^2]; Table[t[m, n], {n, 1, 12}, {m, n, 1, -1}] // Flatten(* Jean-François Alcover, Oct 17 2012 *)

More terms from Naohiro Nomoto, Jul 15 2001

A033567 a(n) = (2n-1)(4*n-1).

Original entry on oeis.org

1, 3, 21, 55, 105, 171, 253, 351, 465, 595, 741, 903, 1081, 1275, 1485, 1711, 1953, 2211, 2485, 2775, 3081, 3403, 3741, 4095, 4465, 4851, 5253, 5671, 6105, 6555, 7021, 7503, 8001, 8515, 9045, 9591, 10153, 10731, 11325, 11935, 12561, 13203, 13861, 14535, 15225

Offset: 0

N. J. A. Sloane

a(n+1) = A005563(1), A061037(3), A061039(5), A061041(7), A061043(9), A061045(11), A061047(13), A061049(15). Lyman, Balmer, Paschen, Brackett, Pfund, Humphreys, Hansen-Strong, ... spectra of hydrogen. - Paul Curtz, Oct 08 2008

Sequence found by reading the segment [1, 3] together with the line from 3, in the direction 3, 21, ..., in the square spiral whose vertices are the triangular numbers A000217. - Omar E. Pol, Sep 03 2011

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

Cf. A045944, A014634.

[(2*n-1)*(4*n-1): n in [0..50]]; // G. C. Greubel, Sep 19 2018

Table[(2*n - 1)*(4*n - 1), {n, 0, 50}] (* G. C. Greubel, Jul 06 2017 *)
LinearRecurrence[{3,-3,1},{1,3,21},50] (* Harvey P. Dale, Aug 25 2019 *)

vector(60, n, n--; (2*n-1)*(4*n-1)) \\ Michel Marcus, Apr 12 2015

a(n) = a(n-1) + 16*n - 14 (with a(0)=1). - Vincenzo Librandi, Nov 17 2010
From G. C. Greubel, Jul 06 2017: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-2).
E.g.f.: (1 + 2*x + 8*x^2)*exp(x).
G.f.: (1 + 15*x^2)/(1 - x)^3. (End)
From Amiram Eldar, Jan 03 2022: (Start)
Sum_{n>=0} 1/a(n) = 1 + Pi/4 - log(2)/2.
Sum_{n>=0} (-1)^n/a(n) = 1 + (sqrt(2)-1)*Pi/4 + log(sqrt(2)-1)/sqrt(2). (End)

More terms from Michel Marcus, Apr 12 2015

A144437 Period 3: repeat [3, 3, 1].

Original entry on oeis.org

3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3, 1, 3, 3

Offset: 1

Paul Curtz, Oct 05 2008

The sequence is generated from numerators in the energy differences of the hydrogen spectrum: A005563(1), A061037(4), A061039(6), A061041(8), A061043(10), A061045(12), A061047(14), A061049(16), ...
Conjecture: a(n) is the separatix. See A045944.
Also the decimal expansion of the constant 3310/999. - R. J. Mathar, May 21 2009
Continued fraction expansion of A171417.
Greatest common divisor of (n+1)^2-1 and (n+1)^2+2. - Bruno Berselli, Mar 08 2017

Index entries for linear recurrences with constant coefficients, signature (0,0,1).

Cf. A000217, A045944, A109007, A164306, A167192, A169609, A256095.

Numerators in the energy differences of the hydrogen spectrum: A005563(1), A061037(4), A061039(6), A061041(8), A061043(10), A061045(12), A061047(14), A061049(16), ...

&cat [[3, 3, 1]^^30]; // Wesley Ivan Hurt, Jul 02 2016

seq(op([3, 3, 1]), n=1..50); # Wesley Ivan Hurt, Jul 02 2016

A144437[n_]:=Denominator[n/3]; Array[A144437,100] (* Enrique Pérez Herrero, Oct 05 2011 *)
CoefficientList[Series[(3 + 3 x + x^2)/(1 - x^3), {x, 0, 120}], x] (* Michael De Vlieger, Jul 02 2016 *)
Table[Mod[2*n^2 + 1, 3,1], {n,1,50}] (* G. C. Greubel, Aug 24 2017 *)

a(n)=if(n%3,3,1) \\ Charles R Greathouse IV, Sep 28 2015

a(n) = (7-4*cos(2*Pi*n/3))/3. - Jaume Oliver Lafont, Nov 23 2008
G.f.: x*(3 + 3*x + x^2)/((1 - x)*(1 + x + x^2)). - R. J. Mathar, May 21 2009
a(n) = 3/gcd(n,3). - Reinhard Zumkeller, Oct 30 2009
a(n) = denominator(n^k/3), where k>0 is an integer. - Enrique Pérez Herrero, Oct 05 2011
a(n) = gcd(T(n+1), T(2)) = A256095(n+1, 2), with the triangular numbers T = A000217, for n >= 1. - Wolfdieter Lang, Mar 17 2015
a(n) = a(n-3) for n>3; a(n) = A169609(n) for n>0. - Wesley Ivan Hurt, Jul 02 2016
E.g.f.: (1/3)*(7*exp(x) - 4*exp(-x/2)*cos(sqrt(3)*x/2) - 3). - G. C. Greubel, Aug 24 2017
From Nicolas Bělohoubek, Nov 11 2021: (Start)
a(n) = 9/(a(n-2)*a(n-1)).
a(n) = 7 - a(n-2) - a(n-1). See also A052901 or A069705. (End)

Edited by R. J. Mathar, May 21 2009

A181318 a(n) = A060819(n)^2.

Original entry on oeis.org

0, 1, 1, 9, 1, 25, 9, 49, 4, 81, 25, 121, 9, 169, 49, 225, 16, 289, 81, 361, 25, 441, 121, 529, 36, 625, 169, 729, 49, 841, 225, 961, 64, 1089, 289, 1225, 81, 1369, 361, 1521, 100, 1681, 441, 1849, 121, 2025, 529, 2209, 144, 2401, 625, 2601, 169, 2809, 729

Offset: 0

Paul Curtz, Jan 26 2011

The first sequence, p=0, of the family A060819(n)*A060819(n+p).
Hence array
p=0:  0, 1, 1,  9,  1, 25,  9,  49,    a(n)=A060819(n)^2,
p=1:  0, 1, 3,  3,  5, 15, 21,  14,     A064038(n),
p=2:  0, 3, 1, 15,  3, 35,  6,  63,     A198148(n),
p=3:  0, 1, 5,  9,  7, 10, 27,  35,     A160050(n),
p=4:  0, 5, 3, 21,  2, 45, 15,  77,     A061037(n),
p=5:  0, 3, 7,  6,  9, 25, 33,  21,     A178242(n),
p=6:  0, 7, 2, 27,  5, 55,  9,  91,     A217366(n),
p=7:  0, 2, 9, 15, 11, 15, 39,  49,     A217367(n),
p=8:  0, 9, 5, 33,  3, 65, 21, 105,     A180082(n).
Compare columns 2, 3 and 5, columns 4 and 7 and columns 6 and 8.
From Peter Bala, Feb 19 2019: (Start)
We make some general remarks about the sequence a(n) = numerator(n^2/(n^2 + k^2)) = (n/gcd(n,k))^2 for k a fixed positive integer (we suppress the dependence of a(n) on k). The present sequence corresponds to the case k = 4.
a(n) is a quasi-polynomial in n. In fact, a(n) = n^2/b(n) where b(n) = gcd(n^2,k^2) is a purely periodic sequence in n.
In addition to being multiplicative these sequences are also strong divisibility sequences, that is, gcd(a(n),a(m)) = a(gcd(n,m)) for n, m >= 1. In particular, it follows that a(n) is a divisibility sequence: if n divides m then a(n) divides a(m).
By the multiplicativeness and strong divisibility property of the sequence a(n) it follows that if gcd(n,m) = 1 then a( a(n)*a(m) ) = a(a(n)) * a(a(m)), a( a(a(n))*a(a(m)) ) = a(a(a(n))) * a(a(a(m))) and so on.
The sequence a(n) has the rational generating function Sum_{d divides k} f(d)*F(x^d), where F(x) = x*(1 + x)/(1 - x)^3 = x + 4*x^2 + 9*x^3 + 16*x^4 + ... is the o.g.f. for the squares A000290, and where f(n) is the Dirichlet inverse of the Jordan totient function J_2(n) - see A007434. The function f(n) is multiplicative and is defined on prime powers p^k by f(p^k) = (1 - p^2). See A046970. Cf. A060819. (End)
a(n-4) is the constant needed to complete the n-polygonal numbers into squares (see A377851); a(-1) = 1, which completes the triangle numbers, is not shown in the data. - Jonathan Dushoff, Nov 12 2024

G. C. Greubel, Table of n, a(n) for n = 0..10000
Wikipedia, Completing the square.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,3,0,0,0,-3,0,0,0,1).

Cf. A181829, A046970, A007434, A060819, A168077, A377851.

[n^2/GCD(n,4)^2: n in [0..100]]; // G. C. Greubel, Sep 19 2018

a:=n->n^2/gcd(n,4)^2: seq(a(n),n=0..60); # Muniru A Asiru, Feb 20 2019

Table[n^2/GCD[n,4]^2, {n,0,100}] (* G. C. Greubel, Sep 19 2018 *)
LinearRecurrence[{0,0,0,3,0,0,0,-3,0,0,0,1},{0,1,1,9,1,25,9,49,4,81,25,121},60] (* Harvey P. Dale, Jan 18 2025 *)

a(n)=n^2/gcd(n,4)^2 \\ Charles R Greathouse IV, Dec 21 2011

[n^2/gcd(n, 4)^2 for n in (0..100)] # G. C. Greubel, Feb 20 2019

a(2*n) = A168077(n), a(2*n+1) = A016754(n).
a(n) = 3*a(n-4) - 3*a(n-8) + a(n-12).
G.f.: x*(1 + x + 9*x^2 + x^3 + 22*x^4 + 6*x^5 + 22*x^6 + x^7 + 9*x^8 + x^9 + x^10)/(1-x^4)^3. - R. J. Mathar, Mar 10 2011
From Peter Bala, Feb 19 2019: (Start)
a(n) = numerator(n^2/(n^2 + 16)) = n^2/(gcd(n^2,16)) = (n/gcd(n,4))^2.
a(n) = n^2/b(n), where b(n) = [1, 4, 1, 16, 1, 4, 1, 16, ...] is a purely periodic sequence of period 4.
a(n) is a quasi-polynomial in n: a(4*n) = n^2; a(4*n + 1) = (4*n + 1)^2; a(4*n + 2) = (2*n + 1)^2; a(4*n + 3) = (4*n + 3)^2.
O.g.f.: Sum_{d divides 4} A046970(d)*x^d*(1 + x^d)/(1 - x^d)^3 = x*(1 + x)/(1 - x)^3 - 3*x^2*(1 + x^2)/(1 - x^2)^3 - 3*x^4*(1 + x^4)/(1 - x^4)^3. (End)
Sum_{n>=1} 1/a(n) = 5*Pi^2/12. - Amiram Eldar, Aug 12 2022
From Amiram Eldar, Nov 25 2022: (Start)
Multiplicative with a(2^e) = 4^max(0, e-2), and a(p^e) = p^(2*e) for p > 2.
Dirichlet g.f.: zeta(s-2)*(1 - 3/2^s - 3/4^s).
Sum_{k=1..n} a(k) ~ (37/192) * n^3. (End)
a(n) = (37 - 27*(-1)^n - 3*(-1)^(n*(n-1)/2) - 3*(-1)^(n*(n+1)/2)) * n^2/64. - Vaclav Kotesovec, Nov 14 2024

Edited by Jean-François Alcover, Oct 01 2012 and Jan 15 2013

More terms from Michel Marcus, Jun 09 2014

cp's OEIS Frontend

A247004 Denominator of (n+4)/gcd(n, 4)^2, a 16-periodic sequence that associates A061037 with A106617.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A144392 Inverse binomial transform of A061037 (read as with offset 0).

Views

Author

Keywords

Formula

Extensions

A002378 Oblong (or promic, pronic, or heteromecic) numbers: a(n) = n*(n+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Scala

Formula

Extensions

A060819 a(n) = n / gcd(n,4).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

Sage

Formula

A061038 Denominator of 1/4 - 1/n^2.

Views

Author

Keywords

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

SageMath

Formula

Extensions

A078371 a(n) = (2*n+5)*(2*n+1).

Views

Author

Keywords

Comments

Links

Crossrefs

A078371 a(n) = (2n+5)(2*n+1).

A033567 a(n) = (2n-1)(4*n-1).