最基础排序算法—归并排序

归并排序算法完全遵循分治模式。直观上操作如下;

1. 分解： 分解待排序的n个元素的序列成各具n/2个元素的两个子序列。
2. 解决： 使用归并排序递归的排序两个子序列
3. 合并： 合并两个已排序的子序列来产生已排序的答案

当待排序的序列长度为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;
}