#include <bits/stdc++.h>

using i64 = long long;
using u64 = unsigned long long;
using u32 = unsigned;
using u128 = unsigned __int128;
i64 mul(i64 a, i64 b, i64 m) {
    i64 c = a * b - i64(1.0L * a * b / m) * m;
    c %= m;
    if (c < 0) {
        c += m;
    }
    return c;
}
i64 power(i64 a, i64 b, i64 m) {
    i64 res = 1 % m;
    for (; b; b >>= 1, a = mul(a, a, m))
        if (b & 1)
            res = mul(res, a, m);
    return res;
}
bool isprime(i64 n) {
    if (n < 2)
        return false;
    static constexpr int A[] = {2, 3, 5, 7, 11, 13, 17, 19, 23};
    int s = __builtin_ctzll(n - 1);
    i64 d = (n - 1) >> s;
    for (auto a : A) {
        if (a == n)
            return true;
        i64 x = power(a, d, n);
        if (x == 1 || x == n - 1)
            continue;
        bool ok = false;
        for (int i = 0; i < s - 1; ++i) {
            x = mul(x, x, n);
            if (x == n - 1) {
                ok = true;
                break;
            }
        }
        if (!ok)
            return false;
    }
    return true;
}
std::vector<std::pair<i64, int>> factorize(i64 n) {
    std::vector<i64> p;
    std::function<void(i64)> f = [&](i64 n) {
        if (n <= 10000) {
            for (int i = 2; i * i <= n; ++i)
                for (; n % i == 0; n /= i)
                    p.push_back(i);
            if (n > 1)
                p.push_back(n);
            return;
        }
        if (isprime(n)) {
            p.push_back(n);
            return;
        }
        auto g = [&](i64 x) {
            return (mul(x, x, n) + 1) % n;
        };
        i64 x0 = 2;
        while (true) {
            i64 x = x0;
            i64 y = x0;
            i64 d = 1;
            i64 power = 1, lam = 0;
            i64 v = 1;
            while (d == 1) {
                y = g(y);
                ++lam;
                v = mul(v, std::abs(x - y), n);
                if (lam % 127 == 0) {
                    d = std::gcd(v, n);
                    v = 1;
                }
                if (power == lam) {
                    x = y;
                    power *= 2;
                    lam = 0;
                    d = std::gcd(v, n);
                    v = 1;
                }
            }
            if (d != n) {
                f(d);
                f(n / d);
                return;
            }
            ++x0;
        }
    };
    f(n);
    std::sort(p.begin(), p.end());
    std::vector<std::pair<i64, int>> ret;
    for (int l = 0, r = 0; l < p.size(); l = r) {
        while (r < p.size() && p[l] == p[r]) {
            r++;
        }
        ret.emplace_back(p[l], r - l);
    }
    return ret;
}
i64 get(i64 n, i64 p) {
    i64 ans = 0;
    for (n /= p; n; n /= p) {
        ans += n;
    }
    return ans;
}
void solve() {
    int N;
    i64 X, Y;
    std::cin >> N >> X >> Y;
    
    std::vector<i64> a(N);
    for (int i = 0; i < N; i++) {
        std::cin >> a[i];
    }
    
    auto ps = factorize(X);
    i64 ans = 4E18;
    for (auto [p, e] : ps) {
        i64 res = get(Y, p);
        for (int i = 0; i < N; i++) {
            res -= get(a[i], p);
        }
        ans = std::min(ans, res / e);
    }
    std::cout << ans << "\n";
}

int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    
    int t;
    std::cin >> t;
    
    while (t--) {
        solve();
    }
    
    return 0;
}