// state expression:

 // dp[i][j] = num of sequences

 // i: considered number i

 // j: there is j inversion pairs

 //

 // find the answer:

 // ans = dp[n][k]

 //

 // transferring:

 // dp[i][j] = sigma dp[i-1][j-i+1 to j]

 //

 // boundary:

 // dp[1][0] = 1

 #include <iostream>

 #include <stdio.h>

 #include <string.h>

 #define MAX_N 1005

 #define MAX_K 1005

 #define MOD 10000

 using namespace std;

 int n,t;

 int dp[MAX_N][MAX_K];

 int sum[MAX_N][MAX_K];

 void read()

 {

     cin>>n>>t;

 }

 void update_sum(int i,int j,int a)

 {

     if(j==) sum[i][j]=a;

     else sum[i][j]=(sum[i][j-]+a)%MOD;

 }

 int query_sum(int i,int x,int y)

 {

     if(x==) return sum[i][y];

     else return ((sum[i][y]-sum[i][x-])%MOD+MOD)%MOD;

 }

 void solve()

 {

     memset(dp,,sizeof(dp));

     memset(sum,,sizeof(sum));

     dp[][]=;

     for(int i=;i<=t;i++)

     {

         sum[][i]=;

     }

     for(int i=;i<=n;i++)

     {

         for(int j=;j<=t;j++)

         {

             dp[i][j]=query_sum(i-,max(,j-i+),j);

             update_sum(i,j,dp[i][j]);

         }

     }

 }

 void print()

 {

     cout<<dp[n][t]<<endl;

 }

 int main()

 {

     read();

     solve();

     print();

 }