N = int(input()) M = int(input()) a = N//1000//M b = N - (a * M) * 1000 c = b//1000 # def perm(n, k): # result = 1 # for i in range(1, k + 1): # result *= (n - i + 1) # result //= i # return result def pascal(n, k): z = n+1 p = [[0 for i in range(z)] for j in range(z)] p[0][0] = 1 for x in range(1, z): for y in range(0, z): p[x][y] = (p[x-1][y-1] + p[x-1][y])%1000000000 return p[n][k] # print(M, c, perm(M, c)) # print(perm(M, c)%1000000000) # print(pascal(M, c)) print(pascal(M, c)%1000000000)