A002144 -id:A002144 - cp's OEIS frontend

A001844 Centered square numbers: a(n) = 2n(n+1)+1. Sums of two consecutive squares. Also, consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z; then sequence gives Z values.

1, 5, 13, 25, 41, 61, 85, 113, 145, 181, 221, 265, 313, 365, 421, 481, 545, 613, 685, 761, 841, 925, 1013, 1105, 1201, 1301, 1405, 1513, 1625, 1741, 1861, 1985, 2113, 2245, 2381, 2521, 2665, 2813, 2965, 3121, 3281, 3445, 3613, 3785, 3961, 4141, 4325, 4513

Offset: 0

N. J. A. Sloane

These are Hogben's central polygonal numbers denoted by
...2...
....P..
...4.n.
Numbers of the form (k^2+1)/2 for k odd.
(y(2x+1))^2 + (y(2x^2+2x))^2 = (y(2x^2+2x+1))^2. E.g., let y = 2, x = 1; (2(2+1))^2 + (2(2+2))^2 = (2(2+2+1))^2, (2(3))^2 + (2(4))^2 = (2(5))^2, 6^2 + 8^2 = 10^2, 36 + 64 = 100. - Glenn B. Cox (igloos_r_us(AT)canada.com), Apr 08 2002
a(n) is also the number of 3 X 3 magic squares with sum 3(n+1). - Sharon Sela (sharonsela(AT)hotmail.com), May 11 2002
For n > 0, a(n) is the smallest k such that zeta(2) - Sum_{i=1..k} 1/i^2 <= zeta(3) - Sum_{i=1..n} 1/i^3. - Benoit Cloitre, May 17 2002
Number of convex polyominoes with a 2 X (n+1) minimal bounding rectangle.
The prime terms are given by A027862. - Lekraj Beedassy, Jul 09 2004
First difference of a(n) is 4n = A008586(n). Any entry k of the sequence is followed by k + 2*(1 + sqrt(2k - 1)). - Lekraj Beedassy, Jun 04 2006
Integers of the form 1 + x + x^2/2 (generating polynomial is Schur's polynomial as in A127876). - Artur Jasinski, Feb 04 2007
If X is an n-set and Y and Z disjoint 2-subsets of X then a(n-4) is equal to the number of 4-subsets of X intersecting both Y and Z. - Milan Janjic, Aug 26 2007
Row sums of triangle A132778. - Gary W. Adamson, Sep 02 2007
Binomial transform of [1, 4, 4, 0, 0, 0, ...]; = inverse binomial transform of A001788: (1, 6, 24, 80, 240, ...). - Gary W. Adamson, Sep 02 2007
Narayana transform (A001263) of [1, 4, 0, 0, 0, ...]. Equals A128064 (unsigned) * [1, 2, 3, ...]. - Gary W. Adamson, Dec 29 2007
k such that the Diophantine equation x^3 - y^3 = x*y + k has a solution with y = x-1. If that solution is (x,y) = (m+1,m) then m^2 + (m+1)^2 = k. Note that this Diophantine equation is an elliptic curve and (m+1,m) is an integer point on it. - James R. Buddenhagen, Aug 12 2008
Numbers k such that (k, k, 2*k-2) are the sides of an isosceles triangle with integer area. Also, k such that 2*k-1 is a square. - James R. Buddenhagen, Oct 17 2008
a(n) is also the least weight of self-conjugate partitions having n+1 different odd parts. - Augustine O. Munagi, Dec 18 2008
Prefaced with a "1": (1, 1, 5, 13, 25, 41, ...) = A153869 * (1, 2, 3, ...). - Gary W. Adamson, Jan 03 2009
Prefaced with a "1": (1, 1, 5, 13, 25, 41, ...) where a(n) = 2n*(n-1)+1, all tuples of square numbers (X-Y, X, X+Y) are produced by ((m*(a(n)-2n))^2, (m*a(n))^2, (m*(a(n)+2n-2))^2) where m is a whole number. - Doug Bell, Feb 27 2009
Equals (1, 2, 3, ...) convolved with (1, 3, 4, 4, 4, ...). E.g., a(3) = 25 = (1, 2, 3, 4) dot (4, 4, 3, 1) = (4 + 8 + 9 + 4). - Gary W. Adamson, May 01 2009
The running sum of squares taken two at a time. - Al Hakanson (hawkuu(AT)gmail.com), May 18 2009
Equals the odd integers convolved with (1, 2, 2, 2, ...). - Gary W. Adamson, May 25 2009
Equals the triangular numbers convolved with [1, 2, 1, 0, 0, 0, ...]. - Gary W. Adamson & Alexander R. Povolotsky, May 29 2009
When the positive integers are written in a square array by diagonals as in A038722, a(n) gives the numbers appearing on the main diagonal. - Joshua Zucker, Jul 07 2009
The finite continued fraction [n,1,1,n] = (2n+1)/(2n^2 + 2n + 1) = (2n+1)/a(n); and the squares of the first two denominators of the convergents = a(n). E.g., the convergents and value of [4,1,1,4] = 1/4, 1/5, 2/9, 9/41 where 4^2 + 5^2 = 41. - Gary W. Adamson, Jul 15 2010
From Keith Tyler, Aug 10 2010: (Start)
Running sum of A008574.
Square open pyramidal number; that is, the number of elements in a square pyramid of height (n) with only surface and no bottom nodes. (End)
For k>0, x^4 + x^2 + k factors over the integers iff sqrt(k) is in this sequence. - James R. Buddenhagen, Aug 15 2010
Create the simple continued fraction from Pythagorean triples to get [2n + 1; 2n^2 + 2n, 2n^2 + 2n + 1]; its value equals the rational number 2n + 1 + a(n) / (4n^4 + 8n^3 + 6n^2 + 2n + 1). - J. M. Bergot, Sep 30 2011
a(n), n >= 1, has in its prime number factorization only primes of the form 4*k+1, i.e., congruent to 1 (mod 4) (see A002144). This follows from the fact that a(n) is a primitive sum of two squares and odd. See Theorem 3.20, p. 164, in the given Niven-Zuckerman-Montgomery reference. E.g., a(3) = 25 = 5^2, a(6) = 85 = 5*17. - Wolfdieter Lang, Mar 08 2012
From Ant King, Jun 15 2012: (Start)
a(n) is congruent to 1 (mod 4) for all n.
The digital roots of the a(n) form a purely periodic palindromic 9-cycle 1, 5, 4, 7, 5, 7, 4, 5, 1.
The units' digits of the a(n) form a purely periodic palindromic 5-cycle 1, 5, 3, 5, 1.
(End)
Number of integer solutions (x,y) of |x| + |y| <= n. Geometrically: number of lattice points inside a square with vertices (n,0), (0,-n), (-n,0), (0,n). - César Eliud Lozada, Sep 18 2012
(a(n)-1)/a(n) = 2*x / (1+x^2) where x = n/(n+1). Note that in this form, this is the velocity-addition formula according to the special theory of relativity (two objects traveling at 1/(n+1) slower than c relative to each other appear to travel at 1/a(n) less than c to a stationary observer). - Christian N. K. Anderson, May 20 2013 [Corrected by Rémi Guillaume, May 22 2025]
A geometric curiosity: the envelope of the circles x^2 + (y-a(n)/2)^2 = ((2n+1)/2)^2 is the parabola y = x^2, the n=0 circle being the osculating circle at the parabola vertex. - Jean-François Alcover, Dec 02 2013
Draw n ellipses in the plane (n>0), any 2 meeting in 4 points; a(n-1) gives the number of internal regions into which the plane is divided (cf. A051890, A046092); a(n-1) = A051890(n) - 1 = A046092(n-1) + 1. - Jaroslav Krizek, Dec 27 2013
a(n) is also, of course, the scalar product of the 2-vector (n, n+1) (or (n+1, n)) with itself. The unique inverse of (n, n+1) as vector in the Clifford algebra over the Euclidean 2-space is (1/a(n))(0, n, n+1, 0) (similarly for the other vector). In general the unique inverse of such a nonzero vector v (odd element in Cl_2) is v^(-1) = (1/|v|^2) v. Note that the inverse with respect to the scalar product is not unique for any nonzero vector. See the P. Lounesto reference, sects. 1.7 - 1.12, pp. 7-14. See also the Oct 15 2014 comment in A147973. - Wolfdieter Lang, Nov 06 2014
Subsequence of A004431, for n >= 1. - Bob Selcoe, Mar 23 2016
Numbers k such that 2k - 1 is a perfect square. - Juri-Stepan Gerasimov, Apr 06 2016
The number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 574", based on the 5-celled von Neumann neighborhood. - Robert Price, May 13 2016
a(n) is the first integer in a sum of (2*n + 1)^2 consecutive integers that equals (2*n + 1)^4. - Patrick J. McNab, Dec 24 2016
Central elements of odd-length rows of the triangular array of positive integers. a(n) is the mean of the numbers in the (2*n + 1)-th row of this triangle. - David James Sycamore, Aug 01 2018
Intersection of A000982 and A080827. - David James Sycamore, Aug 07 2018
An off-diagonal of the array of Delannoy numbers, A008288, (or a row/column when the array is shown as a square). As such, this is one of the crystal ball sequences. - Jack W Grahl, Feb 15 2021 and Shel Kaphan, Jan 18 2023
a(n) appears as a solution to a "Riddler Express" puzzle on the FiveThirtyEight website. The Jan 21 2022 issue (problem) and the Jan 28 2022 issue (solution) present the following puzzle and include a proof. - Fold a square piece of paper in half, obtaining a rectangle. Fold again to obtain a square with 1/4 the size of the original square. Then make n cuts through the folded paper. a(n) is the greatest number of pieces of the unfolded paper after the cutting. - Manfred Boergens, Feb 22 2022
a(n) is (1/6) times the number of 2 X 2 triangles in the n-th order hexagram with 12*n^2 cells. - Donghwi Park, Feb 06 2024
If k is a centered square number, its index in this sequence is n = (sqrt(2k-1)-1)/2. - Rémi Guillaume, Mar 30 2025.
Row sums of the symmetric triangle of odd numbers [1]; [1, 3, 1]; [1, 3, 5, 3, 1]; [1, 3, 5, 7, 5, 3, 1]; .... - Marco Zárate, Jun 15 2025

			G.f.: 1 + 5*x + 13*x^2 + 25*x^3 + 41*x^4 + 61*x^5 + 85*x^6 + 113*x^7 + 145*x^8 + ...
The first few triples are (1,0,1), (3,4,5), (5,12,13), (7,24,25), ...
The first four such partitions, corresponding to n = 0,1,2,3, i.e., to a(n) = 1,5,13,25, are 1, 3+1+1, 5+3+3+1+1, 7+5+5+3+3+1+1. - _Augustine O. Munagi_, Dec 18 2008

T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 3.
A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, p. 125, 1964.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 50.
Pertti Lounesto, Clifford Algebras and Spinors, second edition, Cambridge University Press, 2001.
S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 483.
Ivan Niven, Herbert S. Zuckerman and Hugh L. Montgomery, An Introduction to the Theory Of Numbers, Fifth Edition, John Wiley and Sons, Inc., NY 1991.
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).
Travers et al., The Mysterious Lost Proof, Using Advanced Algebra, (1976), pp. 27.

T. D. Noe, Table of n, a(n) for n = 0..1000
M. Ahmed, J. De Loera and R. Hemmecke, Polyhedral Cones of Magic Cubes and Squares, arXiv:math/0201108 [math.CO], 2002.
U. Alfred, n and n+1 consecutive integers with equal sums of squares, Math. Mag., 35 (1962), 155-164.
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Section 2.3.
Paul Barry, Centered polygon numbers, heptagons and nonagons, and the Robbins numbers, arXiv:2104.01644 [math.CO], 2021.
Matthias Beck, Moshe Cohen, Jessica Cuomo, and Paul Gribelyuk, The number of "magic" squares and hypercubes, arXiv:math/0201013 [math.CO], 2002-2005.
Arthur T. Benjamin and Doron Zeilberger, Pythagorean Primes and Palindromic Continued Fractions, Electronic Journal of Combinatorial Number Theory, 5(1) 2005, #A30.
J. A. De Loera, D. C. Haws and M. Koppe, Ehrhart Polynomials of Matroid Polytopes and Polymatroids, arXiv:0710.4346 [math.CO], 2007; Discrete Comput. Geom., 42 (2009), 670-702.
FiveThirtyEight, "Riddler Express" paper cutting problem and solution, Jan 28 2022.
D. C. Haws, Matroids [Broken link, Oct 30 2017]
D. C. Haws, Matroids [Copy on website of Matthias Koeppe]
D. C. Haws, Matroids [Cached copy, pdf file only]
L. Hogben, Choice and Chance by Cardpack and Chessboard, Vol. 1, Max Parrish and Co, London, 1950, pp. 22 and 36.
Milan Janjic, Two Enumerative Functions. [Broken link; WayBackMachine archive.]
Milan Janjic, On a class of polynomials with integer coefficients, JIS 11 (2008) 08.5.
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.
Clark Kimberling and John E. Brown, Partial Complements and Transposable Dispersions, J. Integer Seqs., Vol. 7, 2004.
Ron Knott, Pythagorean Triples and Online Calculators.
G. Kreweras, Sur les hiérarchies de segments, Cahiers Bureau Universitaire Recherche Opérationnelle, Cahier 20, Inst. Statistiques, Univ. Paris, 1973.
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973). (Annotated scanned copy)
A. O. Munagi, Pairing conjugate partitions by residue classes, Discrete Math., 308 (2008), 2492-2501.
Mitchell Paukner, Lucy Pepin, Manda Riehl, and Jarred Wieser, Pattern Avoidance in Task-Precedence Posets, arXiv:1511.00080 [math.CO], 2015-2016.
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 A. Jr. Rochowicz, Harmonic Numbers: Insights, Approximations and Applications, Spreadsheets in Education (eJSiE), 2015, Vol. 8: Iss. 2, Article 4.
Amelia Carolina Sparavigna, Groupoid of OEIS A001844 Numbers (centered square numbers), Politecnico di Torino, Italy (2019).
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
David James Sycamore, Triangular array
Leo Tavares, Illustration: Diamond Rows
B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985), 4545-4558.
Eric Weisstein's World of Mathematics, Centered Polygonal Number, Centered Square Number, Diamond, Pythagorean Triple, and von Neumann Neighborhood.
Index entries for sequences related to centered polygonal numbers
Index entries for two-way infinite sequences.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
Nicolay Avilov, Graphical representation of the sequence members

X values are A005408; Y values are A046092.
Cf. A008586 (first differences), A005900 (partial sums), A254373 (digital roots).
Cf. A069894, A000217, A000290, A001263, A001788, A001845, A002061, A002144, A003215, A005448, A005891, A005917, A008844 (square terms), A027862 (prime terms), A048395, A051890, A056106, A101096, A127876, A128064, A132778, A147973, A153869, A240876, A251599, A000982, A080827.
Subsequence of A004431.
Right edge of A055096; main diagonal of A069480, A078475, A129312.
Row n=2 (or column k=2) of A008288.
Cf. A016754.

a001844 n = 2 * n * (n + 1) + 1
a001844_list = zipWith (+) a000290_list $ tail a000290_list
-- Reinhard Zumkeller, Dec 04 2012

[2*n^2 + 2*n + 1: n in [0..50]]; // Vincenzo Librandi, Jan 19 2013

[n: n in [0..4400] | IsSquare(2*n-1)]; // Juri-Stepan Gerasimov, Apr 06 2016

A001844:=-(z+1)**2/(z-1)**3; # Simon Plouffe in his 1992 dissertation

Table[2n(n + 1) + 1, {n, 0, 50}]
FoldList[#1 + #2 &, 1, 4 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
maxn := 47; Flatten[Table[SeriesCoefficient[Series[(n + (n - 1)*x)/(1 - x)^2, {x, 0, maxn}], k], {n, maxn}, {k, n - 1, n - 1}]] (* L. Edson Jeffery, Aug 24 2014 *)
CoefficientList[ Series[-(x^2 + 2x + 1)/(x - 1)^3, {x, 0, 48}], x] (* or *)
LinearRecurrence[{3, -3, 1}, {1, 5, 13}, 48] (* Robert G. Wilson v, Aug 01 2018 *)
Total/@Partition[Range[0,50]^2,2,1] (* Harvey P. Dale, Dec 05 2020 *)
Table[ j! Coefficient[Series[Exp[x]*(1 + 4*x + 2*x^2), {x, 0, 20}], x,
j], {j, 0, 20}] (* Nikolaos Pantelidis, Feb 07 2023 *)

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

x='x+O('x^200); Vec((1+x)^2/(1-x)^3) \\ Altug Alkan, Mar 23 2016

print([2*n*(n+1)+1 for n in range(48)]) # Michael S. Branicky, Jan 05 2021

[i**2 + (i + 1)**2 for i in range(46)] # Zerinvary Lajos, Jun 27 2008

a(n) = 2*n^2 + 2*n + 1 = n^2 + (n+1)^2.
a(n) = 1 + 3 + 5 + ... + 2*n-1 + 2*n+1 + 2*n-1 + ... + 3 + 1. - Amarnath Murthy, May 28 2001
a(n) = 1/real(z(n+1)) where z(1)=i, (i^2=-1), z(k+1) = 1/(z(k)+2i). - Benoit Cloitre, Aug 06 2002
Nearest integer to 1/Sum_{k>n} 1/k^3. - Benoit Cloitre, Jun 12 2003
G.f.: (1+x)^2/(1-x)^3.
E.g.f.: exp(x)*(1+4x+2x^2).
a(n) = a(n-1) + 4n.
a(-n) = a(n-1).
a(n) = A064094(n+3, n) (fourth diagonal).
a(n) = 1 + Sum_{j=0..n} 4*j. - Xavier Acloque, Oct 08 2003
a(n) = A046092(n)+1 = (A016754(n)+1)/2. - Lekraj Beedassy, May 25 2004
a(n) = Sum_{k=0..n+1} (-1)^k*binomial(n, k)*Sum_{j=0..n-k+1} binomial(n-k+1, j)*j^2. - Paul Barry, Dec 22 2004
a(n) = ceiling((2n+1)^2/2). - Paul Barry, Jul 16 2006
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0)=1, a(1)=5, a(2)=13. - Jaume Oliver Lafont, Dec 02 2008
a(n)*a(n-1) = 4*n^4 + 1 for n > 0. - Reinhard Zumkeller, Feb 12 2009
Prefaced with a "1" (1, 1, 5, 13, 25, 41, ...): a(n) = 2*n*(n-1)+1. - Doug Bell, Feb 27 2009
a(n) = sqrt((A056220(n)^2 + A056220(n+1)^2) / 2). - Doug Bell, Mar 08 2009
a(n) = floor(2*(n+1)^3/(n+2)). - Gary Detlefs, May 20 2010
a(n) = A000330(n) - A000330(n-2). - Keith Tyler, Aug 10 2010
a(n) = A069894(n)/2. - J. M. Bergot, Jun 11 2012
a(n) = 2*a(n-1) - a(n-2) + 4. - Ant King, Jun 12 2012
Sum_{n>=0} 1/a(n) = (Pi/2)*tanh(Pi/2) = 1.4406595199775... = A228048. - Ant King, Jun 15 2012
a(n) = A209297(2*n+1,n+1). - Reinhard Zumkeller, Jan 19 2013
a(n)^3 = A048395(n)^2 + A048395(-n-1)^2. - Vincenzo Librandi, Jan 19 2013
a(n) = A000217(2n+1) - n. - Ivan N. Ianakiev, Nov 08 2013
a(n) = A251599(3*n+1). - Reinhard Zumkeller, Dec 13 2014
a(n) = A101321(4,n). - R. J. Mathar, Jul 28 2016
From Ilya Gutkovskiy, Jul 30 2016: (Start)
a(n) = Sum_{k=0..n} A008574(k).
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = exp(-1) = A068985. (End)
a(n) = 4 * A000217(n) + 1. - Bruce J. Nicholson, Jul 10 2017
a(n) = A002522(n) + A005563(n) = A002522(n+1) + A005563(n-1). - Bruce J. Nicholson, Aug 05 2017
Sum_{n>=0} a(n)/n! = 7*e. Sum_{n>=0} 1/a(n) = A228048. - Amiram Eldar, Jun 20 2020
a(n) = A000326(n+1) + A000217(n-1). - Charlie Marion, Nov 16 2020
a(n) = Integral_{x=0..2n+2} |1-x| dx. - Pedro Caceres, Dec 29 2020
From Amiram Eldar, Feb 17 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = cosh(sqrt(3)*Pi/2)*sech(Pi/2).
Product_{n>=1} (1 - 1/a(n)) = Pi*csch(Pi)*sinh(Pi/2). (End)
a(n) = A001651(n+1) + 1 - A028242(n). - Charlie Marion, Apr 05 2022
a(n) = A016754(n) - A046092(n). - Leo Tavares, Sep 16 2022
For n>0, a(n) = A101096(n+2) / 30. - Andy Nicol, Feb 06 2025
From Rémi Guillaume, Apr 21 2025: (Start)
a(n) = (2*A003215(n)+1)/3.
a(n) = (4*A005448(n+1)-1)/3.
a(n) + a(n-1) = A001845(n) - A001845(n-1), for n >= 1.
a(n) = (A005917(n+1))/(2n+1). (End)

Partially edited by Joerg Arndt, Mar 11 2010

A002313 Primes congruent to 1 or 2 modulo 4; or, primes of form x^2 + y^2; or, -1 is a square mod p.

Original entry on oeis.org

2, 5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 97, 101, 109, 113, 137, 149, 157, 173, 181, 193, 197, 229, 233, 241, 257, 269, 277, 281, 293, 313, 317, 337, 349, 353, 373, 389, 397, 401, 409, 421, 433, 449, 457, 461, 509, 521, 541, 557, 569, 577, 593, 601, 613, 617

Offset: 1

N. J. A. Sloane

Or, primes p such that x^2 - p*y^2 represents -1.
Primes which are not Gaussian primes (meaning not congruent to 3 mod 4).
Every Fibonacci prime (with the exception of F(4) = 3) is in the sequence. If p = 2n+1 is the prime index of the Fibonacci prime, then F(2n+1) = F(n)^2 + F(n+1)^2 is the unique representation of the prime as sum of two squares. - Sven Simon, Nov 30 2003
Except for 2, primes of the form x^2 + 4y^2. See A140633. - T. D. Noe, May 19 2008
Primes p such that for all p > 2, p XOR 2 = p + 2. - Brad Clardy, Oct 25 2011
Greatest prime divisor of r^2 + 1 for some r. - Michel Lagneau, Sep 30 2012
Empirical result: a(n), as a set, compose the prime factors of the family of sequences produced by A005408(j)^2 + A005408(j+k)^2 = (2j+1)^2 + (2j+2k+1)^2, for j >= 0, and a given k >= 1 for each sequence, with the addition of the prime factors of k if not already in a(n). - Richard R. Forberg, Feb 09 2015
Primes such that when r is a primitive root then p-r is also a primitive root. - Emmanuel Vantieghem, Aug 13 2015
Primes of the form (x^2 + y^2)/2. Note that (x^2 + y^2)/2 = ((x+y)/2)^2 + ((x-y)/2)^2 = a^2 + b^2 with x = a + b and y = a - b. More generally, primes of the form (x^2 + y^2) / A001481(n) for every fixed n > 1. - Thomas Ordowski, Jul 03 2016
Numbers n such that ((n-2)!!)^2 == -1 (mod n). - Thomas Ordowski, Jul 25 2016
Primes p such that (p-1)!! == (p-2)!! (mod p). - Thomas Ordowski, Jul 28 2016
The product of 2 different terms (x^2 + y^2)(z^2 + v^2) = (xz + yv)^2 + (xv - yz)^2 is sum of 2 squares (A000404) because (xv - yz)^2 > 0. If x were equal to yz/v then (x^2 + y^2)/(z^2 + v^2) would be equal to ((yz/v)^2 + y^2)/(z^2 + v^2) = y^2/v^2 which is not possible because (x^2 + y^2) and (z^2 + v^2) are prime numbers. For example, (2^2 + 5^2)(4^2 + 9^2) = (2*4 + 5*9)^2 + (2*9 - 5*4)^2. - Jerzy R Borysowicz, Mar 21 2017

			13 is in the sequence since it is prime and 13 = 4*3 + 1.  Also 13 = 2^2 + 3^2.  And -1 is a square (mod 13): -1 + 2*13 = 25 = 5^2.  Of course, only the first term is congruent to 2 (mod 4). - _Michael B. Porter_, Jul 04 2016

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 872.
David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, p. 219, th. 251, 252.
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).

Zak Seidov, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Dario Alpern, Online program that calculates sum of squares representation
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
J. Todd, A problem on arc tangent relations, Amer. Math. Monthly, 56 (1949), 517-528.
Eric Weisstein's World of Mathematics, Fermat's 4n+1 Theorem
G. Xiao, Two squares
Index entries for Gaussian integers and primes

Apart from initial term, same as A002144. For values of x and y see A002330 and A002331.

Cf. A033203, A038873, A038874, A045331, A008784, A057129, A084163, A084165, A137351.

a002313 n = a002313_list !! (n-1)
a002313_list = filter ((`elem` [1,2]) . (`mod` 4)) a000040_list
-- Reinhard Zumkeller, Feb 04 2014

[p: p in PrimesUpTo(700) | p mod 4 in {1,2}]; // Vincenzo Librandi, Feb 18 2015

with(numtheory): for n from 1 to 300 do if ithprime(n) mod 4 = 1 or ithprime(n) mod 4 = 2 then printf(`%d,`,ithprime(n)) fi; od:
# alternative
A002313 := proc(n)
    option remember ;
    local a;
    if n = 1 then
        2;
    elif n = 2 then
        5;
    else
        for a from procname(n-1)+4 by 4 do
            if isprime(a) then
                return a ;
            end if;
        end do:
    end if;
end proc:
seq(A002313(n),n=1..100) ; # R. J. Mathar, Feb 01 2024

Select[ Prime@ Range@ 115, Mod[#, 4] != 3 &] (* Robert G. Wilson v *)
fQ[n_] := Solve[x^2 + 1 == n*y^2, {x, y}, Integers] == {}; Select[ Prime@ Range@ 115, fQ] (* Robert G. Wilson v, Dec 19 2013 *)

select(p->p%4!=3, primes(1000)) \\ Charles R Greathouse IV, Feb 11 2011

a(n) ~ 2n log n. - Charles R Greathouse IV, Jul 04 2016

a(n) = A002331(n)^2 + A002330(n)^2. See crossrefs. - Wolfdieter Lang, Dec 11 2016

More terms from Henry Bottomley, Aug 10 2000

More terms from James Sellers, Aug 22 2000

A007528 Primes of the form 6k-1.

Original entry on oeis.org

5, 11, 17, 23, 29, 41, 47, 53, 59, 71, 83, 89, 101, 107, 113, 131, 137, 149, 167, 173, 179, 191, 197, 227, 233, 239, 251, 257, 263, 269, 281, 293, 311, 317, 347, 353, 359, 383, 389, 401, 419, 431, 443, 449, 461, 467, 479, 491, 503, 509, 521, 557, 563, 569, 587

Offset: 1

N. J. A. Sloane

For values of k see A024898.
Also primes p such that p^q - 2 is not prime where q is an odd prime. These numbers cannot be prime because the binomial p^q = (6k-1)^q expands to 6h-1 some h. Then p^q-2 = 6h-1-2 is divisible by 3 thus not prime. - Cino Hilliard, Nov 12 2008
a(n) = A211890(3,n-1) for n <= 4. - Reinhard Zumkeller, Jul 13 2012
There exists a polygonal number P_s(3) = 3s - 3 = a(n) + 1. These are the only primes p with P_s(k) = p + 1, s >= 3, k >= 3, since P_s(k) - 1 is composite for k > 3. - Ralf Steiner, May 17 2018
From Bernard Schott, Feb 14 2019: (Start)
A theorem due to Andrzej Mąkowski: every integer greater than 161 is the sum of distinct primes of the form 6k-1. Examples: 162 = 5 + 11 + 17 + 23 + 47 + 59; 163 = 17 + 23 + 29 + 41 + 53. (See Sierpiński and David Wells.)
{2,3} Union A002476 Union {this sequence} = A000040.
Except for 2 and 3, all Sophie Germain primes are of the form 6k-1.
Except for 3, all the lesser of twin primes are also of the form 6k-1.
Dirichlet's theorem on arithmetic progressions states that this sequence is infinite. (End)
For all elements of this sequence p=6*k-1, there are no (x,y) positive integers such that k=6*x*y-x+y. - Pedro Caceres, Apr 06 2019

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 870.
A. Mąkowski, Partitions into unequal primes, Bull. Acad. Polon. Sci. Sér. Sci. Math. Astr. Phys. 8 (1960), 125-126.
Wacław Sierpiński, Elementary Theory of Numbers, p. 144, Warsaw, 1964.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, Revised edition, 1997, p. 127.

T. D. Noe, Table of n, a(n) for n = 1..1000
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
F. S. Carey, On some cases of the Solutions of the Congruence z^p^(n-1)=1, mod p, Proceedings of the London Mathematical Society, Volume s1-33, Issue 1, November 1900, Pages 294-312.
Amelia Carolina Sparavigna, The Pentagonal Numbers and their Link to an Integer Sequence which contains the Primes of Form 6n-1, Politecnico di Torino (Italy, 2021).
Amelia Carolina Sparavigna, Binary operations inspired by generalized entropies applied to figurate numbers, Politecnico di Torino (Italy, 2021).
Wikipedia, Dirichlet's theorem on arithmetic progressions.

Intersection of A016969 and A000040.
Cf. A003627, A010051, A117047, A132231, A214360, A057145.
Prime sequences A# (k,r) of the form k*n+r with 0 <= r <= k-1 (i.e., primes == r (mod k), or primes p with p mod k = r) and gcd(r,k)=1: A000040 (1,0), A065091 (2,1), A002476 (3,1), A003627 (3,2), A002144 (4,1), A002145 (4,3), A030430 (5,1), A045380 (5,2), A030431 (5,3), A030433 (5,4), A002476 (6,1), this sequence (6,5), A140444 (7,1), A045392 (7,2), A045437 (7,3), A045471 (7,4), A045458 (7,5), A045473 (7,6), A007519 (8,1), A007520 (8,3), A007521 (8,5), A007522 (8,7), A061237 (9,1), A061238 (9,2), A061239 (9,4), A061240 (9,5), A061241 (9,7), A061242 (9,8), A030430 (10,1), A030431 (10,3), A030432 (10,7), A030433 (10,9), A141849 (11,1), A090187 (11,2), A141850 (11,3), A141851 (11,4), A141852 (11,5), A141853 (11,6), A141854 (11,7), A141855 (11,8), A141856 (11,9), A141857 (11,10), A068228 (12,1), A040117 (12,5), A068229 (12,7), A068231 (12,11).
Cf. A034694 (smallest prime == 1 (mod n)).
Cf. A038700 (smallest prime == n-1 (mod n)).
Cf. A038026 (largest possible value of smallest prime == r (mod n)).
Cf. A001359 (lesser of twin primes), A005384 (Sophie Germain primes).
Cf. A048265, A324076.

Filtered(List([1..100],n->6*n-1),IsPrime); # Muniru A Asiru, May 19 2018

a007528 n = a007528_list !! (n-1)
a007528_list = [x | k <- [0..], let x = 6 * k + 5, a010051' x == 1]
-- Reinhard Zumkeller, Jul 13 2012

select(isprime,[seq(6*n-1,n=1..100)]); # Muniru A Asiru, May 19 2018

Select[6 Range[100]-1,PrimeQ]  (* Harvey P. Dale, Feb 14 2011 *)

forprime(p=2, 1e3, if(p%6==5, print1(p, ", "))) \\ Charles R Greathouse IV, Jul 15 2011

forprimestep(p=5,1000,6, print1(p", ")) \\ Charles R Greathouse IV, Mar 03 2025

A003627 \ {2}. - R. J. Mathar, Oct 28 2008
Conjecture: Product_{n >= 1} ((a(n) - 1) / (a(n) + 1)) * ((A002476(n) + 1) / (A002476(n) - 1)) = 3/4. - Dimitris Valianatos, Feb 11 2020
From Vaclav Kotesovec, May 02 2020: (Start)
Product_{k>=1} (1 - 1/a(k)^2) = 9*A175646/Pi^2 = 1/1.060548293.... =4/(3*A333240).
Product_{k>=1} (1 + 1/a(k)^2) = A334482.
Product_{k>=1} (1 - 1/a(k)^3) = A334480.
Product_{k>=1} (1 + 1/a(k)^3) = A334479. (End)
Legendre symbol (-3, a(n)) = -1 and (-3, A002476(n)) = +1, for n >= 1. For prime 3 one sets (-3, 3) = 0. - Wolfdieter Lang, Mar 03 2021

A007519 Primes of form 8n+1, that is, primes congruent to 1 mod 8.

Original entry on oeis.org

17, 41, 73, 89, 97, 113, 137, 193, 233, 241, 257, 281, 313, 337, 353, 401, 409, 433, 449, 457, 521, 569, 577, 593, 601, 617, 641, 673, 761, 769, 809, 857, 881, 929, 937, 953, 977, 1009, 1033, 1049, 1097, 1129, 1153, 1193, 1201, 1217, 1249, 1289, 1297, 1321, 1361

Offset: 1

N. J. A. Sloane, Robert G. Wilson v

Discriminant is 32, class is 2. Binary quadratic forms ax^2 + bxy + cy^2 have discriminant d = b^2 - 4ac and gcd(a, b, c) = 1.
Integers n (n > 9) of form 4k + 1 such that binomial(n-1, (n-1)/4) == 1 (mod n) - Benoit Cloitre, Feb 07 2004
Primes of the form x^2 + 8y^2. - T. D. Noe, May 07 2005
Also primes of the form x^2 + 16y^2. See A140633. - T. D. Noe, May 19 2008
Is this the same sequence as A141174?
Being a subset of A001132 and also a subset of A038873, this is also a subset of the primes of the form u^2 - 2v^2. - Tito Piezas III, Dec 28 2008
These primes p are only which possess the property: for every integer m from interval [0, p) with the Hamming distance D(m, p) = 2, there exists an integer h from (m, p) with D(m, h) = 2. - Vladimir Shevelev, Apr 18 2012
Primes p such that p XOR 6 = p + 6. - Brad Clardy, Jul 22 2012
Odd primes p such that -1 is a 4th power mod p. - Eric M. Schmidt, Mar 27 2014
There are infinitely many primes of this form. See Brubaker link. - Alonso del Arte, Jan 12 2017
These primes split in Z[sqrt(2)]. For example, 17 = (-1)(1 - 3*sqrt(2))(1 + 3*sqrt(2)). This is also true of primes of the form 8n - 1. - Alonso del Arte, Jan 26 2017

			a(1) = 17 = 2 * 8 + 1 = (10001)_2. All numbers m from [0, 17) with the Hamming distance D(m, 17) = 2 are 0, 3, 5, 9. For m = 0, we can take h = 3, since 3 is drawn from (0, 17) and D(0, 3) = 2; for m = 3, we can take h = 5, since 5 from (3, 17) and D(3, 5) = 2; for m = 5, we can take h = 6, since 6 from (5, 17) and D(5, 6) = 2; for m = 9, we can take h = 10, since 10 is drawn from (9, 17) and D(9, 10) = 2. - _Vladimir Shevelev_, Apr 18 2012

Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 870.
Z. I. Borevich and I. R. Shafarevich, Number Theory.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 261.

Ray Chandler, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Ben Brubaker, 18.781, Fall 2007 Problem Set 5: Solutions to Selected Problems, MIT (2007).
Peter Luschny, Binary Quadratic Forms
Jorma K. Merikoski, Pentti Haukkanen, and Timo Tossavainen, The congruence x^n = -a^n (mod m): Solvability and related OEIS sequences, Notes. Num. Theor. Disc. Math. (2024) Vol. 30, No. 3, 516-529. See p. 521.
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
D. B. Zagier, Zetafunktionen und quadratische Körper, Springer, 1981.

Subsequence of A017077 and of A038873.
Cf. A139643. Complement in primes of A154264. Cf. A042987.
Cf. A065091, A002144, A094407, A133870, A142925, A208177, A208178, A076339.
Cf. A038872 (d = 5). A038873 (d = 8). A068228, A141123 (d = 12). A038883 (d = 13). A038889 (d = 17). A141111, A141112 (d = 65).
Cf. also A242663.
For a list of sequences giving numbers and/or primes represented by binary quadratic forms, see the "Binary Quadratic Forms and OEIS" link.

a007519 n = a007519_list !! (n-1)
a007519_list = filter ((== 1) . a010051) [1,9..]
-- Reinhard Zumkeller, Mar 06 2012

[p: p in PrimesUpTo(2000) | p mod 8 eq 1 ]; // Vincenzo Librandi, Aug 21 2012

Select[1 + 8 Range@ 170, PrimeQ] (* Robert G. Wilson v *)

forprime(p=2,1e4,if(p%8==1,print1(p", "))) \\ Charles R Greathouse IV, Jun 16 2011

forprimestep(p=17,10^4,8, print1(p", ")) \\ Charles R Greathouse IV, Jul 17 2024

lista(nn)= my(vpr = []); for (x = 0, nn, y = 0; while ((v = x^2+6*x*y+y^2) < nn, if (isprime(v), if (! vecsearch(vpr, v), vpr = concat(vpr, v); vpr = vecsort(vpr););); y++;);); vpr; \\ Michel Marcus, Feb 01 2014

A007519_upto(N, start=1)=select(t->t%8==1,primes([start,N]))
#A7519=A007519_upto(10^5)
A007519(n)={while(#A7519A007519_upto(N*3\2, N+1))); A7519[n]} \\ M. F. Hasler, May 22 2025

# uses[binaryQF]
# The function binaryQF is defined in the link 'Binary Quadratic Forms'.
Q = binaryQF([1, 4, -4])
print(Q.represented_positives(1361, 'prime'))  # Peter Luschny, Jan 26 2017

A031971 a(n) = Sum_{k=1..n} k^n.

Original entry on oeis.org

1, 5, 36, 354, 4425, 67171, 1200304, 24684612, 574304985, 14914341925, 427675990236, 13421957361110, 457593884876401, 16841089312342855, 665478473553144000, 28101527071305611528, 1262899292504270591313, 60182438244917445266889, 3031284048960901518840700

Offset: 1

Chris du Feu (chris(AT)beckingham0.demon.co.uk)

From Alexander Adamchuk, Jul 21 2006: (Start)
p^(3k - 1) divides a(p^k) for prime p > 2 and k = 1, 2, 3, 4, ... or p^2 divides a(p) for prime p > 2. p^5 divides a(p^2) for prime p > 2. p^8 divides a(p^3) for prime p > 2. p^11 divides a(p^4) for prime p > 2.
p^2 divides a(2p) for prime p > 3. p^3 divides a(3p) for prime p > 2. p^2 divides a(4p) for prime p > 5. p^3 divides a(5p) for prime p > 3. p^2 divides a(6p) for prime p > 7.
p divides a(2p - 1) for all prime p. p^3 divides a(2p^2 - 1) for all prime p. p^5 divides a(2p^3 - 1) for all prime p.
p divides a((p - 1)/2) for p = 5, 13, 17, 29, 37, 41, 53, 61, ... = A002144 Pythagorean primes: primes of form 4n + 1.
(End)
If p prime then a(p-1) == -1 (mod p) [see De Koninck & Mercier reference]. Example: for p = 7, a(6) = 67171 = 7 * 9596 - 1. - Bernard Schott, Mar 06 2020

J.-M. De Koninck et A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 327 pp. 48-200, Ellipses, Paris (2004).
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 21.

T. D. Noe, Table of n, a(n) for n = 1..100
J. M. Grau, A. M. Oller-Marcen, and J. Sondow, On the congruence 1^m + 2^m + ... + m^m == n (mod m) with n|m, arXiv 1309.7941 [math.NT], 2013-2014.

A diagonal of array A103438.
Cf. A002144, A014117, A089072, A229303, A256016, A261490.
For a(n) mod n see A182398.

a031971 = sum . a089072_row  -- Reinhard Zumkeller, Mar 18 2013

[&+[(k)^n: k in [0..n]]: n in [1..30]]; // Vincenzo Librandi, Apr 18 2011

a := n->sum('i'^n,'i'=1..n);
# alternative
A031971 := proc(n)
    (bernoulli(n+1,n+1)-bernoulli(n+1))/(n+1) ;
end proc: # R. J. Mathar, May 10 2013

Table[Zeta[-n] - Zeta[-n, n + 1], {n, 25}] (* Alexander Adamchuk, Jul 21 2006 *)
Table[Total[Range[n]^n], {n,25}] (* T. D. Noe, Apr 19 2011 *)
Table[HarmonicNumber[n, -n], {n, 1, 25}] (* Jean-François Alcover, Apr 09 2015 *)

a(n)=sum(k=1,n,k^n) \\ Charles R Greathouse IV, Jun 05 2015

from sympy import harmonic
def A031971(n):
    return harmonic(n,-n) # Chai Wah Wu, Feb 15 2020

a(n) is asymptotic to (e/(e - 1))*n^n. - Benoit Cloitre, Dec 17 2003
a(n) = zeta(-n) - zeta(-n, n + 1), where zeta(s) is the Riemann zeta function and zeta(s, a) is the Hurwitz zeta function, a generalization of the Riemann zeta function. - Alexander Adamchuk, Jul 21 2006
a(n) == 1 (mod n) <==> n is in A014117 = 1, 2, 6, 42, 1806 (see the link "On the congruence ..."). - Jonathan Sondow, Oct 18 2013
a(A054377(n)) = A233045(n). - Jonathan Sondow, Dec 11 2013
a(n) = n! * [x^n] exp(x)*(exp(n*x) - 1)/(exp(x) - 1). - Ilya Gutkovskiy, Apr 07 2018
a(n) ~ ((e*n+1)/((e-1)*(n+1))) * n^n. - N. J. A. Sloane, Oct 13 2018, based on email from Claude F. Leibovici who claims this is slightly better than Cloitre's version when n is small.

A122045 Euler (or secant) numbers E(n).

Original entry on oeis.org

1, 0, -1, 0, 5, 0, -61, 0, 1385, 0, -50521, 0, 2702765, 0, -199360981, 0, 19391512145, 0, -2404879675441, 0, 370371188237525, 0, -69348874393137901, 0, 15514534163557086905, 0, -4087072509293123892361, 0, 1252259641403629865468285, 0, -441543893249023104553682821

Offset: 0

Roger L. Bagula, Sep 13 2006

sign

The convention in the OEIS is that the alternate zeros are normally omitted in such sequences. See A000364 for the official version of this sequence.

Odd primes p such that p | E(p-1) are primes p == 1 (mod 4), A002144. Conjecture: odd composites m such that m | E(m-1) are Carmichael numbers m such that p == 1 (mod 4) for every prime p|m, A265237. - Thomas Ordowski, Feb 06 2020

			G.f. = 1 - x^2 + 5*x^4 - 61*x^6 + 1385*x^8 - 50521*x^10 + 2702765*x^12 + ...

Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 5, page 41.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Fernando Barbero G., Juan Margalef-Bentabol, and Eduardo J. S. Villaseñor, A two-sided Faulhaber-like formula involving Bernoulli polynomials, arXiv:2002.00550 [math.NT], 2020.
Marilena Barnabei, Flavio Bonetti, Niccolò Castronuovo, and Matteo Silimbani, Ascending runs in permutations and valued Dyck paths, Ars Mathematica Contemporanea (2019) Vol. 16, No. 2, 445-463.
Hacène Belbachir and Yassine Otmani, Quadrinomial-Like Versions for Wolstenholme, Morley and Glaisher Congruences, Integers (2023) Vol. 23.
Fatima Zohra Bensaci, Rachid Boumahdi, and Laala Khaldi, Finite Sums Involving Fibonacci and Lucas Numbers, J. Int. Seq. (2024). See p. 2.
Michael D. Hirschhorn, Binomial Identities and Congruences for Euler Numbers, Fibonacci Quart. 53 (2015), no. 4, 319-322.
Guodong Liu, Generating functions and generalized Euler numbers, Proc. Japan Acad. Ser. A Math. Sci., Volume 84, Number 2 (2008), 29-34. See p 32.
F. Luca, A. Pizarro-Madariaga, and C. Pomerance, On the counting function of irregular primes, 2014.
Emanuele Munarini, Two-Parameter Identities for q-Appell Polynomials, Journal of Integer Sequences, Vol. 26 (2023), Article 23.3.1.
A. Niedermaier and J. Remmel, Analogues of up-down permutations for colored permutations, J. Int. Seq. 13 (2010) 10.5.6.
Bruce E. Sagan, Generalized Euler numbers and ordered set partitions, arXiv:2501.07692 [math.NT], 2025. See p. 2.
T. J. Stieltjes, Sur la réduction en fraction continue d'une série procédant suivant les puissances descendantes d'une variable, Ann. Fac. Sc. Toulouse 3 (1889) 1-17.

Cf. A000364, A002144, A028296, A265237, A331978, A350972.

m:=35; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( 1/Cosh(x) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 13 2020

seq(euler(n) , n=0..31); # Zerinvary Lajos, Mar 15 2009
P := proc(n,x) option remember; if n = 0 then 1 else
   (n*x+(1/2)*(1-x))*P(n-1,x)+x*(1-x)*diff(P(n-1,x),x); expand(%) fi end:
A122045 := n -> (-1)^n*subs(x=-1, P(n,x)):
seq(A122045(n), n=0..30);  # Peter Luschny, Mar 07 2014
ptan := proc(n) option remember; if irem(n, 2) = 1 then 0 else
    -add(`if`(k=0, 1, binomial(n, k)*ptan(n - k)), k = 0..n-1,2) fi end:
A122045 := n -> ifelse(n = 0, 1, ptan(n)):
seq(A122045(n), n = 0..30); # Peter Luschny, Jun 06 2022

Table[EulerE[n], {n, 0, 30}]
Range[0, 30]! CoefficientList[ Series[ Sech[x], {x, 0, 30}], x] (* Robert G. Wilson v, Aug 08 2018 *)
a[0] := 1; a[n_?OddQ] := 0; a[n_?EvenQ] := a[n] = -Sum[a[k] Binomial[n, k], {k, 0, n - 1, 2}]; Map[a, Range[0, 30]] (* Oliver Seipel, May 20 2024 *)

a[n]:=if n<2 then 1-n else sum(-a[n-2*k]*binomial(n,2*k),k,1,floor(n/2));
makelist(a[n],n,0,50); /* Tani Akinari, Sep 15 2023 */

x='x+O('x^66); Vec(serlaplace(1/cosh(x))) \\ Joerg Arndt, Mar 10 2014

a(n) = 2^n*2^(n+1)*(subst(bernpol(n+1,x), x, 3/4) - subst(bernpol(n+1,x), x, 1/4))/(n+1); \\ Michel Marcus, May 20 2017

from sympy import bernoulli as B
def a(n): return int(2**n*2**(n + 1)*(B(n + 1, 3/4) - B(n + 1, 1/4))/(n + 1))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 24 2017, after PARI code by Michel Marcus

from functools import cache
from math import comb as binomial
@cache
def ptan(n):   # see also A331978 and A350972.
    return (0 if n % 2 == 1 else
    -sum(binomial(n,k) * ptan(n-k) if k > 0 else 1 for k in range(0, n-1, 2)))
def A122045(n): return 1 if n == 0 else ptan(n)
print([A122045(n) for n in range(31)])  # Peter Luschny, Jun 06 2022

[euler_number(i) for i in range(31)] # Zerinvary Lajos, Mar 15 2009

E.g.f.: sech(x). - Michael Somos, Mar 11 2014
a(n) = (Sum_{k>=0} (-1)^k*(2*k+1)^n)*2. - Gottfried Helms, Mar 09 2012
From Sergei N. Gladkovskii, Oct 14 2012 - Oct 13 2013: (Start)
Continued fractions:
G.f.: 1/U(0) where U(k) = 1 - x + x*(k+1)/(1 - x*(k+1)/U(k+1)).
G.f.: 1/U(0) where U(k) = 1 - x^2 + x^2*(2*k+1)*(2*k+2)/(1 + x^2*(2*k+1)*(2*k+2)/ U(k+1)).
E.g.f.: (1-x)/U(0) where U(k) = 1 - x/(1 - x/(x - (2*k+1)*(2*k+2)/U(k+1))).
E.g.f.: 1 - x^2/U(0) where U(k) = (2*k+1)*(2*k+2) + x^2 - x^2*(2*k+1)*(2*k+2)/U(k+1).
E.g.f.: 1/U(0) where U(k) = 1 + x^2/(2*(2*k+1)*(4*k+1) - 2*x^2*(2*k+1)*(4*k+1)/(x^2 + 4*(4*k+3)*(k+1)/U(k+1))).
E.g.f.: (2 + x^4/(U(0)*(x^2-2) - 2))/(2-x^2) where U(k) = 4*k + 4 + 1/(1 + x^2/(2 - x^2 + (2*k+3)*(2*k+4)/U(k+1))).
G.f.: 1/G(0) where G(k) = 1 + x^2*(2*k+1)^2/(1 + x^2*(2*k+2)^2/G(k+1)) (due to T. J. Stieltjes).
G.f.: 1/S(0) where S(k) = 1 + x^2*(k+1)^2/S(k+1) (due to T. J. Stieltjes).
G.f.: 1 - x/(1+x) + x/(1+x)/Q(0) where Q(k) = 1 - x + x*(k+2)/(1 - x*(k+1)/Q(k+1)).
G.f.: -(1/x)/Q(0) where Q(k) = -1/x + (k+1)^2/Q(k+1) (due to T. J. Stieltjes).
G.f.: (1/(1-x))/Q(0) + 1/(1-x) where Q(k) = 1 - 1/x + (k+1)*(k+2)/Q(k+1).
G.f.: (x/(x-1))/Q(0) + 1/(1-x) where Q(k) = 1 - x + x^2*(k+1)*(k+2)/Q(k+1).
G.f.: 1 - x/(1+x) + (x/(1+x))/Q(0) where Q(k) = 1 + x + (k+1)*(k+2)*x^2/Q(k+1).
E.g.f.: 1 - T(0)*x^2/(2+x^2) where T(k) = 1 - x^2*(2*k+1)*(2*k+2)/(x^2*(2*k+1)*(2*k+2) - ((2*k+1)*(2*k+2) + x^2)*((2*k+3)*(2*k+4) + x^2)/T(k+1)).
G.f.: T(0) where T(k) = 1 - x^2*(k+1)^2/(x^2*(k+1)^2 + 1/T(k+1)). (End)
a(n) = 2^(2*n+1)*(zeta(-n,1/4) - zeta(-n,3/4)), where zeta(a, z) is the generalized Riemann zeta function. - Peter Luschny, Mar 11 2015
a(n) = 2^n*(2^(n+1)/(n+1))*(B(n+1, 3/4) - B(n+1, 1/4)) where B(n,x) is the n-th Bernoulli polynomial. See Liu link. - Michel Marcus, May 20 2017 [This is the same as: a(n) = -4^(n+1)*B(n+1, 1/4)*((n+1) mod 2)/(n+1). Peter Luschny, Oct 30 2020]
a(n) = 2*Im(PolyLog(-n, I)). - Peter Luschny, Sep 29 2020
a(4n) == 5 (mod 60) and a(4n+2) == -1 (mod 60). See Hirschhorn. - Michel Marcus, Jan 11 2022
For n > 1, a(n) = -Sum_{k=1..floor(n/2)} a(n-2*k)*binomial(n,2*k). - Tani Akinari, Sep 15 2023

Edited by N. J. A. Sloane, Sep 17 2006

A009003 Hypotenuse numbers (squares are sums of 2 nonzero squares).

Original entry on oeis.org

5, 10, 13, 15, 17, 20, 25, 26, 29, 30, 34, 35, 37, 39, 40, 41, 45, 50, 51, 52, 53, 55, 58, 60, 61, 65, 68, 70, 73, 74, 75, 78, 80, 82, 85, 87, 89, 90, 91, 95, 97, 100, 101, 102, 104, 105, 106, 109, 110, 111, 113, 115, 116, 117, 119, 120, 122, 123, 125, 130, 135, 136, 137, 140

Offset: 1

David W. Wilson

nonn

Multiples of Pythagorean primes A002144 or of primitive Pythagorean triangles' hypotenuses A008846. - Lekraj Beedassy, Nov 12 2003
This is exactly the sequence of positive integers with at least one prime divisor of the form 4k + 1. Compare A072592. - John W. Layman, Mar 12 2008 and Franklin T. Adams-Watters, Apr 26 2009
Circumradius R of the triangles such that the area, the sides and R are integers. - Michel Lagneau, Mar 03 2012
The 2 squares summing to a(n)^2 cannot be equal because sqrt(2) is not rational. - Jean-Christophe Hervé, Nov 10 2013
Closed under multiplication. The primitive elements are those with exactly one prime divisor of the form 4k + 1 with multiplicity one, which are also those for which there exists a unique integer triangle = A084645. - Jean-Christophe Hervé, Nov 11 2013
a(n) are numbers whose square is the mean of two distinct nonzero squares. This creates 1-to-1 mapping between a Pythagorean triple and a "Mean" triple. If the Pythagorean triple is written, abnormally, as {j, k, h} where j^2 +(j+k)^2 = h^2, and h = a(n), then the corresponding "Mean" triple with the same h is {k, 2j, h} where (k^2 + (k+2j)^2)/2 = h^2. For example for h = 5, the Pythagorean triple is {3, 1, 5} and the Mean triple is {1, 6, 5}. - Richard R. Forberg, Mar 01 2015
Integral side lengths of rhombuses with integral diagonals p and q (therefore also with integral areas A because A = pq/2 is some multiple of 24). No such rhombuses are squares. - Rick L. Shepherd, Apr 09 2017
Conjecture: these are bases n in which exists an n-adic integer x satisfying x^5 = x, and 5 is the smallest k>1 such that x^k =x (so x^2, x^3 and x^4 are not x). Example: the 10-adic integer x = ...499879186432 (A120817) satisfies x^5 = x, and x^2, x^3, and x^4 are not x, so 10 is in this sequence. See also A120817, A210850 and A331548. - Patrick A. Thomas, Mar 01 2020
Didactic comment: When students solve a quadratic equation a*x^2 + b*x + c = 0 (a, b, c: integers) with the solution formula, they often make the mistake of calculating b^2 + 4*a*c instead of b^2 - 4*a*c (especially if a or c is negative). If the root then turns out to be an integer, they feel safe. This sequence lists the absolute values of b for which this error can happen. Reasoning: With p^2 = b^2 - 4*a*c and q^2 = b^2 + 4*a*c it follows by addition immediately that p^2 + q^2 = 2*b^2. If 4*a*c < 0, let p = x + y and q = x - y. If 4*a*c > 0, let p = x - y and q = x + y. In both cases follows that y^2 + x^2 = b^2. So every Pythagorean triple gives an absolute value of b for which this error can occur. Example: From (y, x, b) = (3, 4, 5) follows (q^2, b^2, p^2) = (1, 25, 49) or (p^2, b^2, q^2) = (1, 25, 49) with abs(4*a*c) = 24. - Felix Huber, Jul 22 2023
Conjecture: Numbers m such that the limit: Limit_{s->1} zeta(s)*Sum_{k=1..m} [k|m]*A008683(k)*(i^k)/(k^(s - 1)) exists, which is equivalent to numbers m such that abs(Sum_{k=1..m} [k|m]*A008683(k)*(i^k)) = 0. - Mats Granvik, Jul 06 2024

Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 98-104.

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
R. Chapman, Pythagorean triples and sums of squares
Steven R. Finch, Landau-Ramanujan Constant [Broken link]
Steven R. Finch, Landau-Ramanujan Constant [From the Wayback machine]
Ron Knott, Pythagorean Triples and Online Calculators
J. Pahikkala, On contraharmonic mean and Pythagorean triples, Elemente der Mathematik, 65:2 (2010), 62-67.
Patrick A. Thomas, Solutions of x^5=x up to base 100
Index entries for sequences related to sums of squares

Cf. A000404 (sums of 2 squares), A004431 (sums of 2 distinct squares), A009000 (hypotenuse numbers with repetition), A072592, A004613, A187811.
Complement of A004144. Primes in this sequence give A002144. Same as A146984 (integer contraharmonic means) as sets - see Pahikkala 2010, Theorem 5.
Cf. A083025, A084645 (primitive elements), A084646, A084647, A084648, A084649, A006339.

import Data.List (findIndices)
a009003 n = a009003_list !! (n-1)
a009003_list = map (+ 1) $ findIndices (> 0) a005089_list
-- Reinhard Zumkeller, Jan 07 2013

isA009003 := proc(n)
    local p;
    for p in numtheory[factorset](n) do
        if modp(p,4) = 1 then
            return true;
        end if;
    end do:
    false;
end proc:
for n from 1 to 200 do
    if isA009003(n) then
        printf("%d,",n) ;
    end if;
end do: # R. J. Mathar, Nov 17 2014

f[n_] := Module[{k = 1}, While[(n - k^2)^(1/2) != IntegerPart[(n - k^2)^(1/2)], k++; If[2 * k^2 >= n, k = 0; Break[]]]; k]; A009003 = {}; Do[If[f[n^2] > 0, AppendTo[A009003, n]], {n, 3, 100}]; A009003 (* Vladimir Joseph Stephan Orlovsky, Jun 15 2009 *)
Select[Range[200], Length[PowersRepresentations[#^2, 2, 2]] > 1 &] (* Alonso del Arte, Feb 11 2014 *)

is_A009003(n)=setsearch(Set(factor(n)[,1]%4),1)  \\ M. F. Hasler, May 27 2012

list(lim)=my(v=List(),u=vectorsmall(lim\=1)); forprimestep(p=5,lim,4, forstep(n=p,lim,p, u[n]=1)); for(i=5,lim, if(u[i], listput(v,i))); u=0; Vec(v) \\ Charles R Greathouse IV, Jan 13 2022

from itertools import count, islice
from sympy import primefactors
def A009003_gen(): # generator of terms
    return filter(lambda n:any(map(lambda p: p % 4 == 1,primefactors(n))),count(1))
A009003_list = list(islice(A009003_gen(),20)) # Chai Wah Wu, Jun 22 2022

A005089(a(n)) > 0. - Reinhard Zumkeller, Jan 07 2013
a(n) ~ n. - Charles R Greathouse IV, Jan 13 2022
a(n) = sqrt(n-th square in A000404), where A000404 lists the sums of two nonzero squares. - M. F. Hasler, Jun 20 2025

Definition edited by Jean-Christophe Hervé, Nov 10 2013

A004431 Numbers that are the sum of 2 distinct nonzero squares.

Original entry on oeis.org

5, 10, 13, 17, 20, 25, 26, 29, 34, 37, 40, 41, 45, 50, 52, 53, 58, 61, 65, 68, 73, 74, 80, 82, 85, 89, 90, 97, 100, 101, 104, 106, 109, 113, 116, 117, 122, 125, 130, 136, 137, 145, 146, 148, 149, 153, 157, 160, 164, 169, 170, 173, 178, 180, 181, 185, 193, 194, 197

Offset: 1

N. J. A. Sloane

nonn

Numbers whose prime factorization includes at least one prime congruent to 1 mod 4 and any prime factor congruent to 3 mod 4 has even multiplicity. - Franklin T. Adams-Watters, May 03 2006
Reordering of A055096 by increasing values and without repetition. - Paul Curtz, Sep 08 2008
A063725(a(n)) > 1. - Reinhard Zumkeller, Aug 16 2011
The square of these numbers is also the sum of two nonzero squares, so this sequence is a subsequence of A009003. - Jean-Christophe Hervé, Nov 10 2013
Closed under multiplication. Primitive elements are those with exactly one prime factor congruent to 1 mod 4 with multiplicity one (A230779). - Jean-Christophe Hervé, Nov 10 2013
From Bob Selcoe, Mar 23 2016: (Start)
Numbers c such that there is d < c, d >= 1 where c + d and c - d are square. For example, 53 + 28 = 81, 53 - 28 = 25.
Given a prime p == 1 mod 4, a term appears if and only if it is of the form p^i, p*2^j or p*k^2 {i,j,k >= 1}, or a product of any combination of these forms. Therefore, the products of any terms to any powers also are terms. For example, p(1) = 5 and p(2) = 13 so term 45 appears because 5*3^2 = 45 and term 416 appears because 13*2^5 = 416; therefore 45 * 416 = 18720 appears, as does 45^3 * 416^11 = 18720^3 * 416^8.
Numbers of the form j^2 + 2*j*k + 2*k^2 {j,k >= 1}. (End)
Suppose we have a term t = x^2 + y^2. Then s^2*t = (s*x)^2 + (s*y)^2 is a term for any s > 0. Also 2*t = (y+x)^2 + (x-y)^2 is a term. It follows that q*s^2*t is a term for any s > 0 and q=1 or 2. Examples: 2*7^2*26 = 28^2 + 42^2; 6^2*17 = 6^2 + 24^2. - Jerzy R Borysowicz, Aug 11 2017
To find terms up to some upper bound u, we can search for x^2 + y^2 = t where x is odd and y is even. Then we add all numbers of the form 2^m * t <= u and then remove duplicates. - David A. Corneth, Oct 04 2017
From Bernard Schott, Apr 13 2022: (Start)
The 5th comment "Closed under multiplication" can be proved with Brahmagupta's identity: (a^2+b^2) * (c^2+d^2)  = (ac + bd)^2 + (ad - bc)^2.
The subsequence of primes is A002144. (End)

			53 = 7^2 + 2^2, so 53 is in the sequence.

T. D. Noe, Table of n, a(n) for n = 1..10000.
Wikipedia, Brahmagupta's identity.
Index entries for sequences related to sums of squares.

Complement of A004439.

Cf. A000404, A002144, A007692, A009000, A009003, A009177, A024507, A081324, A081325, A118882, A230779.

import Data.List (findIndices)
a004431 n = a004431_list !! (n-1)
a004431_list = findIndices (> 1) a063725_list
-- Reinhard Zumkeller, Aug 16 2011

isA004431 := proc(n)
    local a,b ;
    for a from 2 do
        if a^2>= n then
            return false;
        end if;
        b := n -a^2 ;
        if b < 1 then
            return false ;
        end if;
        if issqr(b) then
            if ( sqrt(b) <> a ) then
                return true;
            end if;
        end if;
    end do:
    return false;
end proc:
A004431 := proc(n)
    option remember ;
    local a;
    if n = 1 then
        5;
    else
        for a from procname(n-1)+1 do
            if isA004431(a) then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, Jan 29 2013

A004431 = {}; Do[a = 2 m * n; b = m^2 - n^2; c = m^2 + n^2; AppendTo[A004431, c], {m, 100}, {n, m - 1}]; Take[Union@A004431, 63] (* Robert G. Wilson v, May 02 2009 *)
Select[Range@ 200, Length[PowersRepresentations[#, 2, 2] /. {{0, } -> Nothing, {a, b_} /; a == b -> Nothing}] > 0 &] (* Michael De Vlieger, Mar 24 2016 *)

select( isA004431(n)={n>1 && vecmin((n=factor(n)%4)[,1])==1 && ![f[1]>2 && f[2]%2 | f<-n~]}, [1..199]) \\ M. F. Hasler, Feb 06 2009, updated Nov 24 2019

is(n)=if(n<5, return(0)); my(f=factor(n)%4); if(vecmin(f[, 1])>1, return(0)); for(i=1, #f[, 1], if(f[i, 1]==3 && f[i, 2]%2, return(0))); 1
for(n=1, 1e3, if(is(n), print1(n, ", "))) \\ Altug Alkan, Dec 06 2015

upto(n) = {my(res = List(), s); forstep(i=1, sqrtint(n), 2, forstep(j = 2, sqrtint(n - i^2), 2, listput(res, i^2 + j^2))); s = #res; for(i = 1, s, t = res[i]; for(e = 1, logint(n \ res[i], 2), listput(res, t<<=1))); listsort(res, 1); res} \\ David A. Corneth, Oct 04 2017

def aupto(limit):
  s = [i*i for i in range(1, int(limit**.5)+2) if i*i < limit]
  s2 = set(a+b for i, a in enumerate(s) for b in s[i+1:] if a+b <= limit)
  return sorted(s2)
print(aupto(197)) # Michael S. Branicky, May 10 2021

A046080 a(n) is the number of integer-sided right triangles with hypotenuse n.

Original entry on oeis.org

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

Offset: 1

Eric W. Weisstein

nonn

Or number of ways n^2 can be written as the sum of two positive squares: a(5) = 1: 3^2 + 4^2 = 5^2; a(25) = 2: 7^2 + 24^2 = 15^2 + 20^2 = 25^2. - Alois P. Heinz, Aug 01 2019

A. H. Beiler, Recreations in the Theory of Numbers, New York: Dover, pp. 116-117, 1966.

Stanislav Sykora, Table of n, a(n) for n = 1..20000
Ron Knott, Pythagorean Triples and Online Calculators
F. Richman, Pythagorean Triples
A. Tripathi, On Pythagorean triples containing a fixed integer, Fib. Q., 46/47 (2008/2009), 331-340. See Theorem 7.
Eric Weisstein's World of Mathematics, Pythagorean Triple

First differs from A083025 at n=65.
Cf. A000290, A006339, A024362, A046079, A046081, A009000.
A088111 gives records; A088959 gives where records occur.
Cf. A046109, A063725.
Cf. A004144, A084647, A084648, A084649.
Partial sums: A224921.

f:= proc(n) local F,t;
  F:= select(t -> t[1] mod 4 = 1, ifactors(n)[2]);
  1/2*(mul(2*t[2]+1, t=F)-1)
end proc:
map(f, [$1..100]); # Robert Israel, Jul 18 2016

a[1] = 0; a[n_] := With[{fi = Select[ FactorInteger[n], Mod[#[[1]], 4] == 1 & ][[All, 2]]}, (Times @@ (2*fi+1)-1)/2]; Table[a[n], {n, 1, 99}] (* Jean-François Alcover, Feb 06 2012, after first formula *)

a(n)={my(m=0,k=n,n2=n*n,k2,l2);
while(1,k=k-1;k2=k*k;l2=n2-k2;if(l2>k2,break);if(issquare(l2),m++));return(m)} \\ brute force, Stanislav Sykora, Mar 18 2015

{a(n) = if( n<1, 0, sum(k=1, sqrtint(n^2 \ 2), issquare(n^2 - k^2)))}; /* Michael Somos, Mar 29 2015 */

a(n) = {my(f = factor(n/(2^valuation(n, 2)))); (prod(k=1, #f~, if ((f[k,1] % 4) == 1, 2*f[k,2] + 1, 1)) - 1)/2;} \\ Michel Marcus, Mar 08 2016

from math import prod
from sympy import factorint
def A046080(n): return prod((e<<1)+1 for p,e in factorint(n).items() if p&3==1)>>1 # Chai Wah Wu, Sep 06 2022

Let n = 2^e_2 * product_i p_i^f_i * product_j q_j^g_j where p_i == 1 mod 4, q_j == 3 mod 4; then a(n) = (1/2)*(product_i (2*f_i + 1) - 1). - Beiler, corrected
8*a(n) + 4 = A046109(n) for n > 0. - Ralf Stephan, Mar 14 2004
a(n) = 0 for n in A004144. - Lekraj Beedassy, May 14 2004
a(A084645(k)) = 1. - Ruediger Jehn, Jan 14 2022
a(A084646(k)) = 2. - Ruediger Jehn, Jan 14 2022
a(A084647(k)) = 3. - Jean-Christophe Hervé, Dec 01 2013
a(A084648(k)) = 4. - Jean-Christophe Hervé, Dec 01 2013
a(A084649(k)) = 5. - Jean-Christophe Hervé, Dec 01 2013
a(n) = A063725(n^2) / 2. - Michael Somos, Mar 29 2015
a(n) = Sum_{k=1..n} Sum_{i=1..k} [i^2 + k^2 = n^2], where [ ] is the Iverson bracket. - Wesley Ivan Hurt, Dec 10 2021
a(A002144(k)^n) = n. - Ruediger Jehn, Jan 14 2022

A008846 Hypotenuses of primitive Pythagorean triangles.

Original entry on oeis.org

5, 13, 17, 25, 29, 37, 41, 53, 61, 65, 73, 85, 89, 97, 101, 109, 113, 125, 137, 145, 149, 157, 169, 173, 181, 185, 193, 197, 205, 221, 229, 233, 241, 257, 265, 269, 277, 281, 289, 293, 305, 313, 317, 325, 337, 349, 353, 365, 373, 377, 389, 397, 401, 409, 421, 425, 433

Offset: 1

N. J. A. Sloane, Ralph Peterson (RALPHP(AT)LIBRARY.nrl.navy.mil)

Numbers of the form x^2 + y^2 where x is even, y is odd and gcd(x, y)=1. Essentially the same as A004613.
Numbers n for which there is no solution to 4/n = 2/x + 1/y for integers y > x > 0. Related to A073101. - T. D. Noe, Sep 30 2002
Discovered by Frénicle (on Pythagorean triangles): Méthode pour trouver ..., page 14 on 44. First text of Divers ouvrages ... Par Messieurs de l'Académie Royale des Sciences, in-folio, 6+518+1 pp., Paris, 1693. Also A020882 with only one of doubled terms (first: 65). - Paul Curtz, Sep 03 2008
All divisors of terms are of the form 4*k+1 (products of members of A002144). - Zak Seidov, Apr 13 2011
A024362(a(n)) > 0. - Reinhard Zumkeller, Dec 02 2012
Closed under multiplication. Primitive elements are in A002144. - Jean-Christophe Hervé, Nov 10 2013
Not only the square of these numbers is equal to the sum of two nonzero squares, but the numbers themselves also are; this sequence is then a subsequence of A004431. - Jean-Christophe Hervé, Nov 10 2013
Conjecture: numbers p for which sqrt(-1) exists in the p-adic numbering system. For example the 5-adic number ...2431212, when squared, gives ...4444444, which is -1, and 5 is in the sequence. - Thierry Banel, Aug 19 2022
The above conjecture was proven true by George Bergman. 3 known facts: (1) prime factors of a(n) are equal to 1 mod 4, (2) modulo such primes, sqrt(-1) exists, (3) if sqrt(m) exists mod r, r being odd, this extends to sqrt(m) in the r-adic ring. - Thierry Banel, Jul 04 2025

A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, pp. 10, 107.

Zak Seidov, Table of n, a(n) for n = 1..87881 (with a(n) up to 10^6).
Ron Knott, Pythagorean Triples and Online Calculators

Subsequence of A004431 and of A000404 and of A339952; primitive elements: A002144.
Cf. A004613, A020882, A073101.
Cf. A137409 (complement), disjoint union of A024409 and A120960.

a008846 n = a008846_list !! (n-1)
a008846_list = filter f [1..] where
   f n = all ((== 1) . (`mod` 4)) $ filter ((== 0) . (n `mod`)) [1..n]
-- Reinhard Zumkeller, Apr 27 2011

for x from 1 by 2 to 50 do for y from 2 by 2 to 50 do if gcd(x,y) = 1 then print(x^2+y^2); fi; od; od; [ then sort ].

Union[ Map[ Plus@@(#1^2)&, Select[ Flatten[ Array[ {2*#1, 2*#2-1}&, {10, 10} ], 1 ], GCD@@#1 == 1& ] ] ] (* Olivier Gérard, Aug 15 1997 *)
lst = {}; Do[ If[ GCD[m, n] == 1, a = 2 m*n; b = m^2 - n^2; c = m^2 + n^2; AppendTo[lst, c]], {m, 100}, {n, If[ OddQ@m, 2, 1], m - 1, 2}]; Take[ Union@ lst, 57] (* Robert G. Wilson v, May 02 2009 *)
Union[Sqrt[#[[1]]^2+#[[2]]^2]&/@Union[Sort/@({Times@@#,(Last[#]^2-First[#]^2)/2}&/@ (Select[Subsets[Range[1,33,2],{2}],GCD@@#==1&]))]] (* Harvey P. Dale, Aug 26 2012 *)

is(n)=Set(factor(n)[,1]%4)==[1] \\ Charles R Greathouse IV, Nov 06 2015

# for an array from the beginning
from math import gcd, isqrt
hypothenuses_upto = 433
A008846 = set()
for x in range(2, isqrt(hypothenuses_upto)+1):
    for y in range(min(x-1, (yy:=isqrt(hypothenuses_upto-x**2))-(yy%2 == x%2)) , 0, -2):
        if gcd(x,y) == 1: A008846.add(x**2 + y**2)
print(A008846:=sorted(A008846)) # Karl-Heinz Hofmann, Sep 30 2024

# for single k
from sympy import factorint
def A008846_isok(k): return not any([(pf-1) % 4 for pf in factorint(k)]) # Karl-Heinz Hofmann, Oct 01 2024

x^2 + y^2 where x is even, y is odd and gcd(x, y)=1. Essentially the same as A004613.

More terms from T. D. Noe, Sep 30 2002

cp's OEIS Frontend

A001844 Centered square numbers: a(n) = 2*n*(n+1)+1. Sums of two consecutive squares. Also, consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z; then sequence gives Z values.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

PARI

Python

Sage

Formula

Extensions

A002313 Primes congruent to 1 or 2 modulo 4; or, primes of form x^2 + y^2; or, -1 is a square mod p.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Formula

Extensions

A007528 Primes of the form 6k-1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Maple

Mathematica

PARI

PARI

Formula

A007519 Primes of form 8n+1, that is, primes congruent to 1 mod 8.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

PARI

PARI

PARI

SageMath

A031971 a(n) = Sum_{k=1..n} k^n.

Views

Author

Keywords

Comments

References

Links

A001844 Centered square numbers: a(n) = 2n(n+1)+1. Sums of two consecutive squares. Also, consider all Pythagorean triples (X, Y, Z=Y+1) ordered by increasing Z; then sequence gives Z values.