Divide: 把长度为n的输入序列分成两个长度为n/2的子序列。
Conquer: 对这两个子序列分别采用归并排序。
Combine: 将两个排序好的子序列合并成一个最终的排序序列。
2》时间复杂度:
这是一个递推公式(Recurrence),我们需要消去等号右侧的T(n),把T(n)写成n的函数。其实符合一定条件的Recurrence的展开有数学公式可以套,这里我们略去严格的数学证明,只是从直观上看一下这个递推公式的结果。当n=1时可以设T(1)=c1,当n>1时可以设T(n)=2T(n/2)+c2n,我们取c1和c2中较大的一个设为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)。
如果先前取c1和c2中较小的一个设为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;
}