#include <bits/stdc++.h> #define int long long #define mod 998244353 #define fir first #define sec second using namespace std; int t, n, a; map <int, vector <int>> aa; queue <int> q; inline int ksm(int x, int y) { int res = 1; while (y) { if (y & 1) res = res * x % mod; x = x * x % mod, y >>= 1; } return res; } signed main() { //int sttt = clock(); //freopen("inequality.in", "r", stdin); //freopen("inequality.out", "w", stdout); // for (int i = 2; i <= 1000000; i++) // inv[i] = (mod - mod / i) % mod * inv[i - 1] % mod; // inv[i] = ksm(i, mod - 2); // printf("%lld %lld\n", 2 * inv[2] % mod, 6 * inv[3] % mod); scanf("%lld", &t); while (t--) { scanf("%lld", &n); int ans = 0; for (int i = 1; i <= n; i++) { scanf("%lld", &a); for (int j = 2; j * j <= a; j++) { if (a % j == 0) { int cnt = 0; while (a % j == 0) cnt++, a /= j; aa[j].push_back(cnt); } } } for (pair <int, vector <int>> p : aa) sort(p.sec.begin(), p.sec.end()); int cur; for (int i = 1; i <= n; i++) { cur = 1; for (pair <int, vector <int>> p : aa) { if (!(p.sec).empty()) { cur = cur * ksm(p.fir, *(p.sec).rbegin()) % mod; (p.sec).pop_back(); } } ans = (ans + cur) % mod; } printf("%lld\n", ans); } // printf("Run time:%lldms", (clock() - sttt) * 1000 / CLOCKS_PER_SEC); return 0; }