排序---希尔排序及希尔增量与Hibbard增量的比较

前言

个人小记

---

一、比较

希尔（n^2）：n/2,n/4,n/8,…
Hibbard(n^1.5):1,3,7,15,…

二、代码

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>
#define MAX_ARR 100000
#define swap(a,b)\
{\__typeof(a) __c=a;\a=b,b=__c;\
}
#define TEST(func,arr,l,n)\
{\printf("test:%s\n",#func);\int* t=(int *)malloc(sizeof(int )*n);\memcpy(t,arr,n);\long long a=clock();\func(t,l,n);\long long b=clock();\if(check(t,n))printf("OK %lldms\n",(b-a)*1000/CLOCKS_PER_SEC);\else printf("FAIL\n");\free(t);\
}int check(int *t,int n)
{for(int i=1;i<n;i++){if(t[i-1]>t[i])return 0;}return 1;
}int *init_arr(int n)
{int* arr=(int*)malloc(sizeof(int)*n);for(int i=0;i<n;i++)arr[i]=rand()%100000;return arr;
}void UC_insert(int *arr,int l,int r,int step)
{int min=l;for(int i=step+l;i<r;i+=step){if(arr[min]>arr[i])min=i;}while(min>l){swap(arr[min],arr[min-step]);min=min-step;}for(int i=l+2*step;i<r;i+=step){int j=i;while(arr[j]<arr[j-step]){swap(arr[j],arr[j-step]);j=j-step;}}return ;
}void shell_sort(int *arr,int l,int r)
{int k=2,n=r-l,step;do{step=(n/k)==0?1:(n/k);for(int i=l;i<l+step;i++){UC_insert(arr,i,r,step);}k*=2;}while(step!=1);return ;
}void shell_sort_hibbard(int* arr,int l,int r)
{int step=1,n=r-l;while(step<=n/2)step=step*2+1;do{step/=2;for(int i=l;i<step+l;i++){UC_insert(arr,i,r,step);}}while(step>1);return ;
}int main()
{srand((unsigned)time(0));int *arr=init_arr(MAX_ARR);TEST(shell_sort,arr,0,MAX_ARR);TEST(shell_sort_hibbard,arr,0,MAX_ARR);free(arr);return 0;
}

二、测试结果

test:shell_sort
OK 13ms
test:shell_sort_hibbard
OK 12ms