OFFSET
1,1
COMMENTS
The four primes are of the form 4*k + 1.
LINKS
Ray Chandler, Table of n, a(n) for n = 1..10000 (first 1000 terms from Colin Barker)
EXAMPLE
32045 is in the sequence because x^2 + y^2 = 32045 = 5*13*17*29 has solutions (x,y) = (2,179), (19,178), (46,173), (67,166), (74,163), (86,157), (109,142) and (122,131).
PROG
(PARI)
dop(d, nmax) = {
my(L=List(), v=vector(d, m, 1)~, f);
for(n=1, nmax,
f=factorint(n);
if(#f~==d && f[1, 1]>2 && f[, 2]==v && f[, 1]%4==v, listput(L, n))
);
Vec(L)
}
dop(4, 200000)
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Nov 15 2015
STATUS
approved