块状链表

 

 

模拟

enum.exe/enum.in/enum.out

题意描述：

有一列整数，共n个。每次可以对这些整数操作，有两种操作：

1 第i个到第j个整数分别加上数p

2 询问这些数中比t小的数的个数。

 

输入文件(enum.in)：

第一行有两个数，n和m(1<=n<=100000,m<=10,000)，表示数的个数和操作数。

第二行n个整数，表示每个数的初始值。

以后m行，每行开始一个数q。

若q为1,则后面跟三个数i，j (i<=j)，表示两个下标； p(-1000=p<=1000)，表示修改的数。

若q为2，则为询问操作，后面跟一个数t。

 

输出文件(enum.out):

对每个询问操作输出数列中比t小的个数。

 

输入样例：

5 3

1 2 3 4 5

2 0

1 1 5 -10

2 0

 

输出样例：

0

5

 

 

//前几天一直很衰，根据rp守恒定律，今天终于爆发了！

这个题目的变态之处是:操作2适合平衡树，操作1适合线段树，加起来之后？线段树套平衡树？

也许可以，但一定很丑！

 

曾经在学线段树的时候，看到过块状链表，那时候觉得很慢，而对于这道题，就不得不选它了.......

 

将1~ n 分为 若干个长度为trunc（sqrt（n））的“块”，然后利用线段树的lazy思想，对于修改操作，先用一个lazy数组记录，等到

查询时在统一做处理。

 

但这道题即使这样，一次查询仍然是O（n）,于是可以考虑修改过程中维护每一段递增（快排就行了），查询时就可以对每一段进行二分了。

 

于是，这道题的一次修改和查询都是O（sqrt（n）*log（sqrt（n）））;（真难写）

 

只要常数弄得好就可以过了@

 

{$inline on} program lmd; uses math; var a,b,old,link,w,lazy:array[0..100100]of longint; bj,n,m,fen,size,i,j,s,l,r,ll,rr,k,p,q,x,sl,sr,ans,tot:longint; function head(x:longint):longint; inline; begin exit((x-1)*size+1); end; function tail(x:longint):longint; inline; begin if x<fen then exit(x*size); exit(n); end; procedure sort(l,r:longint); inline; var i,j:longint; begin if r-l<10 then begin for i:=l to r-1 do begin bj:=i; for j:=i+1 to r do if b[bj]>b[j] then bj:=j; if bj<>i then begin k:=b[bj];b[bj]:=b[i];b[i]:=k; k:=old[bj];old[bj]:=old[i];old[i]:=k; end; end; exit; end; i:=l;j:=r;x:=b[(l+r)shr 1]; repeat while b[i]<x do inc(i); while b[j]>x do dec(j); if i<=j then begin k:=b[j];b[j]:=b[i];b[i]:=k; k:=old[j];old[j]:=old[i];old[i]:=k; i:=i+1;j:=j-1; end; until i>j; if l<j then sort(l,j); if i<r then sort(i,r); end; function got(x,q:longint):longint; inline; var l,r,z:longint; begin l:=head(x);r:=tail(x); while l<r do begin z:=(l+r+1)shr 1; if b[z]>=q then r:=z-1 else l:=z; end; if b[l]>=q then l:=l-1; exit(l-head(x)+1); end; begin assign(input,'enum.in'); assign(output,'enum.out'); reset(input);rewrite(output); read(n,m); size:=trunc(sqrt(n)); fen:=trunc(n/size); if fen*size<n then fen:=fen+1; bj:=1;tot:=0; for i:=1 to n do begin read(a[i]); w[i]:=bj; tot:=tot+1; if tot=size then begin tot:=0; inc(bj); end; end; for i:=1 to fen do begin l:=head(i); r:=tail(i); for j:=l to r do begin b[j]:=a[j]; old[j]:=j; end; sort(l,r); for j:=l to r do link[old[j]]:=j; end; for i:=1 to m do begin read(k); if k=1 then begin read(l,r,p); if r-l>size then begin sl:=w[l]+1;sr:=w[r]-1; for s:=sl to sr do lazy[s]:=lazy[s]+p; //==========left======= rr:= tail(w[l]); for j:=l to min(rr,r) do b[link[j]]:=b[link[j]]+p; sort(head(w[l]),rr); for j:=head(w[l]) to rr do link[old[j]]:=j; //=========end=========== //=========right========= ll:= head(w[r]); for j:=max(ll,l) to r do b[link[j]]:=b[link[j]]+p; sort(ll,tail(w[r])); for j:=ll to tail(w[r]) do link[old[j]]:=j; //==========end========== end else begin for j:=l to r do b[link[j]]:=b[link[j]]+p; if w[l]=w[r] then begin ll:=head(w[l]);rr:=tail(w[l]); sort(ll,rr); for j:=ll to rr do link[old[j]]:=j; end else begin ll:=head(w[l]);rr:=tail(w[l]); sort(ll,rr); for j:=ll to rr do link[old[j]]:=j; ll:=head(w[r]);rr:=tail(w[r]); sort(ll,rr); for j:=ll to rr do link[old[j]]:=j; end; end; end else begin read(q); ans:=0; for s:=1 to fen do ans:=ans+got(s,q-lazy[s]); writeln(ans); end; end; close(input);close(output); end.