N = int(input()) -1
a = list(map(int,input().split()))

def func(N,prime):
    f = [1] * (N + 1)
    for i in range(1,N + 1):
        f[i] = (f[i -1] * i) % prime
    return f

def inv(x,prime):
    _x = x
    inv = 1
    k = prime -2
    while k > 0:
        if k % 2 == 1:
            inv = (inv * _x) % prime
        _x = (_x * _x) % prime   
        k = int(k /2)
    return inv

def func_inv(func,prime):
    N = len(func)-1
    func_inv = [1] * (N+1)
    func_inv[N] = inv(func[N],prime)
    for i in range(1,N):
        func_inv[N - i] = ((N -i + 1) * func_inv[N -i + 1]) % prime
    return func_inv

def comb(N,k, prime):
    return func[N]*func_inv[k]*func_inv[N - k] % prime

prime = 10 ** 9 + 7
func = func(N,prime)
func_inv = func_inv(func,prime)

print(sum([comb(N,k,prime) * a[k] % prime for k in range(N+1)]) % prime)