#pragma GCC optimize ("Ofast") #include using namespace std; void *wmem; template inline S min_L(S a,T b){ return a<=b?a:b; } template inline void walloc1d(T **arr, int x, void **mem = &wmem){ static int skip[16]={0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; (*mem) = (void*)( ((char*)(*mem)) + skip[((unsigned long long)(*mem)) & 15] ); (*arr)=(T*)(*mem); (*mem)=((*arr)+x); } inline void rd(int &x){ int k, m=0; x=0; for(;;){ k = getchar_unlocked(); if(k=='-'){ m=1; break; } if('0'<=k&&k<='9'){ x=k-'0'; break; } } for(;;){ k = getchar_unlocked(); if(k<'0'||k>'9'){ break; } x=x*10+k-'0'; } if(m){ x=-x; } } inline void wt_L(char a){ putchar_unlocked(a); } inline void wt_L(const char c[]){ int i=0; for(i=0;c[i]!='\0';i++){ putchar_unlocked(c[i]); } } int Prime_L(int N, int res[], void *mem=wmem){ bool *isprime; const int r=23000; int a, b, i, *sf, ss=1, sz=1; walloc1d(&isprime, r, &mem); walloc1d(&sf, r, &mem); isprime = (bool*)mem; sf = (int*)(isprime + r); N /= 2; res[0] = 2; b =min_L(r, N); for(i=(1);i<(b);i++){ isprime[i] = 1; } for(i=(1);i<(b);i++){ if(isprime[i]){ res[sz++] = 2*i+1; sf[ss] = 2*i*(i+1); if(sf[ss] < N){ while(sf[ss] < r){ isprime[sf[ss]] = 0; sf[ss] += res[ss]; } ss++; } } } for(a=r; a 1){ res[i] = 0; } } for(i=0;i<(Q);i++){ if(res[i]){ wt_L("Yes"); wt_L('\n'); } else{ wt_L("NO"); wt_L('\n'); } } return 0; } // cLay varsion 20190820-1 // --- original code --- // int N, A[1d5]; // // int Q, P[1d5], L[1d5], R[1d5]; // int res[1d5]; // // int ok[100001], fc[100001]; // // int ps, prm[1000]; // int cnt[1d5]; // // { // int i, j, k; // // ps = Prime(2000, prm); // // rd(N,A(N),Q,(P,L--,R)(Q)); // rep(i,Q) res[i] = 1; // // rep(i,N) if(A[i]==0) ok[i+1] = 1; // // rep(i,N) ok[i+1] += ok[i]; // rep(i,Q) if(ok[R[i]]-ok[L[i]]){ res[i] = 2; continue; } // // rep(k,ps){ // j = 0; // rep(i,Q) if(res[i]==1){ // cnt[i] = 0; // while(P[i]%prm[k]==0) P[i] /= prm[k], cnt[i]++, j++; // } // if(j==0) continue; // // rep(i,N+1) fc[i] = 0; // rep(i,N){ // if(ok[i+1]-ok[i]) continue; // while(A[i]%prm[k]==0) A[i] /= prm[k], fc[i+1]++; // } // rep(i,N) fc[i+1] += fc[i]; // // rep(i,Q) if(res[i]==1){ // if(fc[R[i]]-fc[L[i]] < cnt[i]) res[i] = 0; // } // } // // rep(i,Q) if(res[i]==1 && P[i] > 1) res[i] = 0; // // rep(i,Q) wt(if[res[i],"Yes","NO"]); // }