#include<bits/stdc++.h>
using namespace std;
constexpr int M=5005,mod=998244353;
int n,ans,now,fac[M],ifac[M],dp[4][M];
int qpow(int x,int n){
	int s=1; while (n){
		if (n&1) s=1ll*s*x%mod;
		x=1ll*x*x%mod; n>>=1;
	} return s;
}
int rdc(int x){return x>=mod?x-mod:x;}
int main(){
    scanf("%d",&n);
	for (int i=fac[0]=1;i<=n;i++) fac[i]=1ll*fac[i-1]*i%mod;
    dp[0][0]=1;
	for (int i=1;i<=n;i++){
        if (now) now--; else now=3;
        dp[now][0]=0;
		for (int j=1;j<=i;j++){
			dp[now][j]=dp[(now+1)&3][j-1];
			if (i>=2) dp[now][j]=rdc(dp[now][j]+dp[(now+2)&3][j-1]);
			if (i>=3) dp[now][j]=rdc(dp[now][j]+dp[(now+3)&3][j-1]);
		}
    }
    // for (int i=0;i<=n;i++) cerr<<dp[now][i]<<", ";cerr<<endl;
	for (int i=1;i<=n;i++){
        const int coef=i==n?fac[n]:((0ll+mod-fac[n-i])%mod*fac[i-1]%mod);
        ans=(ans+1ll*dp[now][i]*coef)%mod;
    }
    if (ans<0) ans+=mod;
	printf("%lld\n",1ll*ans*qpow(fac[n-1],mod-2)%mod);
	return 0;
}