#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int n,a[11],v[11];
ll s,pre[5001]={1},f[5001][5001],g[5001];
const int mod=998244353;
ll qpow(ll a,ll b){
    ll s=1;
    for(;b;b>>=1,a=a*a%mod) if(b&1) s=s*a%mod;
    return s;
}
void dfs(int x){
    if(x==n+1){
        bitset <11> b;
        for(int i=1;i<=n;i++){
            b[a[i]]=b[a[i]%n+1]=b[(a[i]+n-2)%n+1]=1;
            if(b.count()==n) return (void)(s+=i);
        }
    }
    for(int i=1;i<=n;i++) if(!v[i]) a[x]=i,v[i]=1,dfs(x+1),v[i]=0;
}
void dp(){
    for(int i=3;i<n;i++) for(int j=1;j<=i;j++) (f[i][j]+=f[i-1][j-1]+f[i-2][j-1]+f[i-3][j-1])%=mod;
}
int main(){
	scanf("%d",&n);
    for(int i=1;i<=n;i++) pre[i]=pre[i-1]*i%mod;
    if(n<=8) return dfs(1),printf("%lld",s*qpow(pre[n],mod-2)%mod),0;
    memset(f,0,sizeof(f)),f[2][0]=1,dp();
    for(int i=0;i<n;i++) (g[i+1]+=f[n-1][i])%=mod;
    memset(f,0,sizeof(f)),f[2][1]=f[3][1]=f[4][1]=1,dp();
    for(int i=0;i<n;i++) (g[i+1]+=f[n-1][i]+f[n-2][i])%=mod;
    memset(f,0,sizeof(f)),f[2][1]=f[3][1]=1,dp();
    for(int i=0;i<n;i++) (g[i+1]+=f[n-1][i]+f[n-2][i]+f[n-3][i])%=mod;
    memset(f,0,sizeof(f)),f[2][2]=f[3][2]=f[4][2]=1,dp();
    for(int i=0;i<n;i++) (g[i+1]+=f[n-1][i]+f[n-2][i]+f[n-3][i])%=mod;
    for(int i=1;i<=n;i++) (s+=g[i]*pre[i]%mod*pre[n-i])%=mod;
    printf("%lld",(pre[n]*(n+1)-s+mod)%mod*qpow(pre[n],mod-2)%mod);
	return 0;
}