OFFSET
1,1
EXAMPLE
see A165250.
MATHEMATICA
amax=10^5; lst={}; k=0; q=12!; Do[If[(e=((n+1)^2-n^2))>amax, Break[]]; Do[If[GCD[m, n]==1, a=m^2-n^2; b=2*m*n; If[GCD[a, b]==1, If[a>b, {a, b}={b, a}]; If[a>amax, Break[]]; c=m^2+n^2; x=a+b+c; If[PrimeQ[x-1]&&PrimeQ[x+1], k++; AppendTo[lst, b]]]], {m, n+1, 12!, 2}], {n, 1, q, 1}]; Union@lst
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Sep 11 2009
STATUS
approved