初始：d＝549 38 65 97 76 13 27 49* 55 04
49 13
|-------------------|
38 27
|-------------------|
65 49*
|-------------------|
97 55
|-------------------|
76 04
|-------------------|
一趟结果
13 27 49* 55 04 49 38 65 97 76
d=3
13 27 49* 55 04 49 38 65 97 76
13 55 38 76
|------------|------------|------------|
27 04 65
|------------|------------|
49* 49 97
|------------|------------|
二趟结果
13 04 49* 38 27 49 55 65 97 76
d=1
13 04 49* 38 27 49 55 65 97 76
|----|----|----|----|----|----|----|----|----|
三趟结果
04 13 27 38 49* 49 55 65 76 97
--------------------------------------------------------------------------------------------
例如对503,17,512,908,170,7,275,653,462，154，509，612，677，765，703，94排序的C语言算法
================================================
功能：希尔排序
输入：数组名称（也就是数组首地址）、数组中元素个数
================================================
*/
/*
====================================================
算法思想简单描述：
在直接插入排序算法中，每次插入一个数，使有序序列只增加1个节点，
并且对插入下一个数没有提供任何帮助。如果比较相隔较远距离（称为
增量）的数，使得数移动时能跨过多个元素，则进行一次比较就可能消除
多个元素交换。D.L.shell于1959年在以他名字命名的排序算法中实现
了这一思想。算法先将要排序的一组数按某个增量d分成若干组，每组中
记录的下标相差d.对每组中全部元素进行排序，然后再用一个较小的增量
对它进行，在每组中再进行排序。当增量减到1时，整个要排序的数被分成
一组，排序完成。
下面的函数是一个希尔排序算法的一个实现，初次取序列的一半为增量，
以后每次减半，直到增量为1。
希尔排序是不稳定的。
=====================================================
*/
void shell_sort(int *x, int n)
{
int h, j, k, t;
for (h=n/2; h>0; h=h/2) /*控制增量*/
{
for (j=h; j{t = *(x+j);
for (k=j-h; (k>=0 && t<*(x+k)); k-=h)
{
*(x+k+h) = *(x+k);
}
*(x+k+h) = t;
}
}
}
void main()
{
#define MAX 16
int *p, i, a[MAX];
/*录入测试数据*/
/*
p = a;
printf(\"Input %d number for sorting :\\n\
for (i=0; i{scanf(\"%d\
}
*可以自己输入数据
*/
a[] = {503,17,512,908,170,7,275,653,462，154，509，612，677，765，703，94};
printf(\"\\n\");
//503,17,512,908,170,7,275,653,462，154，509，612，677，765，703，94
/*测试排序*/
p = a;
shell_sort(p,MAX);
/**/
for (p=a, i=0; i{printf(\"%d \
}
printf(\"\\n\");
system(\"pause\");
}
pascal算法程序：
program xepx;
const n=7;
type
arr=array[1..n] of integer;
var
a:arr;
i,j,t,d:integer;
bool:boolean;
begin
write('input data:');
for i:=1 to n do read(a[i]); writeln; d:=n; while d>1 do
begin d:=d div 2; for j:=d+1 to n do begin t:=a[j]; i:=j-d; while (i>0) and (a[i]>t) do begin a[i+d]:=a[i]; i:=i-d; end; a[i+d]:=t; end;
end; write('output data:'); for i:=1 to n do write(a:6); writeln; end.