#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#define I inline
#define int long long
using namespace std;

const int mod = 998244353;
const int MAXN = 65600;
const int V = 3010;
const int N = 2010;
int a[N],n,pri[N],cnt,f[2][MAXN],mask[V],mul[MAXN],g[2][MAXN];
int pw[22],id[N];
vector<int> v[V];

int qpow(int x,int y) {
	if(y==0) return 1;
	if(y==1) return x;
	int mk=qpow(x,y/2);
	return y%2?mk*mk%mod*x%mod:mk*mk%mod;
}
/*
3
1 2 3

*/
I void solve() {
	cin>>n;
	for(int i=1;i<=n;++i) cin>>a[i];
	int lst=0; f[0][0]=1;
	for(int i=1;i<=n;++i) {
		
		int now=lst^1,mk=a[i];
		for(int j=1;j<=16;++j) {
			while(mk%pri[j]==0) mk/=pri[j];
		}
		
		if(mk!=1) {
			v[mk].push_back(a[i]);
			continue;
		}
		
		for(int j=0;j<pw[16];++j) f[now][j]=f[lst][j];
		for(int j=0;j<pw[16];++j) {
			int val=mul[mask[a[i]]-(j&mask[a[i]])];
			f[now][j|mask[a[i]]]+=f[lst][j]*val%mod*a[i]%mod;
		}
		lst^=1;
	}
	
	for(int i=1;i<=3000;++i) {
		
		if(!v[i].size()) continue;
		
		for(int j=0;j<pw[16];++j) {
			g[0][j]=f[lst][j];
			g[1][j]=0;
		}
		
		for(auto &t:v[i]) {
			//int fac=i,sfac=t/i;
			for(int j=pw[16];j>=0;--j) {
				int val=mul[mask[t/i]-(j&mask[t/i])];
				g[1][j|mask[t/i]]+=g[1][j]*t%mod*val%mod;
				g[1][j|mask[t/i]]+=g[0][j]*t%mod*qpow(i,mod-2)%mod*(i-1)%mod*val%mod;
				g[1][j|mask[t/i]]%=mod;
			}
		}
		for(int j=0;j<pw[16];++j) {
			f[lst][j]=g[0][j]+g[1][j];
			f[lst][j]%=mod;
		}
	}
	int ans=0;
	for(int i=0;i<pw[16];++i) {
		ans+=f[lst][i];
		ans%=mod;
	}
	cout<<ans<<endl;
}
I void init() {
	for(int i=2;i*i<=3000;++i) {
		bool can=1;
		for(int j=2;j*j<=i;++j) {
			if(i%j==0) can=0;
		}
		if(can) pri[++cnt]=i,id[i]=cnt;
	}
	pw[0]=1;
	for(int i=1;i<=16;++i) pw[i]=pw[i-1]<<1;
	for(int i=1;i<=3000;++i) {
		for(int j=1;j<=16;++j) {
			if(i%pri[j]==0) mask[i]|=pw[j-1];
		} 
	}
	for(int i=0;i<pw[16];++i) {
		mul[i]=1;
		for(int j=1;j<=16;++j) {
			if(i&pw[j-1]) {
				mul[i]*=qpow(pri[j],mod-2)*(pri[j]-1)%mod;
				mul[i]%=mod;
			}
		}
	}
}

signed main() 
{
	init();
	solve();
}