This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A080109 #35 Oct 12 2024 09:05:32
%S A080109 25,169,289,841,1369,1681,2809,3721,5329,7921,9409,10201,11881,12769,
%T A080109 18769,22201,24649,29929,32761,37249,38809,52441,54289,58081,66049,
%U A080109 72361,76729,78961,85849,97969,100489,113569,121801,124609,139129
%N A080109 Square of primes of the form 4k+1 (A002144).
%C A080109 a(n) is the sum of two positive squares in only one way. See the Dickson reference, (B) p. 227.
%C A080109 a(n) is the hypotenuse of two and only two right triangles with integral legs (modulo leg exchange). See the Dickson reference, (A) p. 227.
%C A080109 In 1640 Fermat generalized the 3,4,5 triangle with the theorem: A prime of the form 4n+1 is the hypotenuse of one and only one right triangle with integral arms. The square of a prime of the form 4n+1 is the hypotenuse of two and only two... The cube of three and only three...
%D A080109 L. E. Dickson, History of the Theory of Numbers, Volume II, Diophantine Analysis. Carnegie Institution Publ. No. 256, Vol II, Washington, DC, 1920, p. 227.
%D A080109 Morris Kline, Mathematical Thought from Ancient to Modern Times, 1972, pp. 275-276.
%H A080109 Amiram Eldar, <a href="/A080109/b080109.txt">Table of n, a(n) for n = 1..10000</a>
%H A080109 Jorma K. Merikoski, Pentti Haukkanen, and Timo Tossavainen, <a href="https://doi.org/10.7546/nntdm.2024.30.3.516-529">The congruence x^n = -a^n (mod m): Solvability and related OEIS sequences</a>, Notes. Num. Theor. Disc. Math. (2024) Vol. 30, No. 3, 516-529. See p. 521.
%F A080109 a(n) = A002144(n)^2 = A070079(n)^2 + (4*A070151(n))^2, for n >= 1. - _Wolfdieter Lang_, Jan 13 2015
%F A080109 From _Amiram Eldar_, Dec 02 2022: (Start)
%F A080109 Product_{n>=1} (1 + 1/a(n)) = A243380
%F A080109 Product_{n>=1} (1 - 1/a(n)) = A088539. (End)
%e A080109 a(7) = 2809 is the hypotenuse of triangles 1241, 2520, 2809 and 1484, 2385, 2809, and only of these.
%e A080109 a(7) = 53^2 = 2809 = 45^2 + (4*7)^2, and this is the only way. - _Wolfdieter Lang_, Jan 13 2015
%t A080109 Select[4 Range[96] + 1, PrimeQ]^2 (* _Michael De Vlieger_, Dec 27 2016 *)
%o A080109 (PARI) fermat(n) = { for(x=1,n, y=4*x+1; if(isprime(y),print1(y^2" ")) ) }
%Y A080109 Cf. A002144, A080175. - _Wolfdieter Lang_, Jan 13 2015
%Y A080109 Cf. A070079, A070151, A088539, A243380.
%K A080109 nonn,easy
%O A080109 1,1
%A A080109 _Cino Hilliard_, Mar 16 2003
%E A080109 Edited: Name changed, part of old name as comment. Comments added and changed. Dickson reference added. - _Wolfdieter Lang_, Jan 13 2015