bzoj2506 calc

值域分块，题解网上有。。

根本不会，总是往莫队上去想，但事实上这题询问满足区间可减性，有更简单的方法

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
using namespace std;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
struct Q
{
    int l,p,k,tp,nm;
}q[201000];
bool operator<(const Q &a,const Q &b)
{
    return a.l<b.l;
}
int n,m,nq;
int a[101000],an[101000];
int ma[101][101],ma2[10100];
void add(int p)
{
    for(int i=1;i<=100;i++)
        ma[i][a[p]%i]++;
    ma2[a[p]]++;
}
int main()
{
    int i,j,l,r,p,k,tt;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)    scanf("%d",&a[i]);
    for(i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&l,&r,&p,&k);
        q[++nq].l=l-1;q[nq].p=p;q[nq].k=k;q[nq].tp=-1;q[nq].nm=i;
        q[++nq].l=r;q[nq].p=p;q[nq].k=k;q[nq].tp=1;q[nq].nm=i;
    }
    sort(q+1,q+nq+1);
    for(i=1,j=0;i<=nq;i++)
    {
        while(j+1<=q[i].l)    add(++j);
        if(q[i].p<=100)    an[q[i].nm]+=q[i].tp*ma[q[i].p][q[i].k];
        else
        {
            tt=0;
            for(k=q[i].k;k<=10000;k+=q[i].p)
                tt+=ma2[k];
            an[q[i].nm]+=q[i].tp*tt;
        }
    }
    for(i=1;i<=m;i++)    printf("%d
",an[i]);
    return 0;
}