归并排序算法完全遵循分治模式。直观上操作如下;
- 分解: 分解待排序的n个元素的序列成各具n/2个元素的两个子序列。
- 解决: 使用归并排序递归的排序两个子序列
- 合并: 合并两个已排序的子序列来产生已排序的答案
当待排序的序列长度为1,递归开始回升,不需要做处理,长度为1的序列都是排好序的。
归并排序算法的关键点在于“合并”步骤中两个已排序序列的合并,通过调用一个辅助的函数来合并。
merge(A,p,q,r)
A是一个数组,p,q,r是数组的下标。该过程假设子数组A(p,,,q)和A(q+1,,,,r)都是排好序的,他们合并起来形成单一的已经排序好的子数组代替当前的子数组A(p,,,,r)
过程merge需要耗费n级时间,n是待排序数的个数,n=r-p+1。
算法过程详解
假设桌面有两堆排序好的牌,两堆的顶部都有着自己堆中最小的牌。首先比较两堆中各自最小的牌,选取最小的牌拿出(拿出牌的那一堆将显露一张新牌)并牌面朝下的放置到输出堆中,重复这个步骤,直到一个输出堆为空,这时只剩下一堆没拿的牌,但是也是排好序的,所以依次放到输出堆里。
这里有一个特殊的步骤,以避免在每个基本步骤必须检查是否有堆为空,在每一个堆得底部放一个哨兵牌,用于简化代码,这里用正无穷代替。结果每当显露到哨兵牌时,它不可能是较小的牌,除非两堆都到达哨兵牌。但是,一旦发生这种情况,所有的哨兵牌都已经放置到输出堆去了,因为事先我们知道要有刚好r-p+1张牌将会放到输出堆里去,所以执行完r-p+1步骤后,算法停止了。
合并算法在此
void Merge(int a[100],int p,int q,int r)
{
int i,j,k;
int n1=q-p+1;
int n2=r-q;
int l[100];
int R[100];
for(i=1;i<=n1;i++)
{
l[i]=a[p+i-1];
}
for(j=1;j<=n2;j++)
{
R[j]=a[q+j];
}
l[n1+1]=MAX;
R[n2+1]=MAX;
i=1;
j=1;
for(k=p;k<=r;k++)
{
if(l[i]<=R[j])
{
a[k]=l[i];
i++;
}
else
{
a[k]=R[j];
j++;
}
}
}
现在把merge作为归并算法的一个子程序来用,下面的函数用来排序子数组A【p,,r】中元素,如果p>=r,则子数组最多只有一个元素,所以是排好序的,否则简单计算一个下标q,将数组分解为两个子数组。具体如下
void Mergesort(int a[100],int p,int r)
{
if(p<r)
{
int q=(p+r)/2;
Mergesort(a,p,q);
Mergesort(a,q+1,r);
Merge(a,p,q,r);
}
}
完整代码如下
#define MAX 24125541
#include<stdio.h>
int a[100];
void Merge(int a[100],int p,int q,int r)
{
int i,j,k;
int n1=q-p+1;
int n2=r-q;
int l[100];
int R[100];
for(i=1;i<=n1;i++)
{
l[i]=a[p+i-1];
}
for(j=1;j<=n2;j++)
{
R[j]=a[q+j];
}
l[n1+1]=MAX;
R[n2+1]=MAX;
i=1;
j=1;
for(k=p;k<=r;k++)
{
if(l[i]<=R[j])
{
a[k]=l[i];
i++;
}
else
{
a[k]=R[j];
j++;
}
}
}
void Mergesort(int a[100],int p,int r)
{
if(p<r)
{
int q=(p+r)/2;
Mergesort(a,p,q);
Mergesort(a,q+1,r);
Merge(a,p,q,r);
}
}
int main()
{
int i;
int a[9]={0,5,2,4,7,1,3,2,6};
Mergesort(a,1,8);
for(i=1;i<9;i++)
{
printf("%d ",a[i]);
}
return 0;
}