#include <cstdio>
#include <iostream>
#include <algorithm>
#include <complex>
#include <queue>
#include <map>
using namespace std;

#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define FORR(i,a,b) for (int i=(b)-1;i>=(a);i--)
#define REP(i,n) for (int i=0;i<(n);i++)
#define RREP(i,n) for (int i=(n)-1;i>=0;i--)
#define pb push_back
#define ALL(a) (a).begin(),(a).end()

#define EPS (1e-10)
#define EQ(a,b) (abs((a)-(b)) < EPS)

#define PI 3.1415926535

typedef long long ll;
typedef pair<int, int> P;
//typedef complex<double> C;

const int MAX_M = 10000;

ll mod = 1000000000;
ll N, M;
ll C[MAX_M + 1][MAX_M + 1];

void input() {
	cin >> N >> M;
}

void setPascal(int N) {
	REP(i, N + 1) {
		C[i][0] = C[i][i] = 1;
	}
	FOR(i, 1, N + 1) {
		FOR(j, 1, i) {
			C[i][j] = (C[i-1][j-1] + C[i-1][j]) % mod;
		}
	}
}

void solve() {
	N /= 1000;
	if (N == 0) {
		cout << 1 << endl;
		return;
	}
	int r = N % M;
	setPascal(M);
	cout << C[M][r] << endl;
}

int main() {
	input();
	solve();
}