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!!!仔细斟酌,感觉这两种写法是换汤不换药,故删除后者。