KMP算法

ACM模版

KMP算法

KMP_Pre

/*
 * next[]的含义，x[i - next[i]...i - 1] = x[0...next[i] - 1]
 * next[i]为满足x[i - z...i - 1] = x[0...z - 1]的最大z值（就是x的自身匹配）
 */

void KMP_Pre(char x[], int m, int next[])
{
    int i, j;
    j = next[0] = -1;
    i = 0;
    while (i < m)
    {
        while (-1 != j && x[i] != x[j])
        {
            j = next[j];
        }
        next[++i] = ++j;
    }
    return ;
}

preKMP

/*
 * kmpNext[]的意思：next'[i] = next[next[...[next[i]]]]
 * （直到next'[i] < 0或者x[next'[i]] != x[i]）
 * 这样的预处理可以快一些 
 */

void preKMP(char x[], int m, int kmpNext[])
{
    int i, j;
    j = kmpNext[0] = -1;
    i = 0;
    while (i < m)
    {
        while (-1 != j && x[i] != x[j])
        {
            j = kmpNext[j];
        }
        if (x[++i] == x[++j])
        {
            kmpNext[i] = kmpNext[j];
        }
        else
        {
            kmpNext[i] = j;
        }
    }
    return ;
}

KMP_Count

/*
 * 此函数与上述两个函数中的任意一个搭配使用（即调用上述两个函数中的任意一个）
 * 返回x在y中出现的次数，可以重叠
 */

int next[10010];

int KMP_Count(char x[], int m, char y[], int n)
{
    //  x是模式串，y是主串
    int i, j;
    int ans = 0;
    //  preKMP(x, m, next);
    KMP_Pre(x, m, next);
    i = j = 0;
    while (i < n)
    {
        while (-1 != j && y[i] != x[j])
        {
            j = next[j];
        }
        i++, j++;
        if (j >= m)
        {
            ans++;
            j = next[j];
        }
    }
    return ans;
}

---

修正于2017.3.23 常规KMP indexKMP部分14行修正!!!S[i] == T[i] X => S[i] == T[j]

感谢评论区 xiaonuo911teamo 大神的悉心指导！！！

---

修正于2017.4.20 删除常规KMP！！！仔细斟酌，感觉这两种写法是换汤不换药，故删除后者。