#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e6+7,inf=1e18,mod=998244353,S=15;
vector<int> p={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};
int power(int x,int t)
{
	int b=1;
	while(t)
	{
		if(t&1)
			b=b*x%mod;
		x=x*x%mod; t>>=1;
	}
	return b;
}
void O_o()
{
	int n;
	cin>>n;
	vector<int> a(n+1);
	for(int i=1; i<=n; i++)
	{
		cin>>a[i];
	}
	vector<vector<array<int,2>>> d(n+1);
	for(int i=1; i<=n; i++)
	{
		int t=a[i],s=0;
		for(int j=0; j<S; j++)
		{
			int x=p[j];
			while(t%p[j]==0)
			{
				t/=p[j];
				s|=(1<<j);
			}
		}
		d[t].push_back({a[i],s});
	}
	vector<array<int,2>> dp(vector<array<int,2>>(1<<S,{0,0}));
	dp[0][0]=1;
//	int now=0;
	for(int T=1; T<=n; T++)
	{
		
		if(!d[T].size()) continue;
		int phi=(T==1?1:(T-1)*power(T,mod-2)%mod);
		bool fir=1;
		for(auto [x,t]:d[T])
		{
			vector<array<int,2>> dpn(vector<array<int,2>>(1<<S,{0,0}));
			for(int s=0; s<(1<<S); s++)
			{
				(dpn[s][0]+=dp[s][0])%=mod;
				(dpn[s][fir^1]+=dp[s][1])%=mod;
				(dpn[s|t][1]+=dp[s][0]*x%mod*phi%mod)%=mod;
				int y=(fir?phi:T);
				(dpn[s|t][1]+=dp[s][1]*x%mod*y%mod)%=mod;
			}
			swap(dp,dpn);
		}
	}
	int ans=0;
	for(int i=0; i<(1<<S); i++)
	{
		int t=1;
		for(int j=0; j<S; j++)
		{
			if((i>>j)&1)
			{
				(t*=(p[j]-1)*power(p[j],mod-2)%mod)%=mod;
			}
		}
		(ans+=(dp[i][0]+dp[i][1])%mod*t%mod)%=mod;
	}
	cout<<ans<<"\n";
}
signed main()
{
	ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
	cout<<fixed<<setprecision(12);
	int T=1;
//	cin>>T;
	while(T--)
	{
		O_o();
	}
}
/*
3
1 2 3
*/