bzoj5358 [Lydsy1805月赛]口算训练（分解质因数+主席树/二分）

给定一个序列，求[l,r]的数的乘积是不是d的倍数。
我们考虑给d分解质因数，如果对于每个质因数这个区间内的个数都大于它的次数，那么就是d的倍数。而1e5以内的数最多有6个不同的质因数。
我们考虑给每个数分解质因数，维护区间和即可。我写的主席树。
其实也可以直接对于每种质因数记下它每一次出现的位置及个数，然后每次询问时去二分。
复杂度都是$O(n(\sqrt n+6logn))$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,rt[N],owo=0,prime[N],tot=0,id[N];
bool notprime[N];
struct node{
    int lc,rc,x;
}tr[N*120];
inline void init(){
    notprime[1]=1;
    for(int i=2;i<=1e5;++i){
        if(!notprime[i]) prime[++tot]=i,id[i]=tot;
        for(int j=1;prime[j]*i<=1e5;++j){
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}
inline void ins(int &p,int l,int r,int x,int val){
    tr[++owo]=tr[p];p=owo;if(l==r){tr[p].x+=val;return;}
    int mid=l+r>>1;
    if(x<=mid) ins(tr[p].lc,l,mid,x,val);
    else ins(tr[p].rc,mid+1,r,x,val);
}
inline int ask(int p1,int p2,int l,int r,int x){
    if(l==r) return tr[p2].x-tr[p1].x;int mid=l+r>>1;
    if(x<=mid) return ask(tr[p1].lc,tr[p2].lc,l,mid,x);
    return ask(tr[p1].rc,tr[p2].rc,mid+1,r,x);
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();init();
    while(tst--){
        n=read();m=read();owo=0;
        for(int i=1;i<=n;++i){
            int x=read(),xx=x;rt[i]=rt[i-1];
            for(int j=1;prime[j]*prime[j]<=xx;++j){
                if(x%prime[j]) continue;int res=0;
                while(x%prime[j]==0) x/=prime[j],++res;
                ins(rt[i],1,tot,j,res);
            }if(x!=1) ins(rt[i],1,tot,id[x],1);
        }while(m--){
            int l=read(),r=read(),x=read(),xx=x;bool flag=1;
            for(int j=1;prime[j]*prime[j]<=xx;++j){
                if(x%prime[j]) continue;int res=0;
                while(x%prime[j]==0) x/=prime[j],++res;
                if(ask(rt[l-1],rt[r],1,tot,j)<res){flag=0;break;}
            }if(!flag){puts("No");continue;}
            if(x!=1){if(ask(rt[l-1],rt[r],1,tot,id[x])<1) flag=0;}
            puts(flag?"Yes":"No");
        }
    }return 0;
}