OFFSET
1,1
REFERENCES
"Advanced Number Theory" by Harvey Cohn.
LINKS
PROG
(Python)
from itertools import count, islice
from sympy import primefactors
from sympy.solvers.diophantine.diophantine import diop_DN
def A031399_gen(startvalue=1): # generator of terms >= startvalue
return filter(lambda n: not (any(p&3==3 for p in primefactors(n)) or len(diop_DN(n, -1))), count(max(startvalue, 1)))
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
