MI=lambda:map(int,input().split())
li=lambda:list(MI())
ii=lambda:int(input())
mod=998244353
n=ii()
a=li()
fac=[1]*(2*n)
for i in range(1,2*n):
    fac[i]=fac[i-1]*i%mod
ifac=[1]*(2*n)
ifac[-1]=pow(fac[-1],mod-2,mod)
for i in range(2*n-2,-1,-1):
    ifac[i]=ifac[i+1]*(i+1)%mod
invn=pow(n,mod-2,mod)
ans=0
for k in range(1,n+1):
    t=fac[2*n-k-1]*ifac[n-1]%mod*ifac[n-k]%mod
    c=k*t%mod*invn%mod
    ans=(ans+c*a[k-1])%mod
print(ans)