N = int(input())
A = list(map(int, input().split()))
MOD = 10**9 + 7

def factorial_inverse(n, mod):
    fac = [1] * (n+1)
    inv = [1] * (n+1)
    finv = [1] * (n+1)
    
    for i in range(2, n+1):
        fac[i] = fac[i-1] * i % mod
        inv[i] = mod - inv[mod%i] * (mod//i) % mod
        finv[i] = finv[i-1] * inv[i] % mod
        
    return fac, finv

fac, finv = factorial_inverse(10**5+100, MOD)
p = []
n = N-1
for i in range(N):
    p.append(fac[n] * (finv[i] * finv[n - i] % MOD) % MOD)

ans = 0
for i in range(N):
    ans += (A[i] * p[i]) % MOD
    ans %= MOD
print(ans)