【Hackerrank (70)】—Range Modular Query（莫队+暴力）

传送门

有说保证数据随机或者不充分么。。。
我为什么没看到……

直接调和级数证明暴力复杂度是$O(nsqrt nlogn)$的

#include<bits/stdc++.h>
using namespace std;
const int RLEN=1<<18|1;
inline char gc(){return getchar();
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0,f=1;
	while(!isdigit(ch))ch=gc();
	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
	return res;
}
const int M=10004;
const int N=100005;
int n,m,plc[N],a[N],buc[M],blo,ans[N];
struct ask{
	int l,r,p,k,idx;
	friend inline bool operator <(const ask &a,const ask &b){
		return plc[a.l]==plc[b.l]?(plc[a.l]&1?a.r<b.r:a.r>b.r):a.l<b.l;
	}
}p[N];
inline void update(int p,int k){
	buc[a[p]]+=k;
}
int main(){
	n=read(),m=read(),blo=sqrt(n);
	for(int i=1;i<=n;i++)plc[i]=(i-1)/blo+1;
	for(int i=1;i<=n;i++)a[i]=read();
	for(int i=1;i<=m;i++)
		p[i].l=read(),p[i].r=read(),p[i].p=read(),p[i].k=read(),p[i].idx=i;
	int l=1,r=0;sort(p+1,p+m+1);
	for(int i=1;i<=m;i++){
		while(l<p[i].l)update(l++,-1);
		while(l>p[i].l)update(--l,1);
		while(r<p[i].r)update(++r,1);
		while(r>p[i].r)update(r--,-1);
		int res=0;
		for(int j=p[i].k;j<M;j+=p[i].p)res+=buc[j];
		ans[p[i].idx]=res;
	}
	for(int i=1;i<=m;i++)cout<<ans[i]<<'
';
}