字符串匹配之KMP算法

老规矩，直接上代码：

#include <string>  
#include <stack>  
#include <vector>  
#include <map>  
#include <algorithm>  
#include <iostream>  
using namespace std;  
vector<int> makeNext(string P)  
{  
    int q, k;//q:模版字符串下标；k:最大前后缀长度  
    int m = P.size();//模版字符串长度  
    vector<int> next(m, 0);  
    for (q = 1, k = 0; q < m; ++q)//for循环，从第二个字符开始，依次计算每一个字符对应的next值  
    {  
        while (k > 0 && P[q] != P[k])//递归的求出P[0]···P[q]的最大的相同的前后缀长度k  
            k = next[k - 1];          //不理解没关系看下面的分析，这个while循环是整段代码的精髓所在，确实不好理解    
        if (P[q] == P[k])//如果相等，那么最大相同前后缀长度加1  
        {  
            k++;  
        }  
        next[q] = k;  
    }  
    return next;  
}  
void kmp(string T,string P)  
{  
    int n = T.size();  
    int m = P.size();  
    vector<int> next=makeNext(P);  
    for (int i = 0, q = 0; i < n; ++i)  
    {  
        while (q > 0 && P[q] != T[i])  
            q = next[q - 1];  
        if (P[q] == T[i])  
        {  
            q++;  
        }  
        if (q == m)  
        {  
            printf("Pattern occurs with shift:%d\n", (i - m + 1));  
            q = next[q-1]; //look for next match  
        }          
    }  
}  
  
int main()  
{  
    string T = "ababxbababcdabdfdsabcdabdssabcdabd";  
    string P = "abcdabd";  
    cout << "T:" << T << std::endl;  
    cout << "P:" << P << std::endl;  
    kmp(T, P);  
    return 0;  
}