T4: 求逆序对
A[I]为前缀和
推导 (A[J]-A[I])/(J-I)>=M
A[j]-A[I]>=M(J-I)
A[J]-M*J>=A[I]-M*I
设B[]=A[]-M*();
B[J]>=B[I]
也就是求逆序对;
求逆序对的方法主要有两种:
归并排序;
树状数组;
这里两种方法都学习一下:
1.之前对于树状数组的印象就只有单点修改和区间求和
一直觉得lowbit是一个神奇的东西(至今没有搞懂原理)
上网搜了一下用树状数组求逆序对的方法,发现有一个大神写的很棒....看得很明白(感谢~~)
顺便学了一下离散化,也就是文中的:
for
(i=1;i<=n;i++)
{
scanf
(
"%d"
,&in[i].v);
in[i].order=i;
}
sort(in+1,in+n+1,cmp);
for
(i=1;i<=n;i++) aa[in[i].order]=i;
有点像指针操作。(自己对离散化的初步感受是这样的)
然后就是根据树状数组的原理求逆序对了
其实一直不是很清楚树状数组的原理,翻了翻书才有一点点的领悟吧
1).首先是update中的c数组的理解,上面链接的那个资料中解释的过程会误以为c是统计某个数的个数的
仔细看白书会发现其实并不是这样,c数组是用来维护某一区间的区间和的
getsum的操作也就是把沿途中的长条所表示的区间和加起来,这样可以节省时间(注意理解图)
2).具体的,标程中的c是用来记录在某个区间里有多少个数比i小
3).理解lowbit的过程,同样需要仔细回看白书的解释
无论是i=t-n;i+=lowbit(i)或是i-=lowbit(i)都是一个找父亲的过程
4).其实说白了lowbit到底有什么用呢,其实也就是节省了时间,做到快速简便而已吧
嗯...这样的话,对于树状数组的原理算是搞懂了QAQ感觉这种东西很容易忘记怎么办?
给以后的自己一个忠告呗~先看白书,好好理解图,再看自己的感悟,最后看链接资料;
2.另一种方法就是很传统的用:归并排序求逆序对
这种算法求逆序对也很容易理解,简单的说也无非就是基于快排思想+统计个数
主要的程序:p版的...不要在意这些细节...
procedure mergesort(l,r:longint);
var
mid,i,j,k:longint;
begin
if l=r then exit;
mid:=(l+r) div 2;
mergesort(l,mid);
mergesort(mid+1,r);
i:=l;j:=mid+1;k:=l;
while (i<=mid) and (j<=r) do
if v[i]<v[j] then
begin
push(k,i);
ans:=ans+r-j+1;//在l..mid 和 mid+1 ..r中j处于[mid+1,r]中且这段是升序的,那么如果v[j]>v[i]那么v[j..r]都会大于//v[i]这和求逆序对有异曲同工之妙
end
else
push(k,j);
while i<=mid do push(k,i);
while j<=r do push(k,j);
for i:=l to r do v[i]:=temp[i];
end;
根据推导建立等式,之后输出ans即可
很好理解就不多说了...
晚安民那