#include<bits/stdc++.h>
using namespace std;
const int maxn=4040,p=924844033;
int fac[maxn],dp[maxn][maxn][2],a[maxn],g[maxn],n,k,tot;
int main(){
	cin>>n>>k;
	fac[0]=1;for(int i=1;i<=n;++i)fac[i]=1ll*fac[i-1]*i%p;
	for(int i=1;i<=k;++i){
		for(int j=i;j<=n;j+=k)a[++tot]=j;
        for(int j=i;j<=n;j+=k)a[++tot]=j;
	}
	dp[1][0][0]=1;
	for(int i=2;i<=tot;++i)
		for(int j=0;j<=min(n,i/2);++j){
			dp[i][j][0]=(dp[i-1][j][1]+dp[i-1][j][0])%p;
			if(j>0&&a[i]-a[i-1]==k)dp[i][j][1]=dp[i-1][j-1][0];
		}
	for(int i=0;i<=n;++i)g[i]=1ll*(dp[tot][i][0]+dp[tot][i][1])*fac[n-i]%p;
	int ans=0;
	for(int i=0;i<=n;++i){
		if(i&1)ans=(ans+p-g[i])%p;else ans=(ans+g[i])%p;
	}
	cout<<ans;
}