#include <iostream>
#include <vector>
#include <map>
using namespace std;
const long long MOD = 1000000007;
long long modpow(long long a, long long b){
	long long ans = 1;
	while (b > 0){
		if (b % 2 == 1){
			ans *= a;
			ans %= MOD;
		}
		a *= a;
		a %= MOD;
		b /= 2;
	}
	return ans;
}
long long modinv(long long a){
	return modpow(a, MOD - 2);
}
int main(){
	int N;
	cin >> N;
	vector<int> A(N);
	for (int i = 0; i < N; i++){
		cin >> A[i];
	}
	map<int, int> mp;
	for (int i = 0; i < N; i++){
		int B = A[i];
		for (int j = 2; j * j <= B; j++){
			if (B % j == 0){
				int cnt = 0;
				while (B % j == 0){
					cnt++;
					B /= j;
				}
				mp[j] = max(mp[j], cnt);
			}
		}
		if (B != 1){
			mp[B] = max(mp[B], 1);
		}
	}
	long long lcm = 1;
	for (auto P : mp){
		for (int i = 0; i < P.second; i++){
			lcm *= P.first;
			lcm %= MOD;
		}
	}
	long long ans = 0;
	for (int i = 0; i < N; i++){
		ans += lcm * modinv(A[i]);
		ans %= MOD;
	}
	cout << ans << endl;
}