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 while(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 }