from math import gcd

L = int(input())
S = set()
mod = 1000003
for m in range(10**4):
    for n in range(1, m):
        if 8 * m * (m + n) > L:
            break
        D = [m**2 - n**2, 2 * m * n, m**2 + n**2]
        g = gcd(D[0], gcd(D[1], D[2]))
        D.sort()
        for k in range(3):
            D[k] //= g
        S.add(tuple(D))
        
print(len(S)%mod)