#include <bits/stdc++.h>

using namespace std;

string ltrim(const string &);
string rtrim(const string &);
vector<string> split(const string &);

// Complete the solve function below.

const int maxn = 1e6 + 42;

vector<int> pr[maxn];

int par[maxn];

int get(int v) {
    return v == par[v] ? v : par[v] = get(par[v]);
}

int uni(int a, int b) {
    a = get(a);
    b = get(b);
    if(a == b) {
        return 0;
    } else {
        par[a] = b;
        return 1;
    }
}

const int mod = 1e9 + 7;

int bpow(int x, int n) {
    return n ? n % 2 ? 1LL * x * bpow(x, n - 1) % mod : bpow(1LL * x * x % mod, n / 2) : 1;
}

void init() {
    iota(par, par + maxn, 0);
}

int solve(vector<int> a) {
    init();
    int n = a.size();
    map<int, vector<int>> cnt;
    for(int i = 0; i < n; i++) {
        for(auto jt: pr[a[i]]) {
            cnt[jt].push_back(i);
        }
    }
    int comp = n;
    for(auto it: cnt) {
        for(int i = 1; i < it.second.size(); i++) {
            comp -= uni(it.second[i - 1], it.second[i]);
        }
    }
    return (bpow(2, comp) - 2 + mod) % mod;
}

signed main()
{
    for(int i = 2; i < maxn; i++) {
        if(pr[i].empty()) {
            for(int j = i; j < maxn; j += i) {
                pr[j].push_back(i);
            }
        }
    }
    ofstream fout(getenv("OUTPUT_PATH"));

    string t_temp;
    getline(cin, t_temp);

    int t = stoi(ltrim(rtrim(t_temp)));

    for (int t_itr = 0; t_itr < t; t_itr++) {
        string a_count_temp;
        getline(cin, a_count_temp);

        int a_count = stoi(ltrim(rtrim(a_count_temp)));

        string a_temp_temp;
        getline(cin, a_temp_temp);

        vector<string> a_temp = split(rtrim(a_temp_temp));

        vector<int> a(a_count);

        for (int i = 0; i < a_count; i++) {
            int a_item = stoi(a_temp[i]);

            a[i] = a_item;
        }

        int result = solve(a);

        fout << result << "\n";
    }

    fout.close();

    return 0;
}

string ltrim(const string &str) {
    string s(str);

    s.erase(
        s.begin(),
        find_if(s.begin(), s.end(), not1(ptr_fun<int, int>(isspace)))
    );

    return s;
}

string rtrim(const string &str) {
    string s(str);

    s.erase(
        find_if(s.rbegin(), s.rend(), not1(ptr_fun<int, int>(isspace))).base(),
        s.end()
    );

    return s;
}

vector<string> split(const string &str) {
    vector<string> tokens;

    string::size_type start = 0;
    string::size_type end = 0;

    while ((end = str.find(" ", start)) != string::npos) {
        tokens.push_back(str.substr(start, end - start));

        start = end + 1;
    }

    tokens.push_back(str.substr(start));

    return tokens;
}