from typing import List

class UnionFind:
    """0-indexed"""
    def __init__(self, n):
        self.n = n
        self.parent = [-1] * n
        self.__group_count = n

    def unite(self, x, y) -> bool:
        """xとyをマージ"""
        x = self.root(x)
        y = self.root(y)
        if x == y:
            return False

        self.__group_count -= 1

        if self.parent[x] > self.parent[y]:
            x, y = y, x

        self.parent[x] += self.parent[y]
        self.parent[y] = x
        return True

    def is_same(self, x, y) -> bool:
        """xとyが同じ連結成分か判定"""
        return self.root(x) == self.root(y)

    def root(self, x) -> int:
        """xの根を取得"""
        if self.parent[x] < 0:
            return x
        else:
            # 経路圧縮あり
            # self.parent[x] = self.root(self.parent[x])
            # return self.parent[x]
            # 経路圧縮なし
            return self.root(self.parent[x])

    def size(self, x) -> int:
        """xが属する連結成分のサイズを取得"""
        return -self.parent[self.root(x)]

    def all_sizes(self) -> List[int]:
        """全連結成分のサイズのリストを取得 O(N)"""
        sizes = []
        for i in range(self.n):
            size = self.parent[i]
            if size < 0:
                sizes.append(-size)
        return sizes
    
    def members(self, x) -> List[int]:
        """xが属するグループのリストを返す O(N)"""
        mem = []
        r = self.root(x)
        for i in range(self.n):
            if self.root(i) == r:
                mem.append(i)
        return mem

    def groups(self) -> List[List[int]]:
        """全連結成分の内容のリストを取得 O(N・α(N))"""
        groups = dict()
        for i in range(self.n):
            p = self.root(i)
            if not groups.get(p):
                groups[p] = []
            groups[p].append(i)
        return list(groups.values())

    @property
    def group_count(self) -> int:
        """連結成分の数を取得 O(1)"""
        return self.__group_count

N = int(input())
A = list(map(int, input().split()))
num = [0]*(10**5+1)
id = [list() for _ in range(10**5+1)]
for i,a in enumerate(A):
    num[a] += 1
    id[a].append(i)

uf = UnionFind(N)

res = 0
for i in range(1,10**5+1):
    if num[i]>0:
        for j in range(num[i]-1):
            f = uf.unite(id[i][j],id[i][j+1])
            if f:
                res += i
        for j in range(2*i,10**5+1,i):
            if num[j]>0:
                for k in range(num[j]):
                    f = uf.unite(id[i][0],id[j][k])
                    if f:
                        res += j

group = uf.groups()
if len(group)==1:
    print(res)
else:
    mn = [float('inf')]*len(group)
    for i in range(len(group)):
        for g in group[i]:
            mn[i] = min(mn[i],A[g])
    mm = min(mn)
    res += mm*sum(mn)-mm*mm
    print(res)