#include<bits/stdc++.h>
using namespace std;
constexpr int maxn=2e6+5;
vector<int> d[maxn>>1];
inline void solve(){
    int n;cin>>n;
    unordered_map<int,int> mp;
    for(int i=1,x;i<=n;i++) cin>>x,mp[x]++;
    long long ans=0;
    for(auto [x,y]:mp){
        for(auto z:d[x]){
            int w=(x^z);
            if(mp.count(w)&&gcd(w,x)==z) ans+=y*1ll*mp[w];
        }
    }
    cout<<ans/2<<'\n';
}
int main(){
    for(int i=1;i<=1e6;i++)
    for(int j=i+i;j<=1e6;j+=i) d[j].push_back(i);
    int T;cin>>T;
    while(T--) solve();
}