分治法求一个N个元素数组的逆序数

背景

     逆序数:也就是说,对于n个不同的元素,先规定各元素之间有一个标准次序(例如n个 不同的自然数,可规定从小到大为标准次序),于是在这n个元素的任一排列中,当某两个元素的先后次序与标准次序不同时,就说有1个逆序。一个排列中所有逆序总数叫做这个排列的逆序数。

定义

在一个排列中,如果一对数的前后位置与大小顺序相反,即前面的数大于后面的数,那么它们就称为一个逆序。一个排列中逆序的总数就称为这个排列的逆序数。逆序数为偶数的排列称为偶排列;逆序数为奇数的排列称为奇排列。如2431中,21,43,41,31是逆序,逆序数是4,为偶排列。

问题求解

分治法求逆序数相当于在归并排序的过程中加上相应的逆序数的数目。

假设数组A被划分成前半部分(A[left]->A[mid]),后半部分(A[mid+1],A[right])

假设前半部分与后半部分各自都是从小到大排好序的,若A[left]<A[mid+1]那么由A[mid+1]这个数与数组前半部分造成的逆序数目是mid-left+1,这个数在接下来的逆序数中不用再考虑,反之类同。

使用分治法解决该问题的时间复杂度为O(N*log(N)).

代码如下:

 1 #include<iostream>
 2 #include<conio.h>
 3 using namespace std;
 4 int merge(int *A,int left,int mid,int right){
 5     int *temp=new int[right-left+1];
 6     int num=0;
 7     int i=left;
 8     int j=mid+1;
 9     int k;
10     int index=0;
11     for(;i<=mid&&j<=right;){
12         if(A[i]>A[j]){
13             num+=mid-i+1;
14             temp[index]=A[j];
15             j++;
16         }
17         else{
18             temp[index]=A[i];
19             i++;
20         }
22         index++;
23     }
24     //i=mid的时候,A[i]位置的数还未填充到数组temp中
25     //因此判断条件包含等于号
26     if(i<=mid)
27         for(;i<=mid;i++){
28             temp[index]=A[i];
29             index++;
30         }
32     if(j<=right)
33         for(;j<right;j++){
34             temp[index]=A[j];
35             index++;
36         }
3738     for(k=0;k<right-left+1;k++)
39         A[left+k]=temp[k];
40 
41     delete []temp;
42     return num;
43 }
45 int inversion(int *A,int left,int right){
46     if(left>=right)
47         return 0;
48     int mid=(left+right)/2;
49     int num1=inversion(A,left,mid);
50     int num2=inversion(A,mid+1,right);
51     return num1+num2+merge(A,left,mid,right); 
54 }
55 int main()
56 {
57     int A[5]={9,8,7,6,5};
58     cout<<inversion(A,0,4)<<endl;
59     _getch();
60 }

 

 

posted on 2014-06-20 11:19  熊小熊-chris  阅读(3550)  评论(1编辑  收藏  举报

导航