#include <bits/stdc++.h>
using namespace std;
const int N=1e4+5;
const int MD=1e9+7;
unordered_map<int,int> ma,tmp;
int a[N];
int quick_pow(int x,int y) {
	int ans=1;
	while(y) {
		if(y&1) ans=1LL*ans*x%MD;
		y>>=1;
		x=1LL*x*x%MD;
	}
	return ans;
}
void add(int &x,int y) {
	x+=y;
	if(x>=MD) x-=MD;
}
int main() {
	int n;
	scanf("%d",&n);
	for(int i=0;i<n;i++) {
		scanf("%d",&a[i]);
		int t=a[i];
		for(int j=2;j*j<=t;j++) {
			while(t%j==0) tmp[j]++,t/=j;
		}
		if(t!=1) tmp[t]++;
		for(auto it:tmp) {
			ma[it.first]=max(ma[it.first],it.second);
		}
		tmp.clear();
	}
	int g=1;
	for(auto it:ma) {
		g=1LL*g*quick_pow(it.first,it.second)%MD;
	}
	//printf("%d\n",g);
	int ans=0;
	for(int i=0;i<n;i++) {
		add(ans,1LL*g*quick_pow(a[i],MD-2)%MD);
	}
	printf("%d\n",ans);
	return 0;
}