poi2007--drz。

好开心，一次就ac了。。（其实。。这就是一道10k左右的线段树水题。。）

 

{ps：很囧。。。我想了差不多两小时都没想出来。最后看了一下题解- -|ca。。。波兰文的- -|显然我不是

语言帝我是看不懂的- -|中文题解又有bug- -然后看了一下那些式子- -差不多懂了- -|。。
好吧。。跟上次那道物理题一样。。既然是计算机竞赛题目- -就不要像我一样太往数学方面想，。一开始一直想如何贪心活着数学方法什么的。
}

 

 

看了题解后顿悟这题很水。。竟然没自己想出来，分析没有对路啊。。应该是我没有往那么多种情况考虑- -|

以后再遇到这个类型的题目- -有情况就拉进来分析。。然后大胆的分类讨论。。

 

这种题目就是找出他的表达式，然后进行化简并且分情况讨论。。
当这种求区间最值又有一定大小限制的就是先离散化，然后根据要求按顺序插入。。。
而且当时我一直在想简单的方法。。结果分9种情况。、、啊。。。拿来练手也是不错的- -
我们设vi=|hi+1-hi|+|hi-1 - hi|

然后对调一次i，j的代价就是  |hj+1-hi|+|hj-1 - hi|+|hi+1-hj|+|hi-1 - hj| -vi - vj

然后设

max[i]=max(h[i-1],h[i+1])

min[i]=min(h[i-1],h[i+1])

 

这个就变成了 |maxj-hi|+|minj - hi|+|maxi-hj|+|mini - hj| -vi - vj

 

这些是做这题做的草稿- -|好吧。。我太弱了。。。为了避免出错检查太久，我得一个一个列出来才会

- -然后列了差不多7个后才发现了点规律 。。

当然以上只适用于abs(i-j)>1的和(i,j <>1 or i,j<>n)

所以这些先朴素弄出来。

先拿这个情况分析一下
1.h[i]<min[ j]    h[j]<min[i]

   delta = min_i+max_i+min_j+max_j-2*h_i-2*h_j-v_i-v_j
key[j] = min_j + max_j - 2*h_j - v_j
add[i] = key[i]

我们把 min[i], h[i] ,max[i] 都离散化。从大到小排序。。然后维护一个h值的线段树

按照从大到小插。
遇到min_i 在插入h[i]的位置插入key[i]
遇到max_i 无视
遇到h_i 利用线段树找出 从初始到min_i的最小key值

算了 下划线太不顺手了。。还是用中括号来= =

2. min_j<h_i<max_j    h_j<min_i
delta = max[j]-min[j]+max[i]+min[i]-2*h[j]-v[i]-v[j]
key[j] = max_j - min_j - 2*h_j - v_j
add[i] = min_i + max_i - v_i
按从小到大插
遇到min_i 插入 key[i] 到 h_i
遇到max_i 删除 key[i]
遇到h_i 利用线段树找出 从初始到min_i的最小key值

3.h[i]>max[j]  h[j]<min[i]
   delta = 2*h[i]-max[j]-min[j]+max[i]+min[i]-2*h[j]-v[i]-v[j]
   key[j] = -max[j]-min[j]-2*h[j]-v[j]
   add[i] = 2*h[i]+min[i]+max[i]-v[i];
遇到max[i] 插入 h[i]
min无视
遇到hi利用线段树找出max[i]到终点的最小key值。。。


前三种情况是作为代表性的三种- -|

（1.4.7 ）（2,5,8） （3,6,9)  每组内都是类似的。。 线段树。。同理。。

4.h[i]<min[j]   min[i]<h[j]<max[i]
  delta = min[j]+max[j]-2*h[i] + max[i]-min[i]-v[i]-v[j]
  key[j] = min[j] + max[j] -v[j];
  add[i] = -2*h[i]+max[i]-min[i] -v[i];

类似1

5,min[j]<h[i]<max[j]  min[i]<h[j]<max[i]
  delta = max[i]-min[i]+max[j]-min[j]-v[i]-v[j]
key[j] =max[j]-min[j]-v[j];
add[i] = key[i]
类似二。

6。h[i]>max[j]  min[i]<h[j]<max[i]

7 h[i]<min[j]   h[j]>max[i]
  delta = min[j] + max[j] -2 * h[i]  +2*h[j] - max[i] -min[i]-v[i]-v[j]
key[j] = 2*h[j] +max[j]+min[j] -v[j]
add[i] = -2*h[i]-max[i]-min[i] -v[i]

8 min[j]<h[i]<max[j]  h[j]>max[i]
deta = max[j]-min[j] +2*h[j]-max[i]-min[i] -v[i]-v[j]
key[j] = max[j]-min[j]+2*h[j]-v[j]
add[i] = -max[i]-min[i]-v[i]

9.h[i]>max[j]  h[j]>max[i]

 

program Project1; const inf = 1000000000000000; maxn = 50000; var datal,datar,lsame,rsame,lt,rt,pt,min,max,pos,h,data,v,flag,num :array[0..maxn*4]of int64; use :array[0..maxn]of boolean; ans :array[0..maxn]of int64; p :array[0..maxn*3*4]of int64; n,tot :longint; res :int64; function small(x,y:int64):int64; begin if x<y then exit(x) else exit(y); end; function biger(x,y:int64):int64; begin if x>y then exit(x) else exit(y); end; procedure updata(var x:int64;y:int64); begin if y<x then x:=y; end; procedure swap(var x,y:int64); var tmp:int64; begin tmp:=x;x:=y;y:=tmp; end; procedure init; var i:longint; begin readln(n); for i:=1 to n do read(h[i]); end; procedure qsort(l,r:longint); var mif,i,j:longint; mid:int64; begin i:=l;j:=r; mid:=data[(i+j) shr 1]; mif:=flag[(i+j) shr 1]; repeat while (data[i]<mid)or((data[i]=mid)and(flag[i]<mif)) do inc(i); while (mid<data[j])or((data[j]=mid)and(flag[j]>mif)) do dec(j); if i<=j then begin swap(data[i],data[j]); swap(flag[i],flag[j]); swap(num[i],num[j]); swap(pos[i],pos[j]); inc(i);dec(j); end; until i>j; if l<j then qsort(l,j); if i<r then qsort(i,r); end; procedure predo; var i,now :longint; tmp :int64; begin for i:=1 to n do ans[i]:=inf; res:=0; for i:=1 to n-1 do res:=res+abs(h[i]-h[i+1]); for i:=2 to n-1 do v[i]:=abs(h[i]-h[i-1])+abs(h[i+1]-h[i]); for i:=3 to n-1 do begin tmp:=abs(h[2]-h[i])+abs(h[1]-h[i-1])+abs(h[1]-h[i+1])-abs(h[1]-h[2])-v[i]; updata(ans[1],tmp); updata(ans[i],tmp); end; for i:=2 to n-2 do begin tmp:=abs(h[n-1]-h[i])+abs(h[n]-h[i-1])+abs(h[n]-h[i+1])-abs(h[n]-h[n-1])-v[i]; updata(ans[n],tmp); updata(ans[i],tmp); end; tmp:=abs(h[n-1]-h[1])+abs(h[n]-h[2])-abs(h[n-1]-h[n])-abs(h[1]-h[2]); updata(ans[1],tmp); updata(ans[n],tmp); if n>2 then begin tmp:=abs(h[1]-h[3])-abs(h[2]-h[3]); updata(ans[1],tmp); updata(ans[2],tmp); tmp:=abs(h[n]-h[n-2])-abs(h[n-1]-h[n-2]); updata(ans[n],tmp); updata(ans[n-1],tmp); end; for i:=2 to n-2 do begin tmp:=abs(h[i]-h[i+2])+abs(h[i+1]-h[i-1])+abs(h[i]-h[i+1])-v[i]-abs(h[i+2]-h[i+1]); updata(ans[i],tmp); updata(ans[i+1],tmp); end; for i:=2 to n-1 do begin min[i]:=small(h[i-1],h[i+1]); max[i]:=biger(h[i-1],h[i+1]); inc(tot); flag[tot]:=-1; data[tot]:=min[i]; num[tot]:=i; lt[i]:=tot; inc(Tot); flag[tot]:=0; data[tot]:=h[i]; num[tot]:=i; inc(tot); flag[tot]:=1; data[tot]:=max[i]; num[tot]:=i; rt[i]:=tot; end; for i:=1 to tot do pos[i]:=i; qsort(1,tot); data[0]:=inf; now:=0; lsame[1]:=1; rsame[tot]:=tot; for i:=2 to tot do if data[i]=data[i-1] then lsame[i]:=lsame[i-1] else lsame[i]:=i; for i:=tot-1 downto 1 do if data[i]=data[i+1] then rsame[i]:=rsame[i+1] else rsame[i]:=i; for i:=1 to tot do begin data[pos[i]]:=i; datal[pos[i]]:=lsame[i]; datar[pos[i]]:=rsame[i]; if flag[i]=0 then pt[num[i]]:=i; end; for i:=2 to n-1 do begin lt[i]:=datal[lt[i]]; rt[i]:=datar[rt[i]]; end; end; procedure change(tot,l,r,now:longint;x:int64); begin if l=r then begin p[tot]:=x; exit; end; if now<=(l+r) shr 1 then change(tot*2,l,(l+r) shr 1,now,x) else change(tot*2+1,(l+r) shr 1+1,r,now,x); p[tot]:=small(p[tot*2],p[tot*2+1]); end; function findmin(tot,l,r,ll,rr:longint):int64; var t1,t2:int64; begin if (ll<=l)and(r<=rr) then exit(p[tot]); t1:=inf;t2:=inf; if ll<=(l+r) shr 1 then t1:=findmin(tot*2,l,(l+r) shr 1,ll,rr); if rr>(l+r) shr 1 then t2:=findmin(tot*2+1,(l+r) shr 1+1,r,ll,rr); exit(small(T1,t2)); end; procedure insert(x,k:longint); var tmp:int64; begin case k of 0:tmp:=inf; 1:tmp:=max[x]+min[x]-2*h[x]-v[x]; 4:tmp:=max[x]+min[x]-v[x]; 7:tmp:=max[x]+min[x]+2*h[x]-v[x]; 2:tmp:=max[x]-min[x]-2*h[x]-v[x]; 5:tmp:=max[x]-min[x]-v[x]; 8:tmp:=max[x]-min[x]+2*h[x]-v[x]; 3:tmp:=-max[x]-min[x]-2*h[x]-v[x]; 6:tmp:=-max[x]-min[x]-v[x]; 9:tmp:=-max[x]-min[x]+2*h[x]-v[x]; end; change(1,1,tot,pt[x],tmp); end; procedure work1(k:longint); var i,x:longint; t1,t2:boolean; add,tmp:int64; begin for i:=1 to tot*4 do p[i]:=inf; for i:=1 to n do use[i]:=false; for i:=tot downto 1 do begin x:=num[i]; if flag[i]=0 then begin t1:=false; t2:=false; if use[x-1] then begin t1:=true; insert(x-1,0); end; if use[x+1] then begin t2:=true; insert(x+1,0); end; if k=1 then add:=min[x]+max[x]-2*h[x]-v[x] else if k=2 then add:=-2*h[x]+max[x]-min[x]-v[x] else add:=-2*h[x]-max[x]-min[x]-v[x]; if k=1 then tmp:=findmin(1,1,tot,1,lt[x]) else if k=2 then tmp:=findmin(1,1,tot,lt[x],rt[x]) else tmp:=findmin(1,1,tot,rt[x],tot); updata(ans[x],tmp+add); if t1 then insert(x-1,3*k-2); if t2 then insert(x+1,3*k-2); end; if flag[i]=-1 then begin insert(x,3*k-2); use[x]:=true; end; end; end; procedure work2(k:longint); var i,x:longint; t1,t2:boolean; add,tmp:int64; begin for i:=1 to tot*4 do p[i]:=inf; for i:=1 to n do use[i]:=false; for i:=1 to tot do begin x:=num[i]; if flag[i]=0 then begin t1:=false; t2:=false; if use[x-1] then begin t1:=true; insert(x-1,0); end; if use[x+1] then begin t2:=true; insert(x+1,0); end; if k=1 then add:=max[x]+min[x]-v[x] else if k=2 then add:=max[x]-min[x]-v[x] else add:=-max[x]-min[x]-v[x]; if k=1 then tmp:=findmin(1,1,tot,1,lt[x]) else if k=2 then tmp:=findmin(1,1,tot,lt[x],rt[x]) else tmp:=findmin(1,1,tot,rt[x],tot); updata(ans[x],tmp+add); if t1 then insert(x-1,3*k-1); if t2 then insert(x+1,3*k-1); end; if flag[i]=-1 then begin insert(x,3*k-1); use[x]:=true; end; if flag[i]=1 then begin insert(x,0); use[x]:=false; end; end; end; procedure work3(k:longint); var i,x:longint; t1,t2:boolean; tmp,add:int64; begin for i:=1 to tot*4 do p[i]:=inf; for i:=1 to n do use[i]:=false; for i:=1 to tot do begin x:=num[i]; if flag[i]=0 then begin t1:=false; t2:=false; if use[x-1] then begin t1:=true; insert(x-1,0); end; if use[x+1] then begin t2:=true; insert(x+1,0); end; if k=1 then add:=max[x]+min[x]+2*h[x]-v[x] else if k=2 then add:=max[x]-min[x]+2*h[x]-v[x] else add:=-max[x]-min[x]+2*h[x]-v[x]; if k=1 then tmp:=findmin(1,1,tot,1,lt[x]) else if k=2 then tmp:=findmin(1,1,tot,lt[x],rt[x]) else tmp:=findmin(1,1,tot,rt[x],tot); updata(ans[x],tmp+add); if t1 then insert(x-1,3*k); if t2 then insert(x+1,3*k); end; if flag[i]=1 then begin insert(x,3*k); use[x]:=true; end; end; end; procedure solve; var i:longint; begin for i:=1 to 3 do begin work1(i); work2(i); work3(i); end; for i:=1 to n do if ans[i]>0 then writeln(res) else writeln(res+ans[i]); end; begin assign(input,'1.in');assign(output,'1.out'); reset(input);rewrite(output); init; predo; solve; close(input);close(output); end.