n=int(input())
x=[0]*(5*10**5+1)
a=list(map(int,input().split()));ans=0;v=sum(a)
s=set();z={}
for i in range(n):
	x[a[i]]+=1
	s.add(a[i])
	if a[i] not in z:
		z[a[i]]=0
	z[a[i]]+=1
for i in range(1,10**5*5):
	x[i]+=x[i-1]
for i in list(s):
	c=v*z[i]
	for j in range(10**5*5//i):
		c-=(x[min(5*10**5,(j+1)*i-1)]-x[max(0,j*i-1)])*z[i]*(i*j)
	ans+=c
print(ans)