n = int(input())
a = list(map(int, input().split()))
fact = [0] * (n + 1)
rfact = [0] * (n + 1)
MOD = 10 ** 9 + 7
fact[0] = 1
for i in range(1, n + 1):
    fact[i] = fact[i - 1] * i
    fact[i] %= MOD

rfact[n] = pow(fact[n], MOD - 2, MOD)
for i in range(n - 1, -1, -1):
    rfact[i] = rfact[i + 1] * (i + 1)
    rfact[i] %= MOD


def comb(x, y):
    return fact[x] * rfact[x - y] * rfact[y] % MOD


num = 0
for i in range(n):
    num += comb(n - 1, i) * a[i]
    num %= MOD

print(num)