C语言标准库 qsort bsearch 源码实现
C语言是简洁的强大的,当然也有很多坑。C语言也是有点业界良心的,至少它实现了2个最最常用的算法:快速排序和二分查找。
我们知道,对于C语言标准库 qsort和 bsearch:
a. 它是“泛型”的,可以对任何类型进行排序或二分。
b. 我们使用时必须自定义一个比较函数当作函数指针传入。
c语言要实现泛型,基本上就只有 void指针提供的弱爆了的泛型机制,容易出错。
这篇文章中,我实现了 标准库qsort和bsearch函数,最基本的正确性和泛型当然要保证了。
在这里,不涉及优化(写标准库实现的那帮人恨不得用汇编实现),只展现算法的运行原理和泛型的实现机制。
1.C语言标准库qsort源码实现。我先呈上完整实现,然后具体剖析。
#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <string.h> void swap(const void* a, const void* b, int size)
{
assert(a != NULL && b != NULL);
char tmp = ;
int i = ;
while (size > ) {
tmp = *((char*)a + i);
*((char*)a + i) = *((char*)b + i);
*((char*)b + i) = tmp;
++i;
--size;
}
} void Qsort(void* base, int left, int right, int size, int (*cmp)(const void* a, const void* b))
{
assert(base != NULL && size >= && cmp != NULL); /* left may be < 0 because of the last - 1 */
if (left >= right) return;
char* pleft = (char*)base + left * size;
char* pkey = (char*)base + (left + (right - left) / ) * size;
swap(pleft, pkey, size);
int last = left;
char* plast = (char*)base + last * size;
for (int i = left + ; i <= right; ++i) {
char* pi = (char*)base + i * size;
if (cmp(pi, pleft) < ) {
++last;
plast = (char*)base + last * size;
swap(pi, plast, size);
}
}
swap(pleft, plast, size);
Qsort(base, left, last - , size, cmp);
Qsort(base, last + , right, size, cmp);
} int cmp_string(const void* a, const void* b)
{
assert(a != NULL && b != NULL);
const char** lhs = (const char**)a;
const char** rhs = (const char**)b;
return strcmp(*lhs, *rhs);
} int cmp_int(const void* a, const void* b)
{
assert(a != NULL && b != NULL);
const int* lhs = (const int*)a;
const int* rhs = (const int*)b;
if (*lhs < *rhs) {
return -;
} else if (*lhs == *rhs) {
return ;
} else {
return ;
}
} int main(int argc, char* argv[])
{
int a[] = {-, , , , , , , , , };
int len1 = sizeof(a) / sizeof(a[]);
fprintf(stdout, "before sort:\n");
for (int i = ; i < len1; ++i) {
fprintf(stdout, "%d ", a[i]);
}
fprintf(stdout, "\n"); Qsort(a, , len1 - , sizeof(a[]), cmp_int);
fprintf(stdout, "after sort:\n");
for (int i = ; i < len1; ++i) {
fprintf(stdout, "%d ", a[i]);
}
fprintf(stdout, "\n"); const char* b[] = {"what", "chenwei", "skyline", "wel", "dmr"};
int len2 = sizeof(b) / sizeof(b[]);
fprintf(stdout, "before sort:\n");
for (int i = ; i < len2; ++i) {
fprintf(stdout, "%s-->", b[i]);
}
fprintf(stdout, "\n"); Qsort(b, , len2 - , sizeof(b[]), cmp_string);
fprintf(stdout, "after sort:\n");
for (int i = ; i < len2; ++i) {
fprintf(stdout, "%s-->", b[i]);
}
fprintf(stdout, "\n");
}
qsort基本实现在K&R里有非常详细的描述,我这里重点解释的是 swap函数,这个函数是实现泛型的基石。
首先对qsort函数原型说明几点:
a.Qsort函数原型里面没有标准库qsort的 len, 而是使用 left 和 right 来指明每次待排序的区间,用两个值来表示一个区间非常直观且优雅。
b.对于数据类型大小,我没有使用 size_t 无符号类型。
C语言中无符号类型虽然可以对数组提供负向越界保证和2倍空间,但是由于坑爹的类型提升规则滋生了N多的bug,我是尽量少用这个。
然后就是 swap函数.
首先我们知道 数据元素都是以 比特位 形式暂存在 内存中的,后简化为以 字节 形式暂存在内存中。
我们平常的一个 swap 函数,如果交换的数据是 int 型, 我们就是:
void swap(int* a , int* b) {
int tmp = *a;
*a = *b;
*b = tmp;
}
现在我们不知道元素是什么类型,那我们怎么做到泛型呢? 首先 传入的两个指针肯定是 void 指针,如何确定 void所指元素类型呢? 其实我们这里是不需要知道 元素类型的。
我们这里的目的是什么? 我们就是要 交换两个元素的 字节序列,我们每次交换一个字节,直到交换完为止,这样就达到目的了,这时我们就需要第三个参数:待排元素的所占字节数。
我们知道C语言中结构体是可以直接复制的,比如:
struct test {
int a;
char b;
};
struct test a, b;
a = b;
C语言是怎样支持这种直接复制呢? 其实最终就是 两个元素的字节序列的复制。注意到,如果struct里面有指针,那么这个结构体是引用语义的,两个元素的指针成员指向同一内存。
我们这里的swap函数其实就是 手工模拟了 类似结构体复制 的整个过程。我们每次交换一个 字节,总共交换的次数就是 元素的sizeof大小。
注意两点:
a. 数据的字节序列是有 大端小端之分的,我的这个实现是不能 跨大端小端的。如果排序过程都始终都在同一机器中,那么无需担心。
b. 数据的字节序列可能是有 内存对齐的,主要是在结构体之中。 所以我的这个实现受 机器 内存对齐规则的影响,但是 排序发生在同一机器中的话,这个是不会影响正确性的。
总而言之,代码实现受机器底层 数据大小端表示和 内存对齐策略 的影响。但是实际上 一个排序过程是不可能 一部分进行在这台机,另一部分在另外一台机器上进行的(仅考虑单机),所以我的这个实现是可以很好工作的。
Qsort的快排思想就很简单了,我们最需要注意的就是 每次对 交换元素的首字节地址进行更新,我们都是经由char*转换,因为char*所指元素正好1字节,正好模拟每次一字节的swap.
2.C语言标准库bsearch源码实现
void* Bsearch(void* base, int len, int size, const void* key, int (*cmp)(const void* a, const void* b))
{
assert(base != NULL && len >= && size >= && cmp != NULL);
int low = ;
int high = len - ;
while (low <= high) {
int mid = low + (high - low) / ;
char* pmid = (char*)base + mid * size;
if (cmp(pmid, key) < ) {
low = mid + ;
} else if (cmp(pmid, key) > ) {
high = mid - ;
} else {
return pmid;
} }
return NULL;
} int cmp_int(const void* a, const void* b)
{
assert(a != NULL && b != NULL);
const int* lhs = (const int*)a;
const int* rhs = (const int*)b;
if (*lhs < *rhs) {
return -;
} else if (*lhs == *rhs) {
return ;
} else {
return ;
}
} int cmp_string(const void* a, const void* b)
{
assert(a != NULL && b != NULL);
const char** lhs = (const char**)a;
const char** rhs = (const char**)b;
return strcmp(*lhs, *rhs);
} int main(int argc, char* argv[])
{
int a[] = {-, , , , , , , , , };
int len1 = sizeof(a) / sizeof(a[]); int tmp = ;
int* res1 = (int*)Bsearch(a, len1, sizeof(a[]), &tmp, cmp_int);
if (res1 != NULL) {
fprintf(stdout, "found it\n");
} else {
fprintf(stdout, "Not found\n");
} const char* str[] = {"chenwei", "dmr", "skyline", "wel", "what"};
int len2 = sizeof(str) / sizeof(str[]);
const char* p = "chenwei";
char* res2 = (char*)Bsearch(str, len2, sizeof(str[]), &p, cmp_string);
if (res2 != NULL) {
fprintf(stdout, "found it\n");
} else {
fprintf(stdout, "Not found\n");
}
return ;
}
欢迎大家批评指正,共同学习。转载请注明出处,谢谢。
C语言标准库 qsort bsearch 源码实现的更多相关文章
- python 协程库gevent学习--源码学习(一)
总算还是要来梳理一下这几天深入研究之后学习到的东西了. 这几天一直在看以前跟jd对接的项目写的那个gevent代码.为了查错,基本上深入浅出了一次gevent几个重要部件的实现和其工作的原理. 这里用 ...
- python语言线程标准库threading.local源码解读
本段源码可以学习的地方: 1. 考虑到效率问题,可以通过上下文的机制,在属性被访问的时候临时构建: 2. 可以重写一些魔术方法,比如 __new__ 方法,在调用 object.__new__(cls ...
- requests库核心API源码分析
requests库是python爬虫使用频率最高的库,在网络请求中发挥着重要的作用,这边文章浅析requests的API源码. 该库文件结构如图: 提供的核心接口在__init__文件中,如下: fr ...
- 第三方库Mantle的源码解析
Mantle是一个用于简化Cocoa或Cocoa Touch程序中model层的第三方库.通常我们的应该中都会定义大量的model来表示各种数据结构,而这些model的初始化和编码解码都需要写大量的代 ...
- 2018-08-27 使用JDT核心库解析JDK源码后初步分析API命名
源自术语词典API项目 · Issue #85 · program-in-chinese/overview, 打算先用早先的代码提取JDK API中的类/方法/参数名, 看看有哪些词需要翻译. 源码在 ...
- AndroidStudio导入第三方开源库 --文件夹源码
1 在已打开的项目中 File-New-ImportModule 选择开源项目中的 库所在文件夹比如 library文件夹 然后导入. 2 File-Project Sructure 在Modu ...
- qsort.c源码
/* 版权所有(C) 1991-2019 自由软件资金会. 该文件属于是GUN C语言函数库,由Douglas C. Schmidt(schmidt@ics.uci.edu)所写. GUN C语言函数 ...
- Mysql依赖库Boost的源码安装,linux下boost库的安装
boost‘准标准库’安装过程.安装的是boost_1_60_0. (1)首先去下载最新的boost代码包,网址www.boost.org. (2)进入到自己的目录,解压: bzip2 -d bo ...
- 如何快糙好猛的使用Shiqi.Yu老师的公开人脸检测库(附源码)
前言 本次编写所用的库为于仕祺老师免费提供的人脸检测库.真心好用,识别率和识别速度完全不是Opencv自带的程序能够比拟的.将其配合Opencv的EigenFace算法,基本上可以形成一个小型的毕业设 ...
随机推荐
- Mysql查询指定用户并列排名 类似rank函数
SELECT total.* FROM ( SELECT obj.uid, obj.score, CASE WHEN @rowtotal = obj.score THEN @rownum WHEN @ ...
- python-下拉框处理
在自动化中python对下拉框的处理网上相对实例比较少,其它前辈写的教程中对下拉也仅仅是相对与教程来说的,比如下面: m=driver.find_element_by_id("Shippin ...
- Bootstrap历练实例:标签修饰
您可以使用修饰的 class label-default.label-primary.label-success.label-info.label-warning.label-danger 来改变标签 ...
- OI算法复习汇总
各大排序 图论: spfa floyd dijkstra *拉普拉斯矩阵 hash表 拓扑排序 哈夫曼算法 匈牙利算法 分块法 二分法 费马小定理: a^(p-1) ≡1(mod p) 网络流 二分图 ...
- centos7.4系统部署nodejs前端项目
1.安装nodejs运行环境 wget命令下载Node.js安装包,该安装包是编译好的文件,解压之后,在bin文件夹中就已存在node和npm,无需重复编译 wget https://nodejs.o ...
- Win2008 Server配置PHP环境
Win2008 Server配置PHP环境 阅读目录 创建一个网站 配置PHP环境 配置iis的“处理应用程序映射” 在配置PHP环境之前要先配置好IIS. 传送门-> Win2008 Se ...
- 【php】expose_php 作用
php.ini expose_php = On expose_php = Off
- Java求字符串中出现次数最多的字符
Java求字符串中出现次数最多的字符 [尊重原创,转载请注明出处]http://blog.csdn.net/guyuealian/article/details/51933611 Java ...
- 小谈python里 列表 的几种常用用法
在python中列表的常用方法主要包括增加,删除,查看和修改.下面以举例子的方法具体说明,首先我们创建两个列表,列表是用[ ]表示的,里面的元素用逗号隔开. a=[‘hello’,78,15.6,‘你 ...
- 安装openstack同步数据库时出错解决方法
错误提示:(2003, "Can't connect to MySQL server on 'controller' ([Errno -2] Name or service not know ...