#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const ll N=3003,M=(1<<16)+3,V=(1<<16)-1,H=998244353;
ll n,ans,cnt,a[N],b[N],pri[N],num[N],iv[M],sum[M],F[M],G[M];
bool vis[N];
vector<int>ve[N];
ll Val(ll x){return (x-1)*iv[x]%H;}
void Init()
{
	for(int i=2;i<N;i++)
	{
		if(!vis[i])pri[++cnt]=i,num[i]=cnt;
		for(int j=1;j<=cnt&&pri[j]*i<N;j++)
		{
			vis[pri[j]*i]=1;
			if(i%pri[j]==0)break;
		}
	}
	iv[0]=iv[1]=sum[0]=1;
	for(int i=2;i<N;i++)iv[i]=H-H/i*iv[H%i]%H;
	for(int i=1;i<=V;i++)sum[i]=sum[i-(i&-i)]*Val(pri[__lg(i&-i)+1])%H;
}
void Add(ll &x,ll y){x=(x+y)%H;}
void Sol(ll *f,const vector<int> &now)
{
	for(int i=0;i<=V;i++)f[i]=i==0;
	for(int id:now)
		for(int i=V;i>=0;i--)Add(f[i|b[id]],f[i]*a[id]);
}
int main()
{
	cin>>n;Init();
	for(int i=1,x;i<=n;i++)
	{
		cin>>x,a[i]=x;
		for(int j=1;j<=cnt&&pri[j]<=53;j++)if(x%pri[j]==0)
		{
			b[i]|=1<<(j-1);
			while(x%pri[j]==0)x/=pri[j];
		}
		if(x>1)ve[x].push_back(i);
	}
	vector<int>cur(n,0);iota(cur.begin(),cur.end(),1);
	Sol(F,cur);
	for(int t=1;t<=cnt;t++)if(pri[t]>53&&ve[pri[t]].size())
	{
		Sol(G,ve[pri[t]]);
		for(int i=1;i<=V;i++)F[i]=(F[i]-G[i]+H)%H;
		for(int i=1;i<=V;i++)Add(ans,G[i]*sum[i]%H*Val(pri[t])); 
	}
	for(int i=0;i<=V;i++)Add(ans,F[i]*sum[i]);
	cout<<ans;
}