poj 3368 Frequent values 线段树

题目大意

给一个长度为n的不降序列a1,a2,a3,…,an,有q个询问，每个询问为:i j

询问在子序列ai…aj中出现最多的元素。

数据范围：1 <= n, q <= 100000

 

分析

  注意到题目描述中的“不降序列”，让我们联想到可以使用线段树这一数据结构。

  对于线段树的每一个节点：

  记录left[k]表示当前区间最左边连续的元素有多少个

  记录right[k]表示当前区间最右边连续的元素有多少个

  记录max[k]表示当前区间中出现最多的元素有多少（当然它们是连续的）

  当然，还要储存这些区间是什么数。

  对于每一次询问，则可用类似方法，得到区间[i,j]中最左边连续元素的个数，最右边连续元素的个数，以及最多的连续元素的个数。

自己见代码

代码

  

type
  pnode=^tnode;
  tnode=record
    lc,rc:pnode;
    a,b,c:longint;
    m1,m2,m3:longint;
end;

type
  arr=record
    a,b,c:longint;
    m1,m2,m3:longint;
end;

var
  t:pnode;
  i,j,k:longint;
  x,y:longint;
  n,m:longint;
  ans:arr;
  f:array[0..100002] of longint;

procedure neww(var t:pnode);
begin
  if t=nil then
    begin
      new(t);
      t^.c:=0;
      t^.lc:=nil;
      t^.rc:=nil;
    end;
end;

function max(x,y:longint):longint;
begin
  if x>y then exit(x)
         else exit(y);
end;

procedure insert(var t:pnode; l,r:longint);
var
  i,j,k:longint;
  mid:longint;

begin
  with t^ do
    begin
      {if c=0 then
        begin   }
          mid:=(l+r) div 2;
          if l=r
            then
              begin
                a:=1; b:=1; c:=1;
                m1:=f[l]; m2:=f[l]; m3:=f[l];
                exit;
              end;
          neww(lc);
          neww(rc);
          insert(lc,l,mid);
          insert(rc,mid+1,r);
          c:=max(lc^.c,rc^.c);
          if c=lc^.c then m3:=lc^.m3
                     else m3:=rc^.m3;
          if lc^.m2=rc^.m1
            then
              begin
                c:=max(c,lc^.b+rc^.a);
                m3:=lc^.m2;
              end;
          a:=lc^.a;
          m1:=lc^.m1;
          if lc^.m1=rc^.m1
            then a:=a+rc^.a;
          b:=rc^.b;
          m2:=rc^.m2;
          if lc^.m2=rc^.m2
            then b:=b+lc^.b;
    end;
    {end;}
end;

function find(t:pnode;l,r,x,y:longint):arr;
var
  mid:longint;
  l1,r1:arr;
begin
  if t=nil then exit;
  with t^ do
    begin
      mid:=(l+r) div 2;
      if (l=x) and (r=y)
        then
          begin
            find.a:=a; find.b:=b; find.c:=c;
            find.m1:=m1; find.m2:=m2; find.m3:=m3;
            exit;
          end;
      if (l<=x) and (mid>=y)
        then
          begin
            find:=find(lc,l,mid,x,y);
            exit;
          end;
      if (mid<x) and (r>=y)
        then
          begin
            find:=find(rc,mid+1,r,x,y);
            exit;
          end;
      l1:=find(lc,l,mid,x,mid);
      r1:=find(rc,mid+1,r,mid+1,y);
      with find do begin
        c:=max(l1.c,r1.c);
        if c=l1.c then m3:=l1.m3
                   else m3:=r1.m3;
        if l1.m2=r1.m1
          then
            begin
              c:=max(c,l1.b+r1.a);
              m3:=l1.m2;
            end;
        a:=l1.a;
        m1:=l1.m1;
        if l1.m1=r1.m1
          then a:=a+r1.a;
        b:=r1.b;
        m2:=r1.m2;
        if l1.m2=r1.m2
          then b:=b+l1.b;
      end;
    end;
end;

begin
  while true do
  begin
  read(m);
  if m=0 then break;
  fillchar(f,sizeof(f),0);
  dispose(t);
  readln(n);
  for i:=1 to m do
    read(f[i]);
  readln;
  fillchar(t,sizeof(t),0);
  neww(t);
  insert(t,1,m);
  for i:=1 to n do
    begin
      readln(x,y);
      ans:=find(t,1,m,x,y);
      writeln(ans.c);
    end;
  end;
end.