OFFSET
1,1
REFERENCES
F. Rubin, "Squared" Pythagorean Triples, Solution to problem 2306, J. Recreational Mathematics, Vol. 29, No. 1, 1998, p. 73.
LINKS
Ray Chandler, Table of n, a(n) for n = 1..10000
FORMULA
a(n)=e^2+f^2, where e>f, e=j^2 - jk + k^2/2 and f=jk for coprime pairs (j, k) with k even.
EXAMPLE
9161 is in the sequence because of the triple 5289^2 + 7480^2 = 9161^2 where we have 5289+7480=113^2.
Similarly, 205 belongs to the triple (133,156,205) and 133+156=17^2.
MATHEMATICA
terms = 1000; jmax = 100; kmax = 200;
Reap[Do[If[CoprimeQ[j, k], e = j^2 - j k + k^2/2; f = j k; If[e > f, Sow[e^2 + f^2]]], {j, 1, jmax}, {k, 2, kmax, 2}]][[2, 1]] // Union // Take[#, terms]& (* Jean-François Alcover, Mar 04 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Nov 06 2003
STATUS
approved