p, g = 998244353, 3 invg = pow(g, p-2, p) W = [pow(g, (p - 1) >> i, p) for i in range(24)] iW = [pow(invg, (p - 1) >> i, p) for i in range(24)] def fft(k, f): for l in range(k)[::-1]: d = 1 << l u = 1 for i in range(d): for j in range(i, 1 << k, 2*d): f[j], f[j+d] = (f[j] + f[j+d]) % p, u * (f[j] - f[j+d]) % p u = u * W[l+1] % p def ifft(k, f): for l in range(k): d = 1 << l u = 1 for i in range(d): for j in range(i, 1 << k, 2*d): f[j+d] *= u f[j], f[j+d] = (f[j] + f[j+d]) % p, (f[j] - f[j+d]) % p u = u * iW[l+1] % p def convolve(a, b): n0, n1 = len(a), len(b) k = (max(n0, n1) - 1).bit_length() + 1 n = 1 << k a = a + [0] * (n-n0) b = b + [0] * (n-n1) fft(k, a), fft(k, b) for i in range(n): a[i] = a[i] * b[i] % p ifft(k, a) invn = pow(n, p - 2, p) return [a[i] * invn % p for i in range(n0 + n1 - 1)] n=int(input()) mod=998244353 M=(10**5)*3+1 fac=[1]*M ninv=[1]*M finv=[1]*M for i in range(2,M): fac[i]=fac[i-1]*i%mod ninv[i]=(-(mod//i)*ninv[mod%i])%mod finv[i]=finv[i-1]*ninv[i]%mod def comb(n,k): if k>n: return 0 return (fac[n]*finv[k]%mod)*finv[n-k]%mod a=[0]*n for i in range(n): a[i]=fac[n-i] b=[0]*n for i in range(n): b[i]=finv[i] c=convolve(a, b) #print(c) for i in range(1,n+1): ans=c[n-i]*pow(comb(n,i),mod-2,mod)%mod*fac[n]%mod*(n-i+1)%mod*finv[i-1] print(ans%mod)