#include<bits/stdc++.h>
using namespace std;
#define int long long
const int mo=998244353;
int n;
int dp[5050][5050];
int frac[5050],inv[5050];
int po(int x,int y){
	int z=1;
	while(y){
		if(y%2) z*=x,z%=mo;
		x*=x,x%=mo;
		y/=2;
	}
	return z;
}
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(nullptr);
	cin>>n;
	if(n==3) cout<<1,exit(0);
	if(n==4) cout<<2,exit(0);
	dp[2][1]=1;
	for(int i=3;i<=n;i++){
		for(int j=1;j<i;j++){
			dp[i][j]=(dp[i-1][j-1]+dp[i-2][j-1]+dp[i-3][j-1])%mo;
		}
	}
	frac[0]=1,inv[0]=1;
	for(int i=1;i<=n;i++){
		frac[i]=frac[i-1]*i%mo;
		inv[i]=po(frac[i],mo-2);
	}
	int ans=0;
	for(int i=1;i<=n;i++){
		ans+=dp[n-2][i]*frac[i]%mo*3%mo*n%mo*frac[n-i-1]%mo*inv[n]%mo*(i+1)%mo;
		ans%=mo;
		ans+=dp[n-3][i]*frac[i]%mo*2%mo*n%mo*frac[n-i-1]%mo*inv[n]%mo*(i+1)%mo;
		ans%=mo;
		ans+=dp[n-4][i]*frac[i]%mo*1%mo*n%mo*frac[n-i-1]%mo*inv[n]%mo*(i+1)%mo;
		ans%=mo;
	}
	cout<<ans;
	return 0;
}