OFFSET
1,1
LINKS
Shiv K. Gupta, Areas of Pythagorean Triangles, post to NMBRTHRY list, Dec 23 2008
PROG
(PARI) {Q=Qfb(1, 0, 1); forprime(p=5, 10^5, p%4==1|next; t=qfbsolve(Q, p); A=t[1]*t[2]*(t[1]^2-t[2]^2); fordiv( A, d, if( issquare(d^2+(2*A/d)^2, &q) && isprime(q) && q>p, print1(p", "); break)))}
CROSSREFS
KEYWORD
nonn
AUTHOR
David Broadhurst, Feb 26 2009
STATUS
approved