数学建模社区-数学中国

标题: 透彻理解Shift-And/Shift-Or 算法 字符串匹配/子串查找 [打印本页]

作者: 美人如花    时间: 2015-4-13 10:06
标题: 透彻理解Shift-And/Shift-Or 算法 字符串匹配/子串查找
Shift-And/Shift-Or 算法和KMP 算法一样,也是线性时间复杂度的字符串匹配算法,运行时间上甚至要比KMP 算法快得多。而理解上比KMP 算法更容易一些。Shift-And/Shift-Or 算法设计的非常巧妙,初次接触时同样“吓了一跳”。
; G( G0 x) D3 H% }6 yShift-And 与 Shift-Or 算法的原理完全一样,区别仅在于Shift-Or 对Shift-And 做了一点儿改进。我们先说Shift-And 算法。
9 ~' n9 f: P" e+ Z1 _
+ d: c' k8 P3 @+ t# n7 Z* R3 l0 CS 表示原字符串,T 表示目标串(模式串),我们要在S 中搜索T。* u! o( m0 l: @' @! \
令 S[0..m-1] = abcabcabdabba, T[0..n-1] = abcabd7 s5 Z; N& f% E# c, W
+ k8 V& S; w( C6 U8 D
1,Shift-And 算法思想& S) Y* ]# S# {# e
Shift-And 算法的核心思想是利用掩码D 来记录模式串的前缀匹配情况。(瞧,shift 算法的核心也是前缀匹配)。Shift 算法大量应用了位运算。# ?1 d% Z3 u# j* |: U" @4 K
D 是一个m 位的无符号整数:D[n-1, n-2, ..,1,0] (注意D 并不是一个数组,仅仅是一个整数,D[n-1] 表示其最高位bit)。
3 N. O$ c6 T9 s/ ?1 ]1 w数组索引i 控制S 串的扫描,当扫描的字符S 时,D 的第j 位D[j] = 1 当且仅当T[0..j] 是S[0..i] 的一个后缀。
5 i* [( ~( D7 s: {! }0 ^- o. o) S7 R" H! r% }7 {& `
要使用Shift 算法,需要一个辅助表B。B 是一个字典,key 是问题域字符集中的每个字符,value 是一个n 位无符号整数,记录该字符在模式串T 的哪些位置出现。
/ t$ u& X+ o) F( y. \" [6 a; n+ B例如,字符c 在T[2] 处出现,那么B['c'] = 000100 (对于字符串,低位在左;对于B['c'],低位在右);同理,a 在T[0],T[3] 处出现,B['a'] = 001001.9 H# a+ ?- k3 X+ Y0 M# {
( I. R; l8 z! D, c  o" a: a
假设当前处理到S,需要对D 进行更新。由于D[j] (0<j<n) 标识T[0..j] 是否是S[0..i] 的后缀,所以D[j]=1 当且仅当更新前的D[j-1]=1 并且S==T[j];D[0] 是边界情况,D[0]=1 当且仅当S==T[0]。6 n# b0 d7 y; Y
所以,D = (D << 1 | 1) & B[S[i] ;
& g& i: G, [  j# U2 D2 A显然,当D[n-1]=1 时,表示T[0..n-1] 是S[0..i] 的后缀,此时找到一个T的完全匹配。/ U( g& Q* c& b$ p8 ~- L& P

7 Y' X: a0 u6 q6 A) R# o, h2,Shift-And 算法实现
; x3 [* ]9 Y  q' a" @Shift-And 匹配过程代码:3 f' u( _* x+ f9 o8 t; u6 z3 g

  e" P0 X# h: W4 h9 k& s. q  c" C9 G% R
4 h& M6 g- q. s3 w$ g' P5 E6 g% G$ k由于位运算在计算机中可以并行进行,每次循环的执行是常数时间的,所以上面代码段的复杂度是 O(m)。0 ~: e3 j# |( z0 C; T

$ O$ t- s4 h9 O' S; q3,辅助表 B
1 w& P. |. B& W上面没有提到如何得到辅助表B。很简单,只要获得模式串T 中每个字符出现的位置。8 F& h# L8 u% i( H
, d2 Q+ g; c& }" ]0 ^' d' j
3 u; _0 R& D5 r
显然,上述代码段的复杂度是 O(n)。Shift-And 算法的时间复杂度是O(m+n)。
# I+ F! D" q5 p$ j; v2 z2 ~实际上,shift 算法通常比KMP 算法的匹配速度要快,因为计算机位并行运算是非常高效的。6 n- F' ^8 b1 B& e- ]- O
7 m! \: U& d7 A0 A+ \
注意:数组B 的大小是由字符集决定的,如果字符来自ASCII 码,字符的数值范围是0~127,数组大小是128 即可;否则,可能需要更大的数组B,或者自己构建字符到整数索引之间的散列关系。% P# b0 r4 [5 w7 L3 S$ n/ ?# b* x* u

2 f6 L: l- E  i+ p, [, e4,Shift-Or 算法6 x$ F  ]0 |/ ^' z, }3 |
在Shift-And 中,对掩码D 的更新:D = (D << 1 | 1) & B[S[i] ;
# X+ Y3 F6 x* W  N每次更新D 都需要额外进行D 移位后与"1" 的"或"运算。这是由于我们要保证当字符S 在T[0] 处出现时,D[0] 一定要等于1,而D 向左移位后最低位是0。
6 w. `$ c4 d# p
* z  \! R% I& ~% y0 [如果将Shift-And 中核心的“与” 运算改为“或” 运算,可以节省这一个附加的“或1” 运算。这正是Shift-Or 所改进的地方。6 @3 K  `! |1 a3 ~2 n
Shift-Or 与Shift-And 的唯一区别在于,在Shift-Or 中,“有效位” 是通过0(而不是1)来标识。' n/ o  |5 x2 u4 Q2 ?% j
于是求解辅助表B 和更新掩码D 都会与Shift-And 有一些区别,详见代码。! E/ O- R+ |  p. ^8 ?* Z: K- L
! ]2 U+ M' v8 h7 R1 x7 \2 o0 R
Shift-And 完整代码:C++ 实现  Python 实现) ]/ c* O: u8 g1 J- E7 g; \
Shift-Or 完整代码:C++ 实现  Python 实现2 i7 s# R8 F5 g0 D
8 H+ u# }3 D: m4 F) V% A/ @9 f





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5