#include <time.h>
#include <stdlib.h>
#include <stdio.h>
#define sential RAND_MAX /* 定义哨兵*/
#define SIZE 1000000 /*可以设定size的大小来测试该算法的性能*/
void MergeSort(int *arr,int p ,int r);
void Merge(int arr[],int p ,int q,int r);
void main(void)
{ //分配内存、建立随机数组、打印排序前的数组
int randnumber;
int *pt=(int *)malloc(SIZE*sizeof(int));
if(!pt)
{ puts("Fail to allocate memory");
exit(EXIT_FAILURE);
}
puts("Before Sorting");
for(int i=0;i<SIZE;i++)
{ randnumber=rand()%SIZE;
pt[i]=randnumber;
printf("%7d",randnumber);
if((i+1)%10==0)
putchar('\n');
}
//排序、计时
putchar('\n');
double time_start=time(0);
MergeSort(pt,0,SIZE-1);
double time_end=time(0);
puts("After Sorting");
//打印排序后的数组
for(int i=0;i<SIZE;i++)
{ printf("%7d",pt[i]);
if((1+i)%10==0)
putchar('\n');
}
printf("Consume:%fs\n",time_end-time_start);
putchar('\n');
free(pt);
}
//合并排序函数排序
void MergeSort(int *arr,int p,int r)
{ if(p<r)
{
int q=(p+r)/2;
MergeSort(arr,p,q);
MergeSort(arr,q+1,r);
Merge(arr,p,q,r);
}
}
//辅助合并函数
void Merge(int arr[],int p,int q,int r)
{ int n1=q-p+2;
int n2=r-q+1;
int l_arr[n1];
int r_arr[n2];
for(int i=0;i<n1-1;i++)
l_arr[i]=arr[p+i];
l_arr[n1-1]=sential;
for(int i=0;i<(n2-1);i++)
r_arr[i]=arr[q+1+i];
r_arr[n2-1]=sential;
int k=0;
int j=0;
for(int i=p;i<1+r;i++)
{ if(l_arr[j]<r_arr[k])
{ arr[i]=l_arr[j];
j++;
}
else
{ arr[i]=r_arr[k];
k++;
}
}
}