归并排序

废话不说，直接上代码。

递归版本：

#include<iostream>
template<typename T>
void swap(T *tmp1,T *tmp2);
template<typename T>
void merge(T *tmp,int left,int right){
	int i=left,mid=(left+right)/2,k=0;
	int j = mid+1;
	T *s;
	s = (T *)malloc((right-left+1)*sizeof(T));
	while(i<=mid&&j<=right){
		if(tmp[i]>tmp[j]){
			swap(tmp[i],tmp[j]);
			++i;
		}
		else{
			s[k] = tmp[j];
			++j;
		}
		++k;
	}
	while(i<=mid){
		s[k] = tmp[i];
		++k;++i;
	}
	while(j<=right){
		s[k] = tmp[j];
		++k;++j;
	}
	for(int i=left,k=0;i<=right;++k,++i)tmp[i] = s[k];
	free(s);
}
template<typename T>
void swap(T *tmp1,T *tmp2){
	if(tmp1!=tmp2){
		*tmp1 = *tmp1 + *tmp2;
		*tmp2 = *tmp1 - *tmp2;
		*tmp1 = *tmp1 - *tmp2;
	}
}
void swap()
template<typename T>
void mergesort(T *tmp,int left,int right){
	if(left<right){
		int mid;
		mid = (left+right)/2;
		mergesort(tmp,left,mid);
		mergesort(tmp,mid+1,right);
		merge(tmp,left,right);
	}
}
int main()
{
	using namespace std;
	int A[]= {8,3,2,9,7,1,5,4};
	int size = sizeof(A)/sizeof(*A);
	mergesort(A,0,size-1);
	for(int i=0;i<size;++i)
		cout<<A[i]<<" ";
	cout<<endl;
}

关于递归版本的归并时间复杂度计算：
T(n)=2T(n/2)+n =2(2(T(/4)+N/2))+Cmerge(n)
n=1 Cmerge(n) = 0;
n=2 Cmerge(n) = 1;
............
n Cmerge(n) = n-1
取 n = 2^k，因为n是否去2的幂次方的数对最后的时间复杂度计算没有多大影响，这个可以画2叉树图证明。
显然 T(n) = 2^k+n(log n) n= 2^k; T(n) = n(log n)+n 得出复杂度为O(n logn);
关于最好情况下和平均情况时间复杂度，这里没有说明，感兴趣的朋友私下讨论。
对于将整个数组每次按3份进行划分的情况，相同n的时候k确实会小一些，但是这样会增加归并时候的复杂度，甚至在长度n超过某个数的时候后面的C merge(n)时间复杂度会大过前面。印象中这个界限长度n不是一个大数。
非递归版本： 这里就不做时间复杂度分析。

#include<iostream>
#include<cmath>
int upfoor(int size,int m);
template<typename T>
void mergesort(T *A,int size){
int m=2;
while(m<=size){
int tim =upfoor(size,m);//取上整
for(int i=0;i<tim;++i){
int *tmp = (T *)malloc(m*sizeof(T));
int left=m*i,right=m*(i+1)-1;
if(right>size)right = size;
int mid = (left+right)/2;
int t = left,j=mid+1,k=0;
while(t<=mid&&j<=right){
if(A[t]<=A[j]){
tmp[k] = A[t];
++t;
}
else{
tmp[k] = A[j];
++j;
}
++k;
}
while(t<=mid){
tmp[k] = A[t];
++k;++t;
}
while(j<=right){
tmp[k] = A[j];
++k;++j;
}
for(int i=left,k=0;i<=right;++i,++k)
A[i] = tmp[k];
}
m = m*2;
}
}
int upfoor(int size,int m){
int tmp = size/m;
if(tmp*m==size)return tmp;
else
return tmp+1;
}
int main()
{
using namespace std;
int A[]= {8,3,2,9,7,1,5,4};
int size = sizeof(A)/sizeof(*A);
mergesort(A,size);
for(int i=0;i<size;++i)
cout<<A[i]<<" ";
cout<<endl;
}