bzoj2653

CLJ神牛的可持久化论文的题目,果然厉害
其实第一步能想到后面就还是很简单的
首先是二分答案,转化为判定性问题
然后对于区间内的数,比他大的标为1,小的标为-1
显然,如果存在一个左右端点符合的区间使得这个区间和大于等于0(因为这里中位数是向下取整)
那么中位数一定是大于等于这个数的,这并不难理解
下面我们就是要快速求出左端点在[a,b],右端点[c,d]的最大区间和
显然可以转化为rmax[a,b]+sum[b+1,c-1]+lmax[c,d];
这里我们不难想到对于每个数建立一个线段树
线段树就是以位置为索引的,记录着区间左最大,区间右最大和区间和
考虑到树的形态都是相同的,我们可以先对数进行排序,然后建立主席树
平时我们的主席树都是在每个位置上建立一棵以键值为索引的线段树,而这里刚好相反
二分答案后就对对应数的线段树求区间最大和

 type node=record
       lm,rm,l,r,s:longint;
     end; var tree:array[..*] of node;
    a,c,h:array[..] of longint;
    q:array[..] of longint;
    j,m,t,ans,n,mid,i,l,r,w:longint; function max(a,b:longint):longint;
  begin
    if a>b then exit(a) else exit(b);
  end; procedure swap(var a,b:longint);
  var c:longint;
  begin
    c:=a;
    a:=b;
    b:=c;
  end; procedure deal;
  var i,j:longint;
  begin
    for i:= to do
      for j:=i+ to do
        if q[i]>q[j] then swap(q[i],q[j]);
  end; procedure sort(l,r: longint);
  var i,j,x,y: longint;
  begin
    i:=l;
    j:=r;
    x:=a[(l+r) shr ];
    repeat
      while a[i]<x do inc(i);
      while x<a[j] do dec(j);
      if not(i>j) then
      begin
        swap(c[i],c[j]);
        swap(a[i],a[j]);
        inc(i);
        j:=j-;
      end;
    until i>j;
    if l<j then sort(l,j);
    if i<r then sort(i,r);
  end;
procedure update(var a,b,c:node);
  begin
    a.s:=b.s+c.s;
    a.lm:=max(b.lm,b.s+c.lm);
    a.rm:=max(c.rm,c.s+b.rm);
  end; function build(l,r:longint):longint;
  var m,q:longint;
  begin
    inc(t);
    if l=r then
    begin
      tree[t].s:=;
      tree[t].lm:=;
      tree[t].rm:=;
      exit(t);
    end
    else begin
      m:=(l+r) shr ;
      q:=t;
      tree[q].l:=build(l,m);
      tree[q].r:=build(m+,r);
      update(tree[q],tree[tree[q].l],tree[tree[q].r]);
      exit(q);
    end;
  end; function work(l,r,last,x:longint):longint;
  var m,q:longint;
  begin
    inc(t);
    if l=r then
    begin
      tree[t].s:=-;
      tree[t].lm:=-;
      tree[t].rm:=-;
      exit(t);
    end
    else begin
      m:=(l+r) shr ;
      q:=t;
      if x<=m then
      begin
        tree[q].r:=tree[last].r;
        tree[q].l:=work(l,m,tree[last].l,x);
      end
      else begin
        tree[q].l:=tree[last].l;
        tree[q].r:=work(m+,r,tree[last].r,x);
      end;
      update(tree[q],tree[tree[q].l],tree[tree[q].r]);
      exit(q);
    end;
  end; function ask(l,r,x,ql,qr:longint):node;  //传递整个node的的写法比我以前的写法要简单不少
  var m:longint;
      s,s1,s2:node;
  begin
    s.lm:=;  s.rm:=;    s.s:=;
    if ql>qr then exit(s);
    if (ql<=l) and (qr>=r) then exit(tree[x])
    else begin
      m:=(l+r) shr ;
      if ql>m then exit(ask(m+,r,tree[x].r,ql,qr));
      if qr<=m then exit(ask(l,m,tree[x].l,ql,qr));
      s1:=ask(l,m,tree[x].l,ql,qr);
      s2:=ask(m+,r,tree[x].r,ql,qr);
      update(s,s1,s2);
      exit(s);
    end;
  end; begin
  readln(n);
  for i:= to n do
  begin
    readln(a[i]);
    c[i]:=i;
  end;
  sort(,n);
  h[]:=build(,n);
  for i:= to n do
    h[i]:=work(,n,h[i-],c[i-]);
  readln(m);
  for i:= to m do
  begin
    for j:= to do
    begin
      read(q[j]);
      q[j]:=(q[j]+ans) mod n+;
    end;
    readln;
    deal;
    l:=;
    r:=n;
    while l<=r do  //二分答案
    begin
      mid:=(l+r) shr ;
      w:=ask(,n,h[mid],q[],q[]).rm+ask(,n,h[mid],q[],q[]).lm+ask(,n,h[mid],q[]+,q[]-).s;
      if w>= then
      begin
        ans:=a[mid];
        l:=mid+;
      end
      else r:=mid-;
    end;
    writeln(ans);
  end;
end.
上一篇:BZOJ.1443.[JSOI2009]游戏Game(二分图博弈 匈牙利)


下一篇:【CF860E】Arkady and a Nobody-men 长链剖分