1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
import math

n = int(input())

wynik = 0

nKw = n * n
mozliweKwadratyH = {i * i: i for i in range(1, n + 1)}

for a in range(1, n):
    for b in range(a, n):
        aKwPlusbKw = a * a + b * b

        if aKwPlusbKw > nKw:
            break

        for p in range(int(math.sqrt(aKwPlusbKw)) + 1, n + 1): # p - przekątna
            if p * p - aKwPlusbKw in mozliweKwadratyH:
                wynik += 1

print(wynik)