#include <bits/stdc++.h>
#define int long long
#define N 200007
#pragma GCC optimize(2)
using namespace std;
struct basis
{
    int ba[22],r;
    void init()
    {
        for(int j=0;j<22;j++)
        {
            ba[j]=0;
        }
        r=0;
    }
    void ins(int x)
    {
        for(int j=21;j>=0;j--)
        {
            if((x>>j)&1)
            {
                if(!ba[j])
                {
                    ba[j]=x;
                    r++;
                    break;
                }
                x=x^ba[j];
            }
        }
    }
    int rk()
    {
        return r;
    }
};
basis bs;
int T,n,a[N],c[N];
vector<int>vc[N];
const int p=998244353;
int qpow(int a,int x)
{
    int ans=1;
    while(x)
    {
        if(x&1)
        {
            ans=(ans*a)%p;
        }
        a=(a*a)%p;
        x>>=1;
    }
    return ans;
}
void pre()
{
    for(int i=1;i<N;i++)
    {
        for(int j=i+i;j<N;j+=i)
        {
            vc[j].push_back(i);
        }
    }
}
void solve()
{
    scanf("%lld",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&a[i]);
    }
    for(int i=1;i<=n;i++)
    {
        c[i]=0;
    }
    for(int i=1;i<=n;i++)
    {
        c[a[i]]++;
    }
    int ans=0;
    bs.init();
    for(int i=1;i<=n;i++)
    {
        bs.ins(a[i]);
    }
    ans=(ans+qpow(2,n-bs.rk()))%p;
    ans=((ans-1)%p+p)%p;
    for(int i=1;i<=n;i++)
    {
        if(!c[i])
        {
            continue;
        }
        bs.init();
        int all=0;
        for(int j:vc[i])
        {
            all+=c[j];
            bs.ins(j);
        }
        ans=(ans+qpow(2,all-bs.rk()+c[i]-1))%p;
    }
    printf("%lld\n",ans);
}
signed main()
{
    cin>>T;
    pre();
    while(T--)
    {
        solve();
    }
    return 0;
}