【POI06 TET】二维线段树

终于知道二维线段树怎么写了！

想来也有点好笑，我搞了这么久OI居然连二维线段树都不会写~

以前一直因为这个二维线段树的标记怎么放没搞得清就一直搁着没写，这才终于了却这一“心结”

这个题要求在二维线段树上支持两个操作

1：将一个子矩阵修改成某个值

2：询问某个子矩阵的最大值

整个矩阵中的值都是单调非降的，也就是说可以用一种不用下放标记的线段树写法来处理一维的情况

一维的情况：

线段树上的每个结点记录两个值

cover——被完整覆盖的值

max——子树中cover的最大值（包括自己）

然后询问的时候

如果询问的区间包含该节点的区间，则返回max值，否则用该点的cover更新一下，继续递归

修改的时候对应修改即可

增至二维，与之类似的，第一维没有变化，第二维这个时候就要开两个线段树了，一个cover的线段树，一个max的线段树

询问时，在第一维上做的时候，如果询问的区间包含了该结点，则返回第二维上对应区间的max值

而如果没有包含，就用第二维上的cover值更新一下，继续递归

第二维线段树跟一维的情况完全一样

总的时间复杂度就是O(QlogNlogN)了

ps:我代码中的0域表示cover值，1域表示max值

代码：

{$inline+} program syj; uses math; (* ask2(1, ... ) === max ask2(2, ... ) === label *) const maxd=1005; var p,q,n,a,b,c,x,y:longint; f:array[0..1 shl 11,1..2,0..1 shl 11,0..1]of longint; function ask2(p,u,i,l,r:longint):longint; inline; var m:longint; begin if (y<=l)and(r<=b) then exit(f[p,u,i,1]); ask2:=f[p,u,i,0]; m:=(l+r)>>1; if y<=m then ask2:=max(ask2,ask2(p,u,i*2,l,m)); if b>m then ask2:=max(ask2,ask2(p,u,i*2+1,m+1,r)); end; function ask(i,l,r:longint):longint; inline; var m:longint; begin if (x<=l)and(r<=a) then exit(ask2(i,1,1,1,q)); ask:=ask2(i,2,1,1,q); m:=(l+r)>>1; if x<=m then ask:=max(ask,ask(i*2,l,m)); if a>m then ask:=max(ask,ask(i*2+1,m+1,r)); end; procedure cover2(p,u,i,l,r:longint); inline; var m:longint; begin if (y<=l)and(r<=b) then begin f[p,u,i,0]:=max(f[p,u,i,0],c); f[p,u,i,1]:=max(f[p,u,i,1],c); exit; end; m:=(l+r)>>1; if y<=m then cover2(p,u,i*2,l,m); if b>m then cover2(p,u,i*2+1,m+1,r); f[p,u,i,1]:=max(f[p,u,i,1],c); end; procedure cover(i,l,r:longint); inline; var m:longint; begin cover2(i,1,1,1,q); if (x<=l)and(r<=a) then begin cover2(i,2,1,1,q); exit; end; m:=(l+r)>>1; if x<=m then cover(i*2,l,m); if a>m then cover(i*2+1,m+1,r); end; begin assign(input,'tet.in');reset(input); assign(output,'tet.out');rewrite(output); readln(p,q,n); for n:=1 to n do begin read(a,b,c,x,y); inc(a,x);inc(b,y); inc(x);inc(y); inc(c,ask(1,1,p)); cover(1,1,p); end; writeln(f[1,1,1,1]); close(input);close(output); end.  

 

再说一点，好久没写过blog了，最近搞的东西都懒的发上来，以后会补上的！