[算法导论]归并排序

1》归并排序的步骤如下：

       Divide: 把长度为n的输入序列分成两个长度为n/2的子序列。

       Conquer: 对这两个子序列分别采用归并排序。      

       Combine: 将两个排序好的子序列合并成一个最终的排序序列。

2》时间复杂度：

       这是一个递推公式（Recurrence），我们需要消去等号右侧的T(n)，把T(n)写成n的函数。其实符合一定条件的Recurrence的展开有数学公式可以套，这里我们略去严格的数学证明，只是从直观上看一下这个递推公式的结果。当n=1时可以设T(1)=c₁，当n>1时可以设T(n)=2T(n/2)+c₂n，我们取c₁和c₂中较大的一个设为c，把原来的公式改为：

       这样计算出的结果应该是T(n)的上界。下面我们把T(n/2)展开成2T(n/4)+cn/2（下图中的(c)），然后再把T(n/4)进一步展开，直到最后全部变成T(1)=c（下图中的(d)）：

       把图(d)中所有的项加起来就是总的执行时间。这是一个树状结构，每一层的和都是cn，共有lgn+1层，因此总的执行时间是cnlgn+cn，相比nlgn来说，cn项可以忽略，因此T(n)的上界是Θ(nlgn)。

       如果先前取c₁和c₂中较小的一个设为c，计算出的结果应该是T(n)的下界，然而推导过程一样，结果也是Θ(nlgn)。既然T(n)的上下界都是Θ(nlgn)，显然T(n)就是Θ(nlgn)。

#include <iostream>
#include <cstdlib>

#define Max 65535
using namespace std;

void Merge(int *A, int p, int q, int r)
{
   int n1 = q - p + 1;
   int n2 = r - q;
   //Let L[0.n1] and R[0.n2] be new arrays
   int L[n1], R[n2];
   for (int i = 0; i < n1; i ++)
   {
      L[i] = A[p + i];
   }
   for (int i = 0; i < n2; i ++)
   {
      R[i] = A[q + i + 1];
   }
   L[n1] = Max;
   R[n2] = Max;
   int i = 0, j = 0;
   for (int k = p; k <= r; k ++)
   {
      if (L[i] <= R[j])
      {
         A[k] = L[i];
         i += 1;
      }
      else
      {
         A[k] = R[j];
         j += 1;
      }
   }
}
void MergeSort(int *A, 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 n;
    int *A;
    cout<<"请输入要排序的数目：";
    cin>>n;
    A = (int *)malloc((n + 1) * sizeof(int));
    cout<<"请依次输入要排序的数列：";
    for (int i = 0; i < n; i ++)
    {
       cin>>A[i];
    }
    MergeSort(A, 0, n - 1);
    cout<<"排序后结果为："<<endl;
    for (int i = 0; i < n; i ++)
    {
       cout<<A[i]<<" ";
    }
    cout<<endl;
    return 0;
}