def fac(x):
    for i in range(2,int(x**0.5)+1):
        if x%i==0:
            return i
    return -1




p=[0]*(10**6+15)
#d=[0]*(10**6+15)
for i in range(2,10**6+5):
    if p[i]==0:
        for j in range(i,10**6+5,i):
            if p[j]==0:
                p[j]=i
for i in range(1,10**6+5):
    if p[i]==0:
        p[i]=i
a=int(input())
b=[int(i) for i in input().split()]
#print(p[:100])
ans=[]
for i in b:
    k=i
    an=0
    while(k!=1):
        an+=k
        if k>10**6:
            lk=fac(k)
            if lk==-1:
                k=1
            else:
                k=k//lk

        else:
            k=k//p[k]
    ans.append(an)
print(sum(ans)+len(ans))