n=int(input())


#############################
#############
cnb_max=n+10
mod=998244353
#############

kai=[1]*(cnb_max+1)
rkai=[1]*(cnb_max+1)
for i in range(cnb_max):
    kai[i+1]=kai[i]*(i+1)%mod

rkai[cnb_max]=pow(kai[cnb_max],mod-2,mod)
for i in range(cnb_max):
    rkai[cnb_max-1-i]=rkai[cnb_max-i]*(cnb_max-i)%mod

def cnb(x,y):
    if y>x:
        return 0
    if x<0:return 0
    if y<0:return 0
    return (kai[x]*rkai[y]%mod)*rkai[x-y]%mod


def inv(n):
    return kai[n-1]*rkai[n]%mod

##################################

dp=[0,0,0,0,0]
dp[0]=1
for i in range(1,n+1):
    newdp=[0,0,0,0,0]
    for j in range(4):
        newdp[j]=(dp[j]+inv(i)*dp[j-1])%mod
    dp=newdp[:]
ans=dp[1]+6*dp[2]+6*dp[3]
ans*=kai[n]
print(ans%mod)