loj #6164. 「美团 CodeM 初赛 Round A」数列互质

莫队==每次维护每个数出现的次数，再上链表维护出现次数的次数。

发现不同的出现次数最多根号n个，因为要出现次数有不同的x个，至少要有x*(x-1)/2个数。然后考虑每次询问暴力枚举。

先把k分解成质因数，不超过logk个，分解时要判定大质数，不然GG

对于每个出现次数暴力验证，复杂度

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define rep(j,k,l) for (int j=k;j<=l;++j)
#define red(j,k,l) for (int j=k;j>=l;--j)
#define N 50005

using namespace std;
struct _233{ int l,r,k,id; } q[N];
int n,m,sn,a[N],sum[N],last[N],next[N],Sum[N],ans[N],P[N],cnt,prime[N];
bool isprime[N];

bool cmp(const _233 &x,const _233 &y){ if (x.l/sn==y.l/sn) return x.r<y.r; return x.l<y.l; }
void add(int x){ Sum[x]++; if (Sum[x]==1) next[x]=next[0],next[0]=x,last[x]=0,last[next[x]]=x; }
void del(int x){ Sum[x]--; if (Sum[x]==0) next[last[x]]=next[x],last[next[x]]=last[x]; }

void getprime(){
	
	rep(i,2,n){
		
		if (isprime[i]==0) prime[++cnt]=i;
		for (int j=1;j<=cnt&&prime[j]*i<=n;++j){
			
			isprime[prime[j]*i]=1;
			if (i%prime[j]==0) break;
			
		}
		
	}
	
}

int main(){
	
	scanf("%d%d",&n,&m);sn=(int)sqrt(n);
	rep(i,1,n) scanf("%d",a+i);
	rep(i,1,m) scanf("%d%d%d",&q[i].l,&q[i].r,&q[i].k),q[i].id=i;
	sort(q+1,q+m+1,cmp);
	int l=1,r=1;getprime();
	sum[a[1]]++;next[0]=n+1;last[n+1]=0;add(1);
	rep(i,1,m){
		
		_233 x=q[i];
		while (l>x.l) l--,del(sum[a[l]]),sum[a[l]]++,add(sum[a[l]]);
		while (r<x.r) r++,del(sum[a[r]]),sum[a[r]]++,add(sum[a[r]]);
		while (l<x.l) del(sum[a[l]]),sum[a[l]]--,add(sum[a[l]]),l++;
		while (r>x.r) del(sum[a[r]]),sum[a[r]]--,add(sum[a[r]]),r--;
		int poi=0;cnt=0;
		for (int i=1;isprime[x.k];++i) if (x.k%prime[i]==0){
			
			P[++cnt]=prime[i];
			while (x.k%prime[i]==0) x.k=x.k/prime[i];
			
		}
		if (x.k>1) P[++cnt]=x.k;
		for (int i=next[0];i<=n;i=next[i]){
			
			bool lych=1;
			rep(j,1,cnt) if (i%P[j]==0){ lych=0; break; }
			if (lych) poi+=Sum[i];
			
		}
		ans[x.id]=poi;
		
	}
	rep(i,1,m) printf("%d\n",ans[i]);
	system("pause");
	
}