【混合排序】（STL std::sort | 快速排序 | 插入排序 | 堆排序）

XNB's Not a Beginner

于 2023-12-19 14:09:42 发布

阅读量428

点赞数 9

分类专栏：算法语言特性 Modern Cpp 文章标签： c++ 算法数据结构

本文链接：https://blog.csdn.net/RFoer/article/details/135083452

版权

算法同时被 2 个专栏收录

83 篇文章 0 订阅

订阅专栏

语言特性 Modern Cpp

38 篇文章 0 订阅

订阅专栏

文章介绍了C++STL中的std::sort函数实现，尤其是IntroSort算法，它结合了median-of-threeQuickSort策略，能自动切换到HeapSort以应对性能下降。文中还提及了pivot选择和内存管理的细节。

摘要由CSDN通过智能技术生成

/*
	std::sort 的实现基于 David Musser 1996年提出的一种混合式排序算法 Introspective Sorting，
即IntroSort，其行为几乎与median-of-three QuickSort完全相同，
但是当Partition有恶化的倾向时，能够自我侦测，转而改用HeapSort，
使复杂度维持在HeapSort的 O(NlonN)。
—— 侯捷《STL源码剖析》
std::sort在最后，会在整个数组上做一次插入排序。
ntrosort在GCC源码中涉及了两处 以__unguarded__...... 为名的函数，
它们的名字暗示了它们不会做越界检查。

在快排中 pivot的选择是一件non-trivial的事情，
introsort选用了median-of-three方法，
这又与一个 __unguarded 函数产生了联动。
*/

#define _THRESHOLD_VALUE_ 16 
#define regUInt \
	register unsigned int
void QuickSort(int A[],int beg,int end,int depMon,int N)
{
	if(beg+1+_THRESHOLD_VALUE_ >= end){
		BinaryInsertSort(A,beg,end);//改用二分插入排序插入 
		return;
	}
	if(depMon > log2(N)){
		HeapSort(A,beg,end);
		//递归深度大于log2N阈值 改成堆排序 
		return;
	}
	int mid = beg + ( (end-beg)>>1 ),
		rbeg = end-1;
	if(A[beg] > A[mid])swap(&A[beg],&A[mid]);
	if(A[rbeg] > A[mid])swap(&A[rbeg],&A[mid]);
	if(A[beg] > A[rbeg])swap(&A[beg],&A[rbeg]);
	//三者取中 后置 
	regUInt i = beg,
			_Pivot = rbeg,
		    j = _Pivot-1,
		    k,r;
	while(i < j)
	{
		while((i < j) && A[j] > A[_Pivot])--j;
		while((i < j) && A[i] <= A[_Pivot])++i;
		if(i < j)swap(&A[i],&A[j]);
	}	
	swap(&A[_Pivot],&A[!(i == j)?j:++j]);
	swap(&_Pivot,&j);
	r = _Pivot-1;
	for(r = _Pivot,k = beg;k < _Pivot;++k)
		if(A[k] == A[_Pivot])swap(A[k],A[--r]);
	//元素聚集 
	QuickSort(A,beg,r,depMon+1);
	QuickSort(A,_Pivot+1,end,depMon+1);
}