A-A+

[转]排序算法一览

2008年11月09日 未分类 暂无评论 阅读 1 次

基本概念
排序(Sorting)是计算机程序设计中的一种重要操作,其功能是对一个数据元素集合或序列重新排列成一个按数据元素某个项值有序的序列。作为排序依据的数据项称为“排序码”,也即数据元素的关键码。为了便于查找,通常希望计算机中的数据表是按关键码有序的。如有序表的折半查找,查找效率较高。还有,二叉排序树、B-树和B+树的构造过程就是一个排序过程。若关键码是主关键码,则对于任意待排序序列,经排序后得到的结果是唯一的;若关键码是次关键码,排序结果可能不唯一,这是因为具有相同关键码的数据元素,这些元素在排序结果中,它们之间的的位置关系与排序前不能保持。
若对任意的数据元素序列,使用某个排序方法,对它按关键码进行排序:若相同关键码元素间的位置关系,排序前与排序后保持一致,称此排序方法是稳定的;而不能保持一致的排序方法则称为不稳定的。
排序分为两类:内排序和外排序。
内排序:指待排序列完全存放在内存中所进行的排序过程,适合不太大的元素序列。
外排序:指排序过程中还需访问外存储器,足够大的元素序列,因不能完全放入内存,只能使用外排序。
1.插入排序
1.1直接插入排序
设有n个记录,存放在数组r中,重新安排记录在数组中的存放顺序,使得按关键码有序。即r[1].key≤r[2].key≤……≤r[n].key

先来看看向有序表中插入一个记录的方法:
设1<j≤n,r[1].key≤r[2].key≤……≤r[j-1].key,将r[j]插入,重新安排存放顺序,使得r[1].key≤r[2].key≤……≤r[j].key,得到新的有序表,记录数增1。

【算法1】
① r[0]=r[j]; //r[j]送r[0]中,使r[j]为待插入记录空位
i=j-1; //从第i个记录向前测试插入位置,用r[0]为辅助单元,可免去测试i<1。
②若r[0].key≥r[i].key,转④。 //插入位置确定
③若r[0].key < r[i].key时,
r[i+1]=r[i];i=i-1;转②。 //调整待插入位置
④ r[i+1]=r[0];结束。 //存放待插入记录
【例10.1】向有序表中插入一个记录的过程如下:
r[1] r[2] r[3] r[4] r[5] 存储单元
2 10 18 25 9 将r[5]插入四个记录的有序表中,j=5
r[0]=r[j];i=j-1;初始化,设置待插入位置
2 10 18 25 □ r[i+1]为待插入位置
i=4,r[0] < r[i],r[i+1]=r[i];i--;调整待插入位置
2 10 18 □ 25
i=3,r[0] < r[i],r[i+1]=r[i];i--;调整待插入位置
2 10 □ 18 25
i=2,r[0] < r[i],r[i+1]=r[i];i--;调整待插入位置
2 □ 10 18 25
i=1,r[0] ≥r[i],r[i+1]=r[0];插入位置确定,向空位填入插入记录
2 9 10 18 25 向有序表中插入一个记录的过程结束

直接插入排序方法:仅有一个记录的表总是有序的,因此,对n个记录的表,可从第二个记录开始直到第n个记录,逐个向有序表中进行插入操作,从而得到n个记录按关键码有序的表。

【算法2】
void InsertSort(S_TBL &p)
{ for(i=2;i<=p->length;i++)
if(p->elem[i].key < p->elem[i-1].key) /*小于时,需将elem[i]插入有序表*/
{ p->elem[0].key=p->elem[i].key; /*为统一算法设置监测*/
for(j=i-1;p->elem[0].key < p->elem[j].key;j--)
p->elem[j+1].key=p->elem[j].key; /*记录后移*/
p->elem[j+1].key=p->elem[0].key; /*插入到正确位置*/
}
}
【效率分析】
空间效率:仅用了一个辅助单元。
时间效率:向有序表中逐个插入记录的操作,进行了n-1趟,每趟操作分为比较关键码和移动记录,而比较的次数和移动记录的次数取决于待排序列按关键码的初始排列。
最好情况下:即待排序列已按关键码有序,每趟操作只需1次比较2次移动。
总比较次数=n-1次
总移动次数=2(n-1)次
最坏情况下:即第j趟操作,插入记录需要同前面的j个记录进行j次关键码比较,移动记录的次数为j+2次。
平均情况下:即第j趟操作,插入记录大约同前面的j/2个记录进行关键码比较,移动记录的次数为j/2+2次。
由此,直接插入排序的时间复杂度为O(n2)。是一个稳定的排序方法。

1.2折半插入排序
直接插入排序的基本操作是向有序表中插入一个记录,插入位置的确定通过对有序表中记录按关键码逐个比较得到的。平均情况下总比较次数约为n2/4。既然是在有序表中确定插入位置,可以不断二分有序表来确定插入位置,即一次比较,通过待插入记录与有序表居中的记录按关键码比较,将有序表一分为二,下次比较在其中一个有序子表中进行,将子表又一分为二。这样继续下去,直到要比较的子表中只有一个记录时,比较一次便确定了插入位置。
二分判定有序表插入位置方法:
① low=1;high=j-1;r[0]=r[j]; // 有序表长度为j-1,第j个记录为待插入记录
//设置有序表区间,待插入记录送辅助单元
②若low>high,得到插入位置,转⑤
③ low≤high,m=(low+high)/2; // 取表的中点,并将表一分为二,确定待插入区间*/
④若r[0].key<r[m].key,high=m-1; //插入位置在低半区
否则,low=m+1; // 插入位置在高半区
转②
⑤ high+1即为待插入位置,从j-1到high+1的记录,逐个后移,r[high+1]=r[0];放置待插入记录。

【算法3】
void InsertSort(S_TBL *s)
{ /* 对顺序表s作折半插入排序 */
for(i=2;i<=s->length;i++)
{ s->elem[0]=s->elem[i]; /* 保存待插入元素 */
low=i;high=i-1; /* 设置初始区间 */
while(low<=high) /* 该循环语句完成确定插入位置 */
{ mid=(low+high)/2;
if(s->elem[0].key>s->elem[mid].key)
low=mid+1; /* 插入位置在高半区中 */
else high=mid-1; /* 插入位置在低半区中 */
}/* while */
for(j=i-1;j>=high+1;j--) /* high+1为插入位置 */
s->elem[j+1]=s->elem[j]; /* 后移元素,留出插入空位 */
s->elem[high+1]=s->elem[0]; /* 将元素插入 */
}/* for */
}/* InsertSort */
【时间效率】
确定插入位置所进行的折半查找,关键码的比较次数至多为 ,次,移动记录的次数和直接插入排序相同,故时间复杂度仍为O(n2)。是一个稳定的排序方法。

1.3表插入排序

直接插入排序、折半插入排序均要大量移动记录,时间开销大。若要不移动记录完成排序,需要改变存储结构,进行表插入排序。所谓表插入排序,就是通过链接指针,按关键码的大小,实现从小到大的链接过程,为此需增设一个指针项。操作方法与直接插入排序类似,所不同的是直接插入排序要移动记录,而表插入排序是修改链接指针。用静态链表来说明。
#define SIZE 200
typedef struct{
ElemType elem; /*元素类型*/
int next; /*指针项*/
}NodeType; /*表结点类型*/
typedef struct{
NodeType r[SIZE]; /*静态链表*/
int length; /*表长度*/
}L_TBL; /*静态链表类型*/
假设数据元素已存储在链表中,且0号单元作为头结点,不移动记录而只是改变链指针域,将记录按关键码建为一个有序链表。首先,设置空的循环链表,即头结点指针域置0,并在头结点数据域中存放比所有记录关键码都大的整数。接下来,逐个结点向链表中插入即可。
【例1】表插入排序示例

MAXINT 49 38 65 97 76 13 27 49
0 - - - - - - - -

MAXINT 49 38 65 97 76 13 27 49
1 0 - - - - - - -

MAXINT 49 38 65 97 76 13 27 49
2 0 1 - - - - - -

MAXINT 49 38 65 97 76 13 27 49
2 3 1 0 - - - - -

MAXINT 49 38 65 97 76 13 27 49
2 3 1 4 0 - - - -

MAXINT 49 38 65 97 76 13 27 49
2 3 1 5 0 4 - - -

MAXINT 49 38 65 97 76 13 27 49
6 3 1 5 0 4 2 - -

MAXINT 49 38 65 97 76 13 27 49
6 3 1 5 0 4 7 2 -

MAXINT 49 38 65 97 76 13 27 49
6 8 1 5 0 4 7 2 3

【算法4】
1. j=l->r[0].next;i=1; //指向第一个记录位置,从第一个记录开始调整
2. 若i=l->length时,调整结束;否则,
a. 若i=j,j=l->r[j].next;i++;转(2) //数据元素应在这分量中,不用调整,处理下一个结点
b. 若j>i,l->r[i].elem<-->l->r[j].elem; //交换数据元素
p=l->r[j].next; // 保存下一个结点地址
l->r[j].next=l->[i].next;l->[i].next=j; // 保持后续链表不被中断
j=p;i++;转(2) // 指向下一个处理的结点
c. 若j<i,while(j<i) j=l->r[j].next;//j分量中原记录已移走,沿j的指针域找寻原记录的位置
转到(a)
【例2】对表插入排序结果进行重排示例

MAXINT 49 38 65 97 76 13 27 49
6 8 1 5 0 4 7 2 3

MAXINT 13 38 65 97 76 49 27 49
6 (6) 1 5 0 4 8 2 3

MAXINT 13 27 65 97 76 49 38 49
6 (6) (7) 5 0 4 8 1 3

MAXINT 13 27 38 97 76 49 65 49
6 (6) (7) (7) 0 4 8 5 3

MAXINT 13 27 38 49 76 97 65 49
6 (6) (7) (7) (6) 4 0 5 3

MAXINT 13 27 38 49 49 97 65 76
6 (6) (7) (7) (6) (8) 0 5 4

MAXINT 13 27 38 49 49 65 97 76
6 (6) (7) (7) (6) (8) (7) 0 4

MAXINT 13 27 38 49 49 65 76 97
6 (6) (7) (7) (6) (8) (7) (8) 0

【时效分析】
表插入排序的基本操作是将一个记录插入到已排好序的有序链表中,设有序表长度为i,则需要比较至多i+1次,修改指针两次。因此,总比较次数与直接插入排序相同,修改指针总次数为2n次。所以,时间复杂度仍为O(n2)

1.4希尔排序(Shell’s Sort)

希尔排序又称缩小增量排序,是1959年由D.L.Shell提出来的,较前述几种插入排序方法有较大的改进。
直接插入排序算法简单,在n值较小时,效率比较高,在n值很大时,若序列按关键码基本有序,效率依然较高,其时间效率可提高到O(n)。希尔排序即是从这两点出发,给出插入排序的改进方法。
希尔排序方法:
1. 选择一个步长序列t1,t2,…,tk,其中ti>tj,tk=1;
2. 按步长序列个数k,对序列进行k趟排序;
3. 每趟排序,根据对应的步长ti,将待排序列分割成若干长度为m的子序列,分别对各子表进行直接插入排序。仅步长因子为1时,整个序列作为一个表来处理,表长度即为整个序列的长度。

【例4】待排序列为 39,80,76,41,13,29,50,78,30,11,100,7,41,86。
步长因子分别取5、3、1,则排序过程如下:
p=5 39 80 76 41 13 29 50 78 30 11 100 7 41 86
└─────────┴─────────┘
└─────────┴──────────┘
└─────────┴──────────┘
└─────────┴──────────┘
└─────────┘
子序列分别为{39,29,100},{80,50,7},{76,78,41},{41,30,86},{13,11}。

第一趟排序结果:
p=3 29 7 41 30 11 39 50 76 41 13 100 80 78 86
└─────┴─────┴─────┴──────┘
└─────┴─────┴─────┴──────┘
└─────┴─────┴──────┘
子序列分别为{29,30,50,13,78},{7,11,76,100,86},{41,39,41,80}。
第二趟排序结果:
p=1 13 7 39 29 11 41 30 76 41 50 86 80 78 100
此时,序列基本“有序”,对其进行直接插入排序,得到最终结果:
7 11 13 29 30 39 41 41 50 76 78 80 86 100

【算法5】
void ShellInsert(S_TBL &p,int dk)
{ /*一趟增量为dk的插入排序,dk为步长因子*/
for(i=dk+1;i<=p->length;i++)
if(p->elem[i].key < p->elem[i-dk].key) /*小于时,需elem[i]将插入有序表*/
{ p->elem[0]=p->elem[i]; /*为统一算法设置监测*/
for(j=i-dk;j>0&&p->elem[0].key < p->elem[j].key;j=j-dk)
p->elem[j+dk]=p->elem[j]; /*记录后移*/
p->elem[j+dk]=p->elem[0]; /*插入到正确位置*/
}
}

void ShellSort(S_TBL *p,int dlta[],int t)
{ /*按增量序列dlta[0,1…,t-1]对顺序表*p作希尔排序*/
for(k=0;k<t;t++)
ShellSort(p,dlta[k]); /*一趟增量为dlta[k]的插入排序*/
}

【时效分析】
希尔排序时效分析很难,关键码的比较次数与记录移动次数依赖于步长因子序列的选取,特定情况下可以准确估算出关键码的比较次数和记录的移动次数。目前还没有人给出选取最好的步长因子序列的方法。步长因子序列可以有各种取法,有取奇数的,也有取质数的,但需要注意:步长因子中除1外没有公因子,且最后一个步长因子必须为1。希尔排序方法是一个不稳定的排序方法。

2 交换排序
交换排序主要是通过两两比较待排记录的关键码,若发生与排序要求相逆,则交换之。
2.1冒泡排序(Bubble Sort)
先来看看待排序列一趟冒泡的过程:设1<j≤n,r[1],r[2],···,r[j]为待排序列,
通过两两比较、交换,重新安排存放顺序,使得r[j]是序列中关键码最大的记录。一趟冒泡方法为:
① i=1; //设置从第一个记录开始进行两两比较
②若i≥j,一趟冒泡结束。
③比较r[i].key与r[i+1].key,若r[i].key≤r[i+1].key,不交换,转⑤
④当r[i].key>r[i+1].key时, r[0]=r[i];r[i]=r[i+1];r[i+1]=r[0];
将r[i]与r[i+1]交换
⑤ i=i+1;调整对下两个记录进行两两比较,转②
冒泡排序方法:对n个记录的表,第一趟冒泡得到一个关键码最大的记录r[n],第二趟冒泡对n-1个记录的表,再得到一个关键码最大的记录r[n-1],如此重复,直到n个记录按关键码有序的表。
【算法6】
① j=n; //从n记录的表开始
②若j<2,排序结束
③ i=1; //一趟冒泡,设置从第一个记录开始进行两两比较,
④若i≥j,一趟冒泡结束,j=j-1;冒泡表的记录数-1,转②
⑤比较r[i].key与r[i+1].key,若r[i].key≤r[i+1].key,不交换,转⑤
⑥当r[i].key>r[i+1].key时, r[i]<-->r[i+1];将r[i]与r[i+1]交换
⑦ i=i+1;调整对下两个记录进行两两比较,转④
【效率分析】
空间效率:仅用了一个辅助单元。
时间效率:总共要进行n-1趟冒泡,对j个记录的表进行一趟冒泡需要j-1次关键码比较。
移动次数:
最好情况下:待排序列已有序,不需移动。
2.2快速排序
快速排序是通过比较关键码、交换记录,以某个记录为界(该记录称为支点),将待排序列分成两部分。其中,一部分所有记录的关键码大于等于支点记录的关键码,另一部分所有记录的关键码小于支点记录的关键码。我们将待排序列按关键码以支点记录分成两部分的过程,称为一次划分。对各部分不断划分,直到整个序列按关键码有序。
一次划分方法:
设1≤p<q≤n,r

,r
,...,r[q]为待排序列
① low=p;high=q; //设置两个搜索指针,low是向后搜索指针,high是向前搜索指针
r[0]=r[low]; //取第一个记录为支点记录,low位置暂设为支点空位
②若low=high,支点空位确定,即为low。
r[low]=r[0]; //填入支点记录,一次划分结束
否则,low<high,搜索需要交换的记录,并交换之
③若low<high且r[high].key≥r[0].key //从high所指位置向前搜索,至多到low+1位置
high=high-1;转③ //寻找r[high].key<r[0].key
r[low]=r[high]; //找到r[high].key<r[0].key,设置high为新支点位置,
//小于支点记录关键码的记录前移。
④若low<high且r[low].key<r[0].key //从low所指位置向后搜索,至多到high-1位置
low=low+1;转④ //寻找r[low].key≥r[0].key
r[high]=r[low]; //找到r[low].key≥r[0].key,设置low为新支点位置,
//大于等于支点记录关键码的记录后移。
转② //继续寻找支点空位

【算法7】
int Partition(S_TBL *tbl,int low,int high) /*一趟快排序*/
{ /*交换顺序表tbl中子表tbl->[low…high]的记录,使支点记录到位,并反回其所在位置*/
/*此时,在它之前(后)的记录均不大(小)于它*/
tbl->r[0]=tbl->r[low]; /*以子表的第一个记录作为支点记录*/
pivotkey=tbl->r[low].key; /*取支点记录关键码*/
while(low<higu) /*从表的两端交替地向中间扫描*/
{ while(low<high&&tbl->r[high].key>=pivotkey) high--;
tbl->r[low]=tbl->r[high]; /*将比支点记录小的交换到低端*/
while(low<high&&tbl-g>r[high].key<=pivotkey) low++;
tbl->r[low]=tbl->r[high]; /*将比支点记录大的交换到低端*/
}
tbl->r[low]=tbl->r[0]; /*支点记录到位*/
return low; /*反回支点记录所在位置*/
}
【例3】一趟快排序过程示例
r[1] r[2] r[3] r[4] r[5] r[6] r[7] r[8] r[9] r[10] 存储单元
49 14 38 74 96 65 8 49 55 27 记录中关键码
low=1;high=10;设置两个搜索指针, r[0]=r[low];支点记录送辅助单元,
□ 14 38 74 96 65 8 49 55 27
↑ ↑
low high
第一次搜索交换
从high向前搜索小于r[0].key的记录,得到结果
27 14 38 74 96 65 8 49 55 □
↑ ↑
low high
从low向后搜索大于r[0].key的记录,得到结果
27 14 38 □ 96 65 8 49 55 74
↑ ↑
low high
第二次搜索交换
从high向前搜索小于r[0].key的记录,得到结果
27 14 38 8 96 65 □ 49 55 74
↑ ↑
low high
从low向后搜索大于r[0].key的记录,得到结果
27 14 38 8 □ 65 96 49 55 74
↑ ↑
low high
第三次搜索交换
从high向前搜索小于r[0].key的记录,得到结果
27 14 38 8 □ 65 96 49 55 74
↑↑
low high
从low向后搜索大于r[0].key的记录,得到结果
27 14 38 8 □ 65 96 49 55 74
↑↑
low high
low=high,划分结束,填入支点记录
27 14 38 8 49 65 96 49 55 74

【算法8】
void QSort(S_TBL *tbl,int low,int high) /*递归形式的快排序*/
{ /*对顺序表tbl中的子序列tbl->[low…high]作快排序*/
if(low<high)
{ pivotloc=partition(tbl,low,high); /*将表一分为二*/
QSort(tbl,low,pivotloc-1); /*对低子表递归排序*/
QSort(tbl,pivotloc+1,high); /*对高子表递归排序*/
}
}

【效率分析】
空间效率:快速排序是递归的,每层递归调用时的指针和参数均要用栈来存放,递归调用层次数与上述二叉树的深度一致。因而,存储开销在理想情况下为O(log2n),即树的高度;在最坏情况下,即二叉树是一个单链,为O(n)。
时间效率:在n个记录的待排序列中,一次划分需要约n次关键码比较,时效为O(n),若设T(n)为对n个记录的待排序列进行快速排序所需时间。
理想情况下:每次划分,正好将分成两个等长的子序列,则

T(n)≤cn+2T(n/2) c是一个常数
≤cn+2(cn/2+2T(n/4))=2cn+4T(n/4)
≤2cn+4(cn/4+T(n/8))=3cn+8T(n/8)
······
≤cnlog2n+nT(1)=O(nlog2n)

最坏情况下:即每次划分,只得到一个子序列,时效为O(n2)。
快速排序是通常被认为在同数量级(O(nlog2n))的排序方法中平均性能最好的。但若初始序列按关键码有序或基本有序时,快排序反而蜕化为冒泡排序。为改进之,通常以“三者取中法”来选取支点记录,即将排序区间的两个端点与中点三个记录关键码居中的调整为支点记录。快速排序是一个不稳定的排序方法。

3选择排序
选择排序主要是每一趟从待排序列中选取一个关键码最小的记录,也即第一趟从n个记录中选取关键码最小的记录,第二趟从剩下的n-1个记录中选取关键码最小的记录,直到整个序列的记录选完。这样,由选取记录的顺序,便得到按关键码有序的序列。
3.1 简单选择排序
操作方法:第一趟,从n个记录中找出关键码最小的记录与第一个记录交换;第二趟,从第二个记录开始的n-1个记录中再选出关键码最小的记录与第二个记录交换;如此,第i趟,则从第i个记录开始的n-i+1个记录中选出关键码最小的记录与第i个记录交换,直到整个序列按关键码有序。

【算法9】
void SelectSort(S_TBL *s)
{ for(i=1;i<s->length;i++)
{ /* 作length-1趟选取 */
for(j=i+1,t=i;j<=s->length;j++)
{ /* 在i开始的length-n+1个记录中选关键码最小的记录 */
if(s->elem[t].key>s->elem[j].key)
t=j; /* t中存放关键码最小记录的下标 */
}
s->elem[t]<-->s->elem[i]; /* 关键码最小的记录与第i个记录交换 */
}
}

从程序中可看出,简单选择排序移动记录的次数较少,但关键码的比较次数依然是

3.2树形选择排序
按照锦标赛的思想进行,将n个参赛的选手看成完全二叉树的叶结点,则该完全二叉树有2n-2或2n-1个结点。首先,两两进行比赛(在树中是兄弟的进行,否则轮空,直接进入下一轮),胜出的再兄弟间再两两进行比较,直到产生第一名;接下来,将作为第一名的结点看成最差的,并从该结点开始,沿该结点到根路径上,依次进行各分枝结点子女间的比较,胜出的就是第二名。因为和他比赛的均是刚刚输给第一名的选手。如此,继续进行下去,直到所有选手的名次排定。

【例4】16个选手的比赛(n=24)

从叶结点开始的兄弟间两两比赛,胜者上升到父结点;胜者兄弟间再两两比赛,直到根结点,产生第一名91。比较次数为 23+22+21+20=24-1=n-1。
图10.6中,将第一名的结点置为最差的,与其兄弟比赛,胜者上升到父结点,胜者兄弟间再比赛,直到根结点,产生第二名83。比较次数为4,即log2n次。其后各结点的名次均是这样产生的,所以,对于n个参赛选手来说,即对1,故时间复杂度为O(nlog2n)。该方法占用空间较多,除需输出排序结果的n个单元外,尚需n-1个辅助单元。-n+1)log2n-n个记录进行树形选择排序,总的关键码比较次数至多为(n

3.3 堆排序(Heap Sort)
若以一维数组存储一个堆,则堆对应一棵完全二叉树,且所有非叶结点的值均不大于(或不小于)其子女的值,根结点的值是最小(或最大)的。
设有n个元素,将其按关键码排序。首先将这n个元素按关键码建成堆,将堆顶元素输出,得到n个元素中关键码最小(或最大)的元素。然后,再对剩下的n-1个元素建成堆,输出堆顶元素,得到n个元素中关键码次小(或次大)的元素。如此反复,便得到一个按关键码有序的序列。称这个过程为堆排序。
因此,实现堆排序需解决两个问题:
1. 如何将n个元素的序列按关键码建成堆;
2. 输出堆顶元素后,怎样调整剩余n-1个元素,使其按关键码成为一个新堆。
首先,讨论输出堆顶元素后,对剩余元素重新建成堆的调整过程。
调整方法:设有m个元素的堆,输出堆顶元素后,剩下m-1个元素。将堆底元素送入堆顶,堆被破坏,其原因仅是根结点不满足堆的性质。将根结点与左、右子女中较小(或小大)的进行交换。若与左子女交换,则左子树堆被破坏,且仅左子树的根结点不满足堆的性质;若与右子女交换,则右子树堆被破坏,且仅右子树的根结点不满足堆的性质。继续对不满足堆性质的子树进行上述交换操作,直到叶子结点,堆被建成。称这个自根结点到叶子结点的调整过程为筛选。
【算法10】
void HeapAdjust(S_TBL *h,int s,int m)
{/*r中的记录关键码除r外均满足堆的定义,本函数将对第s个结点为根的子树筛选,使其成为大顶堆*/
rc=h->r
for(j=2*s;j<=m;j=j*2) /* 沿关键码较大的子女结点向下筛选 */
{ if(j<m&&h->r[j].key<h->r[j+1].key)
j=j+1; /* 为关键码较大的元素下标*/
if(rc.key<h->r[j].key) break; /* rc应插入在位置s上*/
h->r=h->r[j]; s=j; /* 使s结点满足堆定义 */
}
h->r=rc; /* 插入 */
}

void HeapSort(S_TBL *h)
{ for(i=h->length/2;i>0;i--) /* 将r[1..length]建成堆 */
HeapAdjust(h,i,h->length);
for(i=h->length;i>1;i--)
{ h->r[1]<-->h->r[i]; /* 堆顶与堆低元素交换 */
HeapAdjust(h,1,i-1); /*将r[1..i-1]重新调整为堆*/
}
}

次,交换记录至多k次。所以,在建好堆后,排序过程中的筛选次数不超过下式:
+ … +û2)-log2(në + û1)-log2(në( 2  )ûlog22 < nlog2n2
而建堆时的比较次数不超过4n次,因此堆排序最坏情况下,时间复杂度也为O(nlog2n)。
4二路归并排序
二路归并排序的基本操作是将两个有序表合并为一个有序表。
设r[u…t]由两个有序子表r[u…v-1]和r[v…t]组成,两个子表长度分别为v-u、t-v+1。合并方法为:
⑴ i=u;j=v;k=u; //置两个子表的起始下标及辅助数组的起始下标
⑵若i>v 或 j>t,转⑷ //其中一个子表已合并完,比较选取结束
⑶ //选取r[i]和r[j]关键码较小的存入辅助数组rf
如果r[i].key<r[j].key,rf[k]=r[i]; i++; k++;转⑵
否则,rf[k]=r[j]; j++; k++;转⑵
⑷ //将尚未处理完的子表中元素存入rf
如果i<v,将r[i…v-1]存入rf[k…t] //前一子表非空
如果j<=t,将r[i…v]存入rf[k…t] //后一子表非空
⑸合并结束。

【算法10.11】
void Merge(ElemType *r,ElemType *rf,int u,int v,int t)
{
for(i=u,j=v,k=u;i<v&&j<=t;k++)
{ if(r[i].key<r[j].key)
{ rf[k]=r[i];i++;}
else
{ rf[k]=r[j];j++;}
}
if(i<v) rf[k…t]=r[i…v-1];
if(j<=t) rf[k…t]=r[j…t];
}

一.两路归并的迭代算法
1个元素的表总是有序的。所以对n个元素的待排序列,每个元素可看成1个有序子

表长度均为2。再进行两两合并,直到生成n个元素按关键码有序的表。
【算法10.12】
void MergeSort(S_TBL *p,ElemType *rf)
{ /*对*p表归并排序,*rf为与*p表等长的辅助数组*/
ElemType *q1,*q2;
q1=rf;q2=p->elem;
for(len=1;len<p->length;len=2*len) /*从q2归并到q1*/
{ for(i=1;i+2*len-1<=p->length;i=i+2*len)
Merge(q2,q1,i,i+len,i+2*len-1); /*对等长的两个子表合并*/
if(i+len-1<p->length)
Merge(q2,q1,i,i+len,p->length); /*对不等长的两个子表合并*/
else if(i<=p->length)
while(i<=p->length) /*若还剩下一个子表,则直接传入*/
q1[i]=q2[i];
q1<-->q2; /*交换,以保证下一趟归并时,仍从q2归并到q1*/
if(q1!=p->elem) /*若最终结果不在*p表中,则传入之*/
for(i=1;i<=p->length;i++)
p->elem[i]=q1[i];
}
}

二.两路归并的递归算法
【算法10.13】
void MSort(ElemType *p,ElemType *p1,int s,int t)
{ /*将p归并排序为p1*/
if(s==t) p1=p
else
{ m=(s+t)/2; /*平分*p表*/
MSort(p,p2,s,m); /*递归地将p归并为有序的p2*/
MSort(p,p2,m+1,t); /*递归地将p[m+1…t]归并为有序的p2[m+1…t]*/
Merge(p2,p1,s,m+1,t); /*将p2和p2[m+1…t]归并到p1*/
}
}

void MergeSort(S_TBL *p)
{ /*对顺序表*p作归并排序*/
MSort(p->elem,p->elem,1,p->length);
}

【效率分析】
需要一个与表等长的辅助元素数组空间,所以空间复杂度为O(n)。
对n个元素的表,将这n个元素看作叶结点,若将两两归并生成的子表看作它们的父结点,则归并过程对应由叶向根生成一棵二叉树的过程。所以归并趟数约等于二叉树的高度-1,即log2n,每趟归并需移动记录n次,故时间复杂度为O(nlog2n)。

5基数排序
基数排序是一种借助于多关键码排序的思想,是将单关键码按基数分成“多关键码”进行排序的方法。

5.1 多关键码排序
扑克牌中52张牌,可按花色和面值分成两个字段,其大小关系为:
花色:梅花 < 方块 < 红心 < 黑心
面值: 2 < 3 < 4 < 5 < 6 < 7 < 8 < 9 < 10 < J < Q < K < A

若对扑克牌按花色、面值进行升序排序,得到如下序列:
梅花2,3,...,A,方块2,3,...,A,红心2,3,...,A,黑心2,3,...,A
即两张牌,若花色不同,不论面值怎样,花色低的那张牌小于花色高的,只有在同花色情况下,大小关系才由面值的大小确定。这就是多关键码排序。
为得到排序结果,我们讨论两种排序方法。
方法1:先对花色排序,将其分为4个组,即梅花组、方块组、红心组、黑心组。再对每个组分别按面值进行排序,最后,将4个组连接起来即可。
方法2:先按13个面值给出13个编号组(2号,3号,...,A号),将牌按面值依次放入对应的编号组,分成13堆。再按花色给出4个编号组(梅花、方块、红心、黑心),将2号组中牌取出分别放入对应花色组,再将3号组中牌取出分别放入对应花色组,……,这样,4个花色组中均按面值有序,然后,将4个花色组依次连接起来即可。
设n个元素的待排序列包含d个关键码{k1,k2,…,kd},则称序列对关键码{k1,k2,…,kd}有序是指:对于序列中任两个记录r[i]和r[j](1≤i≤j≤n)都满足下列有序关系:

其中k1称为最主位关键码,kd称为最次位关键码。
多关键码排序按照从最主位关键码到最次位关键码或从最次位到最主位关键码的顺序
逐次排序,分两种方法:
最高位优先(Most Significant Digit first)法,简称MSD法:先按k1排序分组,同一组中记录,关键码k1相等,再对各组按k2排序分成子组,之后,对后面的关键码继续这样的排序分组,直到按最次位关键码kd对各子组排序后。再将各组连接起来,便得到一个有序序列。扑克牌按花色、面值排序中介绍的方法一即是MSD法。
最低位优先(Least Significant Digit first)法,简称LSD法:先从kd开始排序,再对kd-1进行排序,依次重复,直到对k1排序后便得到一个有序序列。扑克牌按花色、面值排序中介绍的方法二即是LSD法。

10.6.2链式基数排序
将关键码拆分为若干项,每项作为一个“关键码”,则对单关键码的排序可按多关键码排序方法进行。比如,关键码为4位的整数,可以每位对应一项,拆分成4项;又如,关键码由5个字符组成的字符串,可以每个字符作为一个关键码。由于这样拆分后,每个关键码都在相同的范围内(对数字是0~9,字符是'a'~'z'),称这样的关键码可能出现的符号个数为“基”,记作RADIX。上述取数字为关键码的“基”为10;取字符为关键码的“基”为26。基于这一特性,用LSD法排序较为方便。
基数排序:从最低位关键码起,按关键码的不同值将序列中的记录“分配”到RADIX个队列中,然后再“收集”之。如此重复d次即可。链式基数排序是用RADIX个链队列作为分配队列,关键码相同的记录存入同一个链队列中,收集则是将各链队列按关键码大小顺序链接起来。

【例10.8】以静态链表存储待排记录,头结点指向第一个记录。链式基数排序过程如下图。
图(a):初始记录的静态链表。

图(b):第一趟按个位数分配,修改结点指针域,将链表中的记录分配到相应链队列中。
图(c):第一趟收集:将各队列链接起来,形成单链表。

图(d):第二趟按十位数分配,修改结点指针域,将链表中的记录分配到相应链队列中。
图(e):第二趟收集:将各队列链接起来,形成单链表。

图(f):第三趟按百位数分配,修改结点指针域,将链表中的记录分配到相应链队列中。
图(g):第三趟收集:将各队列链接起来,形成单链表。此时,序列已有序。

图10.10
【算法10.14】
#define MAX_KEY_NUM 8 /*关键码项数最大值*/
#define RADIX 10 /*关键码基数,此时为十进制整数的基数*/
#define MAX_SPACE 1000 /*分配的最大可利用存储空间*/
typedef struct{
KeyType keys[MAX_KEY_NUM]; /*关键码字段*/
InfoType otheritems; /*其它字段*/
int next; /*指针字段*/
}NodeType; /*表结点类型*/
typedef struct{
NodeType r[MAX_SPACE]; /*静态链表,r[0]为头结点*/
int keynum; /*关键码个数*/
int length; /*当前表中记录数*/
}L_TBL; /*链表类型*/
typedef int ArrayPtr[radix]; /*数组指针,分别指向各队列*/

void Distribute(NodeType *s,int i,ArrayPtr *f,ArrayPtr *e)
{ /*静态链表ltbl的r域中记录已按(kye[0],keys[1],…,keys[i-1])有序)*/
/*本算法按第i个关键码keys[i]建立RADIX个子表,使同一子表中的记录的keys[i]相同*/
/*f[0…RADIX-1]和e[0…RADIX-1]分别指向各子表的第一个和最后一个记录*/
for(j=0;j<RADIX;j++) f[j]=0; /* 各子表初始化为空表*/
for(p=r[0].next;p;p=r

.next)
{ j=ord(r
.keys[i]); /*ord将记录中第i个关键码映射到[0…RADIX-1]*/
if(!f[j]) f[j]=p;
else r[e[j]].next=p;
e[j]=p; /* 将p所指的结点插入到第j个子表中*/
}
}

void Collect(NodeType *r,int i,ArrayPtr f,ArrayPtr e)
{/*本算法按keys[i]自小到大地将f[0…RADIX-1]所指各子表依次链接成一个链表*e[0…RADIX-1]为各子表的尾指针*/
for(j=0;!f[j];j=succ(j)); /*找第一个非空子表,succ为求后继函数*/
r[0].next=f[j];t=e[j]; /*r[0].next指向第一个非空子表中第一个结点*/
while(j<RADIX)
{ for(j=succ(j);j<RADIX-1&&!f[j];j=succ(j)); /*找下一个非空子表*/
if(f[j]) {r[t].next=f[j];t=e[j];} /*链接两个非空子表*/
}
r[t].next=0; /*t指向最后一个非空子表中的最后一个结点*/
}
void RadixSort(L_TBL *ltbl)
{ /*对ltbl作基数排序,使其成为按关键码升序的静态链表,ltbl->r[0]为头结点*/
for(i=0;i<ltbl->length;i++) ltbl->r[i].next=i+1;
ltbl->r[ltbl->length].next=0; /*将ltbl改为静态链表*/
for(i=0;i<ltbl->keynum;i++) /*按最低位优先依次对各关键码进行分配和收集*/
{ Distribute(ltbl->r,i,f,e); /*第i趟分配*/
Collect(ltbl->r,i,f,e); /*第i趟收集*/
}
}
【效率分析】
时间效率:设待排序列为n个记录,d个关键码,关键码的取值范围为radix,则进行链式基数排序的时间复杂度为O(d(n+radix)),其中,一趟分配时间复杂度为O(n),一趟收集时间复杂度为O(radix),共进行d趟分配和收集。
空间效率:需要2*radix个指向队列的辅助空间,以及用于静态链表的n个指针。

6外排序
6.1外部排序的方法
外部排序基本上由两个相互独立的阶段组成。首先,按可用内存大小,将外存上含n个记录的文件分成若干长度为k的子文件或段(segment),依次读入内存并利用有效的内部排序方法对它们进行排序,并将排序后得到的有序子文件重新写入外存。通常称这些有序子文件为归并段或顺串;然后,对这些归并段进行逐趟归并,使归并段(有序子文件)逐渐由小到大,直至得到整个有序文件为止。
显然,第一阶段的工作已经讨论过。以下主要讨论第二阶段即归并的过程。先从一个例子来看外排序中的归并是如何进行的?

假设有一个含 10000 个记录的
文件,首先通过10次内部排序得到
10个初始归并段 R1~R10 ,其中每
一段都含1000个记录。然后对它们
作如图10.11所示的两两归并,直至
得到一个有序文件为止。

从图10.11可见,由10个初始归并段到一个有序文件,共进行了四趟归并,每一趟

将两个有序段归并成一个有序段的过程,若在内存中进行,则很简单,前面讨论的2-路归并排序中的Merge函数便可实现此归并。但是,在外部排序中实现两两归并时,不仅要调用Merge函数,而且要进行外存的读/写,这是由于我们不可能将两个有序段及归并结果同时放在内存中的缘故。对外存上信息的读/写是以“物理块”为单位。假设在上例中每个物理块可以容纳200个记录,则每一趟归并需进行50次“读”和50次“写”,四趟归并加上内部排序时所需进行的读/写,使得在外排序中总共需进行500次的读/写。
一般情况下,外部排序所需总时间=
内部排序(产生初始归并段)所需时间 m*tis
+外存信息读写的时间 d*tio
+内部归并排序所需时间 s*utmg
其中:tis是为得到一个初始归并段进行的内部排序所需时间的均值;tio是进行一次外存读/写时间的均值;utmg是对u个记录进行内部归并所需时间;m为经过内部排序之后得到的初始归并段的个数;s为归并的趟数;d为总的读/写次数。由此,上例10000个记录利用2-路归并进行排序所需总的时间为:
10*tis+500*tio+4*10000tmg
其中tio取决于所用的外存设备,显然,tio较tmg要大的多。因此,提高排序效率应主要着眼于减少外存信息读写的次数d。
下面来分析d和“归并过程”的关系。若对上例中所得的10个初始归并段进行5-平衡归并(即每一趟将5个或5个以下的有序子文件归并成一个有序子文件),则从下图可见,仅需进行二趟归并,外部排序时总的读/写次数便减少至2×100+100=300,比2-路归并减少了200次的读/写。
R1 R2 R3 R4 R5 R6 R7 R8 R9 R10
└─┴─┼─┴─┘ └─┴─┼─┴─┘
R1' R2'
└────┬────┘
有序文件
图10.12
可见,对同一文件而言,进行外部排序时所需读/写外存的次数和归并的趟数s成正比。而在一般情况下,对m个初始归并段进行k-路平衡归并时,归并的趟数
可见,若增加k或减少m便能减少s。下面分别就这两个方面讨论之。
10.7.2多路平衡归并的实现
从上式可见,增加k可以减少s,从而减少外存读/写的次数。但是,从下面的讨论中又可发现,单纯增加k将导致增加内部归并的时间utmg。那末,如何解决这个矛盾呢?
先看2-路归并。令u个记录分布在两个归并段上,按Merge函数进行归并。每得到归并后的含u个记录的归并段需进行u-1次比较。
再看k-路归并。令u个记录分布在k个归并段上,显然,归并后的第一个记录应是k个归并段中关键码最小的记录,即应从每个归并段的第一个记录的相互比较中选出最小者,这需要进行k-1次比较。同理,每得到归并后的有序段中的一个记录,都要进行k-1次比较。显然,为得到含u个记录的归并段需进行(u-1)(k-1)次比较。由此,对n个记录的文件进行外部排序时,在内部归并过程中进行的总的比较次数为s(k-1)(n-1)。假设所得初始归并段为m个,则可得内部归并过程中进行比较的总的次数为

k而减少外存信息读写时间所得效益,这是我们所不希望的。然而,若在进行k-路归并时利用“败者树”(Tree of Loser),则可使在k个记录中选出关键码最小的记录时仅需进

它不再随k的增长而增长。
何谓“败者树”?它是树形选择排序的一种变型。相对地,我们可称图10.5和图10.6中二叉树为“胜者树”,因为每个非终端结点均表示其左、右子女结点中“胜者”。反之,若在双亲结点中记下刚进行完的这场比赛中的败者,而让胜者去参加更高一层的比赛,便可得到一棵“败者树”。
【例10.9】

(a) (b)
图10.13 实现5-路归并的败者树

图10.13(a)即为一棵实现5-路归并的败者树ls[0…4],图中方形结点表示叶子结点(也可看成是外结点),分别为5个归并段中当前参加归并的待选择记录的关键码;败者树中根结点ls[1]的双亲结点ls[0]为“冠军”,在此指示各归并段中的最小关键码记录为第三段中的记录;结点ls[3]指示b1和b2两个叶子结点中的败者即是b2,而胜者b1和b3(b3是叶子结点b3、b4和b0经过两场比赛后选出的获胜者)进行比较,结点ls[1]则指示它们中的败者为b1。在选得最小关键码的记录之后,只要修改叶子结点b3中的值,使其为同一归并段中的下一个记录的关键码,然后从该结点向上和双亲结点所指的关键码进行比较,败者留在该双亲,胜者继续向上直至树根的双亲。如图10.13(b)所示。当第3个归并段中第2个记录参加归并时,选得最小关键码记录为第一个归并段中的记录。为了防止在归并过程中某个归并段变为空,可以在每个归并段中附加一个关键码为最大的记录。当选出的“冠军”记录的关键码为最大值时,表明此次归并已完成。由于实现k-路归并的败者树

的初始化也容易实现,只要先令所有的非终端结点指向一个含最小关键码的叶子结点,然后从各叶子结点出发调整非终端结点为新的败者即可。

下面程序中简单描述了利用败者树进行k-路归并的过程,为了突出如何利用败者树进行归并,避开了外存信息存取的细节,可以认为归并段已存在。

【算法10.15】
typedef int LoserTree[k]; /*败者树是完全二叉树且不含叶子,可采用顺序存储结构*/

typedef struct{
KeyType key;
}ExNode,External[k]; /*外结点,只存放待归并记录的关键码*/

void K_Merge(LoserTree *ls,External *b) /*k-路归并处理程序*/
{ /*利用败者树ls将编号从0到k-1的k个输入归并段中的记录归并到输出归并段*/
/*b[0]到b[k-1]为败者树上的k个叶子结点,分别存放k个输入归并段中当前记录的关键码*/
for(i=0;i<k;i++) input(b[i].key); /*分别从k个输入归并段读入该段当前第一个记录的*/
/*关键码到外结点*/
CreateLoserTree(ls); /*建败者树ls,选得最小关键码为b[0].key*/
while(b[ls[0]].key!=MAXKEY)
{ q=ls[0]; /*q指示当前最小关键码所在归并段*/
output(q); /*将编号为q的归并段中当前(关键码为b[q].key的记录写至输出归并段)*/
input(b[q].key); /*从编号为q的输入归并段中读入下一个记录的关键码*/
Adjust(ls,q); /*调整败者树,选择新的最小关键码*/
}
output(ls[0]); /*将含最大关键码MAXKEY的记录写至输出归并段*/
}

void Adjust(LoserTree *ls,int s) /*选得最小关键码记录后,从叶到根调整败者树,选下一个最小关键码*/
{ /*沿从叶子结点b到根结点ls[0]的路径调整败者树*/
t=(s+k)/2; /*ls[t]是b的双亲结点*/
while(t>0)
{ if(b.key>b[ls[t]].key) s<-->ls[t]; /*s指示新的胜者*/
t=t/2;
}
ls[0]=s;
}

void CreateLoserTree(LoserTree *ls) /*建立败者树*/
{ /*已知b[0]到b[k-1]为完全二叉树ls的叶子结点存有k个关键码,沿从叶子到根的k条路径*/
/*将ls调整为败者树*/
b[k].key=MINKEY; /*设MINKEY为关键码可能的最小值*/
for(i=0;i<k;i++) ls[i]=k; /*设置ls中“败者”的初值*/
for(i=k-1;k>0;i--) Adjust(ls,i); /*依次从b[k-1],b[k-2],…,b[0]出发调整败者*/
}

最后要提及一点,k值的选择并非越大越好,如何选择合适的k是一个需要综合考虑的问题。

给我留言

Copyright © 浩然东方 保留所有权利.   Theme  Ality 07032740

用户登录