// O(N^2) #pragma GCC target("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include using namespace std; #define rep(i, n) for(long long i=0;i<(long long)(n);i++) #define REP(i,k,n) for(long long i=k;i<(long long)(n);i++) #define all(a) a.begin(),a.end() #define rsort(a) {sort(all(a));reverse(all(a));} #define pb emplace_back #define eb emplace_back #define lb(v,k) (lower_bound(all(v),(k))-v.begin()) #define ub(v,k) (upper_bound(all(v),(k))-v.begin()) #define fi first #define se second #define pi M_PI #define PQ(T) priority_queue #define SPQ(T) priority_queue,greater> #define dame(a) {out(a);return 0;} #define decimal cout< P; typedef tuple PP; typedef tuple PPP; typedef multiset S; using vi=vector; using vvi=vector; using vvvi=vector; using vvvvi=vector; using vp=vector

; using vvp=vector; using vb=vector; using vvb=vector; const ll inf=1001001001001001001; const ll INF=1001001001; const ll mod=998244353; const double eps=1e-10; template bool chmin(T&a,T b){if(a>b){a=b;return true;}return false;} template bool chmax(T&a,T b){if(a void out(T a){cout< void outp(T a){cout<<'('< void outvp(T v){rep(i,v.size())cout<<'('< void outvvp(T v){rep(i,v.size())outvp(v[i]);} template void outv(T v){rep(i,v.size()){if(i)cout<<' ';cout< void outvv(T v){rep(i,v.size())outv(v[i]);} template bool isin(T x,T l,T r){return (l)<=(x)&&(x)<=(r);} template void yesno(T b){if(b)out("yes");else out("no");} template void YesNo(T b){if(b)out("Yes");else out("No");} template void YESNO(T b){if(b)out("YES");else out("NO");} template void noyes(T b){if(b)out("no");else out("yes");} template void NoYes(T b){if(b)out("No");else out("Yes");} template void NOYES(T b){if(b)out("NO");else out("YES");} int main(){ ll n;cin>>n;assert(n<=5000); vp dp(n+1); dp[1]=P(1,0); rep(i,n-1){ rep(j,n+1){ dp[j].se+=dp[j].fi*j%mod; dp[j].se%=mod; } // outvp(dp); vp ndp(n+1); rep(j,n+1)if(dp[j].fi){ if(j>1){ ndp[j-1].fi+=dp[j].fi*j%mod*(j-1)%mod; ndp[j-1].se+=dp[j].se%mod*j%mod*(j-1)%mod; } ndp[j].fi+=dp[j].fi*2*j%mod; ndp[j].se+=dp[j].se%mod*2*j%mod; ndp[j+1].fi+=dp[j].fi; ndp[j+1].se+=dp[j].se; } rep(j,n+1)dp[j]=P(ndp[j].fi%mod,ndp[j].se%mod); // outvp(dp); } out(dp[1].se%mod); }