#include<bits/stdc++.h> using namespace std; const int N = 1e6+10,md=998244353; #define int long long int t,n,a[N],isprime[N],prime[N],c,maxn; vector<int> G[N]; signed main(){ ios::sync_with_stdio(false); cin>>t; for(int i=2;i<=N-10;i++){ if(!isprime[i]){ prime[++c]=i; } for(int j=1;j<=c;j++){ if(i*prime[j]>(N-10)){ break; } isprime[i*prime[j]]=1; if(i%prime[j]) break; } } while(t--){ cin>>n; maxn=0; for(int i=1;i<=c;i++){ G[prime[i]].clear(); } for(int i=1;i<=n;i++){ cin>>a[i]; } for(int i=1;i<=n;i++){ for(int j=1;j<=c;j++){ if(prime[j]>a[i]) break; int tmp=0; while(a[i]%prime[j]==0){ tmp++; a[i]/=prime[j]; } if(tmp){ G[prime[j]].push_back(tmp); maxn=max(maxn,j); } } } for(int i=1;i<=maxn;i++){ sort(G[prime[i]].begin(),G[prime[i]].end(),greater<int>()); } int ans=0; for(int i=1;i<=n;i++){ int tmp=1; for(int j=1;j<=maxn;j++){ if(G[prime[j]].size()){ tmp=tmp%md*pow(prime[j],G[prime[j]][0]); G[prime[j]].erase(G[prime[j]].begin()); } } ans=(ans%md+tmp%md)%md; } cout<<ans<<endl; /* for(int i=1;i<=maxn;i++){ for(int j=0;j<G[prime[i]].size();j++){ cout<<G[prime[i]][j]<<" "; } cout<<endl; } */ } return 0; }