学点字符串匹配——shiftAnd/shiftOr

　　这个算法利用了位运算的优势速度很快，易于实现，缺点是模式串不能太长。据说在适用范围内速度是kmp的至少2倍。算法思想是用一个数字prefix的二进制去表示模式串的前缀，prefix的二进制第i位表示模式串的前缀0～i。如果这个前缀i是当前搜索到的匹配串的后缀时，将第i位至1。算法就是从匹配串的第0个字符到最后一个字符不断更新prefix的值每次更新后看最长的前缀（即模式串）的位是否被至1，如果是说明得到了一个匹配。

　　下面说明下更新的过程。假如求出了比较匹配串第j位时的prefix，现在求比较j+1位的情况：prefix第i位为1当且仅当prefix第i-1位为1（匹配了i-1个字符），并且第i位所在字符和匹配串第j+1个字符相等。按下面实现来说，每次循环将prefix左移一位，第0位至一，在&上一个代表第i个字符是否匹配模式串某位的数字（这个通过初始化得到），就得到更新后的prefix。

　　shiftAnd和shiftOr的思想一致。只是实现时存的是反码，因为prefix左移自然引入第0位的是0，为了减去或运算，才引入反码表示。

上代码～

//return a hashtable 代表字母表中的各个元素是否能成为模式串的后缀，用64位位压缩表示
long long * initAlphbat(char *patternString)
{
    long long *hash = (long long *) malloc(256 * sizeof(long long));
    unsigned int idx = 0;
    int psLen = strlen(patternString);

    if(hash == NULL) return NULL;

    memset(hash, 0, 256 * sizeof(long long));
    for(idx = 0; idx < psLen; ++idx)
    {
        hash[(unsigned int)patternString[idx]] |= ((long long)1 << idx); //表示patternString[idx]代表字符可以是在模式串的哪些位置
    }

    return hash;
}

//the patternString length must be less than 64
MatchingInfo * ShiftAnd(char *patternString, char *dataString, unsigned int dsLen)
{
    unsigned int idxDs = 0, psLen = strlen(patternString);
    long long *hash = initAlphbat(patternString);

    long long prefix = 0;  //每一位表示一个模式串的前缀，如果该位是1则表示当前该前缀被匹配上

    if(hash == NULL) return NULL;

    //匹配信息的初始化
    MatchingInfo *mi = (MatchingInfo *)malloc(sizeof(MatchingInfo));
    if(mi == NULL) return NULL;
    mi->indexOfMatch = (unsigned int *)malloc(sizeof(unsigned int));
    mi->numOfMatch = 0;
    mi->allocLen = 1;

    //下面的循环中动态更新prefix
    for(idxDs = 0; idxDs < dsLen; ++idxDs)
    {
        prefix = (prefix << 1) | 1;
        prefix &= hash[(unsigned int)dataString[idxDs]];
        if(prefix & ((long long)1 << (psLen - 1)))
        {
            //匹配上一个模式串
            if(mi->allocLen > mi->numOfMatch)
            {
                //记录下标的内存够用
                mi->indexOfMatch[mi->numOfMatch++] = idxDs - psLen + 1;
            }
            else
            {
                //扩张匹配下标的内存
                mi->indexOfMatch = (unsigned int *)realloc(mi->indexOfMatch, (mi->allocLen << 1) * sizeof(unsigned int));
                mi->allocLen <<= 1;
                mi->indexOfMatch[mi->numOfMatch++] = idxDs - psLen + 1;
            }
        }
    }

    free(hash);
    hash = NULL;
    return mi;
}