import sys
input=lambda: sys.stdin.readline().rstrip()
import bisect
n=int(input())
m=int(input())
n//=1000
n%=m
ans=1
for i in range(1,n+1):
  ans=ans*(m-n+i)//i
print(ans%(10**9))