#include <stdio.h>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <list>
#include <iterator>
#include <assert.h>
#pragma warning(disable:4996) 
 
typedef long long ll;
#define MIN(a, b) ((a)>(b)? (b): (a))
#define MAX(a, b) ((a)<(b)? (b): (a))
#define LINF 9223300000000000000
#define INF 2140000000
const long long MOD = 1000000007;
using namespace std;

ll gcd(ll a, ll b) {
    if(b == 0) return a;
    return gcd(b,a%b);
}
 
ll lcm(ll a, ll b) {
    ll g = gcd(a,b);
    return a / g * b;   // Be careful not to overflow
}

vector<int> pp;
void preparePrime( int maxp )
{
    int i;
    vector<int> a(maxp+1);
    for(i=2; i<=maxp; i++) {
        int k=i+i;
        while(k<=maxp) {
            a[k]=1;
            k+=i;
        }
    }
    for(i=2; i<=maxp; i++) {
        if(a[i]==0) pp.push_back( i );
    }
    return;
}

void decompositPrime( ll a, vector<pair<ll,ll> >& vv )
{
    //preparePrime( (int)sqrt((double)a)+1 );

    int i;
    for(i=0; i<(int)pp.size(); i++) {
        int count = 0;
        while(a % pp[i] == 0 ) {
            a /= pp[i];
            count ++;
        }
        if(count>0) vv.push_back( make_pair( pp[i], count ) );
        if (a==1) break;
    }
    if(a>1) {
        vv.push_back( make_pair( a, 1 ) );
    }

    return;
}

int s[100005][400];

int main(int argc, char* argv[])
{
    int n;
    scanf("%d", &n);

    preparePrime(2500);
    
    map<int,int> mp;
    int i;
    for(i=0; i<(int)pp.size(); i++) {
        mp[pp[i]]=i;
    }

    vector<pair<ll,ll> > vv;

    for(i=0; i<n; i++) {
        int a;
        scanf("%d", &a);
        vv.clear();
        decompositPrime(a, vv);

        int k;
        for(k=0; k<(int)vv.size(); k++) {
            int id=mp[vv[k].first];
            s[i+1][id]=s[i][id]+vv[k].second;
        }
    }

    int Q;
    scanf("%d", &Q);
    for(i=0; i<Q; i++) {
        int p, l, r;
        scanf("%d %d %d", &p, &l, &r);

        vv.clear();
        decompositPrime(p, vv);
        int k;
        int ng=0;
        for(k=0; k<(int)vv.size(); k++) {
            if( mp.find(vv[k].first) == mp.end() ) {
                ng=1; break;
            }
            int id=mp[vv[k].first];
            if(s[r][id]-s[l-1][id]<vv[k].second) {
                ng=1; break;
            }
        }
        if(ng) {
            printf("NO\n");
        }
        else {
            printf("Yes\n");
        }
    }

    return 0;
}