核心的merge函数:
1 /* L = 左边起始位置, R = 右边起始位置, RightEnd = 右边终点位置*/ 2 void Merge( ElementType A[], ElementType TmpA[], int L, int R, int RightEnd ) 3 { /* 将有序的A[L]~A[R-1]和A[R]~A[RightEnd]归并成一个有序序列 */ 4 int LeftEnd, NumElements, Tmp; 5 int i; 6 7 LeftEnd = R - 1; /* 左边终点位置 */ 8 Tmp = L; /* 有序序列的起始位置 */ 9 NumElements = RightEnd - L + 1; 10 11 while( L <= LeftEnd && R <= RightEnd ) { 12 if ( A[L] <= A[R] ) 13 TmpA[Tmp++] = A[L++]; /* 将左边元素复制到TmpA */ 14 else 15 TmpA[Tmp++] = A[R++]; /* 将右边元素复制到TmpA */ 16 } 17 18 while( L <= LeftEnd ) 19 TmpA[Tmp++] = A[L++]; /* 直接复制左边剩下的 */ 20 while( R <= RightEnd ) 21 TmpA[Tmp++] = A[R++]; /* 直接复制右边剩下的 */ 22 23 for( i = 0; i < NumElements; i++, RightEnd -- ) 24 A[RightEnd] = TmpA[RightEnd]; /* 将有序的TmpA[]复制回A[] */ 25 }
递归算法:
1 void Msort( ElementType A[], ElementType TmpA[], int L, int RightEnd ) 2 { /* 核心递归排序函数 */ 3 int Center; 4 5 if ( L < RightEnd ) { 6 Center = (L+RightEnd) / 2; 7 Msort( A, TmpA, L, Center ); /* 递归解决左边 */ 8 Msort( A, TmpA, Center+1, RightEnd ); /* 递归解决右边 */ 9 Merge( A, TmpA, L, Center+1, RightEnd ); /* 合并两段有序序列 */ 10 } 11 }
非递归算法:
1 /* length = 当前有序子列的长度*/ 2 void Merge_pass( ElementType A[], ElementType TmpA[], int N, int length ) 3 { /* 两两归并相邻有序子列 */ 4 int i, j; 5 6 for ( i=0; i <= N-2*length; i += 2*length ) 7 Merge( A, TmpA, i, i+length, i+2*length-1 ); 8 if ( i+length < N ) /* 归并最后2个子列*/ 9 Merge( A, TmpA, i, i+length, N-1); 10 else /* 最后只剩1个子列*/ 11 for ( j = i; j < N; j++ ) TmpA[j] = A[j]; 12 } 13 14 void Merge_Sort( ElementType A[], int N ) 15 { 16 int length; 17 ElementType *TmpA; 18 19 length = 1; /* 初始化子序列长度*/ 20 TmpA = malloc( N * sizeof( ElementType ) ); 21 if ( TmpA != NULL ) { 22 while( length < N ) { 23 Merge_pass( A, TmpA, N, length ); 24 length *= 2; 25 Merge_pass( TmpA, A, N, length ); 26 length *= 2; 27 } 28 free( TmpA ); 29 } 30 else printf( "空间不足" ); 31 }
要点:
(1)为了节省空间所以两边分别倒,这在主函数里面实现
(2)注意区别奇偶区间