一，kmp算法的原理：
字符串匹配是计算机的基本任务之一。它所执行的任务就是在一个文本（较长的字符串）中查找是否包含着当前指定的模式字符串（较短的字符串），并找到其位置，这就是字符串匹配问题。注：在算法导论里面待搜索的词叫“模式”字符串，并用P来表示，本文有的地方叫他“待搜索词”。如无特殊说明数组下标以0为起始。举例来说，有一个长字符串"BBC ABCDAB ABCDABCDABDE"（即文本串），我想知道里面是否包含另一个模式字符串"ABCDABD"，它的位置在哪呢？
1，朴素匹配算法过程：
1）.

2，KMP匹配算法过程：
****1).

因为 6 - 2 等于4，所以将搜索词向后移动4位（之所以移动4，是因为搜索词的已匹配的前缀串“ABCDAB”的部分匹配值为2，即有“AB”对称，因为当前已匹配字符串的“AB”始终要找到下一个“AB”的开始位置，而这个开始位置恰好就在本字符串中的后缀中，所以直接移动4，我们可以最大减少匹配次数）。
4).

因为空格与Ｃ不匹配，搜索词还要继续往后移。这时，已匹配的字符数为2（"AB"），对应的"部分匹配值"为0。所以，移动位数 = 2 - 0，结果为 2，于是将搜索词向后移2位。
5).

逐位比较，直到发现C与D不匹配。于是，移动位数 = 6 - 2，继续将搜索词向后移动4位。
7).

3,关于前缀函数表
****1）前缀和后缀

3）.前缀函数表的意义

A,首相要搞清楚的是它是模式字符串的所有前缀产生的一张表，保存的值表征了最大对称度，我们一般用next数组来保存其某个前缀的对称值，例如next[6]=2，代表的就是模式串的某个有6个字符的前缀“ABCDAB”，这个前缀的对称度就是2，即“AB”是对称的。
B,模式串的某前缀的对称值的实质是，有时候，搜索词已经部分匹配了（一定是某个前缀），其某个前缀（已经匹配部分）的头部和尾部会有重复。比如，"ABCDAB"之中前缀有"AB"，后缀也有"AB"，那么它的"前缀函数对称值"就是2（"AB"的长度,且“AB”是"ABCDAB"所有前缀和后缀中的最长公有元素）。搜索词移动的时候，第一个"AB"向后移动4位（字符串长度-部分匹配值），就可以来到第二个"AB"的位置，在这个新的位置我们跳过了不必要的比较位置，并且直接就有“AB”匹配。

二，前缀函数表实现原理
　　通过上文完全可以对kmp算法的原理有个清晰的了解，那么下一步就是编程实现了，其中最重要的就是如何根据待匹配的模式字符串求出其所有前缀函数表中的最大对称值，在下文中我们将其存储在next数组中。

1，编程策略：
1）、当前字符的前面所有字符的对称程度为0的时候，只要将当前字符与前面这个子串的第一个字符进行比较。这个很好理解啊，前面所有字符串的对称值都是0，说明都不对称了，如果多加了一个字符，要对称的话只能是当前字符和前面字符串的第一个字符对称。比如“ABCDA”这个里面"ABCD"的最大对称值是0，那么后面的A的对称程度只需要看它是不是等于前面字符串的第一个字符A相等，如果相等就增加1，如果不相等那就保持不变，显然还是为0。

2）、按照这个推理，我们就可以总结一个规律，不仅前面是0呀，如果前面字符串的的最大对称值是1（k），那么我们就把当前字符与前面字符串的第二（k）个字符即P[1]（P[k]）进行比较，因为前面的是1（k），说明前面的字符已经和第一（k）个字符相等了，如果这个又与第二（k+1）个相等了，说明对称程度就是2（k+1）了。有两（k+1）个字符对称了。比如上面“ABCDA”的最大对称值是1，说明它只和第一个A对称了，接着我们就把下一个字符“B”与P[1]（即第二个字符）比较，又相等，自然对称程度就累加了，就是2了。
但是如果不相等呢？那么这个对称值显然要减少，并且我们只能到前面去寻找对称值，而在找的过程我们同时也利用前缀函数表快速搜索找到与当前字符匹配的位置。比如假设是“（AGCTAGC）（AGCTAGC）T”（请无视字符串中的括号，只为方便看出对称），显然最后一个T的前一个位置的对称度是7,说明T的前一个位置的7个字符的后缀必与7个字符的前缀相等，然而T!=P[7]，说明T位置的对称度只能是比7小的长度的前缀，所以递减k值，递减为多少呢？当前字符前一个位置的对称度为k=next[13]=7，显然必须以7为基准减少，即在前缀长度为7以内的范围重新寻找以T结尾的前缀，所以k=next[6]，再接着判断是否相等,

3）、按照上面的推理，我们总是在找当前字符P[q]（q为遍历到的位置下标，见下面程序）通过其前一个位置的对称值判断是否与P[k]相等，如果相等，那么加，如果不相等，那么就减少k值，重新寻找与与P[q]相等的元素位置

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
const int MAXN = 1000;
int next2[MAXN];

void  getnext(string str)
{
       memset(next2, 0, sizeof(next2));

       next2[0] = -1;
       int i = 0;
       int j = -1;
       while (i < str.size())
       {

              if (j == -1 || str[i] == str[j])
              {
                     i++;
                     j++;
                     next2[i] = j;

              }
              else
              {
                     j = next2[j]  ;
              }

       }

}

int kmp(string str, string substr)
{

       int i, j;
       i = -1;
       j = -1;
       int size = str.size();
       while (i < size)
       {
              if (j == -1 || str[i] == substr[j])
              {
                     i++;
                     j++;
                     if (j == substr.size())
                           return i - j;

              }
              else
                     j = next2[j];


       }
       return 0;
}

int main()
{
       string str;
       string substr;
       while (cin >> str >> substr)
       {

              getnext(substr);

              for (int i = 0; i < substr.size(); i++)
                     cout << next2[i] << " ";
                cout<<kmp(str, substr);
       }

       return 0;
}