一、算法原理
归并排序是一种常用的排序算法,属于稳定排序法。 归并排序就是将两个已经分别排好序的数组A和B合并为一个排好序的数组C。 如果数组散乱的数组,则需要将数组元素分别按照长度为d=2^n,n=0,1,2,3,…,进行分组,然后对相邻的两组进行归并排序。具体过程就是首先取长度d=1,即将数组的每个元素作为一个子数组,然后把相邻的两个子数组作为一对进行归并排序,直到整个数组均排序完成。之后进行长度d=2,即把相邻的两个元素作为一个子数组,再对两个相邻的子数组进行归并排序,直到整个数组排序完成。之后再分别进行d=4,8,…的归并排序,直到整个数组已经是排好序的数组结束。 本文给出了归并排序算法的另外一种实现方法,不同于作者在https://blog.csdn.net/sunnyoldman001/article/details/127008922给出的实现方法。 以下Demo演示归并排序的过程。 **Demo:**假设有数据如下表所示: 第一趟归并排序: 每个元素为一个数组,总共有7个,相邻的两个数组依次进行排序,过程如下: 第二趟归并排序: 相邻的2个元素组成一个数组,不足的部分也自成一个数组,总共有4个,相邻的两个数组依次进行排序,过程如下: 第三趟归并排序: 相邻的4个元素组成一个数组,不足的部分也自成一个数组,总共有2个,相邻的两个数组依次进行排序,过程如下: 至此,归并排序结束。 由上述归并排序过程可以看出,数组中的元素在每一趟归并排序中,不一定能刚好分成偶数个数组,可能会存在“尾巴”数据。在编制算法的时候,就需要单独考虑这些“尾巴”的排序问题。当剩余的“尾巴”元素个数小于等于d的长度时,对这些数据不需要处理,如上述Demo的第一趟归并排序,剩余尾巴元素2。当剩余的“尾巴”元素个数大于d的长度小于2d长度时,需要分成2组进行归并排序,如上述Demo的第二趟排序中,剩余的尾巴元素4,7和2,则4,7组成一个数组,2组成一个数组,然后进行归并排序。 从上述Demo可以看出,归并排序的原理非常简单。只是在处理“尾巴”元素时,需要单独考虑而已。
二、归并排序算法之C程序
在https://blog.csdn.net/sunnyoldman001/article/details/127008922中给出了归并排序算法实现的方法一,该方法是把每次归并排序中用到的两个子数组单独拿出来进行排序,然后再还回到原数组中。由于不同趟归并排序中子数组的长度不同,因此每一轮次都需要重新分配子数组的大小,虽然算法对于初学者来说比较容易理解,但是这样也导致了算法实现起来比较麻烦,代码行数比较多。 本文中给出的方法二,只是多使用了一个和原数组长度相同的数组来缓存每趟排序的结果,在每一次归并排序时,从原数组直接截取子数组进行归并排序,排序后缓存到新数组的相应位置上,这样就省去了每一趟都分配子空间的麻烦。
1.两个子数组的归并排序
void Merge( int a[], int b[], int m, int n, int start )
{
int i, j, k;
i = start;
j = start + m;
k = start;
while( i < start + m && j < start + m + n )
{
if( a[i] < a[j] )
{
b[k++] = a[i++];
}
else
{
b[k++] = a[j++];
}
}
while( i < start + m )
{
b[k++] = a[i++];
}
while( j < start + m + n )
{
b[k++] = a[j++];
}
}
2.对数组进行归并排序
void MergeSort( int arr[], int len )
{
int k, n, m, x, xx;
unsigned int d;
x = int ( log(len) / log(2) );
int *b = new int[len];
if( int( pow( 2, x ) ) == len )
{
xx = x - 1;
}
else
{
xx = x;
}
for( int t = 0; t <= xx; t++ )
{
d = 1 << t;
int newLen = int( len / ( 2*d ) ) * 2 * d;
for( k = 0; k < newLen; k += 2*d )
{
Merge( arr, b, d, d, k );
}
int tail = len % (2*d);
if( tail != 0 )
{
if( tail <= d )
{
for( n = newLen; n < len; n++ )
{
b[n] = arr[n];
}
}
else if( tail > d )
{
Merge( arr, b, d, tail - d, newLen );
}
}
for( n = 0; n < len; n++ )
{
arr[n] = b[n];
}
printf( " No. %d time sort: ", t+1 );
for( n = 0; n < len; n++ )
{
printf( "%5d", arr[n] );
}
printf( "\n" );
}
delete[] b;
}
3.完整的代码
#include "stdio.h"
#include"math.h"
void MergeSort( int arr[], int len );
void Merge( int a[], int b[], int m, int n, int start );
int main()
{
int i;
int arr[] = { 3, 1, 5, 6, 7, 4, 2 };
int len = 7;
printf( " Initial Array: " );
for( i = 0; i < len; i++ )
{
printf( "%5d", arr[i] );
}
printf( "\n" );
MergeSort( arr, len );
return 0;
}
void Merge( int a[], int b[], int m, int n, int start )
{
int i, j, k;
i = start;
j = start + m;
k = start;
while( i < start + m && j < start + m + n )
{
if( a[i] < a[j] )
{
b[k++] = a[i++];
}
else
{
b[k++] = a[j++];
}
}
while( i < start + m )
{
b[k++] = a[i++];
}
while( j < start + m + n )
{
b[k++] = a[j++];
}
}
void MergeSort( int arr[], int len )
{
int k, n, m, x, xx;
unsigned int d;
x = int ( log(len) / log(2) );
int *b = new int[len];
if( int( pow( 2, x ) ) == len )
{
xx = x - 1;
}
else
{
xx = x;
}
for( int t = 0; t <= xx; t++ )
{
d = 1 << t;
int newLen = int( len / ( 2*d ) ) * 2 * d;
for( k = 0; k < newLen; k += 2*d )
{
Merge( arr, b, d, d, k );
}
int tail = len % (2*d);
if( tail != 0 )
{
if( tail <= d )
{
for( n = newLen; n < len; n++ )
{
b[n] = arr[n];
}
}
else if( tail > d )
{
Merge( arr, b, d, tail - d, newLen );
}
}
for( n = 0; n < len; n++ )
{
arr[n] = b[n];
}
printf( " No. %d time sort: ", t+1 );
for( n = 0; n < len; n++ )
{
printf( "%5d", arr[n] );
}
printf( "\n" );
}
delete[] b;
}
4.测试用例 测试用例一 测试用例二 测试用例三 测试用例四
|