目录
有序向量:Fibonacci查找(仅使用加法减法实现二分查找)
接口与实现
ADT接口
操作实例
Vector模板类
#pragma once
typedef int Rank; //秩
#define DEFAULT_CAPACITY 3 //默认的初始容量(实际应用中可设置为更大)
template <typename T> class Vector { //向量模板类
protected:
Rank _size; int _capacity; T* _elem; //规模、容量、数据区
void copyFrom ( T const* A, Rank lo, Rank hi ); //复制数组区间A[lo, hi)
void expand(); //空间不足时扩容
void shrink(); //装填因子过小时压缩
bool bubble ( Rank lo, Rank hi ); //扫描交换
void bubbleSort ( Rank lo, Rank hi ); //起泡排序算法
Rank max ( Rank lo, Rank hi ); //选取最大元素
void selectionSort ( Rank lo, Rank hi ); //选择排序算法
void merge ( Rank lo, Rank mi, Rank hi ); //归并算法
void mergeSort ( Rank lo, Rank hi ); //归并排序算法
Rank partition ( Rank lo, Rank hi ); //轴点构造算法
void quickSort ( Rank lo, Rank hi ); //快速排序算法
void heapSort ( Rank lo, Rank hi ); //堆排序(稍后结合完全堆讲解)
public:
// 构造函数
Vector ( int c = DEFAULT_CAPACITY, int s = 0, T v = 0 ) //容量为c、规模为s、所有元素初始为v
{ _elem = new T[_capacity = c]; for ( _size = 0; _size < s; _elem[_size++] = v ); } //s<=c
Vector ( T const* A, Rank n ) { copyFrom ( A, 0, n ); } //数组整体复制
Vector ( T const* A, Rank lo, Rank hi ) { copyFrom ( A, lo, hi ); } //区间
Vector ( Vector<T> const& V ) { copyFrom ( V._elem, 0, V._size ); } //向量整体复制
Vector ( Vector<T> const& V, Rank lo, Rank hi ) { copyFrom ( V._elem, lo, hi ); } //区间
// 析构函数
~Vector() { delete [] _elem; } //释放内部空间
// 只读访问接口
Rank size() const { return _size; } //规模
bool empty() const { return !_size; } //判空
int disordered() const; //判断向量是否已排序
Rank find ( T const& e ) const { return find ( e, 0, _size ); } //无序向量整体查找
Rank find ( T const& e, Rank lo, Rank hi ) const; //无序向量区间查找
Rank search ( T const& e ) const //有序向量整体查找
{ return ( 0 >= _size ) ? -1 : search ( e, 0, _size ); }
Rank search ( T const& e, Rank lo, Rank hi ) const; //有序向量区间查找
// 可写访问接口
T& operator[] ( Rank r ) const; //重载下标操作符,可以类似于数组形式引用各元素
Vector<T> & operator= ( Vector<T> const& ); //重载赋值操作符,以便直接克隆向量
T remove ( Rank r ); //删除秩为r的元素
int remove ( Rank lo, Rank hi ); //删除秩在区间[lo, hi)之内的元素
Rank insert ( Rank r, T const& e ); //插入元素
Rank insert ( T const& e ) { return insert ( _size, e ); } //默认作为末元素插入
void sort ( Rank lo, Rank hi ); //对[lo, hi)排序
void sort() { sort ( 0, _size ); } //整体排序
void unsort ( Rank lo, Rank hi ); //对[lo, hi)置乱
void unsort() { unsort ( 0, _size ); } //整体置乱
int deduplicate(); //无序去重
int uniquify(); //有序去重
// 遍历
void traverse ( void (* ) ( T& ) ); //遍历(使用函数指针,只读或局部性修改)
template <typename VST> void traverse ( VST& ); //遍历(使用函数对象,可全局性修改)
}; //Vector
#include "vector_implementation.h"
可扩充向量
动态空间管理
扩容算法expand()
template <typename T> void Vector<T>::expand() { //向量空间不足时扩容 if ( _size < _capacity ) return; //尚未满员时,不必扩容 if ( _capacity < DEFAULT_CAPACITY ) _capacity = DEFAULT_CAPACITY; //不低于最小容量 T* oldElem = _elem; _elem = new T[_capacity <<= 1]; //容量加倍 for ( int i = 0; i < _size; i++ ) _elem[i] = oldElem[i]; //复制原向量内容(T为基本类型,或已重载赋值操作符'=') /*DSA*/ //printf("\n_ELEM [%x]*%d/%d expanded and shift to [%x]*%d/%d\n", oldElem, _size, _capacity/2, _elem, _size, _capacity); delete [] oldElem; //释放原空间 }
无序向量
有序向量:唯一化
唯一化,uniquify( )
逆序对是指相邻元素之间前者大于后者,即逆序对。
低效版 O(n^2)
高效版O(n):希望能成批的删除,而不是每次remove一个。
有序向量:二分查找
接口
template <typename T> //在有序向量的区间[lo, hi)内,确定不大于e的最后一个节点的秩 Rank Vector<T>::search ( T const& e, Rank lo, Rank hi ) const { //assert: 0 <= lo < hi <= _size return ( rand() % 2 ) ? //按各50%的概率随机使用二分查找或Fibonacci查找 binSearch ( _elem, e, lo, hi ) : fibSearch ( _elem, e, lo, hi ); }
版本A:实现
有序向量:Fibonacci查找(仅使用加法减法实现二分查找)
思路和原理
实现
有序向量:二分查找(改进)
版本B:改进思路
版本B:实现
// 二分查找算法(版本B):在有序向量的区间[lo, hi)内查找元素e,0 <= lo <= hi <= _size template <typename T> static Rank binSearch ( T* A, T const& e, Rank lo, Rank hi ) { /*DSA*/printf ( "BIN search (B)\n" ); while ( 1 < hi - lo ) { //每步迭代仅需做一次比较判断,有两个分支;成功查找不能提前终止 /*DSA*/ for ( int i = 0; i < lo; i++ ) printf ( " " ); if ( lo >= 0 ) for ( int i = lo; i < hi; i++ ) printf ( "....^" ); printf ( "\n" ); Rank mi = ( lo + hi ) >> 1; //以中点为轴点 ( e < A[mi] ) ? hi = mi : lo = mi; //经比较后确定深入[lo, mi)或[mi, hi) } //出口时hi = lo + 1,查找区间仅含一个元素A[lo] /*DSA*/ for ( int i = 0; i < lo; i++ ) printf ( " " ); if ( lo >= 0 ) printf ( "....|\n" ); else printf ( "<<<<|\n" ); return ( e == A[lo] ) ? lo : -1 ; //查找成功时返回对应的秩;否则统一返回-1 } //有多个命中元素时,不能保证返回秩最大者;查找失败时,简单地返回-1,而不能指示失败的位置
版本C:实现
// 二分查找算法(版本C):在有序向量的区间[lo, hi)内查找元素e,0 <= lo <= hi <= _size template <typename T> static Rank binSearch ( T* A, T const& e, Rank lo, Rank hi ) { /*DSA*/printf ( "BIN search (C)\n" ); while ( lo < hi ) { //每步迭代仅需做一次比较判断,有两个分支 /*DSA*/ for ( int i = 0; i < lo; i++ ) printf ( " " ); if ( lo >= 0 ) for ( int i = lo; i < hi; i++ ) printf ( "....^" ); printf ( "\n" ); Rank mi = ( lo + hi ) >> 1; //以中点为轴点 ( e < A[mi] ) ? hi = mi : lo = mi + 1; //经比较后确定深入[lo, mi)或(mi, hi) } //成功查找不能提前终止 /*DSA*/ for ( int i = 0; i < lo - 1; i++ ) printf ( " " ); if ( lo > 0 ) printf ( "....|\n" ); else printf ( "<<<<|\n" ); return --lo; //循环结束时,lo为大于e的元素的最小秩,故lo - 1即不大于e的元素的最大秩 } //有多个命中元素时,总能保证返回秩最大者;查找失败时,能够返回失败的位置
有序向量:插值查找 (按比例查找法)
原理:
实例
性能
时间复杂度:平均情况O(loglog(n)),最坏O(log(n))
实现
int Insert_Search(int *a, int n, int key) { int low, high, mid; low = 0; high = n - 1; while (low <= high) { mid = low + (key - a[low]) / (a[high] - a[low])*(high - low); if (a[mid] < key) low = mid + 1; else if (a[mid]>key) high = mid - 1; else return mid; } return -1; }
起泡排序
排序器:统一接口
template <typename T> void Vector<T>::sort ( Rank lo, Rank hi ) { //向量区间[lo, hi)排序 switch ( rand() % 5 ) { //随机选取排序算法。可根据具体问题的特点灵活选取或扩充 case 1: bubbleSort ( lo, hi ); break; //起泡排序 case 2: selectionSort ( lo, hi ); break; //选择排序(习题) case 3: mergeSort ( lo, hi ); break; //归并排序 case 4: heapSort ( lo, hi ); break; //堆排序(稍后介绍) default: quickSort ( lo, hi ); break; //快速排序(稍后介绍) } /*DSA*/ //selectionSort(lo, hi); }
实现
template <typename T> //向量的起泡排序 void Vector<T>::bubbleSort ( Rank lo, Rank hi ) //assert: 0 <= lo < hi <= size { while ( !bubble ( lo, hi-- ) ); } //逐趟做扫描交换,直至全序 template <typename T> bool Vector<T>::bubble ( Rank lo, Rank hi ) { //一趟扫描交换 bool sorted = true; //整体有序标志 while ( ++lo < hi ) //自左向右,逐一检查各对相邻元素 if ( _elem[lo - 1] > _elem[lo] ) { //若逆序,则 sorted = false; //意味着尚未整体有序,并需要 swap ( _elem[lo - 1], _elem[lo] ); //通过交换使局部有序 } return sorted; //返回有序标志 }
改进
template <typename T> Rank Vector<T>::bubble ( Rank lo, Rank hi ) { //一趟扫描交换 Rank last = lo; //最右侧的逆序对初始化为[lo - 1, lo] while ( ++lo < hi ) //自左向右,逐一检查各对相邻元素 if ( _elem[lo - 1] > _elem[lo] ) { //若逆序,则 last = lo; //更新最右侧逆序对位置记录,并 swap ( _elem[lo - 1], _elem[lo] ); //通过交换使局部有序 } return last; //返回最右侧的逆序对位置 }
评价
归并排序
原理
实现
template <typename T> //向量归并排序 void Vector<T>::mergeSort ( Rank lo, Rank hi ) { //0 <= lo < hi <= size /*DSA*/printf ( "\tMERGEsort [%3d, %3d)\n", lo , hi ); if ( hi - lo < 2 ) return; //单元素区间自然有序,否则... int mi = ( lo + hi ) / 2; //以中点为界 mergeSort ( lo, mi ); mergeSort ( mi, hi ); //分别排序 merge ( lo, mi, hi ); //归并 }
二路归并:原理
实现
template <typename T> //有序向量的归并 void Vector<T>::merge ( Rank lo, Rank mi, Rank hi ) { //各自有序的子向量[lo, mi)和[mi, hi) T* A = _elem + lo; //合并后的向量A[0, hi - lo) = _elem[lo, hi) int lb = mi - lo; T* B = new T[lb]; //前子向量B[0, lb) = _elem[lo, mi) for ( Rank i = 0; i < lb; B[i] = A[i++] ); //复制前子向量 int lc = hi - mi; T* C = _elem + mi; //后子向量C[0, lc) = _elem[mi, hi) for ( Rank i = 0, j = 0, k = 0; ( j < lb ) || ( k < lc ); ) { //B[j]和C[k]中的小者续至A末尾 if ( ( j < lb ) && ( ! ( k < lc ) || ( B[j] <= C[k] ) ) ) A[i++] = B[j++]; if ( ( k < lc ) && ( ! ( j < lb ) || ( C[k] < B[j] ) ) ) A[i++] = C[k++]; } delete [] B; //释放临时空间B } //归并后得到完整的有序向量[lo, hi)
精简实现
/*DSA*/ //原基本实现易于理解、可移植性强,但针对向量归并排序的情况,可优化如下 for ( Rank i = 0, j = 0, k = 0; j < lb; ) { //将B[j]和C[k]中的小者续至A末尾 if ( ( k < lc ) && ( C[k] < B[j] ) ) A[i++] = C[k++]; if ( ( lc <= k ) || ( B[j] <= C[k] ) ) A[i++] = B[j++]; } /*DSA*/ //交换循环体内两句的次序,删除冗余逻辑
评价