import sys
input = lambda :sys.stdin.readline()[:-1]
ni = lambda :int(input())
na = lambda :list(map(int,input().split()))
yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")
no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")
#######################################################################

from collections import defaultdict
 
class UnionFind():
    def __init__(self, n):
        self.n = n
        self.parents = [-1] * n
 
    def find(self, x):
        if self.parents[x] < 0:
            return x
        else:
            self.parents[x] = self.find(self.parents[x])
            return self.parents[x]
 
    def union(self, x, y):
        x = self.find(x)
        y = self.find(y)
 
        if x == y:
            return
 
        if self.parents[x] > self.parents[y]:
            x, y = y, x
 
        self.parents[x] += self.parents[y]
        self.parents[y] = x
 
    def size(self, x):
        return -self.parents[self.find(x)]
 
    def same(self, x, y):
        return self.find(x) == self.find(y)
 
    def members(self, x):
        root = self.find(x)
        return [i for i in range(self.n) if self.find(i) == root]
 
    def roots(self):
        return [i for i, x in enumerate(self.parents) if x < 0]
 
    def group_count(self):
        return len(self.roots())
 
    def all_group_members(self):
        group_members = defaultdict(list)
        for member in range(self.n):
            group_members[self.find(member)].append(member)
        return group_members
 
    def __str__(self):
        return '\n'.join(f'{r}: {m}' for r, m in self.all_group_members().items())

N = 2*10**5+10
ff = [0]*N
primes = []
for i in range(2, N):
    if ff[i] == 0:
        primes.append(i)
        for j in range(i, N, i):
            if ff[j] == 0:
                ff[j] = i

def factorize(n):
    res = []
    while n > 1:
        z = ff[n]
        res.append(z)
        #print(ff[:20], n, z)

        while n % z == 0:
            n //= z
            
    return res



n = ni()
a = na()

if all(x == 1 for x in a):
    print(2 * n)
    exit()

m = len(primes)
d = {p: i for i, p in enumerate(primes)}
uf = UnionFind(m+N)

z = 0
f = [0] * m

for i in range(n):
    if a[i] == 1:
        z += 1
    else:
        #print(factorize(a[i])   )
        for x in factorize(a[i]):
            f[d[x]] += 1
            uf.union(d[x], m+a[i])

ans = 10**18
s = set()
for i in range(m):
    if f[i] > 0:
        s.add(uf.find(i))
if f[0]> 0:
    ans = min(ans, (z+len(s)-1)*2)
else:
    ans = min(ans, (z+len(s))*2)
    x = min([i for i in range(m) if f[i] > 0])
    ans = min(ans, (z + len(s) - 1) * primes[x])

print(ans)