#include<bits/stdc++.h>
using namespace std;
const int mod=998244353;
const int N=3005;
int primes[16],cnt;
bool check(int x){
    if(x<=1)return false;
    for(int i=2;i*i<=x;i++){
        if(x%i==0)return false;
    }
    return true;
}
void init(){
    int idx=2;
    while(cnt<16){
        while(!check(idx))idx++;
        primes[cnt++]=idx;
        idx++;
    }
}
struct node{
    int v,s;
    node(int v,int s):v(v),s(s){}
};
int n,a[N];
vector<node> g[N];
int val[1<<16];
int f[2][1<<16];
int pos[N],idx;
int ksm(int a,int b){
    int res=1;
    while(b){
        if(b&1){
            res=1ll*res*a%mod;
        }
        a=1ll*a*a%mod;
        b>>=1;
    }
    return res;
}
int calc(int p){
    if(p==1)return 1;
    return (1-ksm(p,mod-2)+mod)%mod;
}
int tmp[1<<16];
int main(){
    init();
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        int x=a[i];
        int s=0;
        for(int j=0;j<16;j++){
            if(x%primes[j]==0){
                s|=1<<j;
                while(x%primes[j]==0)x/=primes[j];
            }
        }
        g[x].push_back(node(a[i],s));
    }
    for(int i=1;i<=n;i++){
        if(g[i].size()>0){
            idx++;
            pos[idx]=i;
        }
    }
    for(int j=0;j<(1<<16);j++){
        int res=1;
        for(int i=0;i<16;i++){
            if((j>>i)&1){
                res=1ll*res*calc(primes[i])%mod;
            }
        }
        val[j]=res;
    }
    f[0][0]=1;
    for(int p=1,cur=1;p<=idx;p++,cur^=1){
        for(int j=0;j<(1<<16);j++){
            f[cur][j]=f[cur^1][j];
            tmp[j]=f[cur^1][j];
        }
        int at=pos[p];
        for(auto x:g[at]){
            for(int j=(1<<16)-1;j>=0;j--){
                tmp[j|x.s]=(tmp[j|x.s]+1ll*tmp[j]*x.v%mod)%mod;
            }
        }
        int t=calc(at);
        for(int j=0;j<(1<<16);j++){
            tmp[j]-=f[cur^1][j];
            tmp[j]=(tmp[j]+mod)%mod;
            f[cur][j]=(f[cur][j]+1ll*tmp[j]*t%mod)%mod;
        }
    }
    int ans=0;
    for(int j=0;j<(1<<16);j++){
        ans=(ans+1ll*f[idx&1][j]*val[j]%mod)%mod;
    }
    cout<<ans<<endl;
    return 0;
}