OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..315
EXAMPLE
a(1) = 1799820 is a term because it is the product of the sides (99, 20, 101) and inradius 9 of a primitive Pythagorean triple and 2*1799820+-1 = 3599639 and 3599641 and 3*1799820+-1 = 8999099 and 8999101 are all prime.
MAPLE
N:= 10^23: # for terms <= N
R:= {}:
for m from 1 while 2*(m-1)^2*(m+1)*m*(m^2+1) <= N do
for n from 1 + (m mod 2) to m-1 by 2 do
v:= 2*(m-n)^2 * (m+n)*m*n^2*(m^2+n^2);
if v/(m-n)^2 > N then break fi;
if v <= N and igcd(m, n) = 1 and andmap(isprime, [2*v-1, 2*v+1, 5*v-1, 5*v+1]) then
R:= R union {v};
fi od od:
sort(convert(R, list));
CROSSREFS
KEYWORD
nonn
AUTHOR
Will Gosnell and Robert Israel, Apr 19 2026
STATUS
approved
