#include <bits/stdc++.h>
using namespace std;

#define int long long

#define rep(i,s,n) for(int i = s;i<n;i++)
#define repe(i,s,n) for(int i = s;i<=n;i++)
#define rrep(i,s,n) for(int i = (n)-1;i>=(s);i--)
#define all(v) (v).begin(),(v).end()
#define pb push_back
#define fi first
#define se second
#define chmin(a,b) a=min((a),(b))
#define chmax(a,b) a=max((a),(b))
typedef long long ll;
typedef pair<int,int>pint;
typedef vector<int>vint;
typedef vector<pint>vpint;
typedef pair<pint,int> P1;
typedef pair<int,pint> P2;
typedef pair<pint,pint>PP;
static const ll maxLL = (ll)1 << 62;
const ll MOD=1000000007;
const ll INF=1e18;
int dx[]={0,0,1,-1};
int dy[]={-1,1,0,0};

int N,M;
int mod=1e9;
vint a[10010];

signed main(){
    cin.tie(0);
	ios::sync_with_stdio(false);
    cin>>N>>M;
    N%=1000*M;
    N/=1000;
    rep(i,0,10001){
        a[i].resize(i+2,0);
        rep(j,0,i+1){
            if(j==0)a[i][j]=1;
            else a[i][j]=(a[i-1][j]+a[i-1][j-1])%mod;
        }
    }
    cout<<a[M][N]<<endl;
































    return 0;
}