BZOJ 2506 分块

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int N=;
int a[N],n,m,f[][],g[N],tmp=;
struct Node{int l,r,p,k,ans,id;}ask[N];
bool cmp(Node a,Node b){return a.l<b.l;}
bool cmp2(Node a,Node b){return a.id<b.id;}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++)scanf("%d",&a[i]);
    for(int i=;i<=m;i++)
        scanf("%d%d%d%d",&ask[i].l,&ask[i].r,&ask[i].p,&ask[i].k),ask[i].id=i,ask[i].r++;
    sort(ask+,ask++m,cmp);
    for(int i=;i<=m;i++){
        while(tmp<ask[i].l){
            for(int i=;i<=;i++)f[i][a[tmp]%i]++;
            g[a[tmp]]++;
            tmp++;
        }
        if(ask[i].p<=)ask[i].ans-=f[ask[i].p][ask[i].k];
        else{
            for(int i=;i*ask[i].p+ask[i].k<=;i++)ask[i].ans-=g[i*ask[i].p+ask[i].k];
        }
    }tmp=;
    memset(f,,sizeof(f)),memset(g,,sizeof(g));
    for(int i=;i<=m;i++){
        while(tmp<ask[i].r){
            for(int i=;i<=;i++)f[i][a[tmp]%i]++;
            g[a[tmp]]++;
            tmp++;
        }
        if(ask[i].p<=)ask[i].ans+=f[ask[i].p][ask[i].k];
        else{
            for(int i=;i*ask[i].p+ask[i].k<=;i++)ask[i].ans+=g[i*ask[i].p+ask[i].k];
        }
    }
    sort(ask+,ask++m,cmp2);
    for(int i=;i<=m;i++)printf("%d\n",ask[i].ans);
}