#include<bits/stdc++.h>
using namespace std;typedef long long ll;const int MAX=3e5+9,MOD=1e9+7;vector<int>g[MAX];ll ans,fact[MAX],inv[MAX],mul[MAX],n;ll pw(ll a,ll b){return b?pw(a*a%MOD,b>>1)*(b&1?a:1)%MOD:1;}ll c(ll n,ll r){return fact[n]*inv[r]%MOD*inv[n-r]%MOD;}int main(){fact[0]=inv[0]=1;for(int i=1;i<MAX;i++)fact[i]=fact[i-1]*i%MOD,inv[i]=pw(fact[i],MOD-2);cin>>n;for(int i=0,x;i<n;i++){cin>>x;for(int j=2;j*j<=x;j++)if(x%j==0){int t=0;while(x%j==0)t++,x/=j;g[j].push_back(t);}if(x>1)g[x].push_back(1);}mul[0]=(pw(2,n-1)-1+MOD)%MOD;for(int i=1;i<n;i++)mul[i]=((mul[i-1]-c(n-1,i)-c(n-1,i-1))%MOD+MOD)%MOD;for(int i=1;i<MAX;i++){sort(g[i].begin(),g[i].end(),greater<int>());for(int j=0;j<g[i].size();j++)ans=(ans+g[i][j]*mul[j])%MOD;}cout<<ans;}