It is known that if $x^2 + y^2 = z^2$ is a primitive Pythagorean triplet then $z$ is not divisible by any prime of the form $4k-1$. The following is a generalization of this classical result which shows that the source of this property is not the hypotenuse $z$, but the two orthogonal side $x$ and $y$:
Conjecture: Let $f(x,y) = a_0x^n + a_1 x^{n-1}y + a_2x^{n-2}y^2 + \cdots + a_ny^n$, $n \ge 2, a_0a_n \ne 0$. Then there are infinitely many primes of the form $8k+3$ which do not divide $f(x,y)$ for any primitive Pythagorean triplet $x^2 + y^2 = z^2$.
Questions:
- Is the conjecture known? I am looking for a proof or disproof of the conjecture.
- Or can we prove the simpler case for $x^n +y^n$?
Update 1: If has been proved for the special case $x+y$ (in the answer below) and is already known to be true for $x^2 + y^2$. Experimental data shows that $x^3 + y^3$ is not divisible by infinitely many primes of the form $8k+3$ while $x^4 + y^4$ is not divisible by infinitely many primes of the form $8k+3, 8k+5$ and $8k+7$.
Update 2: Posted in MO since the general case of the conjecture is open
Sagemath Code
r = 2
fac = prime_factors(1)
while r <= 200:
s = 1 + r%2
while(s < r):
if gcd(s,r)== 1:
b = r^2 - s^2
c = 2rs
# t = 5b^3 +7b^1c^2 + 5b^2c^1 + 2c^3
# t = b^4 - 3b^1c^3 - b^2c^2 - 1c^4
# t = b^2 - 11bc - c^2
t = b + c
fac = fac + prime_factors(t)
fac = list(dict.fromkeys(fac))
s = s + 2
r = r + 1
fac = sorted(fac)
fac2 = fac
fac = fac[:floor(0.5*len(fac))]
P = Primes()
prime = P[:prime_pi(max(fac))]
diff = list(set(prime) - set(fac))
diff = sorted(diff)
print diff