码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj3110: [Zjoi2013]K大数查询 【树套树,标记永久化】

时间:2015-04-01 23:27:59      阅读:216      评论:0      收藏:0      [点我收藏+]

标签:

好久没写题解了。

但是这题太神了然后做法太神了于是写一下。

这题做法很多,比如黄学长hzw的权值线段树套线段树,比如学长云的bit套主席树(其实是写法更神然后我不会用)。

然后看到hzhwcmhf大神题解。

http://tieba.baidu.com/p/2246783535

震惊了。

 

好了开说说做法。建一颗朴素的线段树,树的每个点表示每个区间,然后每个区间建两棵树,一棵是mark树,一棵是all树,两棵都是权值线段树。

“mark表示该区间每个点上都会加上mark线段树里的元素
 all表示该区间所有点的元素集合”

这道题麻烦的地方就在与标记,如果标记要下传的话,那么每次可能要下传很多个点。

于是标记永久化,就是不要下传。

对于修改:如果修改区间是当前点区间,那么就加入mark中。如果不是,那么就加入all。操作都是单点修改。加入mark时为c+1,加入all是为c+len

对于询问:首先包含在查询区间内的区间的所有的all,然后mark要去掉没有重合在一起的地方,乘回差值。

(太久没写题解都不知道怎么说了)

 

速度还是挺快的

技术分享
type
  arr=record
    all,mark:longint;
  end;
const
  maxn=10000000;
var
  tree:array[0..maxn]of arr;
  size,lson,rson,root1,root2:array[0..maxn]of longint;
  tot,sum,n,m,i,j,k,l:longint;
 
function add:longint;
begin
  inc(tot);
  exit(tot);
end;
 
procedure sinsert(var x:longint;y,z:longint);
var
  left,right,mid,old:longint;
begin
  if x=0 then x:=add;
  old:=x;
  left:=1;
  right:=n;
  while left<=right do begin
    mid:=(left+right)>>1;
    inc(size[x],z);
    if left=right then break;
    if y<=mid then begin
      right:=mid;
      if lson[x]=0 then lson[x]:=add;
      x:=lson[x];
    end
    else begin
      left:=mid+1;
      if rson[x]=0 then rson[x]:=add;
      x:=rson[x];
    end;
  end;
  x:=old;
end;
 
procedure binsert(x,l,r,tl,tr,tc:longint);
var
  mid:longint;
begin
  if (tl=l) and (tr=r) then begin
    sinsert(tree[x].mark,tc,1);
    exit;
  end;
  mid:=(l+r)>>1;
  sinsert(tree[x].all,tc,tr-tl+1);
  if tl>mid then binsert(x<<1+1,mid+1,r,tl,tr,tc)
  else
    if tr<=mid then binsert(x<<1,l,mid,tl,tr,tc)
    else begin
      binsert(x<<1,l,mid,tl,mid,tc);
      binsert(x<<1+1,mid+1,r,mid+1,tr,tc);
    end;
end;
 
procedure before(x,l,r,tl,tr:longint);
var
  mid:longint;
begin
  inc(sum);
  root1[sum]:=tree[x].mark;
  root2[sum]:=tr-tl+1;
  if (tl=l) and (tr=r) then begin
    inc(sum);
    root1[sum]:=tree[x].all;
    root2[sum]:=1;
    exit;
  end;
  mid:=(l+r)>>1;
  if tl>mid then before(x<<1+1,mid+1,r,tl,tr)
  else
    if tr<=mid then before(x<<1,l,mid,tl,tr)
    else begin
      before(x<<1,l,mid,tl,mid);
      before(x<<1+1,mid+1,r,mid+1,tr);
    end;
end;
 
procedure query(x,y,z:longint);
var
  left,right,now,mid:longint;
begin
  sum:=0;
  before(1,1,n,x,y);
  left:=1;
  right:=n;
  while left<right do begin
    now:=0;
    mid:=(left+right)>>1;
    for i:=1 to sum do
      now:=size[rson[root1[i]]]*root2[i]+now;
    if now<z then begin
      dec(z,now);
      right:=mid;
      for i:=1 to sum do
        root1[i]:=lson[root1[i]];
    end
    else begin
      left:=mid+1;
      for i:=1 to sum do
        root1[i]:=rson[root1[i]];
    end;
  end;
  writeln(left);
end;
 
begin
  readln(n,m);
  while m>0 do begin
    dec(m);
    readln(i,j,k,l);
    if i=1 then binsert(1,1,n,j,k,l)
      else query(j,k,l);
  end;
  readln;
  readln;
end.
View Code

数组开小竟然wa了两次。

 

bzoj3110: [Zjoi2013]K大数查询 【树套树,标记永久化】

标签:

原文地址:http://www.cnblogs.com/Macaulish/p/4385489.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!