Day2:T4求逆序对(树状数组+归并排序)

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即可

很好理解就不多说了...

晚安民那

原文地址:https://www.cnblogs.com/polebug/p/3961907.html