#include <bits/stdc++.h>
using namespace std;
constexpr int N = 2e5 + 5, mod = 1e9 + 7;
struct Base{
    int a[20];
    void init() {
        for (int i = 0; i < 20; i++)
            a[i] = 0;
    }
    void insert(int x) {
        for (int i = 19; i >= 0; i--) {
            if (!((x >> i) & 1)) {
                continue;
            }
            if (!a[i]) {
                a[i] = x;
                break;
            }
            x ^= a[i];
        }
    }
    int count() {
        int cnt = 0;
        for (int i = 0; i < 20; i++)
            if (a[i])
                cnt++;
        return cnt;
    }
};
int n;
int pw[N];
int a[N];
int cnt[N];
vector<int> d[N];
void solve() {
    scanf("%d", &n);
    for (int i = 1; i <= n; i++) {
        scanf("%d", &a[i]);
        cnt[a[i]]++;
    }
    for (int i = 1; i <= n; i++)
        for (int j = 2 * i; j <= n; j += i)
            d[j].push_back(i);
    Base alls;
    alls.init();
    for (int i = 1; i <= n; i++)
        alls.insert(a[i]);
    int ans = pw[n - alls.count()] - 1;
    for (int i = 1; i <= n; i++) {
        if (!cnt[i]) continue;
        Base B;
        B.init();
        int sum = 0;
        for (int x : d[i]) {
            if (cnt[x]) {
                B.insert(x);
                sum += cnt[x];
            }
        }
        ans = (ans + 1ll * cnt[i] * pw[sum - B.count()]) % mod;
    }
    printf("%d\n", ans);
}
void clear() {
    for (int i = 0; i <= n; i++)
        cnt[i] = 0, d[i].clear();
}
int main() {
    pw[0] = 1;
    for (int i = 1; i < N; i++)
        pw[i] = 2ll * pw[i - 1] % mod;
    int T;
    scanf("%d", &T);
    while (T--) {
        solve();
        clear();
    }
}