y = int(input())
ar = [int(t) for t in input().split()]
def fact(n):
	return [x for x in range(2,(n//2)+1) if n % x == 0]

def sol(n,z):
    if n == 2 or n % 2 == 1:
        return z
    q = 0
    for i in fact(n):
        k = z//i
        q = max(q,k+sol(i,z))
    return q
m = 0
for j in ar:
    if j == 1:
        m += 1
    elif j % 2 == 1:
        m += j+1
    else:
        m += sol(j,j) + 1
print(m)