import sys

sys.setrecursionlimit(10**7)
def I(): return int(sys.stdin.readline().rstrip())
def MI(): return map(int,sys.stdin.readline().rstrip().split())
def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))
def LI2(): return list(map(int,sys.stdin.readline().rstrip()))
def S(): return sys.stdin.readline().rstrip()
def LS(): return list(sys.stdin.readline().rstrip().split())
def LS2(): return list(sys.stdin.readline().rstrip())


def factorial(N,p):  # 0!~N!(mod p) のリスト
    res = [0]*(N+1)
    res[0] = 1
    for i in range(1,N+1):
        res[i] = (res[i-1]*i) % p
    return res


def factorial_inverse(N,p):  # 0!~N!(mod p) の逆元のリスト
    fac = factorial(N,p)
    res = [0]*(N+1)
    a = pow(fac[-1],p-2,p)
    for i in range(N,-1,-1):
        res[i] = a
        a *= i
        a %= p
    return fac,res


N = I()
A = LI()
mod = 10**9+7

fac,fac_inv = factorial_inverse(N-1,mod)

ans = 0
for i in range(N):
    ans += fac[N-1]*fac_inv[i]*fac_inv[N-1-i]*A[i]
    ans %= mod

print(ans)