#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef double db;
#define int long long
const int N = 2e5 + 10;
const double eps = 1e-10;
#define ull unsigned long long
#define endl '\n'
using i64 = int64_t;
const int mod =998244353;
int ans = 0;
int ksm(int a, int b) {
    int ans = 1;
    while (b) {
        if (b % 2 == 1) {
            ans = ans * a % mod;
        }
        b /= 2;
        a = a * a % mod;
    }
    return ans;
}
struct node {
    int h[32];

    //线性基增加了return true
    bool in(int x) {
        for (int i = 31; i >= 0; i--) {
            if ((1ll << i) & x) {
                if (h[i] == 0) {
                    h[i] = x;
                    return true;
                } else {
                    x ^= h[i];
                }
            }
        }
        return false;
    }

    int ma(int x,int y) {
        for (int i = 31; i >= 0; i--) {
            if(ans>max(x^h[i],y^h[i])){
                ans=max(x^h[i],y^h[i]);
                x^=h[i];
                y^=h[i];
            }
        }
        return x;
    }
}t1;

void solve() {
    int n, m;
    cin >> n ;
    int ans=0;
    int p=0;
    for(int i=31;i>=0;i--){
        t1.h[i]=0;
    }
    vector<int>a(n+5),num(n+5);
    for(int i=1;i<=n;i++){
        cin>>a[i];
        num[a[i]]++;
        if(!t1.in(a[i])){
            p++;
        }
    }
    ans+=ksm(2,p)-1;
   // ans%=mod;
  //  ans+=mod;
    ans%=mod;
   // cout<<ans<<" ";
    for(int i=1;i<=n;i++){
        if(num[i]==0)continue;
        int c1=ksm(2,num[i]-1),c2=1;
        for(int j=1;j<=sqrt(i);j++){
            if(i%j==0){
                if(j!=i&&num[j]) {
                    c2 *= ksm(2, num[j] - 1);
                    c2 %= mod;
                }
                int f2=i/j;
                if(f2!=i&&f2!=j&&num[f2]){
                    c2*=ksm(2,num[f2]-1);
                    c2%=mod;
                }
            }
        }
        ans+=c1*c2%mod;
        ans%=mod;
        //cout<<ans<<" ";
    }
    cout <<ans << '\n';
}

signed main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    // cout<<fixed<<setprecision(2);
    int T = 1;
    cin>>T;
    while (T--) solve();
    return 0;
}