from functools import reduce
from math import gcd

n = int(input())
a = [int(k) for k in input().split()]

sm = sum(a)
mx = max(a)
a.sort()
a = [mx - k for k in a[:-1]]

z = reduce(lambda u, v: gcd(u, v), a)
y = sum([k // z for k in a])
print(y, z)
