QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3369|回复: 0
打印 上一主题 下一主题

[其他资源] 透彻理解Shift-And/Shift-Or 算法 字符串匹配/子串查找

[复制链接]
字体大小: 正常 放大

14

主题

10

听众

43

积分

升级  40%

  • TA的每日心情
    慵懒
    2015-5-5 09:46
  • 签到天数: 10 天

    [LV.3]偶尔看看II

    自我介绍
    跳转到指定楼层
    1#
    发表于 2015-4-13 10:06 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    Shift-And/Shift-Or 算法和KMP 算法一样,也是线性时间复杂度的字符串匹配算法,运行时间上甚至要比KMP 算法快得多。而理解上比KMP 算法更容易一些。Shift-And/Shift-Or 算法设计的非常巧妙,初次接触时同样“吓了一跳”。4 A) w2 {* A& [% Q! S
    Shift-And 与 Shift-Or 算法的原理完全一样,区别仅在于Shift-Or 对Shift-And 做了一点儿改进。我们先说Shift-And 算法。
    + U% D# `* [# k6 v, i/ E+ Q  A3 ?
    % {. X* R' S  U2 OS 表示原字符串,T 表示目标串(模式串),我们要在S 中搜索T。
    2 X, E' ?0 B! [( M, u( q6 [9 G令 S[0..m-1] = abcabcabdabba, T[0..n-1] = abcabd
    , c) r% [) o3 O- Z, K" T: ]
    7 {- ^2 u$ k/ m; y; W1,Shift-And 算法思想9 A' `, Z) F3 [" V$ I9 l+ K- F
    Shift-And 算法的核心思想是利用掩码D 来记录模式串的前缀匹配情况。(瞧,shift 算法的核心也是前缀匹配)。Shift 算法大量应用了位运算。
    1 t6 c$ a: E( _& ]D 是一个m 位的无符号整数:D[n-1, n-2, ..,1,0] (注意D 并不是一个数组,仅仅是一个整数,D[n-1] 表示其最高位bit)。% X) f, f! ?2 V7 |9 ^6 X' s- C
    数组索引i 控制S 串的扫描,当扫描的字符S 时,D 的第j 位D[j] = 1 当且仅当T[0..j] 是S[0..i] 的一个后缀。
    / ^4 o* F3 j& R* V, U
    0 E* E/ q% V$ c4 h" q要使用Shift 算法,需要一个辅助表B。B 是一个字典,key 是问题域字符集中的每个字符,value 是一个n 位无符号整数,记录该字符在模式串T 的哪些位置出现。
    ! B) t  m. M, A; U/ l例如,字符c 在T[2] 处出现,那么B['c'] = 000100 (对于字符串,低位在左;对于B['c'],低位在右);同理,a 在T[0],T[3] 处出现,B['a'] = 001001.( G' h8 J; Z& J1 u
    9 R6 A, ^: v* A3 E' L8 J% c4 s5 R8 b
    假设当前处理到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]。
    ' |; t. z8 M. X! p" f4 R所以,D = (D << 1 | 1) & B[S[i] ;
    4 B6 G. }" @7 n+ j% l显然,当D[n-1]=1 时,表示T[0..n-1] 是S[0..i] 的后缀,此时找到一个T的完全匹配。  T; `$ n% h( F1 n* T
    0 j1 ]: C+ z! f% [& {% [
    2,Shift-And 算法实现
    % p" o0 k+ M0 {8 _3 nShift-And 匹配过程代码:
    & I' r" k& L  F/ i/ A( K, V$ D
    + b, E2 [/ F# g& T
    / W# V& p1 q, B; b" m2 P2 q+ P
    由于位运算在计算机中可以并行进行,每次循环的执行是常数时间的,所以上面代码段的复杂度是 O(m)。
    ' ^2 ]6 z* ~. b( j% D$ ?4 r- C* F, y: u8 z! T. W5 k! U
    3,辅助表 B
    ) n! T9 E% L. X  ?  ?6 O上面没有提到如何得到辅助表B。很简单,只要获得模式串T 中每个字符出现的位置。
    & M# s9 J* E6 t' o$ Q

    5 ^7 a" u: |$ s5 P' V( \8 ?" G) R$ h' D$ z* y$ p
    显然,上述代码段的复杂度是 O(n)。Shift-And 算法的时间复杂度是O(m+n)。5 K& M) _" a' G! X/ |( Q3 f  `
    实际上,shift 算法通常比KMP 算法的匹配速度要快,因为计算机位并行运算是非常高效的。
    : a2 W/ ~8 m3 l% g" L  p; t
    % i& S$ @& W/ F. a, n; @, |注意:数组B 的大小是由字符集决定的,如果字符来自ASCII 码,字符的数值范围是0~127,数组大小是128 即可;否则,可能需要更大的数组B,或者自己构建字符到整数索引之间的散列关系。
    6 N6 O& R7 M: h* z3 Q, V( ^" g- i! \  a" C. n5 ?/ j) ^9 Y
    4,Shift-Or 算法5 s# y$ i6 n* {4 q9 l+ a9 V9 j% u
    在Shift-And 中,对掩码D 的更新:D = (D << 1 | 1) & B[S[i] ;
    ; D! }4 g, L9 I每次更新D 都需要额外进行D 移位后与"1" 的"或"运算。这是由于我们要保证当字符S 在T[0] 处出现时,D[0] 一定要等于1,而D 向左移位后最低位是0。
      m; q+ K8 X' ^' w) z  Q, E) D8 r$ ^* y
    如果将Shift-And 中核心的“与” 运算改为“或” 运算,可以节省这一个附加的“或1” 运算。这正是Shift-Or 所改进的地方。0 h9 B4 t0 v, j1 x' D
    Shift-Or 与Shift-And 的唯一区别在于,在Shift-Or 中,“有效位” 是通过0(而不是1)来标识。
    4 {8 q# R/ M: }1 o5 v6 V) `于是求解辅助表B 和更新掩码D 都会与Shift-And 有一些区别,详见代码。' t$ X4 a* F( r/ j! d4 Q
    ) s: {( L. g6 j* j4 G" g
    Shift-And 完整代码:C++ 实现  Python 实现3 E1 g3 p$ \6 H# J  w' s
    Shift-Or 完整代码:C++ 实现  Python 实现
    & Y. d$ z7 Z- z" [$ G! ^
    2 F. N; e# L2 a( K: ]/ b) ?! b6 L
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-15 12:03 , Processed in 0.572573 second(s), 54 queries .

    回顶部