OFFSET
1,2
COMMENTS
If only one side of a PPT is divisible by primes of the form 4k+1 that side is the hypotenuse. Furthermore this hypotenuse is divisible by 5 because all PPT's must have at least one side divisible by 5. However this is only a necessary condition. It is not sufficient. For instance, PPT's can share the same hypotenuse divisible by 5, some will have no other side divisible by primes of the form 4k+1 and others will. E.g. PPT's (16,63,65) and (33,56,65) so 13 is in the sequence.
Also an identical sequence is generated from the superset of primitive Heronian triangles if only one side is divisible by primes of the form 4k+1 and this side is divided by 5 (see Yiu link).
LINKS
EXAMPLE
a(1)=1 as (3,4,5) is the first PPT where only the hypotenuse is divisible by primes of the form 4k+1, whereas the PPT (5,12,13) has two sides and the PPT (39,80,89) has three sides divisible by primes of the form 4k+1.
MATHEMATICA
pyprimeQ[n0_] := If[Length@Select[FactorInteger[n0], Mod[#[[1]], 4]==1 &] > 0, 1, 0]; lst1 = {}; Do[If[GCD[m, n]==1&&m<n&&OddQ[m+n], AppendTo[lst1, {n^2-m^2, 2m*n, n^2+m^2}]], {n, 1, 100}, {m, 1, n}]; SetAttributes[pyprimeQ, Listable]; lst2=Select[lst1, Total[pyprimeQ[#]]==1 &]; Union@Table[lst2[[k]][[3]]/5, {k, 1, Length[lst2]}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Frank M Jackson, Oct 18 2014
STATUS
approved