OFFSET
1,1
COMMENTS
a(22) > 10^8.
Conjecture: All terms can be expressed in the form of 4*m^2 + 1 or 12*m^2 - 1.
LINKS
XiaoYang Zhang, How many prime numbers P are currently known to satisfy the equation y^2 = x^3 +- P^2? (in Chinese).
EXAMPLE
Prime 101 is in the sequence because both Mordell equations y^2 = x^3 + 101^2 and y^2 = x^3 - 101^2 have integer solutions as {{24, 155}} and {{101, 1010}}.
PROG
(SageMath)
def f(p):
for k in primes(p, p + 1000):
try:
E1 = EllipticCurve([0, 0, 0, 0, -k^2])
E2 = EllipticCurve([0, 0, 0, 0, k^2])
pts1 = E1.integral_points()
pts2 = E2.integral_points()
if any(pt[0] != 0 for pt in pts1) and any(pt[0] != 0 for pt in pts2):
print(k)
except:
pass
f(1)
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Zhining Yang, Jun 11 2026
STATUS
approved
