noi2005 维护序列。 splay

挑战自己写了一下这道noi历史上最变态的数据结构，

noi2005的维护序列，非常变态非常变态，在经过4小时的debug后一次ac了= =

如果再赛场上以我现在水平一定搞不出来。。。

吃一堑长一智，这题对平衡树的操作更加深入、

由于oj上的内存限制只有64mb。。。所以为了内存，加了个回收栈。

后来又加了个回收栈，

回收栈指的是被删除的节点进入栈中，需要从栈中拿出，栈为空那个则再新增节点。。

从回收栈取出节点的时候要清除节点的标记。。

节约很大的空间（尽管noi时期不用管它）

 

 

这题改天要重新写写。。

 

 

之前一直不敢写是因为有翻转操作不敢写，翻转则是给根结点一个标记，左右儿子交换，值传递，

其余的与普通的splay值传递是一样的。。

注意给给前后各加一个点，max【0】=-inf。这样进行初始化。。。

 

 

 

 

 

const inf=maxlongint shr 3; maxn=1000000; var b,p,bj1,bj2,st,sum,a,s,rt,lt,max,lx,rx:array[0..maxn]of longint; root,t,top,tot,n,tmp:longint; procedure cover(z,x:longint); begin sum[z]:=s[z]*x; if x>0 then begin max[z]:=sum[z]; lx[z]:=sum[z]; rx[z]:=sum[z]; end else begin max[z]:=x; lx[z]:=0; rx[z]:=0; end; end; function bigger(x,y:longint):longint; begin if x>y then exit(x) else exit(y); end; procedure swap(X:longint); var tmp:longint; begin tmp:=rt[x];rt[x]:=lt[x];lt[x]:=tmp; tmp:=rx[x];rx[x]:=lx[x];lx[x]:=tmp; end; procedure updata(X:longint); begin s[x]:=s[lt[x]]+s[rt[x]]+1; sum[x]:=sum[lt[x]]+sum[rt[x]]+a[x]; max[x]:=bigger(max[lt[x]],max[rt[x]]); max[x]:=bigger(max[x],rx[lt[x]]+lx[rt[x]]+a[x]); lx[x]:=bigger(lx[lt[x]],sum[lt[x]]+a[x]+lx[rt[x]]); rx[x]:=bigger(rx[rt[x]],sum[rt[x]]+a[x]+rx[lt[x]]); end; procedure mark(X:longint); begin if bj2[x]<>0 then begin if lt[x]<>0 then begin bj2[lt[x]]:=bj2[lt[x]] xor 1; swap(lt[x]); end; if rt[x]<>0 then begin bj2[rt[x]]:=bj2[rt[x]] xor 1; swap(rt[x]); end; bj2[x]:=0; end; if bj1[x]<>-inf then begin if lt[x]<>0 then begin bj1[lt[x]]:=bj1[x]; a[lt[x]]:=bj1[x]; COVER(lt[x],bj1[x]); end; if rt[x]<>0 then begin bj1[rt[x]]:=bj1[x]; a[rt[x]]:=bj1[x]; cover(rt[x],bj1[x]); end; bj1[x]:=-inf; end; end; procedure same(x:longint); var y,z:longint; begin y:=lt[root]; z:=rt[y]; mark(z); a[z]:=x; bj1[z]:=x; cover(z,bj1[z]); updata(y); updata(root); end; procedure rev; var y,z:longint; begin y:=lt[root]; z:=rt[y]; mark(z); swap(z); bj2[z]:=bj2[z] xor 1; updata(y); updata(root); end; procedure right(var X:longint); var tmp:longint; begin tmp:=lt[x]; mark(x); mark(tmp); lt[x]:=rt[tmp]; rt[tmp]:=x; updata(x); updata(tmp); x:=tmp; end; procedure left(var X:longint); var tmp:longint; begin tmp:=rt[x]; mark(x); mark(tmp); rt[x]:=lt[tmp]; lt[tmp]:=x; updata(x); updata(tmp); x:=tmp; end; procedure splay(var root:longint;k:longint); begin mark(root); if s[lt[root]]+1=k then exit; if s[lt[root]]>=k then begin if s[lt[lt[root]]]+1>k then begin right(root); splay(lt[root],k); end else if s[lt[lt[root]]]+1<k then splay(lt[root],k); right(root); end else begin if s[lt[root]]+s[lt[rt[root]]]+2<k then begin left(root); splay(rt[root],k-s[lt[root]]-1); end else if s[lt[root]]+s[lt[rt[root]]]+2>k then splay(rt[root],k-s[lt[root]]-1); left(root); end; end; procedure build(var now:longint;l,r:longint); begin if l>r then begin now:=0; exit;end; if l=r then begin now:=b[l]; sum[b[l]]:=a[b[l]]; max[b[l]]:=a[b[l]]; if a[b[l]]>0 then begin lx[b[l]]:=a[b[l]]; rx[b[l]]:=a[b[l]]; end else begin lx[b[l]]:=0; rx[b[l]]:=0; end; s[b[l]]:=1; exit; end; now:=(l+r) shr 1; build(lt[b[now]],l,now-1); build(rt[b[now]],now+1,r); now:=b[now]; updata(now); end; procedure clear(var x:longint); begin if x=0 then exit; inc(top); p[top]:=x; clear(lt[x]); clear(rt[x]); x:=0; end; procedure insert(n:longint); var now,t:longint; begin for t:=1 to n do begin if top>0 then begin now:=p[top];bj1[now]:=-inf; bj2[now]:=0; dec(top);end else begin inc(tot);now:=tot;end; b[t]:=now; read(a[now]); end; readln; build(lt[b[n]],1,n-1); s[b[n]]:=n; tmp:=rt[root]; rt[b[n]]:=tmp; updata(b[n]); rt[root]:=b[n]; updata(root); end; procedure predo; begin root:=1;tot:=2; a[1]:=-inf;a[2]:=-inf; sum[1]:=-inf*2;sum[2]:=-inf; max[1]:=-inf;max[2]:=-inf;max[0]:=-inf; s[1]:=2;s[2]:=1; rt[1]:=2; end; procedure main; var i,m,a,b,c:longint; ch:char; procedure readspace; begin while ch<>' ' do read(ch); end; begin readln(n,m); for i:=1 to maxn do begin bj1[i]:=-inf; bj2[i]:=0; end; insert(n); for m:=1 to m do begin read(ch,ch,ch); case ch of 'S': begin readspace; read(a,n); splay(root,a+1); insert(n); end; 'L':begin readspace; readln(a,b); b:=a+b+1; splay(root,b);splay(lt[root],a); clear(rt[lt[root]]); updata(lt[root]); updata(root); end; 'K':begin readspace; readln(a,b,c); b:=a+b+1; splay(root,b);splay(lt[root],a); same(C); end; 'V':begin readspace; readln(a,b); b:=a+b+1; splay(root,b);splay(lt[root],a); rev; end; 'T':begin readspace; readln(a,b); b:=a+b+1; splay(root,b); splay(lt[root],a); writeln(sum[rt[lt[root]]]); end; 'X':begin for i:=1 to 4 do read(Ch); readln; writeln(max[root]); end; end; end; end; begin predo; main; end.