这个算法实现起来花了一些时间,主要错误出在临近问题上,有运行时错误调试了不短的时间。
这个算法对比之前的几类算法,运行复杂度是都有所降低,如果去g=3*g+1,这个算法的复杂度维持在O(N^1.25)左右,对比之下,这是一个高效的算法,十分值得学习与借鉴。
有利于对数据进行灵活的排序处理
#include <iostream>
#include <cstdio>
#include <vector>
using namespace std;
const int max_n=1e6;
int a[max_n];
int n,m,cnt;
vector<int> G;
void Insert_sort(int g)
{
for(int i=g;i<n;++i)
{
int tmp=a[i];
int j;
for(j=i-g;j>=0 && tmp<a[j];j-=g)
{
a[j+g]=a[j];
++cnt;
}
a[j+g]=tmp;
}
}
void Shell_sort()
{
// 首先,要生成数列G.
int g=1;
m=0;
while(g<=n)
{
G.push_back(g);
++m;
g=3*g+1;
}
cnt=0;
for(int i=G.size()-1;i>=0;--i)
{
Insert_sort(G[i]);
}
}
int main()
{
// 输入数据
scanf("%d",&n);
for(int i=0;i<n;++i)
{
scanf("%d",&a[i]);
}
// 希尔排序
Shell_sort();
printf("%d\n",m);
for(int i=m-1;i>0;--i)
{
printf("%d ",G[i]);
}
printf("%d\n",G[0]);
printf("%d\n",cnt);
for(int i=0;i<n;++i)
{
printf("%d\n",a[i]);
}
return 0;
}