#include<bits/stdc++.h>
using namespace std;
#define pii pair<int,int>
#define ll long long
#define inf 0x3f3f3f3f
#define ull unsigned long long
#define pll pair<ll,ll>
#define ls d*2
#define rs d*2+1
#define mid (l+r)/2
#define lowbit(x) (x&(-x))
//#define endl "\n"
#define all(x) x.begin(),x.end()
#define int long long
//mt19937 seed;
//uniform_int_distribution<int>num(0,2e9);
const int N=2e5+10;
const int M=33;
int mod=998244353;
int dp[1<<17];
int prime[20]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,1};
int kp(int x,int y){
int re=1;
while(y){
if(y&1)re=re*x%mod;
x=x*x%mod;
y>>=1;
}
return re;
}
void solve() {
map<int, vector<pair<int, int>>> mp;
cin >> n;
for (int i = 1; i <= n; i++) {
int x;
cin >> x;
int nd = 0, sum = 1;
for (int i = 0; i < 16; i++) {
while (x % prime[i] == 0) {
x /= prime[i];
sum = (sum * prime[i]) % mod;
nd |= (1 << i);
}
}
mp[x].push_back({nd, sum});
}
vector<int> dp(1 << 17);
dp[0] = 1;
for (auto [x, vec] : mp) {
prime[16] = x;
for (auto [nd, sum] : vec) {
if (x != 1)
nd |= 1 << 16;
for (int i = (1 << 17) - 1; i >= 0; i--) {
int nsum = sum, nnd = nd | i;
if ((i & (1 << 16)) != 0)
nsum = (nsum * x) % mod;
else
nsum = (nsum * max(x - 1, 1ll)) % mod;
dp[nnd] = (dp[nnd] + dp[i] * nsum % mod) % mod;
}
}
for (int i = (1 << 16); i < (1 << 17); i++) {
int ni = i ^ (1 << 16);
dp[ni] = (dp[ni] + dp[i]) % mod;
dp[i] = 0;
}
}
int sum = 0;
for (int i = 0; i < (1 << 17); i++) {
for (int j = 0; j < 16; j++) {
if ((i & (1 << j)) != 0) {
dp[i] = (dp[i] * kp(prime[j], mod - 2)) % mod;
dp[i] = (dp[i] * max(prime[j] - 1, 1ll)) % mod;
}
}
sum = (sum + dp[i]) % mod;
}
cout << sum << '\n';
return;
}
signed main(){
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
int _=1;//cin>>_;
while(_--){
// int n;
// cin>>n;
// map<int,vector<pair<int,int>>>mp;
// for(int i=0;i<n;i++){
// int x;
// cin>>x;
// int num=0,sum=1;
// for(int j=0;j<16;j++){
// while(x%pri[j]==0)num|=(1<<j),sum*=pri[j],x/=pri[j];
// }
// mp[x].push_back({num,sum});
// }
// dp[0]=1;
// for(auto [x,vec]:mp){
// pri[16]=x;
// for(auto [num,sum]:vec){
// if(x!=1)num|=(1<<16);
// for(int i=(1<<17)-1;i>=0;i--){
// if(i>>16&1)dp[i|num]=(dp[i|num]+dp[i]*sum%mod*x%mod)%mod;
// else dp[i|num]=(dp[i|num]+dp[i]*sum%mod*max(x-1,1ll)%mod)%mod;
// }
// }
// for(int i=(1<<16);i<(1<<17);i++){
// dp[i-(1<<16)]=(dp[i-(1<<16)]+dp[i])%mod;
// dp[i]=0;
// }
// }
// int ans=0;
// for(int i=0;i<(1<<17);i++){
// for(int j=0;j<16;j++){
// if(i>>j&1){
// dp[i]=dp[i]*(pri[j]-1)%mod;
// dp[i]=dp[i]*qp(pri[j],mod-2)%mod;
// }
// }
// ans=(ans+dp[i])%mod;
// }
// cout<<ans<<'\n';
solve();
}
}