[bzoj2506] calc

　　莫队+暴力。。

　　对于>100的模数p，每次暴力查k,p+k,2p+k.....就好了。

　　对于<=100的模数，就在莫队的时候处理出来就好了。。。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100233;
struct zs{int r,id,p,k;}a[maxn<<1];
int i,j,k,n,m,mx,mx1;
int mp[maxn],ans[maxn<<1];
int now[103][103],sm[10023];
 
int ra;char rx;
inline int read(){
    rx=getchar(),ra=0;
    while(rx<'0'||rx>'9')rx=getchar();
    while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra;
}
bool cmp(zs a,zs b){return a.r<b.r;}
int main(){
    n=read(),m=read();
    for(i=1;i<=n;i++)mp[i]=read(),mx=max(mx,mp[i]);mx1=min(70,mx);
    for(i=1;i<=m;i++)
        a[i+m].r=read()-1,a[i].r=read(),a[i+m].p=a[i].p=read(),a[i+m].k=a[i].k=read(),
        a[i+m].id=i+m,a[i].id=i;
    sort(a+1,a+1+m+m,cmp);int top=1,m1=m<<1;register int j;
    while(top<=m1&&a[top].r==0)top++;
    for(i=1;i<=n;i++){
        sm[mp[i]]++;
        for(j=1;j<=mx1;j++)now[j][mp[i]%j]++;
        while(top<=m1&&a[top].r==i){
            int tmp=0;
            if(a[top].p<=mx1)tmp=now[a[top].p][a[top].k];
            else for(j=a[top].k;j<=mx;j+=a[top].p)tmp+=sm[j];
            ans[a[top].id]=tmp;
            top++;
        }
        if(top>m1)break;
    }
    for(i=1;i<=m;i++)printf("%d\n",ans[i]-ans[i+m]);
}