#include <bits/stdc++.h>
using namespace std;

const int N=5e3+5,H=998244353;
int n,f[N][N][3],ans=0,fac[N],ifac[N];

int adc(int a,int b){return a+b>=H?a+b-H:a+b;}
int dec(int a,int b){return a<b?a-b+H:a-b;}
int mul(int a,int b){return 1ll*a*b%H;}
void add(int &a,int b){a=adc(a,b);}
void del(int &a,int b){a=dec(a,b);}

int qpow(int a,int b=H-2){
  int res=1;
  while(b){if(b&1) res=mul(res,a);a=mul(a,a),b>>=1;}
  return res;
}

void init(){
  fac[0]=1;
  for(int i=1;i<N;i++) fac[i]=mul(fac[i-1],i);
  ifac[N-1]=qpow(fac[N-1]);
  for(int i=N-1;i>=1;i--) ifac[i-1]=mul(ifac[i],i);
}

int binom(int n,int m){
  if(m<0||n<m) return 0;
  return mul(fac[n],mul(ifac[m],ifac[n-m]));
}


int main(){
  ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
  cin>>n;init();

  f[1][1][0]=1;
  for(int i=2;i<=n;i++){
	for(int j=1;j<=i;j++) for(int k=0;k<3;k++) add(f[i][j][0],f[i-1][j-1][k]);
	for(int j=0;j<i;j++) for(int k=1;k<3;k++) add(f[i][j][k],f[i-1][j][k-1]);
  }

  for(int i=0;i<=n;i++){
	int res=0;
    for(int j=0;j<3;j++) add(res,f[n][i][j]);
	for(int j=0;j<2;j++) add(res,f[n-1][i][j]);
	add(res,f[n-2][i][0]);
    add(ans,mul(qpow(binom(n,i)),dec(binom(n,i),res)));
  }
  cout<<ans;
  return 0;
}