# ------------------------------------INPUT STARTS------------------------------------------------
q = int(input(""))
data = input("")
data = data.split(" ")
for i in range(q) :
	data[i] = int(data[i])
# ------------------------------------INPUT ENDS--------------------------------------------------
def magic(n,f) :
	ori = n
	if n == 2 or (2**(n-1))%n == 1 :
		return (1+n)
	elif n == 1 :
		return (1)
	else :
		pass
	split = 0
	coef = {}
	temp = 0
	i = 2 # divisor
	while(n!=1):
		if n%i == 0:
			if temp == i : # i.e i has not changed
				coef[i] += 1
				n = n/i
			else :
				coef[i] = 1
				n = n/i
			temp = i
		else :
			i += 1
	n = ori # resetting n for further use
	while(n!=1):
		m = max(coef)
		for i in range(1,coef[m]+1):
			split += f
			n = int(n/m)
			f = f*int(m)
		del coef[m]
	
	return(split+ori)
sums = 0
for i in range(q) :
	sums += magic(data[i],1)
print(sums)