t=int(input())
l=map(int,raw_input().split())
ar=[]
def prm(m,n):    
    k=int(n**0.5)
    a=list(range(k+1))
    a[1]=0
    for i in range(2,int(k**0.5)+1):
        if a[i]:
            for j in range(i**2,k+1,i):
                a[j]=0
    if m<2:
        m=2
    b=list(range(m,n+1))
    for i in a:
        if i:
            if m%i==0:
                p=0
            else:
                p=i-m%i
            if b[p] in a:
                p+=i
            for j in range(p,n+1-m,i):
                b[j]=0
    for i in b:
        if i:
            ar.append(i)
prm(2,max(l))
an=0
for i in range(t):
    n=l[i]
    an+=1
    k=0
    if n>1:
        while True:
            if n in ar:
                an+=n
                break
            else:
                if n%ar[k]==0:
                    an+=n
                    n/=ar[k]
                    k=0
                else:
                    k+=1
print an