【模板】KMP字符串匹配

题目描述

如题，给出两个字符串s1和s2，其中s2为s1的子串，求出s2在s1中所有出现的位置。

为了减少骗分的情况，接下来还要输出子串的前缀数组next。

（如果你不知道这是什么意思也不要问，去百度搜[kmp算法]学习一下就知道了。）

输入输出格式

输入格式：

第一行为一个字符串，即为s1（仅包含大写字母）

第二行为一个字符串，即为s2（仅包含大写字母）

输出格式：

若干行，每行包含一个整数，表示s2在s1中出现的位置

接下来1行，包括length(s2)个整数，表示前缀数组next[i]的值。

输入输出样例

输入样例#1：

ABABABC
ABA

输出样例#1：

1
3
0 0 1

说明

时空限制：1000ms,128M

数据规模：

设s1长度为N，s2长度为M

对于30%的数据：N<=15，M<=5

对于70%的数据：N<=10000，M<=100

对于100%的数据：N<=1000000，M<=1000

样例说明：

所以两个匹配位置为1和3，输出1、3

http://www.cnblogs.com/SYCstudio/p/7194315.html

http://blog.csdn.net/no1_terminator/article/details/52925547

http://blog.csdn.net/qq_30974369/article/details/74276186

KMP算法思维性很强，重点是理解next数组，提供以上3篇文章

就是要利用前面的比较，不要每一次都从开始判断，找到一个可以继续匹配的位置

根据条件next[j]=max(k)   (p₁p₂......p_k-1==p_j-k+1.....p_j-1)

就变成了自己匹配自己的问题

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 using namespace std;
 int nxt[],l1,l2;
 char s1[],s2[],a[],b[];
 void get_next()
 {int i,j;
     nxt[]=;
     j=;
      for (i=;i<=l2;i++)
      {
         while (j&&s2[i]!=s2[j+]) j=nxt[j];
         if (s2[i]==s2[j+]) j++;
         nxt[i]=j;
      }
 }
 void KMP()
 {int i,j;
     j=;
     for (i=;i<=l1;i++)
     {
         while (j&&s1[i]!=s2[j+]) j=nxt[j];
         if (s1[i]==s2[j+]) j++;
         if (j==l2)
         {
             printf("%d\n",i-l2+);
             j=nxt[j];
         }
     }
 }
 int main()
 {int i;
     scanf("%s%s",a,b);
     l1=strlen(a);l2=strlen(b);
     for (i=;i<=l1;i++)
     s1[i]=a[i-];
     for (i=;i<=l2;i++)
     s2[i]=b[i-];
      get_next();
      KMP();
     for (i=;i<=l2;i++)
     printf("%d ",nxt[i]);
     cout<<endl;
 }