'''input
3
1 7 24
'''

Pr=[1]*(10**6+1)

for i in range(4,10**6+1,2):
    Pr[i]=2
for i in range(3,1001,2):
    for j in range(i*i,10**6+1,2*i):
        Pr[j]=i
P=[]
for i in range(2,10**6+1):
    if Pr[i]==1:
        P.append(i)
mem=[0]*(10**6+1)
for i in range(2,1000):
    mem[i]=max(i+1,mem[i])
    for j in range(2,i+1):
        #print i*j,mem[i],mem[j],i,j,max(mem[i*j],mem[j]*i+1,mem[i]*j+1)
        mem[i*j]=max(mem[i*j],mem[j]*i+1,mem[i]*j+1)

M={}
def pf(a):
    pf=[]
    sq=(a**0.5)
    v=a
    for i in P:
        if i>sq:
            break
        while v%i==0:
            pf.append(i)
            v/=i

    if v>1:
        pf.append(v)
    return pf
mem[1]=1
def solve(a):
    prfac=pf(a)
    if a==1:
        return 1
    if len(prfac)>1:
        ans=(mem[prfac[0]])
    else:
        return prfac[0]+1
    for i in range(0,len(prfac)-1):

        ans=ans*prfac[i+1]+1
    return ans
n=input()
print sum([solve(int(x)) for x in raw_input().split()])