再讲KMP算法

       真正理解KMP  

    今天我想再仔细讲一下这个算法，因为通过不断的做题，我发现，这个算法真的很有用。因为我是主要搞图论的，所以一般碰到字符串的题目都是把题意讲给队友听，然后去看其他题，可是这也凸显了一个弊端，就是有些水题我做不出来，原因很简单，我自己不会写KMP。自己博客有个字符串匹配的专题，我看了一下是之前转发的别人的，而且觉得讲的也不是很好，所以趁着自己搞懂了，今天我就自己写一下对KMP的理解，希望那些不知道的同学也可以有所收获。

    字符串匹配，假设我们有两个串s和t,长度分别是n和m，判断t在s中首次出现的位置，也就是匹配了。最暴力的方法就是一个一个比较，遇到不匹配的字符，整个t串就往前移一位接着喝s串匹配，这种是效率很低的，复杂度是O(n*m)。这种方法有个缺陷那就是，当匹配到t串中的第i位时，没有利用到t串中第i为之前的字符都已经匹配完了这个条件，我们完全可以利用这个条件来减少算法的复杂度的。怎么减少呢？很简单就是尽可能的往前多移动几位呗，怎么实现这个功能，假设，当t串中匹配到第i位不能再和s串中相应的字符匹配了，但是在第i位之前的t的子串中，有一个最长的前缀和最长的后缀是一样的，那么我们是不是可以直接把最长的后缀作为下次匹配的前缀？假设最长的公共子串（相同的最长前缀和最长后缀）的长度是k。现在匹配到t串的第i-1位，第i位那么我们可以往前移动i-k位，串的下坐标是从0开始的。怎么样，是不是很神奇！KMP的复杂度是O(n+m)的，因此，效率会很高。

    当然，这只是空讲，懂了理论，当然还是要写出代码，这才叫完美。在贴出代码之前我们约定，一个next数组，next[i]表示，在串t的第i位之前（也就是从0到i-1）这个子串中，前缀和后缀的公共部分最长的长度。那么怎么算next[i+1]呢？很简单，如果t[i]和next[i]相等的话（t数组的下坐标是从0开始的）,那么next[i+1] = next[i]+1。如果不相等，那么就判断t[i]是不是和next[next[i]]相等，直到next[k]等于0。

    现在贴出代码：

void make_next(int len)//len是字符串s的长度
{
    int next1 = 0;
    next[0] = next[1] = 0;

    for(int i = 1; i < len; i++)
    {
        while((next1>0)&&(s[i]!=s[next1])) next1 = next[next1];
        if(s[i]==s[next1]) next1++;
        next[i+1] = next1;//更新next[i+1]
    }

    return ;
}

    那么解决了next数组之后，怎么用这个next数组去匹配母串呢？其实和求next数组一样的道理，如果和母串的第i为匹配的话，那么i++，匹配串的位置next1也++，如果不匹配呢，那也好办，就看s[next[next1]]是否匹配，如果也不匹配就继续分，直到next1等于0，就开始从匹配串的初始位置匹配了。贴出代码：

//寻找s是否是s1的子串
void Search()
{
    int len,len1; int next1 = 0;
    len = strlen(s); len1 = strlen(s1);

    for(int i = 0; i < len1; i++)
    {
        while((s1[i]!=s[next1]) && (next1>0)) next1 = next[next1];
        if(s1[i]==s[next1]) next1++;
        if(next1 == len)
        {
            num++;
            cout<<"在第"<<i-len+1<<"个位置发现了匹配子串"<<endl;
            next1 = next[next1];
        }
    }
    cout<<"共"<<num<<"个子串"<<endl;

    return ;
}