OFFSET
1,1
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 3091 terms from Robert Israel)
EXAMPLE
1848 = 463^2 - 461^2 = 233^2 - 229^2 = 157^2 - 151^2.
MAPLE
N:= 10^6: # to get all terms <= N
V:= Vector(N/4):
p:= 3:
while p < N/2 do
q:= p;
p:= nextprime(p);
r:= (p^2-q^2)/4;
if r <= N/4 then
V[r]:= V[r]+1
fi
od:
map(`*`, select(t -> V[t]=3, [$1..N/4]), 4); # Robert Israel, Aug 13 2018
PROG
(PARI) is(n) = my(i=0, v=[]); forprime(p=5, n, v=[precprime(p-1), p]; if(v[2]^2-v[1]^2==n, i++)); i==3 \\ Felix Fröhlich, Aug 13 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Ray G. Opao, Feb 08 2004
EXTENSIONS
More terms from Ray Chandler, Feb 11 2004
STATUS
approved