#include using namespace std; typedef long long ll; #define rep(i,n) for(int i = 0;i<((int)(n));i++) #define reg(i,a,b) for(int i = ((int)(a));i<=((int)(b));i++) #define irep(i,n) for(int i = ((int)(n)-1);i>=0;i--) #define ireg(i,a,b) for(int i = ((int)(b));i>=((int)(a));i--) /* 何故かエラー出る 解決。配列範囲には気をつけよう 倍数かどうかだけ見れば良い。 */ int n,m,a[100010],q,p[100010]; int l[100010],r[100010],fact[2010][305],dp[100010][305],zero[100010]; vector prime,pri2; void prime_list(ll n,vector& prime_){ ll visited[300010]={}; reg(i,2,n){ if(visited[i]==0){ prime_.push_back(i); reg(j,2,n){ if(j*i>n)break; visited[i*j]=1; } } } } void factoring(ll n,vector>& v,vector& p){ rep(i,p.size()){ int t=0; while(n%p[i]==0){ n/=p[i]; t++; } if(t>0)v.push_back({i,t}); } if(n>1)v.push_back({10000,1}); } ll modpow(__int128_t a, ll n,ll mo) { __int128_t r=1; while(n) r=r*((n%2)?a:1)%mo,a=a*a%mo,n>>=1; return r; } void init(){ cin>>n; reg(i,1,n)cin>>a[i]; cin>>q; rep(i,q)cin>>p[i]>>l[i]>>r[i]; prime_list(2000,prime); m = prime.size(); rep(j,m)fact[0][j]=0; reg(i,1,2000){ ll t = i; rep(j,m){ fact[i][j]=0; while(t%prime[j]==0){ // if(prime[j]>=45) t/=prime[j]; fact[i][j]++; } } } // 12ms zero[0]=0; rep(j,m)dp[0][j]=0; for(int i=1;i<=n;++i){ zero[i]=zero[i-1]; if(a[i]==0)zero[i]++; for(int j=0;j> v; factoring(p[i],v,prime); if(zero[r[i]]-zero[l[i]-1]==0){ irep(j,v.size()){ if(v[j].first>m){ ok=false; break; } if(v[j].second > dp[r[i]][v[j].first] - dp[l[i]-1][v[j].first]){ ok=false; break; } } } cout<<(ok?"Yes":"NO")<<"\n"; } return 0; }