http://bailian.openjudge.cn/practice/1961/
总时间限制: 3000ms 内存限制: 65536kB
描述
一个字符串的前缀是从第一个字符开始的连续若干个字符,例如”abaab”共有5个前缀,分别是a, ab, aba, abaa, abaab。
我们希望知道一个N位字符串S的前缀是否具有循环节。换言之,对于每一个从头开始的长度为 i (i 大于1)的前缀,是否由重复出现的子串A组成,即 AAA…A (A重复出现K次,K 大于 1)。如果存在,请找出最短的循环节对应的K值(也就是这个前缀串的所有可能重复节中,最大的K值)。
输入
输入包括多组测试数据。每组测试数据包括两行。
第一行包括字符串S的长度N(2 <= N <= 1 000 000)。
第二行包括字符串S。
输入数据以只包括一个0的行作为结尾。
输出
对于每组测试数据,第一行输出 “Test case #“ 和测试数据的编号。
接下来的每一行,输出前缀长度i和重复测数K,中间用一个空格隔开。前缀长度需要升序排列。
在每组测试数据的最后输出一个空行。
样例输入
3
aaa
12
aabaabaabaab
0
样例输出
Test case #1
2 2
3 3
Test case #2
2 2
6 2
9 3
12 4
思路:
运用KMP算法中的next[]
(1) 判断是否循环: k*2>j
(2) 循环节的长度 = j-k, 看j是否是(j-k)的整数倍
#include <iostream>
#include <memory.h>
using namespace std;
int len;
char s[10001010];
int KVal[10001010];//KVal用来储存next数组
int j, k;
int n=1;
int main(){
while (cin>>len){
if (len==0) break;
memset(KVal,0,sizeof(KVal));
cin>>s;
KVal[0]=-1; KVal[1]=0;
j=0; k=-1;
cout<<"Test case #"<<n<<endl;
while(j<len){
while (k>=0 && s[p2]!= s[p1]) {k = KVal[k];}
j++; k++;
KVal[j]=k;
// j-k 为循环节的长度; j%(j-k)==0表明前缀长度是循环节长度的整数倍
// k是前缀j中,最长相等的首尾缀的长度,k必须大于j的一半,才说明是循环的
if (j%(j-k)==0 && k*2>=j) cout<<j<<" "<<j/(j-KVal[j])<<endl;
}
cout<<endl; n++;
}
return 0;
}