kmp算法详解

 1 #include<iostream>//***28&&16行相当于递归算法****** 
 2 #include<cstdio>
 3 #include<cstring>
 4 using namespace std;
 5 const int MAXN=1000001;
 6 char a[MAXN],b[MAXN];
 7 int la,lb;
 8 int p[MAXN];
 9 void makep()
10 {
11     int j=0;// 记录最长公共前后缀的长度；
12     p[0]=0;// 第一个数的前缀和后缀都是空集，故为0 
13     for(int i=1;i<lb;i++)// 第一个是0，所以从1开始，半闭半开区间 
14     {
15         while(j>0&&b[j]!=b[i])//当出现不匹配的情况的的话，就让j从p数组中取一个能继续匹配的元素 
16         j=p[j-1];// 因为对j的最后一次处理是j++，而p是从0开始，故j=p[j-1]
17         if(b[j]==b[i])//如果相同的话就让j++ 自己可以举例子试试 
18         j++;
19         p[i]=j;// 将每一个元素的最长公共前后缀的长度记录下来，方便以后处理 
20     }
21 }
22 void KMP()
23 {
24     int j=0;// 已经匹配的最长的长度 
25     for(int i=0;i<la;i++)
26     {
27         if(j>0&&b[j]!=a[i])// 当将要匹配的元素不想同且找到相同的匹配元素的时候 
28         j=p[j-1];// j==0的时候p[-1]不存在，所以j>0 
29         if(b[j]==a[i])// 同理，当将要匹配的元素相同的时候，j++ 
30         j++;
31         if(j==lb)// 找到一个匹配 
32             printf("%d\n",i+1-lb+1);
33         // 因为i，lb都是从0开始计数的，so+1 
34     }
35 }
36 int main()
37 {
38     scanf("%s %s",a,b);
39     la=strlen(a);lb=strlen(b);//统计两个字符串的长度 
40     makep();// 构造p数组 
41     KMP();//
42     for(int i=0;i<lb;i++)printf("%d ",p[i]);
43     return 0;
44 }