#include <bits/stdc++.h>
using namespace std;
#define int long long
#define pii pair<int,int>
#define maxn 3005
int a[maxn];
vector<int>prime;
int id[maxn];
int dp[1<<17][2],lst[1<<17][2];
vector<pii>vt[maxn];
const int p=998244353;
int func[maxn];

int qpow(int a,int b){
    int ans=1;
    int tmp=a;
    while(b){
        if(b&1){
            ans=ans*tmp%p;
        }
        tmp=tmp*tmp%p;
        b>>=1;
    }
    return ans;
}

void solve(){
    int n;cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
    }
    for(int i=2;i<=3000;i++){
        bool flag=1;
        for(int j=2;j*j<=i;j++){
            if(i%j==0){
                flag=0;
                break;
            }
        }
        if(flag){
            id[i]=prime.size();
            prime.push_back(i);
            func[i]=(i-1)*qpow(i,p-2)%p;
        }
    }

    int up=sqrt(3000);
    int ans=0;

    dp[0][0]=1;

    for(int i=1;i<=n;i++){
        int sit=0,flag=0;
        for(auto v:prime){
            if(a[i]%v!=0)continue;
            if(v>up){
                vt[v].emplace_back(sit,a[i]);
                flag=1;
                break;
            }
            else{
                sit|=(1<<id[v]);
            }
        }
        
        if(flag)continue;
        for(int j=0;j<(1<<17);j++){
            if(!dp[j][0])continue;
            int tmp=a[i];
            for(int k=0;k<17;k++){
                if((sit>>k&1)&&!(j>>k&1)){
                    tmp=tmp*func[prime[k]]%p;
                    
                }
            }
            lst[j|sit][0]+=dp[j][0]*tmp%p;
            if(lst[j|sit][0]>=p)lst[j|sit][0]-=p;
        }
        for(int j=0;j<(1<<17);j++){
            dp[j][0]+=lst[j][0];
            if(dp[j][0]>=p)dp[j][0]-=p;
            lst[j][0]=0;
        }
    }

    for(int i=1;i<=3000;i++){
        if(vt[i].empty())continue;
        for(auto [sit,x]:vt[i]){
            for(int j=0;j<(1<<17);j++){
                int tmp=x;
                for(int k=0;k<17;k++){
                    if((sit>>k&1)&&!(j>>k&1)){
                        tmp=tmp*func[prime[k]]%p;
                    }
                }
                (lst[j|sit][1]+=dp[j][0]*tmp%p*func[i]%p+dp[j][1]*tmp%p)%=p;
            }
            for(int j=0;j<(1<<17);j++){
                dp[j][0]+=lst[j][0];
                dp[j][1]+=lst[j][1];
                if(dp[j][0]>=p)dp[j][0]-=p;
                if(dp[j][1]>=p)dp[j][1]-=p;
                lst[j][0]=lst[j][1]=0;
            }
        }
        for(int j=0;j<(1<<17);j++){
            dp[j][0]+=dp[j][1];
            dp[j][1]=0;
            if(dp[j][0]>=p)dp[j][0]-=p;
        }
    }

    for(int i=0;i<(1<<17);i++){
        ans+=dp[i][0];
        if(ans>=p)ans-=p;
    }
    cout<<ans<<'\n';

}

signed main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    int TT=1;//cin>>TT;
    while(TT--){
        solve();
    }
    return 0;
}