海量数据处理算法—Bloom Filter

1. Bloom-Filter算法简介

Bloom-Filter,即布隆过滤器,1970年由Bloom中提出。它可以用于检索一个元素是否在一个集合中。

Bloom Filter(BF)是一种空间效率很高的随机数据结构,它利用位数组很简洁地表示一个集合,并能判断一个元素是否属于这个集合。它是一个判断元素是否存在集合的快速的概率算法。Bloom Filter有可能会出现错误判断,但不会漏掉判断。也就是Bloom Filter判断元素不再集合,那肯定不在。如果判断元素存在集合中,有一定的概率判断错误。因此,Bloom Filter不适合那些“零错误”的应用场合。而在能容忍低错误率的应用场合下,Bloom Filter比其他常见的算法(如hash,折半查找)极大节省了空间。

它的优点是空间效率和查询时间都远远超过一般的算法,缺点是有一定的误识别率和删除困难。

Bloom Filter的详细介绍:Bloom Filter

2、 Bloom-Filter的基本思想

Bloom-Filter算法的核心思想就是利用多个不同的Hash函数来解决“冲突”。

计算某元素x是否在一个集合中,首先能想到的方法就是将所有的已知元素保存起来构成一个集合R,然后用元素x跟这些R中的元素一一比较来判断是否存在于集合R中;我们可以采用链表等数据结构来实现。但是,随着集合R中元素的增加,其占用的内存将越来越大。试想,如果有几千万个不同网页需要下载,所需的内存将足以占用掉整个进程的内存地址空间。即使用MD5,UUID这些方法将URL转成固定的短小的字符串,内存占用也是相当巨大的。

于是,我们会想到用Hash table的数据结构,运用一个足够好的Hash函数将一个URL映射到二进制位数组(位图数组)中的某一位。如果该位已经被置为1,那么表示该URL已经存在。

Hash存在一个冲突(碰撞)的问题,用同一个Hash得到的两个URL的值有可能相同。为了减少冲突,我们可以多引入几个Hash,如果通过其中的一个Hash值我们得出某元素不在集合中,那么该元素肯定不在集合中。只有在所有的Hash函数告诉我们该元素在集合中时,才能确定该元素存在于集合中。这便是Bloom-Filter的基本思想。


原理要点:一是位数组, 而是k个独立hash函数。

1)位数组:

假设Bloom Filter使用一个m比特的数组来保存信息,初始状态时,Bloom Filter是一个包含m位的位数组,每一位都置为0,即BF整个数组的元素都设置为0。


2)添加元素,k个独立hash函数

为了表达S={x1, x2,…,xn}这样一个n个元素的集合,Bloom Filter使用k个相互独立的哈希函数(Hash Function),它们分别将集合中的每个元素映射到{1,…,m}的范围中。

当我们往Bloom Filter中增加任意一个元素x时候,我们使用k个哈希函数得到k个哈希值,然后将数组中对应的比特位设置为1。即第i个哈希函数映射的位置hashi(x)就会被置为1(1≤i≤k)。

注意,如果一个位置多次被置为1,那么只有第一次会起作用,后面几次将没有任何效果。在下图中,k=3,且有两个哈希函数选中同一个位置(从左边数第五位,即第二个“1“处)。


3)判断元素是否存在集合

在判断y是否属于这个集合时,我们只需要对y使用k个哈希函数得到k个哈希值,如果所有hashi(y)的位置都是11ik),即k个位置都被设置为1了,那么我们就认为y是集合中的元素,否则就认为y不是集合中的元素。下图中y1就不是集合中的元素(因为y1有一处指向了“0”位)。y2或者属于这个集合,或者刚好是一个false positive



显然这 个判断并不保证查找的结果是100%正确的。

Bloom Filter的缺点:

1)Bloom Filter无法从Bloom Filter集合中删除一个元素因为该元素对应的位会牵动到其他的元素。所以一个简单的改进就是 counting Bloom filter,用一个counter数组代替位数组,就可以支持删除了。此外,Bloom Filter的hash函数选择会影响算法的效果。

2)还有一个比较重要的问题,如何根据输入元素个数n,确定位数组m的大小及hash函数个数hash函数选择会影响算法的效果当hash函数个数k=(ln2)*(m/n)时错误率最小。在错误率不大于E的情况 下,m至少要等于n*lg(1/E)才能表示任意n个元素的集合。但m还应该更大些,因为还要保证bit数组里至少一半为0,则m应 该>=nlg(1/E)*lge ,大概就是nlg(1/E)1.44(lg表示以2为底的对数)。

举个例子我们假设错误率为0.01,则此时m应大概是n的13倍。这样k大概是8个。

注意:

这里m与n的单位不同,m是bit为单位,而n则是以元素个数为单位(准确的说是不同元素的个数)。通常单个元素的长度都是有很多bit的。所以使用bloom filter内存上通常都是节省的。

一般BF可以与一些key-value的数据库一起使用,来加快查询。由于BF所用的空间非常小,所有BF可以常驻内存。这样子的话,对于大部分不存在的元素,我们只需要访问内存中的BF就可以判断出来了,只有一小部分,我们需要访问在硬盘上的key-value数据库。从而大大地提高了效率。


一个Bloom Filter有以下参数:


m bit数组的宽度(bit数)
n 加入其中的key的数量
k 使用的hash函数的个数
f False Positive的比率

Bloom Filter的f满足下列公式:


在给定m和n时,能够使f最小化的k值为:

此时给出的f为:

根据以上公式,对于任意给定的f,我们有:


n = m ln(0.6185) / ln(f) <wbr><wbr>[1]</wbr></wbr>

同时,我们需要k个hash来达成这个目标:

k = - ln(f) / ln(2) <wbr><wbr><wbr><wbr><wbr><wbr><wbr><wbr>[2]</wbr></wbr></wbr></wbr></wbr></wbr></wbr></wbr>

由于k必须取整数,我们在Bloom Filter的程序实现中,还应该使用上面的公式来求得实际的f:

f = (1 – e -kn/m) k<wbr><wbr><wbr><wbr><wbr><wbr><wbr><wbr><wbr>[3]</wbr></wbr></wbr></wbr></wbr></wbr></wbr></wbr></wbr>

以上3个公式是程序实现Bloom Filter的关键公式。

3、 扩展CounterBloom Filter

CounterBloom Filter

BloomFilter有个缺点,就是不支持删除操作,因为它不知道某一个位从属于哪些向量。那我们可以给Bloom Filter加上计数器,添加时增加计数器,删除时减少计数器。

但这样的Filter需要考虑附加的计数器大小,假如同个元素多次插入的话,计数器位数较少的情况下,就会出现溢出问题。如果对计数器设置上限值的话,会导致Cache Miss,但对某些应用来说,这并不是什么问题,如Web Sharing。

Compressed Bloom Filter

为了能在服务器之间更快地通过网络传输Bloom Filter,我们有方法能在已完成Bloom Filter之后,得到一些实际参数的情况下进行压缩。

将元素全部添加入Bloom Filter后,我们能得到真实的空间使用率,用这个值代入公式计算出一个比m小的值,重新构造Bloom Filter,对原先的哈希值进行求余处理,在误判率不变的情况下,使得其内存大小更合适。


4、Bloom-Filter的应用

Bloom-Filter一般用于在大数据量的集合中判定某元素是否存在。例如邮件服务器中的垃圾邮件过滤器。在搜索引擎领域,Bloom-Filter最常用于网络蜘蛛(Spider)的URL过滤,网络蜘蛛通常有一个URL列表,保存着将要下载和已经下载的网页的URL,网络蜘蛛下载了一个网页,从网页中提取到新的URL后,需要判断该URL是否已经存在于列表中。此时,Bloom-Filter算法是最好的选择。

1.key-value 加快查询

一般Bloom-Filter可以与一些key-value的数据库一起使用,来加快查询。

一般key-value存储系统的values存在硬盘,查询就是件费时的事。Storage的数据都插入Filter,在Filter中查询都不存在时,那就不需要去Storage查询了。False Position出现时,只是会导致一次多余的Storage查询。

由于Bloom-Filter所用的空间非常小,所有BF可以常驻内存。这样子的话,对于大部分不存在的元素,我们只需要访问内存中的Bloom-Filter就可以判断出来了,只有一小部分,我们需要访问在硬盘上的key-value数据库。从而大大地提高了效率。如图:

海量数据处理算法—Bloom Filter


2 .GoogleBigTable

Google的BigTable也使用了Bloom Filter,以减少不存在的行或列在磁盘上的查询,大大提高了数据库的查询操作的性能。

3.Proxy-Cache

在Internet Cache Protocol中的Proxy-Cache很多都是使用Bloom Filter存储URLs,除了高效的查询外,还能很方便得传输交换Cache信息。

4.网络应用

1)P2P网络中查找资源操作,可以对每条网络通路保存Bloom Filter,当命中时,则选择该通路访问。

2)广播消息时,可以检测某个IP是否已发包。

3)检测广播消息包的环路,将Bloom Filter保存在包里,每个节点将自己添加入Bloom Filter。

4)信息队列管理,使用Counter Bloom Filter管理信息流量。

5. 垃圾邮件地址过滤

像网易,QQ这样的公众电子邮件(email)提供商,总是需要过滤来自发送垃圾邮件的人(spamer)的垃圾邮件。

一个办法就是记录下那些发垃圾邮件的email地址。由于那些发送者不停地在注册新的地址,全世界少说也有几十亿个发垃圾邮件的地址,将他们都存起来则需要大量的网络服务器。

如果用哈希表,每存储一亿个email地址,就需要1.6GB的内存(用哈希表实现的具体办法是将每一个email地址对应成一个八字节的信息指纹,然后将这些信息指纹存入哈希表,由于哈希表的存储效率一般只有50%,因此一个email地址需要占用十六个字节。一亿个地址大约要1.6GB,即十六亿字节的内存)。因此存贮几十亿个邮件地址可能需要上百GB的内存。

而Bloom Filter只需要哈希表1/8到1/4的大小就能解决同样的问题。

BloomFilter决不会漏掉任何一个在黑名单中的可疑地址。而至于误判问题,常见的补救办法是在建立一个小的白名单,存储那些可能别误判的邮件地址。



5、Bloom-Filter的具体实现

c语言实现:

stdafx.h:

  1. #pragmaonce
  2. #include<stdio.h>
  3. #include"stdlib.h"
  4. #include<iostream>
  5. #include<time.h>
  6. usingnamespacestd;
  1. #include"stdafx.h"
  2. #defineARRAY_SIZE256/*wegetthe256charsofeachline*/
  3. #defineSIZE48000000/*sizeshouldbe1/8ofmax*/
  4. #defineMAX384000000/*themaxbitspace*/
  5. #defineSETBIT(ch,n)ch[n/8]|=1<<(7-n%8)
  6. #defineGETBIT(ch,n)(ch[n/8]&1<<(7-n%8))>>(7-n%8)
  7. unsignedintlen(char*ch);/*functionstocalculatethelengthoftheurl*/
  8. unsignedintRSHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  9. unsignedintJSHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  10. unsignedintPJWHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  11. unsignedintELFHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  12. unsignedintBKDRHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  13. unsignedintSDBMHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  14. unsignedintDJBHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  15. unsignedintDEKHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  16. unsignedintBPHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  17. unsignedintFNVHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  18. unsignedintAPHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  19. unsignedintHFLPHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  20. unsignedintHFHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  21. unsignedintStrHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  22. unsignedintTianlHash(char*str,unsignedintlen);/*functionstocalculatethehashvalueoftheurl*/
  23. intmain()
  24. {
  25. inti,num,num2=0;/*thenumbertorecordtherepeatedurlsandthetotalofit*/
  26. unsignedinttt=0;
  27. intflag;/*ithelpstocheckweathertheurlhasalreadyexisted*/
  28. charbuf[257];/*ithelpstoprintthestarttimeoftheprogram*/
  29. time_ttmp=time(NULL);
  30. charfile1[100],file2[100];
  31. FILE*fp1,*fp2;/*pointertothefile*/
  32. charch[ARRAY_SIZE];
  33. char*vector;/*thebitspace*/
  34. vector=(char*)calloc(SIZE,sizeof(char));
  35. printf("Pleaseenterthefilewithrepeatedurls:\n");
  36. scanf("%s",&file1);
  37. if((fp1=fopen(file1,"rb"))==NULL){/*openthegoalfile*/
  38. printf("Connotopenthefile%s!\n",file1);
  39. }
  40. printf("Pleaseenterthefileyouwanttosaveto:\n");
  41. scanf("%s",&file2);
  42. if((fp2=fopen(file2,"w"))==NULL){
  43. printf("Connotopenthefile%s\n",file2);
  44. }
  45. strftime(buf,32,"%Y-%m-%d%H:%M:%S",localtime(&tmp));
  46. printf("%s\n",buf);/*printthesystemtime*/
  47. for(i=0;i<SIZE;i++){
  48. vector[i]=0;/*set0*/
  49. }
  50. while(!feof(fp1)){/*thecheckprocess*/
  51. fgets(ch,ARRAY_SIZE,fp1);
  52. flag=0;
  53. tt++;
  54. if(GETBIT(vector,HFLPHash(ch,len(ch))%MAX)){
  55. flag++;
  56. }else{
  57. SETBIT(vector,HFLPHash(ch,len(ch))%MAX);
  58. }
  59. if(GETBIT(vector,StrHash(ch,len(ch))%MAX)){
  60. flag++;
  61. }else{
  62. SETBIT(vector,StrHash(ch,len(ch))%MAX);
  63. }
  64. if(GETBIT(vector,HFHash(ch,len(ch))%MAX)){
  65. flag++;
  66. }else{
  67. SETBIT(vector,HFHash(ch,len(ch))%MAX);
  68. }
  69. if(GETBIT(vector,DEKHash(ch,len(ch))%MAX)){
  70. flag++;
  71. }else{
  72. SETBIT(vector,DEKHash(ch,len(ch))%MAX);
  73. }
  74. if(GETBIT(vector,TianlHash(ch,len(ch))%MAX)){
  75. flag++;
  76. }else{
  77. SETBIT(vector,TianlHash(ch,len(ch))%MAX);
  78. }
  79. if(GETBIT(vector,SDBMHash(ch,len(ch))%MAX)){
  80. flag++;
  81. }else{
  82. SETBIT(vector,SDBMHash(ch,len(ch))%MAX);
  83. }
  84. if(flag<6)
  85. num2++;
  86. else
  87. fputs(ch,fp2);
  88. /*printf("%d",flag);*/
  89. }
  90. /*theresult*/
  91. printf("\nThereare%durls!\n",tt);
  92. printf("\nThereare%dnotrepeatedurls!\n",num2);
  93. printf("Thereare%drepeatedurls!\n",tt-num2);
  94. fclose(fp1);
  95. fclose(fp2);
  96. return0;
  97. }
  98. /*functionsmaybeusedinthemain*/
  99. unsignedintlen(char*ch)
  100. {
  101. intm=0;
  102. while(ch[m]!='\0'){
  103. m++;
  104. }
  105. returnm;
  106. }
  107. unsignedintRSHash(char*str,unsignedintlen){
  108. unsignedintb=378551;
  109. unsignedinta=63689;
  110. unsignedinthash=0;
  111. unsignedinti=0;
  112. for(i=0;i<len;str++,i++){
  113. hash=hash*a+(*str);
  114. a=a*b;
  115. }
  116. returnhash;
  117. }
  118. /*EndOfRSHashFunction*/
  119. unsignedintJSHash(char*str,unsignedintlen)
  120. {
  121. unsignedinthash=1315423911;
  122. unsignedinti=0;
  123. for(i=0;i<len;str++,i++){
  124. hash^=((hash<<5)+(*str)+(hash>>2));
  125. }
  126. returnhash;
  127. }
  128. /*EndOfJSHashFunction*/
  129. unsignedintPJWHash(char*str,unsignedintlen)
  130. {
  131. constunsignedintBitsInUnsignedInt=(unsignedint)(sizeof(unsignedint)*8);
  132. constunsignedintThreeQuarters=(unsignedint)((BitsInUnsignedInt*3)/4);
  133. constunsignedintOneEighth=(unsignedint)(BitsInUnsignedInt/8);
  134. constunsignedintHighBits=(unsignedint)(0xFFFFFFFF)<<(BitsInUnsignedInt-OneEighth);
  135. unsignedinthash=0;
  136. unsignedinttest=0;
  137. unsignedinti=0;
  138. for(i=0;i<len;str++,i++){
  139. hash=(hash<<OneEighth)+(*str);
  140. if((test=hash&HighBits)!=0){
  141. hash=((hash^(test>>ThreeQuarters))&(~HighBits));
  142. }
  143. }
  144. returnhash;
  145. }
  146. /*EndOfP.J.WeinbergerHashFunction*/
  147. unsignedintELFHash(char*str,unsignedintlen)
  148. {
  149. unsignedinthash=0;
  150. unsignedintx=0;
  151. unsignedinti=0;
  152. for(i=0;i<len;str++,i++){
  153. hash=(hash<<4)+(*str);
  154. if((x=hash&0xF0000000L)!=0){
  155. hash^=(x>>24);
  156. }
  157. hash&=~x;
  158. }
  159. returnhash;
  160. }
  161. /*EndOfELFHashFunction*/
  162. unsignedintBKDRHash(char*str,unsignedintlen)
  163. {
  164. unsignedintseed=131;/*31131131313131131313etc..*/
  165. unsignedinthash=0;
  166. unsignedinti=0;
  167. for(i=0;i<len;str++,i++)
  168. {
  169. hash=(hash*seed)+(*str);
  170. }
  171. returnhash;
  172. }
  173. /*EndOfBKDRHashFunction*/
  174. unsignedintSDBMHash(char*str,unsignedintlen)
  175. {
  176. unsignedinthash=0;
  177. unsignedinti=0;
  178. for(i=0;i<len;str++,i++){
  179. hash=(*str)+(hash<<6)+(hash<<16)-hash;
  180. }
  181. returnhash;
  182. }
  183. /*EndOfSDBMHashFunction*/
  184. unsignedintDJBHash(char*str,unsignedintlen)
  185. {
  186. unsignedinthash=5381;
  187. unsignedinti=0;
  188. for(i=0;i<len;str++,i++){
  189. hash=((hash<<5)+hash)+(*str);
  190. }
  191. returnhash;
  192. }
  193. /*EndOfDJBHashFunction*/
  194. unsignedintDEKHash(char*str,unsignedintlen)
  195. {
  196. unsignedinthash=len;
  197. unsignedinti=0;
  198. for(i=0;i<len;str++,i++){
  199. hash=((hash<<5)^(hash>>27))^(*str);
  200. }
  201. returnhash;
  202. }
  203. /*EndOfDEKHashFunction*/
  204. unsignedintBPHash(char*str,unsignedintlen)
  205. {
  206. unsignedinthash=0;
  207. unsignedinti=0;
  208. for(i=0;i<len;str++,i++){
  209. hash=hash<<7^(*str);
  210. }
  211. returnhash;
  212. }
  213. /*EndOfBPHashFunction*/
  214. unsignedintFNVHash(char*str,unsignedintlen)
  215. {
  216. constunsignedintfnv_prime=0x811C9DC5;
  217. unsignedinthash=0;
  218. unsignedinti=0;
  219. for(i=0;i<len;str++,i++){
  220. hash*=fnv_prime;
  221. hash^=(*str);
  222. }
  223. returnhash;
  224. }
  225. /*EndOfFNVHashFunction*/
  226. unsignedintAPHash(char*str,unsignedintlen)
  227. {
  228. unsignedinthash=0xAAAAAAAA;
  229. unsignedinti=0;
  230. for(i=0;i<len;str++,i++){
  231. hash^=((i&1)==0)?((hash<<7)^(*str)*(hash>>3)):
  232. (~((hash<<11)+(*str)^(hash>>5)));
  233. }
  234. returnhash;
  235. }
  236. /*EndOfAPHashFunction*/
  237. unsignedintHFLPHash(char*str,unsignedintlen)
  238. {
  239. unsignedintn=0;
  240. inti;
  241. char*b=(char*)&n;
  242. for(i=0;i<strlen(str);++i){
  243. b[i%4]^=str[i];
  244. }
  245. returnn%len;
  246. }
  247. /*EndOfHFLPHashFunction*/
  248. unsignedintHFHash(char*str,unsignedintlen)
  249. {
  250. intresult=0;
  251. char*ptr=str;
  252. intc;
  253. inti=0;
  254. for(i=1;c=*ptr++;i++)
  255. result+=c*3*i;
  256. if(result<0)
  257. result=-result;
  258. returnresult%len;
  259. }
  260. /*EndOfHKHashFunction*/
  261. unsignedintStrHash(char*str,unsignedintlen)
  262. {
  263. registerunsignedinth;
  264. registerunsignedchar*p;
  265. for(h=0,p=(unsignedchar*)str;*p;p++){
  266. h=31*h+*p;
  267. }
  268. returnh;
  269. }
  270. /*EndOfStrHashFunction*/
  271. unsignedintTianlHash(char*str,unsignedintlen)
  272. {
  273. unsignedlongurlHashValue=0;
  274. intilength=strlen(str);
  275. inti;
  276. unsignedcharucChar;
  277. if(!ilength){
  278. return0;
  279. }
  280. if(ilength<=256){
  281. urlHashValue=16777216*(ilength-1);
  282. }else{
  283. urlHashValue=42781900080;
  284. }
  285. if(ilength<=96){
  286. for(i=1;i<=ilength;i++){
  287. ucChar=str[i-1];
  288. if(ucChar<='Z'&&ucChar>='A'){
  289. ucChar=ucChar+32;
  290. }
  291. urlHashValue+=(3*i*ucChar*ucChar+5*i*ucChar+7*i+11*ucChar)%1677216;
  292. }
  293. }else{
  294. for(i=1;i<=96;i++)
  295. {
  296. ucChar=str[i+ilength-96-1];
  297. if(ucChar<='Z'&&ucChar>='A')
  298. {
  299. ucChar=ucChar+32;
  300. }
  301. urlHashValue+=(3*i*ucChar*ucChar+5*i*ucChar+7*i+11*ucChar)%1677216;
  302. }
  303. }
  304. returnurlHashValue;
  305. }
  306. /*EndOfTianlHashFunction*/


网上找到的php简单实现:

  1. <?php
  2. /**
  3. *ImplementsaBloomFilter
  4. */
  5. classBloomFilter{
  6. /**
  7. *Sizeofthebitarray
  8. *
  9. *@varint
  10. */
  11. protected$m;
  12. /**
  13. *Numberofhashfunctions
  14. *
  15. *@varint
  16. */
  17. protected$k;
  18. /**
  19. *Numberofelementsinthefilter
  20. *
  21. *@varint
  22. */
  23. protected$n;
  24. /**
  25. *Thebitsetholdingthefilterinformation
  26. *
  27. *@vararray
  28. */
  29. protected$bitset;
  30. /**
  31. *计算最优的hash函数个数:当hash函数个数k=(ln2)*(m/n)时错误率最小
  32. *
  33. *@paramint$mbit数组的宽度(bit数)
  34. *@paramint$n加入布隆过滤器的key的数量
  35. *@returnint
  36. */
  37. publicstaticfunctiongetHashCount($m,$n){
  38. returnceil(($m/$n)*log(2));
  39. }
  40. /**
  41. *ConstructaninstanceoftheBloomfilter
  42. *
  43. *@paramint$mbit数组的宽度(bit数)Sizeofthebitarray
  44. *@paramint$khash函数的个数Numberofdifferenthashfunctionstouse
  45. */
  46. publicfunction__construct($m,$k){
  47. $this->m=$m;
  48. $this->k=$k;
  49. $this->n=0;
  50. /*Initializethebitset*/
  51. $this->bitset=array_fill(0,$this->m-1,false);
  52. }
  53. /**
  54. *FalsePositive的比率:f=(1–e-kn/m)k
  55. *Returnstheprobabilityforafalsepositivetooccur,giventhecurrentnumberofitemsinthefilter
  56. *
  57. *@returndouble
  58. */
  59. publicfunctiongetFalsePositiveProbability(){
  60. $exp=(-1*$this->k*$this->n)/$this->m;
  61. returnpow(1-exp($exp),$this->k);
  62. }
  63. /**
  64. *Addsanewitemtothefilter
  65. *
  66. *@parammixedEitherastringholdingasingleitemoranarrayof
  67. *stringholdingmultipleitems.Inthelattercase,all
  68. *itemsareaddedonebyoneinternally.
  69. */
  70. publicfunctionadd($key){
  71. if(is_array($key)){
  72. foreach($keyas$k){
  73. $this->add($k);
  74. }
  75. return;
  76. }
  77. $this->n++;
  78. foreach($this->getSlots($key)as$slot){
  79. $this->bitset[$slot]=true;
  80. }
  81. }
  82. /**
  83. *QueriestheBloomfilterforanelement
  84. *
  85. *IfthismethodreturnFALSE,itis100%certainthattheelementhas
  86. *notbeenaddedtothefilterbefore.Incontrast,ifTRUEisreturned,
  87. *theelement*may*havebeenaddedtothefilterpreviously.Howeverwith
  88. *aprobabilityindicatedbygetFalsePositiveProbability()theelementhas
  89. *notbeenaddedtothefilterwithcontains()stillreturningTRUE.
  90. *
  91. *@parammixedEitherastringholdingasingleitemoranarrayof
  92. *stringsholdingmultipleitems.Inthelattercasethe
  93. *methodreturnsTRUEifthefiltercontainsallitems.
  94. *@returnboolean
  95. */
  96. publicfunctioncontains($key){
  97. if(is_array($key)){
  98. foreach($keyas$k){
  99. if($this->contains($k)==false){
  100. returnfalse;
  101. }
  102. }
  103. returntrue;
  104. }
  105. foreach($this->getSlots($key)as$slot){
  106. if($this->bitset[$slot]==false){
  107. returnfalse;
  108. }
  109. }
  110. returntrue;
  111. }
  112. /**
  113. *Hashestheargumenttoanumberofpositionsinthebitsetandreturnsthepositions
  114. *
  115. *@paramstringItem
  116. *@returnarrayPositions
  117. */
  118. protectedfunctiongetSlots($key){
  119. $slots=array();
  120. $hash=self::getHashCode($key);
  121. mt_srand($hash);
  122. for($i=0;$i<$this->k;$i++){
  123. $slots[]=mt_rand(0,$this->m-1);
  124. }
  125. return$slots;
  126. }
  127. /**
  128. *使用CRC32产生一个32bit(位)的校验值。
  129. *由于CRC32产生校验值时源数据块的每一bit(位)都会被计算,所以数据块中即使只有一位发生了变化,也会得到不同的CRC32值。
  130. *Generatesanumerichashforthegivenstring
  131. *
  132. *RightnowtheCRC-32algorithmisused.Alternativelyonecoulde.g.
  133. *useAdlerdigestsormimickthebehaviourofJava'shashCode()method.
  134. *
  135. *@paramstringInputforwhichthehashshouldbecreated
  136. *@returnintNumerichash
  137. */
  138. protectedstaticfunctiongetHashCode($string){
  139. returncrc32($string);
  140. }
  141. }
  142. $items=array("firstitem","seconditem","thirditem");
  143. /*Addallitemswithonecalltoadd()andmakesurecontains()finds
  144. *themall.
  145. */
  146. $filter=newBloomFilter(100,BloomFilter::getHashCount(100,3));
  147. $filter->add($items);
  148. //var_dump($filter);exit;
  149. $items=array("firsttem","seconditem","thirditem");
  150. foreach($itemsas$item){
  151. var_dump(($filter->contains($item)));
  152. }
  153. /*Addallitemswithmultiplecallstoadd()andmakesurecontains()
  154. *findsthemall.
  155. */
  156. $filter=newBloomFilter(100,BloomFilter::getHashCount(100,3));
  157. foreach($itemsas$item){
  158. $filter->add($item);
  159. }
  160. $items=array("firsttem","seconditem","thirditem");
  161. foreach($itemsas$item){
  162. var_dump(($filter->contains($item)));
  163. }


问题实例】 给你A,B两个文件,各存放50亿条URL,每条URL占用64字节,内存限制是4G,让你找出A,B文件共同的URL。如果是三个乃至n个文件呢?

根据这个问题我们来计算下内存的占用,4G=2^32大概是40亿*8大概是340亿bit,n=50亿,如果按出错率0.01算需要的大概是650亿个bit。 现在可用的是340亿,相差并不多,这样可能会使出错率上升些。另外如果这些urlip是一一对应的,就可以转换成ip,则大大简单了。


http://blog.csdn.net/hguisu/article/details/7866173

你可能感兴趣的:(filter,bloom)