OFFSET
1,2
COMMENTS
EXAMPLE
PROG
(Python)
from math import isqrt, prod
from itertools import count, islice
from sympy import primefactors
def A185615_gen(startvalue=1): # generator of terms >= startvalue
return filter(lambda n: not (n+isqrt(5*n**2)>>1)%prod(primefactors(n)), count(max(startvalue, 1)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna and Sean A. Irvine, Jan 31 2011
STATUS
approved