#include<cstdio>
const int N=4007,P=924844033;
int f[N][N][2],vis[N];
void inc(int&a,int b){a+=b-P,a+=a>>31&P;}
void dec(int&a,int b){a-=b,a+=a>>31&P;}
int mul(int a,int b){return 1ll*a*b%P;}
int main()
{
    int n,k,x,y=0,ans=0;
    scanf("%d%d",&n,&k),x=(n-1)/k+1,f[0][0][1]=1;
    for(int i=((n-1)%k+1)*2;i;--i) vis[y+=x]=1;
    for(;y<2*n;) vis[y+=x-1]=1;
    for(int i=1;i<=2*n;++i)
	for(int j=0;j<=n;++j)
	{
	    f[i][j][0]=f[i][j][1]=f[i-1][j][1];
	    if(!vis[i-1]&&j) inc(f[i][j][1],f[i-1][j-1][0]);
	}
    for(int i=0,k=1;i<=n;k=mul(k,++i)) ((n-i)&1? dec:inc)(ans,mul(f[2*n][n-i][1],k));
    printf("%d",ans);
}
