#include <bits/stdc++.h>
using namespace std;
const int N=(1<<16)+7,p=998244353;
int n,k,f[N],w[N],dp[3][N]; int g[N]; vector<int>vt[N];
long long pows(long long u,int v){
  long long ans=1;
  while(v>0){
    if(v&1) ans=ans*u%p; u=u*u%p,v=v>>1;
  }
  return ans;
}
int main() {
  std::ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  for(int i=2;i<=54;i++){
    g[i]=++k; w[k]=i;
    for(int j=2;j*j<=i;j++){
      if(i%j==0){
        k--,g[i]=0; break;
      }
    }
  }
  cin>>n;
  while(n--){
    int x; cin>>x; int num=x;
    for(int i=2;i<=54;i++) while(x%i==0) x/=i;
    vt[x].push_back(num);
  }
  dp[0][0]=1;
  for(int i=1;i<=2000;i++){
    if(vt[i].size()){
      for(int x:vt[i]){
        int D=0;
        for(int c=2;c<=54;c++){
          if(!g[c]) continue;
          if(x%c==0) D|=(1<<(g[c]-1));
        }
        for(int z=0;z<(1<<16);z++){
          dp[2][z|D]=(dp[2][z|D]+1ll*(dp[1][z]+dp[0][z])*x)%p;
        }
        for(int z=0;z<(1<<16);z++)
          dp[1][z]=(dp[1][z]+dp[2][z])%p,dp[2][z]=0;
      }
      long long inv=pows(i,p-2);
      for(int z=0;z<(1<<16);z++){
        if(i>1) dp[0][z]=(dp[0][z]+dp[1][z]*inv%p*(i-1))%p;
        else dp[0][z]=(dp[0][z]+dp[1][z])%p;
        dp[1][z]=0;
      }
    }
  }
  int ans=0;
  for(int z=0;z<(1<<16);z++){
    for(int x=1;x<=k;x++)
      if((z>>(x-1))&1)
        dp[0][z]=dp[0][z]*pows(w[x],p-2)%p*(w[x]-1)%p;
    ans=(ans+dp[0][z])%p;
  }
  cout<<ans;

}
/*
3
1 2 3
*/