Poj KMP 总结

 题目列表：

Poj1226 http://poj.org/problem?id=1226

Poj1961 http://poj.org/problem?id=1961

Poj2185 http://poj.org/problem?id=2185

Poj2406 http://poj.org/problem?id=2406

Poj2541 http://poj.org/problem?id=2541

Poj2752 http://poj.org/problem?id=2752

Poj3080 http://poj.org/problem?id=3080

Poj3167 http://poj.org/problem?id=3167

Poj3450 http://poj.org/problem?id=3450

Poj3461 http://poj.org/problem?id=3461

 

这里转有一篇介绍KMP算法的比较好的文章：

http://blog.csdn.net/dooder_daodao/article/details/6299504

 

下面就各题进行分析：

1226 求最长公共子串，并且这个公共子串可以是某个字符串反串的子串，可以找最小的一个字符串，

         然后枚举其所有的子串，并在其它串中查找，查找时用KMP或者strstr函数；

 

1961 这题考察对next数组的理解 next数组中的j本质上维护了从距离i最近的上一个重复串中与i对应的位置

          因此如果i % (i - j) == 0 && i / (i - j ) > 1则这个刚好扫描到了重复的串 

 

2185  这题就不多说了，各种错误的代码都AC了，不过相较而言，Discuss里 Iamjw 讲了一段比较好：

           比如：2 8

           aaabcaaa

           abababab

           这个思路会给出16吧？但是答案是12。

           因为横排的计算出错了~其实是因为第一个串符合题意的长度不止5，还有6,7,8。所以事实上应该把每个串符合

           题意的长度都记录下来，然后选一个公共最小的。比如第一个串是 5,6,7,8第二个串是2,4,6,8所以横排应该选6。

           这里只需解决两个问题：

           1) 如何找到串s的所有符合题意的长度？ 只需要对s进行一次kmp特征值计算，如果发

           现后面有多余的后缀，就把后缀再拿来做一次kmp，不断重复，直到没有后缀。这个过程就可以记录所有符合题

           意的串了。

           比如aaabcaaa，先求一遍kmp，得到5，再对后缀aaa求一遍，得到1(即6,7,8)。

           又比如aabaabaabaa，先求一遍kmp，得到3(即3,6,9)，再对后缀aa求一遍，得到1，(即10,11)。

           2)

           如何找到所有串公共的最短符合题意的长度？ 直接用(列数)c个桶来统计一下就行了~ 如果某个长度出现了(行数)r次，

           那么它就是公共的。最后从小到大遍历一遍就求出了最小长度。

 

2406  还是next数组的性质,if(len%(len-next[len])==0),则重复子串的长度为 len-next[len].

2541 这一题可直接用KMP来做的，反方向取后十三个字符做为模式串，向前进行匹配

 

2752 同 2406

 

3080 枚举子串，逐个匹配。

 

3167 比较不同的一个KMP题目，而且比较好的一个题，参照了下面这篇博客，讲的很细致：

          同样对模式串求next数组，并据next进行匹配，不同处在于使用了不同的比较方式。

          http://www.cppblog.com/zxb/archive/2010/10/06/128782.aspx?opt=admin

 

3450 同3080，甚至就修改几行代码

 

3461 最显而易见的KMP题目，下面给出此题的代码：

 

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

#define N 10008
#define M 1000008

char s[M];
char t[N];
int nxt[N];

void GetNext(char *t,int lt)
{
	int i,j;
	nxt[0]=-1;
	for(i=0,j=-1;i<lt;i++){
		while(j!=-1&&t[i]!=t[j])
			j=nxt[j];
		nxt[i+1]=++j;
	}
}
int KMP(char *s,int ls,char *t,int lt)
{
	int i,j,r=0;
	for(i=j=0;i<ls;i++){
		while(j!=-1&&s[i]!=t[j])
			j=nxt[j];
		if(++j==lt) r++;
	}
	return r;
}

int main()
{
	int T;
	int lt;

	scanf("%d",&T);

	while(T--){
		scanf("%s%s",t,s);
		GetNext(t,lt=strlen(t));
		printf("%d\n",KMP(s,strlen(s),t,lt));
	}
	return 0;
}