散列表(hashtable)——算法导论

王朝学院·作者佚名  2016-08-27
窄屏简体版  字體: 小  |  中  |  大  |  超大  

1. 引言

许多应用都需要动态集合结构,它至少需要支持Insert,search和delete字典操作。散列表(hash table)是实现字典操作的一种有效的数据结构。

2. 直接寻址表

在介绍散列表之前,我们先介绍直接寻址表。

当关键字的全域U(关键字的范围)比较小时,直接寻址是一种简单而有效的技术。我们假设某应用要用到一个动态集合,其中每个元素的关键字都是取自于全域U={0,1,…,m-1},其中m不是一个很大的数。另外,假设每个元素的关键字都不同。

为表示动态集合,我们用一个数组,或称为直接寻址表(direct-address table),记为T[0~m-1],其中每一个位置(slot,槽)对应全域U中的一个关键字,对应规则是,槽k指向集合中关键字为k的元素,如果集合中没有关键字为k的元素,则T[k]=NIL。

几种字典操作实现起来非常简单:

上述的每一个操作的时间均为O(1)时间。

在某些应用中,我们其实可以把对象作为元素直接保存在寻址表的槽中,而不需要像上图所示使用指针指向该对象,这样可以节省空间。

3. 散列表

(1) 直接寻址的缺点

我们可以看出,直接寻址技术有几个明显的缺点:如果全域U很大,那么表T 将要申请一段非常长的空间,很可能会申请失败;对于全域较大,但是元素却十分稀疏的情况,使用这种存储方式将浪费大量的存储空间。

(2) 散列函数

为了克服直接寻址技术的缺点,而又保持其快速字典操作的优势,我们可以利用散列函数(hash function)

h:U→{0,1,2,…,m-1}

来计算关键字k所在的的位置,简单的讲,散列函数h(k)的作用是将范围较大的关键字映射到一个范围较小的集合中。这时我们可以说,一个具有关键字k的元素被散列到槽h(k)上,或者说h(k)是关键字k的散列值。

示意图如下:

这时会产生一个问题:两个关键字可能映射到同一槽中(我们称之为冲突(collision)),并且不管你如何优化h(k)函数,这种情况都会发生(因为|U|>m)。

因此我们现在面临两个问题,一是遇到冲突时如何解决;二是要找出一个的函数h(k)能够尽量的减少冲突;

(3) 通过链表法解决冲突

我们先来解决第一个问题。

解决办法就是,我们把同时散列到同一槽中的元素以链表的形式“串联”起来,而该槽中保存的是指向该链表的指针。如下图所示:

采用该解决办法后,我们可以通过如下的操作方式来进行字典操作:

下面我们来分析上图各操作的性能。

首先是插入操作,很明显时间为O(1)。

然后分析删除操作,其花费的时间相当于从链表中删除一个元素的时间:如果链表T[h(k)]是双链表,花费的时间为O(1);如果链表T[h(k)]是单链表,则花费的时间和查找操作的渐进运行时间相同。

下面我们重点分析查找运行时间:

首先,我们假定任何一个给定元素都等可能地散列在散列表T的任何一个槽位中,且与其他元素被散列在T的哪个位置无关。我们称这个假设为简单均匀散列(simple uniform hashing)。

不失一般性,我们设散列表T的m个槽位散列了n个元素,则平均每个槽位散列了α = n/m个元素,我们称α为T的装载因子(load factor)。我们记位于槽位j的链表为T[j](j=1,2,…,m-1),而nj表示链表T[j]的长度,于是有

n = n0+n1+…+nm-1,

且E[nj] = α = n / m。

现在我们分查找成功和查找不成功两种情况讨论。

① 查找不成功

在查找不成功的情况下,我们需要遍历链表T[j]的每一个元素,而链表T[j]的长度是α,因此需要时间O(α),加上索引到T(j)的时间O(1),总时间为θ(1 + α)。

② 查找成功

在查找成功的情况下,我们无法准确知道遍历到链表T[j]的何处停止,因此我们只能讨论平均情况。

我们设xi是散列表T的第i个元素(假设我们按插入顺序对散列表T中的n个元素进行了1~n的编号),ki表示xi.key,其中i = 1,2,…,n,再定义随机变量Xij=I{h(ki)=h(kj)},即:

在简单均匀散列的假设下有

P{h(ki)=h(kj)} = 1 / m,

E[Xij] = 1 / m。

则所需检查的元素的数目的期望是:

因此,一次成功的检查所需要的时间是O(2 + α / 2 –α / 2n) = θ(1 + α)。

综合上面的分析,在平均下,全部的字典操作都可以在O(1)时间下完成。

4. 散列函数

现在我们来解决第二个问题:如何构造一个好的散列函数。

一个好的散列函数应(近似地)满足简单均匀散列:每个关键字都等可能的被散列到各个槽位,并与其他关键字散列到哪一个槽位无关(但很遗憾,我们一般无法检验这一条件是否成立)。

在实际应用中,常常可以可以运用启发式方法来构造性能好的散列函数。设计过程中,可以利用关键字分布的有用信息。一个好的方法导出的散列值,在某种程度上应独立于数据可能存在的任何模式。

下面给出两种基本的构造散列函数的方法:

(1) 除法散列法

除法散列法的做法很简单,就是让关键字k去除以一个数m,取余数,这样就将k映射到m个槽位中的某一个,即散列函数是:

h(k) = k mod m ,

由于只做一次除法运算,该方法的速度是非常快的。但应当注意的是,我们在选取m的值时,应当避免一些选取一些值。例如,m不应是2的整数幂,因为如果m = 2 ^ p,则h(k)就是k的p个最低位数字。除非我们已经知道各种最低p位的排列是等可能的,否则我们最好慎重的选择m。而一个不太接近2的整数幂的素数,往往是较好的选择。

(2) 乘法散列法

该方法包含两个步骤。第一步:用关键字k乘以A(0 < A < 1),并提取kA的小数部分;第二步:用m乘以这个值,在向下取整,即散列函数是:

h(k) = [m (kA mod 1)],

这里“kA mod 1”的是取kA小数部分的意思,即kA –[kA]。

乘法散列法的一个优点是,一般我们对m的选择不是特别的关键,一般选择它为2的整数幂即可。虽然这个方法对任意的A都适用,但Knuth认为,A ≈ (√5 - 1)/ 2 = 0.618033988…是一个比较理想的值

(5) 布隆过滤器

布隆过滤器(Bloom Filter)是一种常被用来检验一个元素是否在一个集合里面的算法(从这里我们可以看出,这个集合只需要保存比对元素的“指纹”即可,而不需要保存比对元素的全部信息),由一个很长的二进制向量和一系列随机映射函数组成。相较于其他算法,它具有空间利用率高,检测速度快等优点。

在介绍布隆过滤器之前,我们先假设这样一种场景:某公司致力于解决用户常常遭遇骚扰电话的问题。该公司打算建立一个骚扰电话号码的黑名单,即把所有骚扰电话的号码保存到一张hash表中。当用户接到某个陌生电话时,服务器会立即将该号码与黑名单进行比对,若比对成功,则对该号码进行拦截。

他们当然不会直接将骚扰电话号码保存在hash表中,而是对每一个号码利用某种算法进行数据压缩,最终得到一个8字节的信息指纹,然后将其存入表中。但即便如此,问题还是来了:由于hash表的空间利用率大约只有50%,等价换算过来,储存一个号码将要花费16字节的空间。按照这样计算,储存1亿个号码将要花费大约1.6G的空间,储存几十亿的号码可能需要上百G的空间。那么有没有更好解决办法呢?

这时,布隆过滤器就派上用场了。假设我们有1亿条骚扰电话号码需要记录,我们的做法是,首先建立一个2亿字节(即16亿位,并假设我们对这16亿位以1~16亿的顺序进行了编号)的向量,将每位都置为0。当要插入某个电话号码x时,我们使用某种算法(该算法可以做到每个位被映射的概率是一样的,且某个映射的分布与其他的映射分布无关)让号码x映射到1~16亿中的8个位上,然后把这8个位设为1。当查找时,利用同样的方法将号码映射到8个位上,若这8个位都为1,则说明该号码在黑名单中,否则就不在。

我们可以发现,布隆过滤器的做法在思想上和hash函数将关键字映射到hash表的做法很相似,因此布隆过滤器也会遇到冲突问题,这会导致将一个“好”的号码误判为骚扰号码(但绝对不会将骚扰号码误判为一个“好”的号码)。下面我们通过计算来证明,在大多数情况和场景中,这种误判我们是可以忍受的。

假设某布隆过滤器共有的m个槽位,我们要把n个号码添加到其中,而每个号码会映射k个槽位。那么,添加这n个号码将会产生kn次映射。因为这m个槽位中,每个槽被映射到的概率是相等的。因此,

在一次映射中,某个槽位被映射到的概率(即该槽位值为1的概率)为

该槽位值为0的概率为

经过kn次映射后,某个槽值为0的概率为

为1的概率为

所以,误判(k个槽位均为1)的概率就为

利用,上式可化为:

这时我们注意到,当k=1时,情况就就变成了hash table的情况,

根据自变量的不同我们分以下两种方式讨论:

① 我们把误判率p看作关于装载因子α的函数(k看作常数),这时我们从函数的函数图像

中可以得出一下结论:

随着装载因子α(α = n / m)的增大,误判率(或者是产生冲突的概率)也将增大,但增长速度逐渐减慢。

要使误判率小于0.5,装载因子必须小于0.7。这也从某种程度上解释了为什么JDK HashMap的装载因子默认是0.75。

② 我们把误判率p看作关于k的函数(α作为常数),通过对p求导分析,我们发现,当k=ln2 / α时,误判率p取得最小值。此时,p = 2^(-k)(或者k = – ln p / ln 2),这个结论让我们能够根据可以忍受的误判率计算出最为合适的k值。

下面给出一个BloomFilter的java实现代码(来自:https://github.com/MagnusS/Java-BloomFilter,只是把其中的变量和方法名换成了上文提及的):

publicclassBloomFilter<E>implementsSerializable {PRivatestaticfinallongserialVersionUID = -9077350041930475408L;privateBitSet bitset;//二进制向量privateintslotSize;//二进制向量的总位(槽)数(文中的m)privatedoubleloadFactor;//装载因子 (文中的α)privateintcapacity;//布隆过滤器的容量(文中的n)privateintsize;//装载的数目privateintk;//一个元素对应的位数(文中的k)staticfinalCharset charset = Charset.forName("UTF-8");staticfinalString hashName = "md5";//默认采用MD5算法,也可改为SHA1staticfinalMessageDigest digestFunction;static{

MessageDigest tmp;try{

tmp=java.security.MessageDigest.getInstance(hashName);

}catch(NoSuchAlgorithmException e) {

tmp=null;

}

digestFunction=tmp;

}publicBloomFilter(intslotSize,intcapacity) {this(slotSize / (double) capacity, capacity, (int) Math.round((slotSize / (double) capacity) * Math.log(2.0)));

}publicBloomFilter(doublefalsePositiveProbability,intcapacity) {this(Math.log(2) / (Math.ceil(-(Math.log(falsePositiveProbability) / Math.log(2)))),//loadFactor = ln2 / k;capacity,//(int) Math.ceil(-(Math.log(falsePositiveProbability) / Math.log(2))));//k = -ln p / ln2}publicBloomFilter(intslotSize,intcapacity,intsize, BitSet filterData) {this(slotSize, capacity);this.bitset =filterData;this.size =size;

}publicBloomFilter(doubleloadFactor,intcapacity,intk) {

size= 0;this.loadFactor =loadFactor;this.capacity =capacity;this.k =k;this.slotSize = (int) Math.ceil(capacity *loadFactor);this.bitset =newBitSet(slotSize);

}publicstaticintcreateHash(String val, Charset charset) {returncreateHash(val.getBytes(charset));

}publicstaticintcreateHash(String val) {returncreateHash(val, charset);

}publicstaticintcreateHash(byte[] data) {returncreateHashes(data, 1)[0];

}publicstaticint[] createHashes(byte[] data,inthashes) {int[] result =newint[hashes];intk = 0;bytesalt = 0;while(k <hashes) {byte[] digest;synchronized(digestFunction) {

digestFunction.update(salt);

salt++;

digest=digestFunction.digest(data);

}for(inti = 0; i < digest.length / 4 && k < hashes; i++) {inth = 0;for(intj = (i * 4); j < (i * 4) + 4; j++) {

h<<= 8;

h|= ((int) digest[j]) & 0xFF;

}

result[k]=h;

k++;

}

}returnresult;

}/*** Compares the contents of two instances to see if they are equal.

*

*@paramobj

* is the object to compare to.

*@returnTrue if the contents of the objects are equal.*/@Overridepublicbooleanequals(Object obj) {if(obj ==null) {returnfalse;

}if(getClass() !=obj.getClass()) {returnfalse;

}finalBloomFilter<E> other = (BloomFilter<E>) obj;if(this.capacity !=other.capacity) {returnfalse;

}if(this.k !=other.k) {returnfalse;

}if(this.slotSize !=other.slotSize) {returnfalse;

}if(this.bitset != other.bitset && (this.bitset ==null|| !this.bitset.equals(other.bitset))) {returnfalse;

}returntrue;

}/*** Calculates a hash code for this class.

*

*@returnhash code representing the contents of an instance of this class.*/@OverridepublicinthashCode() {inthash = 7;

hash= 61 * hash + (this.bitset !=null?this.bitset.hashCode() : 0);

hash= 61 * hash +this.capacity;

hash= 61 * hash +this.slotSize;

hash= 61 * hash +this.k;returnhash;

}/*** Calculates the expected probability of false positives based on the

* number of expected filter elements and the size of the Bloom filter.

* <br />

* <br />

* The value returned by this method is the <i>expected</i> rate of false

* positives, assuming the number of inserted elements equals the number of

* expected elements. If the number of elements in the Bloom filter is less

* than the expected value, the true probability of false positives will be

* lower.

*

*@returnexpected probability of false positives.*/publicdoubleexpectedFalsePositiveProbability() {returngetFalsePositiveProbability(capacity);

}/*** Calculate the probability of a false positive given the specified number

* of inserted elements.

*

*@paramnumberOfElements

* number of inserted elements.

*@returnprobability of a false positive.*/publicdoublegetFalsePositiveProbability(doublenumberOfElements) {//(1 - e^(-k * n / m)) ^ kreturnMath.pow((1 - Math.exp(-k * (double) numberOfElements / (double) slotSize)), k);

}/*** Get the current probability of a false positive. The probability is

* calculated from the size of the Bloom filter and the current number of

* elements added to it.

*

*@returnprobability of false positives.*/publicdoublegetFalsePositiveProbability() {returngetFalsePositiveProbability(size);

}/*** Returns the value chosen for K.<br />

* <br />

* K is the optimal number of hash functions based on the size of the Bloom

* filter and the expected number of inserted elements.

*

*@returnoptimal k.*/publicintgetK() {returnk;

}/*** Sets all bits to false in the Bloom filter.*/publicvoidclear() {

bitset.clear();

size= 0;

}/*** Adds an object to the Bloom filter. The output from the object's

* toString() method is used as input to the hash functions.

*

*@paramelement

* is an element to register in the Bloom filter.*/publicvoidadd(E element) {

add(element.toString().getBytes(charset));

}/*** Adds an array of bytes to the Bloom filter.

*

*@parambytes

* array of bytes to add to the Bloom filter.*/publicvoidadd(byte[] bytes) {int[] hashes =createHashes(bytes, k);for(inthash : hashes)

bitset.set(Math.abs(hash%slotSize));

size++;

}/*** Adds all elements from a Collection to the Bloom filter.

*

*@paramc

* Collection of elements.*/publicvoidaddAll(Collection<?extendsE>c) {for(E element : c)

add(element);

}/*** Returns true if the element could have been inserted into the Bloom

* filter. Use getFalsePositiveProbability() to calculate the probability of

* this being correct.

*

*@paramelement

* element to check.

*@returntrue if the element could have been inserted into the Bloom

* filter.*/publicbooleancontains(E element) {returncontains(element.toString().getBytes(charset));

}/*** Returns true if the array of bytes could have been inserted into the

* Bloom filter. Use getFalsePositiveProbability() to calculate the

* probability of this being correct.

*

*@parambytes

* array of bytes to check.

*@returntrue if the array could have been inserted into the Bloom filter.*/publicbooleancontains(byte[] bytes) {int[] hashes =createHashes(bytes, k);for(inthash : hashes) {if(!bitset.get(Math.abs(hash %slotSize))) {returnfalse;

}

}returntrue;

}/*** Returns true if all the elements of a Collection could have been inserted

* into the Bloom filter. Use getFalsePositiveProbability() to calculate the

* probability of this being correct.

*

*@paramc

* elements to check.

*@returntrue if all the elements in c could have been inserted into the

* Bloom filter.*/publicbooleancontainsAll(Collection<?extendsE>c) {for(E element : c)if(!contains(element))returnfalse;returntrue;

}/*** Read a single bit from the Bloom filter.

*

*@parambit

* the bit to read.

*@returntrue if the bit is set, false if it is not.*/publicbooleangetBit(intbit) {returnbitset.get(bit);

}/*** Set a single bit in the Bloom filter.

*

*@parambit

* is the bit to set.

*@paramvalue

* If true, the bit is set. If false, the bit is cleared.*/publicvoidsetBit(intbit,booleanvalue) {

bitset.set(bit, value);

}/*** Return the bit set used to store the Bloom filter.

*

*@returnbit set representing the Bloom filter.*/publicBitSet getBitSet() {returnbitset;

}/*** Returns the number of bits in the Bloom filter. Use count() to retrieve

* the number of inserted elements.

*

*@returnthe size of the bitset used by the Bloom filter.*/publicintslotSize() {returnslotSize;

}/*** Returns the number of elements added to the Bloom filter after it was

* constructed or after clear() was called.

*

*@returnnumber of elements added to the Bloom filter.*/publicintsize() {returnsize;

}/*** Returns the expected number of elements to be inserted into the filter.

* This value is the same value as the one passed to the constructor.

*

*@returnexpected number of elements.*/publicintcapacity() {returncapacity;

}/*** Get expected number of bits per element when the Bloom filter is full.

* This value is set by the constructor when the Bloom filter is created.

* See also getBitsPerElement().

*

*@returnexpected number of bits per element.*/publicdoublegetLoadFactor() {returnthis.loadFactor;

}

}

 
 
免责声明:本文为网络用户发布,其观点仅代表作者个人观点,与本站无关,本站仅提供信息存储服务。文中陈述内容未经本站证实,其真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
 
© 2005- 王朝網路 版權所有 導航