bzoj2506

不知道该给这题贴什么标签了，但是这是一道好题
1.注意这道题的询问是满足区间减法的，我们可以把一个询问拆成两个询问离线搞搞，这个思想在后面会见到
2.分类的思想，这道题的求相同余数看似很烦，实际上我们发现每个数只有10^4，不难发现
当p>=100时，我们可以完全不虚的穷举k,k+p*1,k+p*2……，计算每个数出现的次数
当p<100时,我们完全可以直接记录f[p,k]代表到现在出现的次数
因此这题做法就很明显了

type node=record
x,w,p,k:longint;
end;

var q:array[0..200010] of node;
a,f2:array[0..100010] of longint;
f1:array[0..100,0..100] of longint;
ans:array[0..100010,0..1] of longint;
mx,j,l,r,p,k,n,m,i,t:longint;

procedure swap(var a,b:node);
var c:node;
begin
c:=a;
a:=b;
b:=c;
end;

procedure sort(l,r:longint);
var i,j,x:longint;
begin
i:=l;
j:=r;
x:=q[(l+r) shr 1].x;
repeat
while q[i].x<x do inc(i);
while x<q[j].x do dec(j);
if not(i>j) then
begin
swap(q[i],q[j]);
inc(i);
dec(j);
end;
until i>j;
if l<j then sort(l,j);
if i<r then sort(i,r);
end;

begin
readln(n,m);
for i:=1 to n do
begin
read(a[i]);
if mx<a[i] then mx:=a[i];
end;
for i:=1 to m do
begin
readln(l,r,p,k);
inc(t); q[t].p:=p; q[t].k:=k; q[t].x:=r; q[t].w:=i*2;
inc(t); q[t].p:=p; q[t].k:=k; q[t].x:=l-1; q[t].w:=i*2-1;
end;
sort(1,2*m);
t:=1;
while q[t].x=0 do inc(t);
for i:=1 to n do
begin
inc(f2[a[i]]);
for j:=1 to 100 do
inc(f1[j,a[i] mod j]);

while q[t].x=i do
begin
if q[t].p>100 then
begin
for j:=0 to mx div q[t].p do
inc(ans[q[t].w shr 1+q[t].w mod 2,q[t].w mod 2],f2[q[t].p*j+q[t].k]);
end
else ans[q[t].w shr 1+q[t].w mod 2,q[t].w mod 2]:=f1[q[t].p,q[t].k];
inc(t);
end;
end;
for i:=1 to m do
writeln(ans[i,0]-ans[i,1]);
end.

View Code