kmp算法

关于字符串匹配问题

比如A,B串，现在要在A中找B

扫描A串，并更新能匹配到B的什么位置

关于kmp算法如何工作，有经典说法  神犇解说

---

摘选

我们用两个指针i和j分别表示，A[i-j+ 1..i]与B[1..j]完全相等。也就是说，i是不断增加的，随着i的增加j相应地变化，且j满足以A[i]结尾的长度为j的字符串正好匹配B串的前 j个字符（j当然越大越好），现在需要检验A[i+1]和B[j+1]的关系。当A[i+1]=B[j+1]时，i和j各加一；什么时候j=m了，我们就说B是A的子串（B串已经整完了），并且可以根据这时的i值算出匹配的位置。当A[i+1]<>B[j+1]，KMP的策略是调整j的位置（减小j值）使得A[i-j+1..i]与B[1..j]保持匹配且新的B[j+1]恰好与A[i+1]匹配（从而使得i和j能继续增加）。

#include <iostream>
#include <cstring>
using namespace std;
 const int N=1e6+6;
 int n,m;
 int pos[N];
 char a[N],b[N];
 
 void init(){
 	int i,j=0;
 	pos[1]=0; 
 	for(i=1;i<m;i++){
 		while(j>0&&b[i+1]!=b[j+1]) j=pos[j];
 		
 		if(b[i+1]==b[j+1]) j++;
		pos[i+1]=j;
	 }
 }
 void kmp(){
 	int i,j=0;
 	
 	for(i=0;i<n;i++){
 		while(j>0&&a[i+1]!=b[j+1]) j=pos[j];
 		
 		if(a[i+1]==b[j+1]) j++;
 		if(j==m){
 			printf("%d\n",i-j+2); 
 			j=pos[j];
		 }
	 }
 }
 int main(){
 	cin>>a+1>>b+1; 
 	n=strlen(a+1); m=strlen(b+1);
 	init();
 	kmp();
 	for(int i=1;i<=m;i++) cout<<pos[i]<<' ';
 }