OFFSET
1,3
FORMULA
k is a term iff each prime factor congruent to 5 or 7 (mod 8) of 8*k!+1 comes in even degree.
PROG
(Python)
from math import factorial
from itertools import count, islice
from sympy.solvers.diophantine.diophantine import diop_quadratic
from sympy.abc import x, y
def A394730_gen(startvalue=0): # generator of terms >= startvalue
return filter(lambda k:bool(len(diop_quadratic(2*(factorial(k)-x**2)-y*(y+1)))), count(max(startvalue, 0)))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Max Alekseyev, Mar 30 2026
STATUS
approved
