#include <bits/stdc++.h>
#define endl '\n'
#define F first
#define S second
using namespace std;
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace __gnu_pbds;
#define ll long long
#define pb push_back
#define mod 1000000007
// #define int long long
#define all(x)      x.begin(),x.end()
#define allr(x)     x.rbegin(),x.rend()
#define CheckBit(x,k)   (x & (1LL << k))
#define SetBit(x,k)     (x |= (1LL << k))
#define ClearBit(x,k)   (x &= ~(1LL << k))
#define LSB(mask)       __builtin_ctzll(mask)
#define MSB(mask)       63-__builtin_clzll(mask) 
#define print(x)    cout << #x << " : " << x << endl
#define error1(x)   cerr << #x << " = " << (x) <<endl
#define coutall(v)  for(auto &it: v) cout<<it<<' '; cout<<endl
#define Abid_52     ios::sync_with_stdio(false);cin.tie(0),cin.tie(0)
#define error2(a,b) cerr<<"( "<<#a<<" , "<<#b<<" ) = ( "<<(a)<<" , "<<(b)<<" )\n"
#define UNIQUE(x)   sort(all(x)), x.erase(unique(all(x)), x.end()), x.shrink_to_fit()
template <typename T, typename U> T ceil(T x, U y) {return (x > 0 ? (x + y - 1) / y : x / y);}
template <typename T, typename U> T floor(T x, U y) {return (x > 0 ? x / y : (x - y + 1) / y);}

const int N = 2e5 + 10;

void solve()
{
    int n, m;
    cin >> n >> m;
    vector<int> v(n);
    gp_hash_table<int, vector<int>>mp;
    for (int i = 0; i < n; i++)
    {
        cin >> v[i];
        for (int j = 2; j * j <= v[i]; j++)
        {
            if (v[i] % j == 0)
            {
                while (v[i] % j == 0)
                {
                    mp[j].push_back(i + 1);
                    v[i] /= j;
                }
            }
        }
        if (v[i] > 1)
            mp[v[i]].push_back(i + 1);
    }
    gp_hash_table<int, int> gp;
    for (int j = 2; j * j <= m; j++)
    {
        if (m % j == 0)
        {
            while (m % j == 0)
            {
                gp[j]++;
                m /= j;
            }
        }
    }
    if (m > 1)
        gp[m]++;
    
    vector<pair<int,int>>vec;
    for(auto it : gp){
        vec.push_back({it.first, it.second});
    }
    
    // for(auto it : mp){
    //     cout << it.first << "-->";
    //     for(auto i : it.second){
    //         cout << i << " ";
    //     }
    //     cout << endl;
    // }

    int q, l, r;
    cin >> q;
    while(q--){
        cin >> l >> r;
        bool ok = 0;
        for(auto it : vec)
        {
            auto first = lower_bound(all(mp[it.first]), l);
            if(first == mp[it.first].end()){
                ok = 1;
                break;
            }
            auto sec = upper_bound(all(mp[it.first]), r);
            sec--;
            int left = first - mp[it.first].begin();
            int right = sec - mp[it.first].begin();
            if(right < 0){
                ok = 1;
                break;
            }
            
            if(right < left){
                ok = 1;
                break;
            }
            int ase = right - left + 1;
            if(ase < it.second){
                ok = 1;
                break;
            }
        }
        if(ok){
            cout << "No" << endl;
        }
        else{
            cout << "Yes" << endl;
        }

    }

}

int32_t main()
{
    Abid_52;
    int t = 1;
    cin >> t;
    for (int tc = 1; tc <= t; ++tc)
    {
        // cout << "Case " << tc << ": ";
        solve();
    }
}