hdu6287可持久线段树+线性筛

传送门

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
int n,m;
int pri[maxn],vis[maxn],pos[maxn],c=0;//质数表，最大质因子及其序号 
int tot=-1,ver[maxn];
struct st{
	int l,r,s;//s为质因子出现次数 
}tre[maxn<<8];
void init(){//线性筛 
	memset(pri,0,sizeof(pri));
	memset(vis,0,sizeof(vis));
	for(int i=2;i<maxn;i++){
		if(!vis[i]){
			pri[++c]=i;
			vis[i]=i;
			pos[i]=c;
		}
		for(int j=1;j<=c;j++){
			if(i*pri[j]>=maxn||pri[j]>vis[i]) break;
			vis[i*pri[j]]=pri[j];
			pos[i*pri[j]]=j;
		}
	}
}
int build(int l,int r){
	int root=++tot;
	tre[root].l =root;
	tre[root].r =root;
	tre[root].s =0;
	return root;
}
int update(int node,int l,int r,int p,int v){
	int root=++tot;
	st &rt=tre[root];
	rt=tre[node];
	if(l==r){
		rt.s +=v;
		return root;
	}
	int mid=(l+r)>>1;
	if(p<=mid) rt.l =update(rt.l ,l,mid,p,v);
	else rt.r =update(rt.r ,mid+1,r,p,v);
	rt.s =tre[rt.l ].s +tre[rt.r ].s ;
	return root;
}
int query(int t1,int t2,int l,int r,int p){
	st &s1=tre[t1];
	st &s2=tre[t2];
	if(l==r) return s2.s -s1.s ;//找到叶子结点(即质因子)，就返回该质因子出现的次数
	int mid=(l+r)>>1;
	if(p<=mid) return query(s1.l ,s2.l ,l,mid,p);
	else return query(s1.r ,s2.r ,mid+1,r,p); 
}
int main()
{
	int t,x,y,z,p,ok,l,r;
	init();
	scanf("%d",&t);
	while(t--){
		scanf("%d%d",&n,&m);
		memset(ver,0,sizeof(ver));
		ver[0]=build(1,c);//建质因子树 
		for(int i=1;i<=n;i++){
			scanf("%d",&x);
			ver[i]=ver[i-1];//x为1时，没有质数 
			while(x>1){
				y=0;
				z=vis[x];
				p=pos[x];
				while(x%z==0){
					y++;
					x=x/z;
				}
				ver[i]=update(ver[i],1,c,p,y);//p是质因子的位置，y是质因子出现的次数 
			}
		}
		for(int i=1;i<=m;i++){
			scanf("%d%d%d",&l,&r,&x);
			ok=1;
			while(x>1){
				y=0;
				z=vis[x];
				p=pos[x];//p标记第几个质因子 
				while(x%z==0){
					y++;
					x=x/z;
				}
				if(y>query(ver[l-1],ver[r],1,c,p)){//分母的该质因子数>分子的该质因子数，不能整除 
					ok=0;
					x=1;
				}
			}
			if(ok==1) printf("Yes\n");
			else printf("No\n");
		}
	}
	return 0;
}