m = 1000000007
N = int(input())
a = list(map(int, input().split()))
f = [1]
for i in range(1, N):
    f.append((f[-1] * i) % m)
c = []
for i in range(N):
    c.append((f[N-1]*pow(f[i],m-2,m)*pow(f[N-1-i],m-2,m))%m)
ans = 0
for i in range(N):
    ans = (ans + a[i] * c[i]) % m
print(ans)