OFFSET
1,1
COMMENTS
No primitive Pythagorean triangle (PPT) can have all its integer sides squarefree since at least one side must be divisible by 4. However it is possible to find PPT's where none of the integer sides and the semiperimeter are squarefree. a(n) is the ordered occurrences of such semiperimeters.
EXAMPLE
a(5)=3225 as the PPT (825, 2752, 2873) has a semiperimeter of 3225, no member of (825, 2752, 2873, 3225) is squarefree and it is the 5th occurrence of such a semiperimeter.
MATHEMATICA
getpairs[k_] := Reverse[Select[IntegerPartitions[k, {2}], GCD[#[[1]], #[[2]]]==1 &]]; getlist[j_] := (newlist=getpairs[j]; Table[{newlist[[m]][[1]]^2-newlist[[m]][[2]]^2, 2newlist[[m]][[1]]*newlist[[m]][[2]], newlist[[m]][[1]]^2+newlist[[m]][[2]]^2, newlist[[m]][[1]]*j}, {m, 1, Length[newlist]}]); p=100; quads={}; Do[AppendTo[quads, getlist[2n+1]], {n, 1, p}]; sqquads=Select[Flatten[quads, 1], Union@Table[SquareFreeQ@#[[r]], {r, 1, 4}]=={False} &]; lst=Table[sqquads[[k]][[4]], {k, 1, Length@sqquads}]; Select[Sort@lst, #<2 p^2 &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Frank M Jackson, Feb 10 2014
STATUS
approved