Я решил эту задачу но она слишком долго проходить компиляцию и поэтому задачу не засчитывает, можете помочь ускорить её?
n, x, y = map(int, input().split())
solutions = 0
for range1 in range(n+1):
for range2 in range(n+1):
for range3 in range(n+1):
for range4 in range(n+1):
if x*range1**2+y*range2**2-x*range3**2-y*range4**2 == 0:
solutions += 1
else:
solutions += 0
print(solutions)