KMP算法/POJ3461(求匹配字符子串的个数)

1.Next数组

next[j]={

j=1时，0

j!=1时，最大子串长度+1，即P1P2….Pk-1==Pj-k+1….Pj-1,1<k<j,(看来j=2时是没有最大子串的,因此next[2]一定为1)

其它情况，1

}

next[j]=0说明模式串与匹配串都要后移一位

next[j]！=0说明匹配串的当前字符与next[j]匹配

next数组的下标从1开始，因此我们在模式串（T）的前面添加一个字符‘#’，初始化next[1]=0

求每个next[j](j>1)时实际上利用到上next[j-1]的值

如求next[5]时，则查看next[4]的值，如果T[4]=T[next[4]]，则next[5]=next[4]+1

如果T[4]!=T[next[4]]，T[4]再与next[T[next[4]]]比，如果相等，则next[5]=next[T[next[4]]]+1,如此继续，如果它的next值为0，则跳到其它情况

因此，标准算法如下：

void get_next()
{
	unsigned int length = strlen(t);
	next[0]=0;
	next[1]=0;
	unsigned int i=1,j=next[1];
	while(i<length-1){
		if(j==0){
			i++;
			j=1;
			next[i]=1;
			continue;
		}
		if(t[i]==t[j]){
				i++;
				j++;
				next[i]=j;
		}else 
			j = next[j];
	}
	
}

这段代码优化一下：

scanf("%s",t+1);

t[0]=’#’;

void get_next()
{
	unsigned int length = strlen(t);
	//next[0]=0;
	next[1]=0;
	unsigned int i=1,j=next[1];
	while(i<length-1){
		
		if(j==0||t[i]==t[j]){
				i++;
				j++;
				next[i]=j;
		}else 
			j = next[j];
	}
	
}

对于POJ3461,next数组要多求一位

要求：求出上串在下串中出现的次数

Sample Input
3
BAPC
BAPC
AZA
AZAZAZA
VERDI
AVERDXIVYERDIAN

Sample Output
1
3
0

#include <iostream>
using namespace std;
#define MAX 1000003
unsigned int next[MAX];
char s[MAX];
char t[MAX];
void get_next()
{
	unsigned int length = strlen(t);
	next[0]=0;
	next[1]=0;
	unsigned int i=1,j=next[1];
	while(i<length){
		
		if(j==0||t[i]==t[j]){
				i++;
				j++;
				next[i]=j;
		}else 
			j = next[j];
	}
	
}
unsigned int kmp(){
	unsigned int length = strlen(s)-1;
	unsigned int matLen = strlen(t)-1;
	if(length<matLen)return 0;
	unsigned int i=1,j=1,count = 0;
	while(i<=length){
		
		if(j==0||s[i]==t[j]){
			i++;
			j++;
			
		}
		else j = next[j];
		if(j==matLen+1){
			count++;
			j = next[j]; 
			}

	}
	return count;
}

int main()
{
	freopen("i://in.txt","r",stdin);
	int count;
	scanf("%d",&count);
	t[0]='#';
	s[0]='#';
	for(int i=0;i<count;i++){
	scanf("%s",t+1);
	scanf("%s",s+1);
	get_next();
	printf("%d\n",kmp());
	
	}
	return 0;
}