CH Round #45 能量释放

能量释放 CH Round #45 - alan有一些陷阱 III

题目描述

alan得到一块由个能量晶体构成的矿石，对于矿石中的每一个能量晶体，如果用化学物质刺激某一个能量晶体，就能使它释放能量。

它的能量释放强度与晶体本身的能量值以及能量晶体的位置有关。

为了方便研究，alan做了如下的定义。

能量集：一块矿石中的第个能量晶体到第个能量晶体（包含和，）构成的集合。

能量储存点：对于一块矿石中的能量晶体和,若有，则称是能量储存点。

能量释放点：在一个能量集中，若存在一个能量晶体，使得除它之外的所有能量晶体都是它的能量储存点，则称这个能量晶体是该能量集的能量释放点。

能量释放强度：对于一个能量集中的能量释放点来说，刺激这个能量晶体，该能量集中其余能量晶体会释放能量。该能量集的能量释放强度就等于该能量集中其余某个能量晶体的能量值与能量释放点的能量值的差的最大值。而第个能量晶体的能量释放强度则是所在的所有能量集的能量释放强度的最大值。

alan将给出矿石中能量晶体的个数，以及每个能量晶体的能量值，要求你求出每个能量晶体的能量释放强度。

输入格式

第一行1个整数，表示矿石中能量晶体的个数。

第二行个整数，第个整数表示第个能量晶体的能量值。

输出格式

仅包含一行个整数，第个整数表示第个能量晶体的能量释放强度。

样例输入

5
3 2 1 2 3

样例输出

0 1 2 1 0

数据范围与约定

• 对于20%的数据:。
• 对于60%的数据:。
• 对于100%的数据:。
• 温馨提醒：本题数据量较大，请尽量使用输入输出优化。

题解：
在考场上我写的是O（n）的单调队列（其实也相当于是单调栈。。。）+O（nlogn）的rmq
自己为能AC，结果只有一半
后来才知道如果数据规模达到百万以上，就一定得用O（n）的算法了，我还是太年轻啊。。。
今天准备做NOI2005瑰丽的华尔兹时，又用到了单调队列，忽然对这题有了想法，就又来做

这次我的算法，均摊应该是O（n）的吧。。。我想是这样。。。

其实考虑到a[i]左边的元素如果比a[j]小的话，那 j 管辖的范围 i 也也一定可以，(j 的初始值为 i-1）

所以我们直接跳到 j 管辖的范围的左边即可，即 j=l[j]-1 同时用 它来更新f[i]

最后跳不动的时候，j+1就是 l[i]

右边类似。。。

不知到复杂度是多少？怎么估计？求大神指教

自测稍有超时，后几个点都是1.2s左右的样子

代码：

1.考场 跑完所有数据需要 22s

 var  i,j,n,h,t:longint;
      f:array[..+,..] of longint;
      a,q,l,r:array[..+] of longint;
   procedure init;
    begin
      readln(n);
      for i:= to n do read(a[i]);
    end;
  procedure queue;
    begin
      a[]:=-;a[n+]:=a[];
      fillchar(q,sizeof(q),);
      h:=;t:=;
      for i:= to n+ do
       begin
         while (h<t) and (a[i]<a[q[t]]) do
            begin
              r[q[t]]:=i-;
              dec(t);
            end;
         inc(t);q[t]:=i;
       end;
      fillchar(q,sizeof(q),);
      h:=;t:=;
      for i:=n downto  do
        begin
          while (h<t) and (a[i]<a[q[t]]) do
            begin
              l[q[t]]:=i+;
              dec(t);
            end;
          inc(t);q[t]:=i;
        end;
      end;
 function max(x,y:longint):longint;
  begin
    if x>y then exit(x) else exit(y);
  end;

 procedure rmq;
   begin
     fillchar(f,sizeof(f),);
     for i:= to n do f[i,]:=a[i];
     for i:= to  do
      begin
        if (<<i)>n then break;
        for j:= to n-(<<i)+ do
         f[j,i]:=max(f[j,i-],f[j+<<(i-),i-]);
      end;
   end;
 function query(x,y:longint):int64;
  var k:longint;
  begin
   k:=trunc(ln(y-x+1.0)/ln(2.0));
   exit(max(f[x,k],f[y-(<<k)+,k]));
  end;

 procedure main;
   begin
     queue;
     rmq;
   //  for i:= to n do writeln(l[i],' ',r[i]);
     for i:= to n do
       write(query(l[i],r[i])-a[i],' ');
   end;
 begin
   assign(input,'input.txt');assign(output,'output.txt');
   reset(input);rewrite(output);
   init;
   main;
   close(input);close(output);
 end.
                                    

2.20140806 跑完所有数据需要 6s

 {$inline on}
 const maxn=+;
 var l,r,a,f:array[..maxn] of longint;
     i,n,j:longint;
     function max(x,y:longint):longint;inline;
      begin
        if x>y then exit(x) else exit(y);
      end;

 procedure init;
  begin
    readln(n);
    for i:= to n do begin read(a[i]);f[i]:=a[i];end;
  end;
 procedure main;
  begin
    a[]:=-maxlongint;a[n+]:=-maxlongint;
    l[]:=;
    for i:= to n do
     begin
     j:=i-;
     while a[i]<=a[j] do
      begin
      if f[j]>f[i] then f[i]:=f[j];
      j:=l[j]-;
      end;
     l[i]:=j+;
     end;
    r[n]:=n;
    for i:=n- downto  do
     begin
     j:=i+;
     while a[i]<=a[j] do
      begin
      if f[j]>f[i] then f[i]:=f[j];
      j:=r[j]+;
      end;
     r[i]:=j-;
     end;
    for i:= to n- do write(f[i]-a[i],' ');writeln(f[n]-a[n]);
   end;
 begin
   assign(input,'input.txt');assign(output,'output.txt');
   reset(input);rewrite(output);
   init;
   main;
   close(input);close(output);
 end.