from math import gcd

K=int(input())
C=list(map(int,input().split()))

mod=10**9+7

FACT=[1]
for i in range(1,2*10**6+1):
    FACT.append(FACT[-1]*i%mod)

FACT_INV=[pow(FACT[-1],mod-2,mod)]
for i in range(2*10**6,0,-1):
    FACT_INV.append(FACT_INV[-1]*i%mod)

FACT_INV.reverse()

def Combi(a,b):
    if 0<=b<=a:
        return FACT[a]*FACT_INV[b]%mod*FACT_INV[a-b]%mod
    else:
        return 0

S=sum(C)
LIST=[]
for i in range(1,S+1):
    if S%i==0:
        LIST.append(i)

D=dict()

for x in LIST:
    D[x]=0

for i in range(1,S+1):
    x=gcd(i,S)
    D[x]+=1

ANS=0

for x in D:
    C2=[]

    waru=S//x
    flag=1

    for c in C:
        if c%waru!=0:
            flag=0
            break
        else:
            C2.append(c//waru)

    if flag==0:
        continue

    #print(C2)

    SUM=sum(C2)

    score=1

    for c2 in C2:
        score=score*Combi(SUM,c2)%mod
        SUM-=c2

    ANS+=score*D[x]

LANS=ANS*pow(S,mod-2,mod)

print(LANS%mod)