数据存在?-‘布隆过滤器’
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了数据存在?-‘布隆过滤器’相关的知识,希望对你有一定的参考价值。
布隆过滤器是一种能够在大量数据中判断数据是否存在的算法。它实际上是一个很长的二进制向量和一系列随机映射函数。布隆过滤器可以用于检索一个元素是否在一个集合中。它的优点是空间效率和查询时间都远远超过一般的算法,缺点是有一定的误识别率和删除困难。在介绍‘布隆过滤器’之前,先介绍一下‘位图’的思想:
这里有这样一个问题:给40亿个没有排序、不重复的无符号整数,如何快速的判断一个数据是否在这40亿个数据之中?
--对于40亿个数据,如果我们将这40亿个数据都放入内存中,我们需要多大的存储空间呢?假设每个数据都是char类型的,这样消耗多少的空间?如果是int类型呢?或者是更多的数据呢?不难知道对于大量的数据,如果采用将数据放入内存中,这种方式是很不理智的。这里介绍一种方法—‘位图’。
位图:主要算法思想就是充分的利用bit位,假设数据都是int类型,每个int类型都占32个bit位。将一个int类型数据的32个bit用来表示32个数据是否存在, 0表示不存在,1表示存在(能够极大地缩小空间)。先计算出数据在哪一个int类型的空间中,然后计算在这个int类型的第几个bit位上,然后将此位置更改为1,表明这个位置上存在数据。
下面是‘位图’的实现:
class BitMap { public: BitMap(size_t size = 0) //构造 :_size(0) { _a.resize((size >> 5) + 1); //resize开辟空间(int类型的个数),并进行初始化 } void set(size_t x) //插入数据 { size_t index = x >> 5; //index表示的是数据存在哪一个int类型的位置上 size_t num = x % 32; //num表示数据存在32bit的具体位置 if (!(_a[index] & (1 << num))) //1<<num表示数据x的位置下标,&结果为0,表示此位置上没有数据 { ++_size; _a[index] |= (1 << num); //利用按位或关系将位置更改为1,表示此位置上现在存在数据 } } void Reset(size_t x) //删除数据 { size_t index = x >> 5; size_t num = index % 32; _a[index] &= (~(1 << num)); --_size; } bool Test(size_t x) //判断数据是否在40亿数据中 { size_t index = x >> 5; size_t num = x % 32; if (_a[index] & (1 << num)) return true; return false; } size_t size() //求数据的有效个数 { return _size; } void Resize(size_t size) //开辟空间 { _a.resize((size >> 5) + 1); } protected: vector<size_t> _a; size_t _size; };
‘布隆过滤器’也是利用位图的思想,它有一个m个比特个数的空间,每一个bit位都初始化为0,通过k种不同的hash函数,每个函数都确定出元素所在的不同位置,将这k个位置的bit位置为1,则将这个元素添加到m个bit的空间中。当需要对数据进行查找时,将k中hash函数得到的k个位置的bit位进行检查,若k个位置都为1,则数据存在,否则数据不存在。布隆过滤器是不允许进行删除数据的,因为那样会将k个位置置为0,可能会影响其他数据的存在性,从而存在错误。
下面是‘布隆过滤器’的实现:
//实现布隆过滤器 template <class K> //使用搜索到的5种Hash函数 struct _HashFunc1 { size_t DJBHash(const char *str) { if (!*str) return 0; register size_t hash = 5381; while (size_t ch = (size_t)*str++) { hash += (hash << 5) + ch; } return hash; } size_t operator()(const K& str) { return DJBHash(str.c_str()); } }; template <class K> struct _HashFunc2 { size_t SDBMHash(const char *str) { register size_t hash = 0; while (size_t ch = (size_t)*str++) { hash = 65599 * hash + ch; } return hash; } size_t operator()(const K& str) { return SDBMHash(str.c_str()); } }; template <class K> struct _HashFunc3 { size_t RSHash(const char *str) { register size_t hash = 0; size_t magic = 63689; while (size_t ch = (size_t)*str++) { hash = hash * magic + ch; magic *= 378551; } return hash; } size_t operator()(const K& str) { return RSHash(str.c_str()); } }; template <class K> struct _HashFunc4 { size_t APHash(const char *str) { register size_t hash = 0; size_t ch; for (long i = 0; ch = (size_t)*str++; i++) { if ((i & 1) == 0) { hash ^= ((hash << 7) ^ ch ^ (hash >> 3)); } else { hash ^= (~((hash << 11) ^ ch ^ (hash >> 5))); } } return hash; } size_t operator()(const K& str) { return APHash(str.c_str()); } }; template <class K> struct _HashFunc5 { size_t JSHash(const char *str) { if (!*str) return 0; register size_t hash = 1315423911; while (size_t ch = (size_t)*str++) { hash ^= ((hash << 5) + ch + (hash >> 2)); } return hash; } size_t operator()(const K& str) { return JSHash(str.c_str()); } }; size_t GetPrimeSize(size_t size) //求大于等于size的最小素数 { static const int _prime = 28; static const unsigned long _PrimeList[_prime] = { 53ul, 97ul, 193ul, 389ul, 769ul, 1543ul, 3079ul, 6151ul, 12289ul, 24593ul, 49157ul, 98317ul, 196613ul, 393241ul, 786433ul, 1572869ul, 3145739ul, 6291469ul, 12582917ul, 25165843ul, 50331653ul, 100663319ul, 201326611ul, 402653189ul, 805306457ul, 1610612741ul, 3221225473ul, 4294967291ul }; for (size_t i = 0; i < _prime; i++) { if (_PrimeList[i] >= size) { return _PrimeList[i]; } } return _PrimeList[_prime - 1]; } template <class K = string, class HashFunc1 = _HashFunc1<K>, class HashFunc2 = _HashFunc2<K>, class HashFunc3 = _HashFunc3<K>, class HashFunc4 = _HashFunc4<K>, class HashFunc5 = _HashFunc5<K>> class BloomFilter { public: BloomFilter(size_t size = 0) //构造 { _capacity = GetPrimeSize(size); _bm.Resize(_capacity); } void set(const K& key) { size_t index1 = HashFunc1()(key); size_t index2 = HashFunc2()(key); size_t index3 = HashFunc3()(key); size_t index4 = HashFunc4()(key); size_t index5 = HashFunc5()(key); _bm.set((index1) % _capacity); _bm.set((index2) % _capacity); _bm.set((index3) % _capacity); _bm.set((index4) % _capacity); _bm.set((index5) % _capacity); } bool Test(const K& key) //测试数据是否存在 { size_t index1 = HashFunc1()(key); if (!_bm.Test((index1) % _capacity)) { return false; } size_t index2 = HashFunc2()(key); if (!_bm.Test((index2) % _capacity)) { return false; } size_t index3 = HashFunc3()(key); if (!_bm.Test((index3) % _capacity)) { return false; } size_t index4 = HashFunc4()(key); if (!_bm.Test((index4) % _capacity)) { return false; } size_t index5 = HashFunc5()(key); if (!_bm.Test((index5) % _capacity)) { return false; } return true; } protected: BitMap _bm; size_t _capacity; };
本文出自 “无心的执着” 博客,转载请与作者联系!
以上是关于数据存在?-‘布隆过滤器’的主要内容,如果未能解决你的问题,请参考以下文章