深圳幻海软件技术有限公司 欢迎您!

一篇带给你索引技术之位图

2023-02-28

要点位图基本算法及其应用场景。位图算法的优化实现。概述位图算法,是指使用一个bit位来表示数据状态。通常应用于海量数据去重、海量数据计算及判断海量数据中是否存在某个数据的场景中。以海量数据中是否存在某个数据的应用场景为例,假设用16个bit位,分别表示数字0-15。bit位的值,表示该数字是否存在,

要点

  • 位图基本算法及其应用场景。
  • 位图算法的优化实现。

概述

位图算法,是指使用一个bit位来表示数据状态。通常应用于海量数据去重、海量数据计算及判断海量数据中是否存在某个数据的场景中。

以海量数据中是否存在某个数据的应用场景为例,假设用16个bit位,分别表示数字0-15。bit位的值,表示该数字是否存在,0表示不存在,1表示存在。如上图所示,在该数据集合中,存在的元素有1、2、6、10、11和13。

可以发现,在数据比较稠密的情况下,位图算法能够节约存储空间,如图中,使用2个字节便可以表示16个数字,同时可以在O(1)的时间复杂度下,判断是否存在某个数字,大大提高了计算速度。

但是,在数据稀疏时,存储空间会存在一定程度的浪费。由于位图算法中,位图空间的大小是一定的,并不会根据存储数据量的大小而改变。因此,当位图空间中存储的数据量很小时,大量地位图空间是空闲的,存在大量的浪费。

算法实现

位图算法在主流开发语言中,都有对应的实现。基本操作主要有写入、查询、删除、交集、并集等。下面通过一个示例来了解一下,位图算法的实现。

位图结构定义例子使用char类型数组来存储位图信息(通常的实现中,会使用长整型数组),一个char类型有8个bit位。定义结构如下:

// 为了简化问题,LEN必须定义为CHAR_SIZE的倍数
#define LEN 16
#define CHAR_SIZE 8
typedef char BitSet[LEN/CHAR_SIZE];
  • 1.
  • 2.
  • 3.
  • 4.

写入在某个bit位写入数据时,首先通过整除,计算出该bit位在数组的哪个下标,然后,用取余计算出char元素中的哪个bit上。最后通过或运算将对应位设置为1。

// 置bit位
void set(BitSet& bits, int pos) {
    // 查找对应数组下标
    int unit = pos / CHAR_SIZE;
    // 查找在字节中的bit位
    int p = pos % CHAR_SIZE;
    // 通过与运算实现对应bit位置1
    bits[unit] = bits[unit] | (0x1 << p);
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.
  • 7.
  • 8.

查询同写入操作,先计算出bit位置,查找到对应的bit位,然后返回该位置的数值。

// 查询bit位
int get(BitSet& bits, int pos) {
    // 查找对应数组下标
    int unit = pos / CHAR_SIZE;
    // 查找在字节中的bit位
    int p = pos % CHAR_SIZE;
    // 通过与运算实现对应bit位置1
    return  (bits[unit] & (0x1 << p)) > 0 ? 1 : 0;
}
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.
  • 7.
  • 8.
  • 9.

删除首先查找到对应的位置,然后通过与运算将该位置清空。

// 清空bit位
void clear(BitSet& bits, int pos) {
    // 查找对应数组下标
    int unit = pos / CHAR_SIZE;
    // 查找在字节中的bit位
    int p = pos % CHAR_SIZE;
    // 通过与运算实现对应bit位置1
    bits[unit] = bits[unit] & (~(0x1 << p));
}
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.
  • 7.
  • 8.
  • 9.

交集对数组逐个元素进行或运算。

// 求位图b1和b2的并集
void unionn(const BitSet& b1, const BitSet& b2, BitSet& res) {
    for (auto i = 0; i < (LEN/CHAR_SIZE); ++i) {
        res[i] = b1[i] | b2[i];
    }
}
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.

并集对数组逐元素进行与运算。

// 求位图b1和b2的交集
void inter(const BitSet& b1, const BitSet& b2, BitSet& res) {
    for (auto i = 0; i < (LEN/CHAR_SIZE); ++i) {
        res[i] = b1[i] & b2[i];
    }
}
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.

在生产实现时,可能会进行一些优化:

  • 使用CPU指令优化,如SSE等,一次能进行128位的运算,可以提高计算速度。
  • 某些业务场景下,一个数据状态可能有大于2个,可以使用多个bit位来表示一个数据状态。

扩展

为了解决位图稀疏时,位图低效率的问题,工业界,有多种位图压缩算法,其中,最经典的是RoaringBitmap。

RoaringBitmap的核心思想是,将整数进行分桶,高16位的值作为其桶的索引,每个桶对应一个容器。如下图所示:

roaring bitmap

容器的结构有三种类型:有序数组、未压缩位图、和行程长度编码。

  • 有序数组:当低16位中,元素个数小于4096时,采用有序数组的结构进行存储。在查找元素时,使用二分查找方法。取值4096的原因是,存储4096个16位的整数,所占用的存储空间:
  • 未压缩位图:未压缩位图的存储结果就是本文所描述的位图存储结构,使用一个固定的连续内存块实现。
  • 行程长度编码(run-length encoding):行程长度编码是一种无损数据压缩技术,其原理是,将连续出现的数据存储为起始值和计算两部分。比如,数据列表[1,2,3,4,5,6]存储为[1,5],表示以1开始,后面连续递增5个数值。在很多实现中,行程长度编码容器,需要手动调用,才能转换为该容器。

在进行插入和删除操作之后,需要根据元素个数进行容器转换。插入元素时,若元素个数达到4096,则需要转换为未压缩位图进行存储。删除元素时,若元素个数小于4096时,则需要转换为有序数组存储。

参考

  • Better bitmap performance with Roaring bitmaps。
  • Consistently faster and smaller compressed bitmaps with Roaring。
  • https://github.com/RoaringBitmap/CRoaring.git。