逆序数算法

在网上看了很多用mergesort求逆序数的代码，发觉很多都不太对，主要是在merge时求本次merge是前后两部分的逆序数过程，其实merge时由于前后两部分都是已经排好序的， 假设merge的两个数组分别为a[n1],b[n2],合并时的逆序数就应该=b中所有小于a[1]的个数，小于a[2]的个数，小于a[3]的个数...小于a[n1]的个数之和=a中所有大于b[1]的个数，大于b[2]的个数，大于b[3]的个数...大于b[n2-1]的个数之和。因此对于此次合并逆序数的个数有两种计算方法，一是以a中元素作为判断来计算，一是以b中元素作为判断来计算。以a中元素作为为例：计算过程为，每当把a中某一个元素a[i]插入合并后的数组中时，计算b中已经插入合并后的数组中的元素个数k_i，把这些所有的k_i加起来即为本次合并的逆序数，具体的代码如下其中的iversion为逆序数。

 int inversion=0;
void merge(int* a,int p,int q,int r){
 int n1=r-p+1;
 int n2=q-r;
 int i,j,k;
 k=0;
 int *tmp=new int[n1+n2];
 for(i=p,j=r+1;i<=r&&j<=q;){
  if(a[i]<=a[j]){
   tmp[k]=a[i];
   inversion+=j-r;
   i++;
   k++;
  }else{
   tmp[k]=a[j];
   j++;
   k++;
  }
 }
 while(i<=r){
  tmp[k++]=a[i++];
  inversion+=q-r;
 }
 while(j<=q){
  tmp[k++]=a[j++];
 }
 for(i=p;i<=q;i++)a[i]=tmp[i-p];
 
 delete [] tmp;
}