OFFSET
1,1
EXAMPLE
See A165236.
MATHEMATICA
amax=6*10^4; 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; If[PrimeQ[2*a+1], b=2*m*n; If[PrimeQ[2*b+1], If[GCD[a, b]==1, If[a>b, {a, b}={b, a}]; If[a>amax, Break[]]; c=m^2+n^2; If[PrimeQ[2*c+1], k++; AppendTo[lst, b]]]]]]; If[a>amax, Break[]], {m, n+1, 12!, 2}], {n, 1, q, 1}]; Union@lst
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Sep 09 2009
STATUS
approved