字符串匹配之KMP算法

KMP算法使用前缀函数来模拟有限自动机的后缀函数，前缀函数通过计算模式与其自身的偏移匹配的信息，本身的证明很复杂，关键在于弄懂其核心思想，下面就不赘述了，仅仅贴出代码：

#include <iostream>
#include <string>
using namespace std;

int const  MAX_N = 1000;
int pi[MAX_N];

void COMPUTE_PREFIX_FUNCTION(string& P)
{
    int m = P.length() - 1;
    pi[1] = 0;
    int k = 0;
    for(int q = 2; q <= m; q++)
    {
        while(k > 0 && P[k + 1] != P[q])
            k = pi[k];
        if(P[k + 1] == P[q])
            k++;
        pi[q] = k;
    }
}

void KMP_MATCHER(string& T, string& P)
{
    int n = T.length() - 1;
    int m = P.length() - 1;
    int q = 0;
    for(int i = 1; i <= n; i++)
    {
        while(q > 0 && P[q + 1] != T[i])
            q = pi[q];
        if(P[q + 1] == T[i])
            q++;
        if(q == m)
        {
            cout<<"Pattern occurs with shift "<<i - m<<endl;
            q = pi[q];
        }
    }
}

int main()
{
    //为了和书上保持一致，这里在字符串前面添加了一个空格，使得真正有效的字符串的起始下标为1
    string T = " abababacaba";
    string P = " ababaca";
    COMPUTE_PREFIX_FUNCTION(P);
    KMP_MATCHER(T,P);
    return 0;
}