#include <bits/stdc++.h>
#define ll long long
#define poly vector<int>
using namespace std;
const int N = 1e6 + 10, P = 998244353;
inline int ksm(int a, int b, int mod)
{
    int z = 1;
    while (b)
    {
        if (b & 1)
            z = 1ll * z * a % mod;
        a = 1ll * a * a % mod;
        b >>= 1;
    }
    return z;
}
namespace Poly
{
    const int mod = 998244353, g = 3, invg = 332748118;
    int lim, len, rev[N], invlim;
    inline void init(int l1, int l2)
    {
        lim = 1, len = 0;
        while (lim <= l1 + l2)
            lim <<= 1, len++;
        for (int i = 0; i < lim; i++)
            rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (len - 1));
        invlim = ksm(lim, mod - 2, mod);
    }
    inline void NTT(poly &f, int tgpe)
    {
        for (int i = 0; i < lim; i++)
            if (i < rev[i])
                swap(f[i], f[rev[i]]);
        for (int m = 2; m <= lim; m <<= 1)
        {
            int wn = ksm(tgpe ? g : invg, (mod - 1) / m, mod);
            for (int i = 0; i < lim; i += m)
            {
                int w = 1;
                for (int j = 0; j < m / 2; j++)
                {
                    int u = f[i + j], v = 1ll * w * f[i + j + m / 2] % mod;
                    f[i + j] = (u + v) % mod, f[i + j + m / 2] = (u - v + mod) % mod;
                    w = 1ll * wn * w % mod;
                }
            }
        }
        if (!tgpe)
        {
            for (int i = 0; i < lim; i++)
                f[i] = 1ll * f[i] * invlim % mod;
        }
    }
    inline void mul(poly &f, poly g)
    {
        int lf = f.size(), lg = g.size();
        init(lf, lg);
        f.resize(lim), g.resize(lim);
        NTT(f, 1);
        NTT(g, 1);
        for (int i = 0; i < lim; i++)
            f[i] = 1ll * f[i] * g[i] % mod;
        NTT(f, 0);
    }
}
int n;
poly f, g;
poly A, B;
int fac[N], ifac[N], inv[N];
void cdq(int l, int r)
{
    if (l == r)
    {
        if (!l)
            g[l] = 1;
        else if (l == 1)
            g[l] = P - 2;
        if (l && l < n)
        {
            g[l + 1] = (g[l + 1] + 1ll * (P - g[l]) * inv[l + 1]) % P;
        }
        return;
    }
    int mid = l + r >> 1;
    cdq(l, mid);
    if (!l)
    {
        A.resize(mid - l + 1);
        B.resize(mid - l + 1);
        for (int i = l; i <= mid; i++)
            A[i - l] = 1ll * g[i] * i % P;
        for (int i = 0; i <= mid; i++)
            B[i] = g[i];
        Poly::mul(A, B);
        for (int i = mid + 1; i <= r; i++)
            g[i] = (g[i] + 1ll * (P - A[i - l]) * inv[i]) % P;
    }
    else
    {
        A.resize(mid - l + 1);
        B.resize(r - l + 1);
        for (int i = l; i <= mid; i++)
            A[i - l] = 1ll * g[i] * i % P;
        for (int i = 0; i <= r - l; i++)
            B[i] = g[i];
        Poly::mul(A, B);
        for (int i = mid + 1; i <= r; i++)
            g[i] = (g[i] + 1ll * (P - A[i - l]) * inv[i]) % P;
        A.resize(mid - l + 1);
        B.resize(r - l + 1);
        for (int i = l; i <= mid; i++)
            A[i - l] = g[i];
        for (int i = 0; i <= r - l; i++)
            B[i] = 1ll * i * g[i] % P;
        Poly::mul(A, B);
        for (int i = mid + 1; i <= r; i++)
            g[i] = (g[i] + 1ll * (P - A[i - l]) * inv[i]) % P;
    }
    cdq(mid + 1, r);
}
void sol(int sz)
{
    if (sz == 1)
    {
        f[0] = 1;
        return;
    }
    int mid = sz + 1 >> 1;
    sol(mid);
    A.resize(sz);
    B.resize(sz);
    for (int i = 0; i < sz; i++)
    {
        A[i] = f[i];
        B[i] = g[i];
    }
    Poly::mul(B, A);
    Poly::mul(B, A);
    for (int i = 0; i < sz; i++)
    {
        f[i] = (2ll * f[i] + P - B[i]) % P;
    }
}
int main()
{
    fac[0] = fac[1] = inv[1] = ifac[0] = ifac[1] = 1;
    for (int i = 2; i < N; i++)
    {
        fac[i] = 1ll * fac[i - 1] * i % P;
        inv[i] = 1ll * (P - P / i) * inv[P % i] % P;
        ifac[i] = 1ll * ifac[i - 1] * inv[i] % P;
    }
    scanf("%d", &n);
    g.resize(n + 1), f.resize(n + 1);
    cdq(0, n);
    g[1] = P - 1;
    sol(n + 1);
    printf("%d", 1ll * f[n] * fac[n] % P);
}