from functools import reduce
import math


def main():
    M = 10**9+7
    def fact(x): return reduce(lambda a, b: a*b % M, range(1, x+1), 1)

    int(input())  # K
    *C, = map(int, input().split())
    g = reduce(math.gcd, C, 0)
    s = sum(C)
    cnt = [0 for i in range(g+1)]

    for i in range(s):
        x = s // math.gcd(s, i)
        if g % x == 0:
            cnt[x] += 1

    ans = 0
    for i in range(1, g+1):
        if cnt[i]:
            a = 1
            for c in C:
                a = a*fact(c // i) % M
            ans = (ans + cnt[i] * fact(s//i) * pow(a, -1, M)) % M

    ans = ans * pow(s, -1, M) % M
    print(ans)

    return


if __name__ == '__main__':
    main()