package main import ( "bufio" "fmt" "os" "sort" "strconv" ) func out(x ...interface{}) { fmt.Println(x...) } var sc = bufio.NewScanner(os.Stdin) func getInt() int { sc.Scan() i, e := strconv.Atoi(sc.Text()) if e != nil { panic(e) } return i } func getInts(N int) []int { ret := make([]int, N) for i := 0; i < N; i++ { ret[i] = getInt() } return ret } func getString() string { sc.Scan() return sc.Text() } // min, max, asub, absなど基本関数 func max(a, b int) int { if a > b { return a } return b } func min(a, b int) int { if a < b { return a } return b } func asub(a, b int) int { if a > b { return a - b } return b - a } func abs(a int) int { if a >= 0 { return a } return -a } func lowerBound(a []int, x int) int { idx := sort.Search(len(a), func(i int) bool { return a[i] >= x }) return idx } func upperBound(a []int, x int) int { idx := sort.Search(len(a), func(i int) bool { return a[i] > x }) return idx } const mod = int(1e9) func nCk(n, k int) int { if k == 0 { return 1 } x := 1 y := 1 for i := 1; i <= k; i++ { // out(n-i+1, i) x *= (n - i + 1) x %= mod y *= i y %= mod } return x / y } var C [10001][10001]int32 // C[n][k] -> nCk func combTable(N int) { C[0][0] = 1 for i := 1; i <= N; i++ { for j := 0; j <= N; j++ { if j == 0 || j == i { C[i][j] = 1 } else { C[i][j] = (C[i-1][j-1] + C[i-1][j]) % int32(mod) } } } } func main() { sc.Split(bufio.ScanWords) N, M := getInt(), getInt() N /= 1000 // x := N / M y := N % M combTable(10000) out(C[M][y]) }