from collections import * from heapq import * import sys sys.setrecursionlimit(10 ** 6) int1 = lambda x: int(x) - 1 p2D = lambda x: print(*x, sep="\n") def II(): return int(sys.stdin.readline()) def MI(): return map(int, sys.stdin.readline().split()) def MF(): return map(float, sys.stdin.readline().split()) def LI(): return list(map(int, sys.stdin.readline().split())) def LF(): return list(map(float, sys.stdin.readline().split())) def LLI(rows_number): return [LI() for _ in range(rows_number)] def gcd(a,b): while b>0:a,b=b,a%b return a def main(): md=10**6+3 l=II() ans=0 for m in range(1,l): for n in range(1,m): a=m**2-n**2 b=2*m*n c=m**2+n**2 if 4*(a+b+c)>l:break if gcd(gcd(a,b),c)>1:continue ans+=1 print(ans%md) main()