#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10,mod=998244353;
int a[N],n,ans=1,b[30],vis[N],cnt[N];
vector<int>g[N];
int qpow(int x){
	int y=2,ret=1;
	while(x){
		if(x&1)ret=1ll*ret*y%mod;
		y=1ll*y*y%mod;
		x>>=1;
	}
	return ret;
}
void solve(){
	ans=1;
	cin>>n;
	memset(b,0,sizeof(b));
	for(int i=1;i<=n;i++)cnt[i]=0,g[i].clear();
	for(int i=1;i<=n;i++){
		cin>>a[i];cnt[a[i]]++;
		int x=a[i];
		for(int j=20;j>=0;j--){
			if(x&(1<<j)){
				if(b[j])x^=b[j];
				else{
					b[j]=x;
					break;
				}
			}
		}
		if(!x)ans=1ll*ans*2%mod;
	}ans--;
	for(int i=1;i<=n;i++)
		if(cnt[i])
			for(int j=a[i];j<=n;j+=a[i])
				if(cnt[j])g[j].push_back(a[i]);
	for(int i=1;i<=n;i++){
		if(cnt[i]){
			int Cnt=0;
			memset(b,0,sizeof(b));
			for(int j:g[i]){
				Cnt+=cnt[j];
					for(int k=20;k>=0;k--){
						if(j&(1<<k)){
						if(b[k])j^=b[k];
						else{
							b[k]=j;
							break;
						}
					}
				}
				if(j)Cnt--;
			}
			ans=(1ll*ans+qpow(Cnt))%mod;
		}
	}
	cout<<ans<<"\n";
}
int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);
	int T;cin>>T;
	while(T--)solve(); 
}