MOD = 10**9+7
N = int(input())
F = list(map(int,input().split()))
A = [1]*(N+1)
for i in range(2,N+1):
    A[i] = (A[i-1]*i)%MOD
B = [1]*(N+1)
B[N] = pow(A[N],MOD-2,MOD)
for i in range(N-1,1,-1):
    B[i] = (B[i+1]*(i+1))%MOD
def comb(n,k):
    if k>n or k<0:
        return 0
    return (A[n]*B[k]*B[n-k])%MOD
ans = 0
for i in range(N):
    ans = (ans+comb(N-1,i)*F[i])%MOD
print(ans)