在网上看了很多用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中已经插入合并后的数组中的元素个数ki,把这些所有的ki加起来即为本次合并的逆序数,具体的代码如下其中的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;
}