首先,SimHash算法主要是用于文本去重的。文本去重的第一步就是判断文本的相似度,如果两个文本的相似度很
高,那么我们可以认为它们是相同的文本。
对于文本相似度的计算,传统的方法是使用向量空间模型,即Vector Space Model,即VSM,VSM计算文本相似
度的方法是这样的:先对文本进行分词,提取出特征词,然后建立文本向量,把相似度的计算转化成某种特征向量
距离的计算,比如余弦夹角,欧氏距离的计算等等。这样做的结果就是,如果文本向量的维度很高,那么计算的代
价将会很大,所以必须另辟蹊径。
可以看出,VSM计算文本的相似度是两两文本进行计算比较,对于小数量文本处理是可以的。但对于百度,Google
这样的搜索引擎,爬虫每天爬取的网页数目大得惊人,对于新爬取的网页,为了防止重复收录网页,需要对网页进行
判重处理,对于这样的海量数据,VSM无能为力。所以针对这种海量文本的去重,Google提出了SimHash算法,它
把高维度的文本向量映射成一个指纹,进行降维,从而减少计算量。
SimHash算法的详细描述如下
(1)输入一个N维的文本特征向量V,每个特征具有一定权重
(2)初始化一个C维的向量Q为0,C位的二进制签名S为0
(3)对于向量V的每一个特征,使用传统Hash算法计算出一个C位的散列值H。对1<=i<=C,如果H的第i位为1,
则Q的第i个元素加上该特征的权重;否则,Q的第i个元素减去该特征的权重。
(4)如果Q的第i个元素大于0,则S的第i位为1,否则为0
(5)返回这个C位的二进制签名S
下面以N = 3为例进行如图说明
对每篇文档根据SimHash算出签名S后,再计算两个签名的海明距离,海明距离等于这两个二进制数异或后1的个
数。根据经验,对于64位的SimHash,海明距离在3以内的都可以认为相似度比较高。
可以看出经过SimHash算法得到的是一个文档的指纹,我们根据这个指纹来判断文档的相似度,前面说过,如果海明
距离在3以内,那么就认为这两个文档基本是相同的。以64位的指纹来说,把它以每16位划分块,如果两个指纹海明
距离小于等于3,那说明至少有一个16位块完全相同(组合数学里的鸽巢原理)。但是我们不知道具体是哪两个块相
同,所以要进行枚举,这里只需要枚举4次。
以8位的数字指纹01100011来说明
我们划分为4块,每块两位,即01 10 00 11,然后在所有的8位二进制签名中分别查找以01,10,00,11开始的
签名。当然在这所有的8位二进制签名数据中,我们可以以前两位进行索引,比如把01000000和01111111放在一个
簇中,如果找不到以01,10,00,11开始的二进制签名,说明海明距离肯定大于3,这样就直接可以判断这两个文
本不相似,否则再比较后面的部分。
以上是8位二进制签名的比较方法,64位的签名类似方法处理。
假如现在有10亿个签名,约为2^34,那么每个块最多可能有2^(34-16)=262144个结果,那么4个块最终比较次数
为4*262144,大约为100万,这样原本需要比较10亿次,现在只需要比较100万次,可以看出计算量大大减少。
问题:有10亿个不重复的64位的01串,再另外给一个64位的01串str,如何快速从这10亿个串中找出与str的海明
距离小于等于3的字符串?
SimHash的Python实现
#coding:utf-8 class SimHash: def __init__(self, tokens='', hashbits=128): self.hashbits = hashbits self.hash = self.simhash(tokens) def __str__(self): return str(self.hash) def simhash(self, tokens): v = [0] * self.hashbits #等价于定义一个长度为hashbits的列表,列表中的元素全为0 for t in [self.string_hash(x) for x in tokens]: for i in range(self.hashbits): bitmask = 1 << i if t & bitmask: v[i] += 1 else: v[i] -= 1 fingerprint = 0 for i in range(self.hashbits): if v[i] >= 0: fingerprint += 1 << i return fingerprint def hamming_distance(self,other): x = (self.hash^other.hash ) & ((1 << self.hashbits) - 1) cnt = 0 while x: cnt += 1 x &= x - 1 return cnt def similarity(self,other): a = float(self.hash) b = float(other.hash) if a > b: return b / a else: return a / b def string_hash(self,source): if source == '': return 0 else: x = ord(source[0]) << 7 m = 1000003 mask = 2 ** self.hashbits - 1 for c in source: x = ((x * m) ^ ord(c)) & mask x ^= len(source) if x == -1: x = -2 return x if __name__ == '__main__': s = 'This is a test string for testing' hash1 = SimHash(s.split()) s = 'This is a test string for testing also' hash2 = SimHash(s.split()) s = 'Hello ,World!!!' hash3 = SimHash(s.split()) print hash1.hamming_distance(hash2),' ',hash1.similarity(hash2) print hash1.hamming_distance(hash3),' ',hash1.similarity(hash3)
原文地址:http://blog.csdn.net/acdreamers/article/details/27964475