BZOJ2506 calc

首先我们发现这个区间特别的蛋疼……所以我们离线，把询问拆成两个，然后从前往后扫整个序列

对于p<=100的，拿个数组记录所有模p得k的，到时候直接出解，对于p>100的，开个数组记录值为多少的都有多少个，然后暴力#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<map>
#include<set>
#include<bitset>
#include<queue>
#include<stack>
using namespace std;
#define MAXN 200010
#define MAXM 110
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
#define ll long long
struct que{
int x;
int p;
int k;
int bel;
int f;
que(){

}
que(int _x,int _p,int _k,int _bel,int _f){
x=_x;
p=_p;
k=_k;
bel=_bel;
f=_f;
}
friend bool operator <(que x,que y){
return x.x<y.x;
}
};
int n,m;
int a[MAXN];
que q[MAXN];
int s[MAXM][MAXM];
int c[MAXN];
int ans[MAXN];
int main(){
freopen("data.txt","r",stdin);
freopen("dui.txt","w",stdout);
int i,j;
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",&q[i*2-1].x,&q[i*2].x,&q[i*2-1].p,&q[i*2-1].k);
q[i*2].p=q[i*2-1].p;
q[i*2].k=q[i*2-1].k;
q[i*2].bel=q[i*2-1].bel=i;
q[i*2-1].x--;
q[i*2-1].f=-1;
q[i*2].f=1;
}
sort(q+1,q+2*m+1);
int wzh=1;
while(q[wzh].x==0){
wzh++;
}
for(i=1;i<=n;i++){
for(j=1;j<=100;j++){
s[j][a[i]%j]++;
}
c[a[i]]++;
while(q[wzh].x==i){
if(q[wzh].p<=100){
ans[q[wzh].bel]+=q[wzh].f*s[q[wzh].p][q[wzh].k];
}else{
for(j=q[wzh].k;j<=10000;j+=q[wzh].p){
ans[q[wzh].bel]+=q[wzh].f*c[j];
}
}
wzh++;
}
}
for(i=1;i<=m;i++){
printf("%d\n",ans[i]);
}
return 0;
}

/*
5 2
1 5 2 3 7
1 3 2 1
2 5 3 0
*/