【BZOJ2506】calc

Description

     给一个长度为n的非负整数序列A1,A2,…,An。现有m个询问，每次询问给出l,r,p,k，问满足l<=i<=r且Ai mod p = k的值i的个数。

Input
第一行两个正整数n和m。
第二行n个数，表示A1,A2,…,An。
以下m行，每行四个数分别表示l,r,p,k。满足1<=l<=r<=n。

Output
对于每个询问，输出一行，表示可行值i的个数。

Sample Input
5 2
1 5 2 3 7
1 3 2 1
2 5 3 0

Sample Output
2
1
HINT

数据范围：

     0<n,m<=10^5，任意1<=i<=n满足Ai<=10^4，0<p<=10^4，0<=k<p。

Source

2011福建集训

乱搞题>_<
考虑把询问拆分成两个然后用容斥原理得到答案
然后离线按坐标排序
把询问按p大小划分成两类:p<=100,p>100
p<=100的我们可以预处理出每个数对1-100每个数取模结果个数
p>100爆枚,因为a不大于10000所以最多枚100次

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 100010
#define MAXINT 0x7fffffff
#define GET (ch>='0'&&ch<='9')
using namespace std;
int n,m,top;
int a[MAXN],ans[MAXN][2];
int c[110][110],Count[MAXN];
struct Query
{
    int opt,pos,p,k,id;
    bool operator <(const Query& a)const    {return pos<a.pos;}
}q[MAXN<<1];
void in(int &x)
{
    char ch=getchar();x=0;
    while (!GET)    ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();
}
void insert(int x)
{
    for (int i=1;i<=100;i++)    c[i][x%i]++;
    Count[x]++;
}
int main()
{
    in(n);in(m);int l,r,p,k;
    for (int i=1;i<=n;i++)  in(a[i]);
    for (int i=1;i<=m;i++)
    {
        in(l);in(r);in(p);in(k);
        q[++top].pos=l-1;q[top].p=p;q[top].k=k;q[top].id=i;q[top].opt=0;
        q[++top].pos=r;q[top].p=p;q[top].k=k;q[top].id=i;q[top].opt=1;
    }
    sort(q+1,q+top+1);int now=0;
    for (int i=1;i<=top;i++)
    {
        while (now<q[i].pos)    insert(a[++now]);
        if (q[i].p<=100)    ans[q[i].id][q[i].opt]+=c[q[i].p][q[i].k];
        else
            for (int j=q[i].k;j<=10000;j+=q[i].p)   ans[q[i].id][q[i].opt]+=Count[j];
    }
    for (int i=1;i<=m;i++)  printf("%d\n",ans[i][1]-ans[i][0]);
}