BZOJ2506: calc

2506: calc

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 268  Solved: 144
[Submit][Status]

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福建集训

题解：
一看题目感觉就是不可做神题，看题解发现真是。。。
首先要把ans(l,r)=ans(1,r)-ans(1,l-1)，然后所有的询问都是从1开始的了，然后我们离线处理。
按终点排序，然后如果p<=100，我们直接暴力开数组记录这样的值有多少个就可以，否则我们枚举mod p==k的数（这样的数很少），加上它到目前为止出现的次数。
代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 200000+5

#define maxm 500+100

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define mod 1000000007

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
int n,m,tot,mx,ans[maxn][2],a[maxn],f[105][105],g[maxn];
struct rec{int x,p,q,id,f;}b[maxn];
inline void add(int x)
{
    for1(i,100)f[i][x%i]++;
    g[x]++;
}
inline bool cmp(rec a,rec b){return a.x<b.x;}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    n=read();m=read();
    for1(i,n)a[i]=read(),mx=max(mx,a[i]);
    for1(i,m)
    {
        int x=read(),y=read(),p=read(),q=read();
        b[++tot]=(rec){x-1,p,q,i,0};
        b[++tot]=(rec){y,p,q,i,1};
    }
    sort(b+1,b+tot+1,cmp);
    int now=0;
    for1(i,tot)
    {
        while(now<b[i].x){now++;add(a[now]);}
        int p=b[i].p,q=b[i].q;
        if(p<=100)ans[b[i].id][b[i].f]=f[p][q];
        else for(int j=q;j<=mx;j+=p)ans[b[i].id][b[i].f]+=g[j];
    }
    for1(i,m)printf("%d\n",ans[i][1]-ans[i][0]);

    return 0;

}