算法模板——splay区间反转 2

实现功能：同splay区间反转 1（基于BZOJ3223 文艺平衡树）
这次改用了一个全新的模板（HansBug：琢磨了我大半天啊有木有），大大简化了程序，同时对于splay的功能也有所完善
这里面没有像一般二叉排序树那样子用一个参量进行排序，而是直接以中序遍历来构建了一个普通的二叉树（当然也可以把每个点的中序遍历排名视作参量），然后插入的时候就是指定位置插入（这个就比较像是文本插入了）
总之得到了较大的提升，代码优美程度也提高不少
  1 var
  2    i,j,k,l,m,n,head,tot,ll:longint;
  3    a,lef,rig,b,c:array[0..100000] of longint;
  4 procedure swap(var x,y:longint);
  5           var z:longint;
  6           begin
  7                z:=x;x:=y;y:=z;
  8           end;
  9 procedure ext(x:longint);  //经典的下推操作，有木有很像线段树呢= =
 10           begin
 11                if x=0 then exit;
 12                if c[x]=0 then exit;
 13                if lef[x]<>0 then c[lef[x]]:=1-c[lef[x]];
 14                if rig[x]<>0 then c[rig[x]]:=1-c[rig[x]];
 15                swap(lef[x],rig[x]);
 16                c[x]:=0;
 17           end;
 18 procedure rt(var x:longint);  //lt、rt函数和Treap的极其相似!!!
 19           var f,l:longint;
 20           begin
 21                ext(x);ext(lef[x]);
 22                if (x=0) or (lef[x]=0) then exit;
 23                f:=x;l:=lef[x];
 24                b[lef[x]]:=b[x];
 25                b[x]:=1+b[rig[x]]+b[rig[lef[x]]];
 26                lef[f]:=rig[l];
 27                rig[l]:=f;
 28                x:=l;
 29           end;
 30 procedure lt(var x:longint);
 31           var f,r:longint;
 32           begin
 33                ext(x);ext(rig[x]);
 34                if (x=0) or (rig[x]=0) then exit;
 35                f:=x;r:=rig[x];
 36                b[rig[x]]:=b[x];
 37                b[x]:=1+b[lef[x]]+b[lef[rig[x]]];
 38                rig[f]:=lef[r];
 39                lef[r]:=f;
 40                x:=r;
 41           end;
 42 procedure splay(var head:longint;x:longint);   //萌萌哒伸展——而且还可以支持伸展的任意位置，只要设定head的位置，比如像后面的splay(rig[head],x);
 43           begin
 44                if head=0 then exit;
 45                ext(head);
 46                if (b[lef[head]]+1)=x then exit;
 47                if x<(b[lef[head]]+1) then
 48                   begin
 49                        ext(lef[head]);
 50                        if x=(b[lef[lef[head]]]+1) then rt(head) else
 51                           if x<(b[lef[lef[head]]]+1) then
 52                              begin
 53                                   splay(lef[lef[head]],x);
 54                                   rt(head);rt(head);
 55                              end
 56                           else
 57                               begin
 58                                    splay(rig[lef[head]],x-1-b[lef[lef[head]]]);
 59                                    lt(lef[head]);rt(head);
 60                               end;
 61                   end
 62                else
 63                    begin
 64                         ext(rig[head]);
 65                         x:=x-1-b[lef[head]];
 66                         if x=(b[lef[rig[head]]]+1) then lt(head) else
 67                            if x<(b[lef[rig[head]]]+1) then
 68                               begin
 69                                    splay(lef[rig[head]],x);
 70                                    rt(rig[head]);lt(head);
 71                               end
 72                            else
 73                                begin
 74                                     splay(rig[rig[head]],x-1-b[lef[rig[head]]]);
 75                                     lt(head);lt(head);
 76                                end;
 77                    end;
 78           end;
 79 procedure ins(x,y:longint);  //这里的ins操作已经可以支持指定位置插入子伸展树了
 80           begin
 81                if head=0 then
 82                   begin
 83                        head:=y;
 84                        exit;
 85                   end;
 86                if x=b[head] then
 87                   begin
 88                        splay(head,b[head]);
 89                        rig[head]:=y;
 90                        inc(b[head],b[y]);
 91                   end
 92                else if x=0 then
 93                     begin
 94                          splay(head,1);
 95                          lef[head]:=y;
 96                          inc(b[head],b[y]);
 97                     end
 98                else begin
 99                     splay(head,x);
100                     splay(rig[head],x+1);
101                     lef[rig[head]]:=y;
102                     inc(b[rig[head]],b[y]);
103                     inc(b[head],b[y]);
104                end;
105           end;
106 procedure putinvalue(x,y:longint);  //加入单个值
107          begin
108               inc(tot);
109               A[TOT]:=x;b[tot]:=1;
110               lef[tot]:=0;rig[tot]:=0;
111               ins(y,tot);
112          end;
113  
114 procedure turnover(x,y:longint);  //翻转，核心思想还是打标记
115           begin
116                if x=y then exit;
117                if (x=1) and (y=n) then
118                   begin
119                        c[head]:=1-c[head];
120                        exit;
121                   end;
122                if (x=1) then
123                   begin
124                        splay(head,y+1);
125                        c[lef[head]]:=1-c[lef[head]];
126                   end
127                else if (y=n) then
128                     begin
129                          splay(head,x-1);
130                          c[rig[head]]:=1-c[rig[head]];
131                     end
132                else begin
133                     splay(head,x-1);
134                     splay(rig[head],y-x+2);
135                     c[lef[rig[head]]]:=1-c[lef[rig[head]]];
136                end;
137           end;
138 function showtree(head:longint):ansistring;
139          var s1,s2,s3,s4:ansistring;
140          begin
141               if head=0 then exit('');
142               str(a[head],s1);
143               s2:=showtree(lef[head]);
144               s3:=showtree(rig[head]);
145               if c[head]=1 then
146                  begin
147                       s4:=s2;
148                       s2:=s3;
149                       s3:=s4;
150                  end;
151               if s3<>'' then s3:=','+s3;
152               s2:=s2+s3;
153               if s2<>'' then s2:='('+s2+')';
154               exit(s1+s2);
155          end;
156 procedure putout(x:longint);
157           begin
158                if x=0 then exit;
159                ext(x);
160                putout(lef[x]);
161                write(a[x],' ');
162                putout(rig[x]);
163           end;
164 begin
165      readln(n,m);head:=0;tot:=0;
166      for i:=1 to n do putinvalue(i,i-1);
167      fillchar(c,sizeof(c),0);
168      for i:=1 to m do
169          begin
170               readln(j,k);
171               turnover(j,k);
172               l:=0;
173  
174          end;
175      putout(head);
176      writeln;
177      readln;
178 end.