Based on experimental data for primes $< 1.4 \times 10^{10}$, I observed that
Every natural number $x$ which is a solution of $x^2 = y^2 - z^2$ where $y$ and $z$ are primes $> 5000$ has a prime factor greater than $17$.
Is this true in general or can we have a counter example?
Note: Posted in MO since it is unanswered in MSE
Code: Generates all all solutions in which the largest prime factor of $x$ is less than $101$.
s = 5
i = 1
f = 1
target = set = 10^6
q_max = 0
while True:
if s*(s+1)%30 == 0:
q = 2*s + 1
p = 2*s^2 + q
n = p - 1
if is_prime(p) and is_prime(q):
i = i + 1
F = prime_divisors(n)
if F[-1] <= 101:
f = f + 1
q_max = q
print (i,s,f,n,p,q_max, F[-1])
if s > target:
print "Reached", target, f,q_max
target = target + set
s = s + 1
Where we have used the fact that the conditions imply that $x$ is even.
Since $z=(m+n)(m-n)$ is prime we must have $m=n+1$. Thus $$x=2n(n+1)\quad y=2n^2+2n+1\quad z=2n+1$$
In particular, $x$ is divisible by the product of two consecutive integers, $n, n+1$. It is very difficult for two, large, consecutive numbers to both be divisible only by small primes.
– lulu Dec 20 '19 at 13:15