OFFSET
1,1
COMMENTS
Each term is a prime power of the form p^e where p is in A002144 and e>1.
LINKS
Ray Chandler, Table of n, a(n) for n = 1..10000
MATHEMATICA
lst1 = {1, 1}; lst2 = {}; Do[ If[ GCD[m, n] == 1, a = 2m*n; b = m^2 - n^2; c = m^2 + n^2; If[ !PrimeQ@c, AppendTo[lst1, c]]], {m, 3, 1000}, {n, If[OddQ@m, 2, 1], m - 1, 2}]; lst1 = Sort@ lst1; Do[ If[ lst1[[n - 1]] != lst1[[n]] && lst1[[n]] != lst1[[n + 1]], AppendTo[lst2, lst1[[n]]]], {n, 2, Length@ lst1 - 1}]; Take[lst2, 50] (* Robert G. Wilson v, May 02 2009 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
John Harrison (harrison_uk_2000(AT)yahoo.co.uk), Apr 20 2009
EXTENSIONS
a(7) corrected by and a(17) and further terms from Robert G. Wilson v, May 02 2009
Minor edits to comments. - Ray Chandler, Nov 27 2019
STATUS
approved