Description
Solution
比赛的时候冲着60分而来,所以打了20棵线段树,可惜忘了开int64,但是分数位思考为正解提供了一条很好的思路。
对于每一个数位都开一个树状数组,a[num,x]表示第num个树状数组中,模
2num
后值为x,且第num位上的数为1的数的个数,当x=0时,询问区间为
[2i−1,2i−1]
,而当x>0时,区间为
[2i−1−x,2i−1−x]
,因为对于
[2i−1,2i−1]
这个区间是加上了加数后的范围。因为有可能会减成负数,所以我们可以对加数进行处理,先把它改成小于
2i
的数字,处理出正数部分的答案。对于负数部分可以发现,当一些数加上加数后会大于
2i
,但是在取模后会出现在前面的部分
画画就发现了。
Code
var
a:array[1..100000] of longint;
tr:array[1..20,1..2100000] of longint;
c:array[0..25] of longint;
n,m,i,j,x,y,cz,sum:longint;
ans:int64;
function lowbit(x:longint):longint;
begin
exit(x and(-x));
end;
procedure change(bh,x,sum:longint);
begin
while x<c[bh+1] do
begin
tr[bh,x]:=tr[bh,x]+sum;
x:=x+lowbit(x);
end;
end;
function sigma(bh,x:longint):int64;
begin
sigma:=0;
while x>0 do
begin
sigma:=sigma+tr[bh,x];
x:=x-lowbit(x);
end;
exit(sigma);
end;
begin
c[0]:=1;
for i:=1 to 25 do c[i]:=c[i-1]*2;
readln(n,m);
for i:=1 to n do
begin
read(a[i]);
for j:=1 to 20 do change(j,a[i] mod c[j]+1,1);
end;
for i:=1 to m do
begin
readln(cz,x,y);
case cz of
1:
begin
for j:=1 to 20 do change(j,a[x] mod c[j]+1,-1);
a[x]:=y;
for j:=1 to 20 do change(j,a[x] mod c[j]+1,1);
end;
2:
begin
ans:=0;
for j:=1 to 20 do
if (y and c[j-1])=c[j-1] then
begin
sum:=x mod c[j];
ans:=ans+sigma(j,c[j]-sum)*c[j-1];
if c[j-1]-sum>0 then
ans:=ans-sigma(j,c[j-1]-sum)*c[j-1]
else
ans:=ans+(n-sigma(j,c[j-1]+c[j]-sum))*c[j-1];
end;
writeln(ans);
end;
end;
end;
end.