Submission #3456310

---

Source Code Expand

#include<bits/stdc++.h>
#define IL inline
#define RG register
#define LL long long
#define Upd(a) (a)>=mod?(a)-=mod:1
using namespace std;
int n, k;
const int _ =2001, mod = 924844033;
bool Use[_<<1];
LL dp[_<<1][_][2], pre=0, pg=0;
LL jc[_], g[_];
IL int Add (RG int now) {
	if (now<=n) {
		if (now+k>n) return -1;
		return now+k+n;
	}
	else {
		if (now+k-n>n) return -1;
		return now+k-n;
	}
}
IL void work (RG int now) {
	for (RG int i=0; i<=n; ++i) dp[pre+1][i][0]=dp[pre][i][0]+dp[pre][i][1], Upd(dp[pre+1][i][0]);
	++pre;
	Use[now]=1;
	for (RG int i=Add(now); ~i; i=Add(i)) {
		Use[i]=1;
		for (RG int j=1; j<=n; ++j) {
			dp[pre+1][j][1]=dp[pre][j-1][0];
		}
		for (RG int j=0; j<=n; ++j) {
			dp[pre+1][j][0]=dp[pre][j][1]+dp[pre][j][0];
			Upd(dp[pre+1][j][0]);
		}
		++pre;
	}
}
int main () {
	cin>>n>>k;
	dp[0][0][0]=1;
	for (RG int i=1; i<=n; ++i) {
		if (!Use[i]) work (i);
		if (!Use[i+n]) work (i+n);
	}
	jc[0]=1;
	for (RG int i=1; i<=n; ++i) jc[i]=1LL*jc[i-1]*i%mod;
	for (RG int i=0; i<=n; ++i) g[i]=1LL*(dp[n*2][i][0]+dp[n*2][i][1])%mod*jc[n-i]%mod;
	for (RG int i=1, j=1; i<=n; ++i, j=-j) g[0]-=1LL*g[i]*j%mod, g[0]=(g[0]+mod)%mod;
	cout<<g[0];
}

Submission Info

Judge Result