归并排序-求逆序数算法

 

归并排序一直以为比较简单，不算特别重视，今天试着写了一下，很尴尬，失败了几次．最后还是决定写博记录一下．废话不多说，直接进正题．

 

算法：归并排序是分治法(分而治之)的一种典型应用，应用递归的思想，自顶向下思考：先假定MergeSort()可以将一个乱序的数组排好序，因此就可以开始"分"(将一个数组平均分成两部分），再"治"(分别对前后部分调用MergeSort()使它们有序)，最后再写一个合并子函数Merge()，它可以将两个有序的数组合并，Merge()实现起来比较容易．只需要管理两个指针，分别指向两个子数组的开头，开辟新内存保存中间结果，遍历完两个数组就可以完成，时间是Θ(n).

 

假定n个元素的数组调用MergeSort()需要时间T(n).因此，T(n)=2T(n/2)+Θ(n)．由主定理可知：T(n)=Θ(nlogn).归并排序算法的时间复杂度是Θ(nlogn)，空间复杂度是Θ(n).

 

     由归并排序扩展可以求逆序数．

 

　　逆序数的定义：如果 i < j 且A[i] > A[j].则A[i]和A[j]即为逆序数对．逆序数对的个数就叫逆序数．因此求逆序数可以通过管理两个指针，两次扫描数组，蛮力法求出，显然时间复杂度是Θ(n^2).那么有更快的办法吗？答案是肯定的，利用归并排序法，稍做改进即可．在Merge()中，合并两个已经有序的数组A,B.因为A.B有序,所以,A,B各自的逆序数是0,所以AB的逆序数等于A,B之间的逆序数.

 

举个例子: A=1,4,6,7,9 B=2,3,5,10,13,21.在Merge中发现当前i号元素4比2大,那么4的逆序数需要+1,又因6,7,9都排在4后面,那么6,7,9的逆序数也应该+1,所以总体的逆序数应该加上last-i+1.(如果i号元素比B[j]小（比如4比5小）,无法确定逆序数的变化,不作任何修改).

代码:

 

#include<iostream>
using namespace std;
int count=0;//新增
void Merge(int* A,int left,int mid,int right,int* C)
{	//Merge可以将两个有序的数组排好序,时间复杂度:o(n)
	int i=left;
	int j=mid+1;
	int k=left;
	while(i <= mid && j <= right)
	{
		if(A[i] <= A[j])
			C[k++]=A[i++];
		else
		{
			C[k++]=A[j++];
			count += mid-i+1;//新增
		}
	}
	while(i <= mid)
		C[k++]=A[i++];
	while(j <= right)
		C[k++]=A[j++];
	//C[]已经有序，将C[]中数据复制回原数组A[]
	for(int i=left;i <= right;++i)
		A[i]=C[i];
}	
void MergeSort(int* A,int left,int right,int* C)//假定MergeSort能将一个乱序数组A排好序．
{
	if(left < right)
	{
		int mid=(left+right)/2;
		MergeSort(A,left,mid,C);//排好一个数组1
		MergeSort(A,mid+1,right,C);//排好一个数组2
		Merge(A,left,mid,right,C); //合并两个有序的数组
	}
}

void main()
{
	int A[]={2,1,3,6,4,0,11,3,5};
	int len=sizeof(A)/sizeof(A[0]);
	int *C=new int[len];
	MergeSort(A,0,len-1,C);
	for(int i=0;i<len;++i)
		cout<<A[i]<<' ';
	cout<<endl;
	cout<<count<<endl;
	delete[] C;
}