POJ2406 Power Strings POJ1961 Period
http://poj.org/problem?id=2406
http://poj.org/problem?id=1961
几乎是一个题
1961是对于第 \(i\) 位,求 \([1,i]\) 能不能由一段字符循环一次以上组成,输出 \(i\) 和这样的长度最小的循环节循环次数
2406是只对第 \(n\) 位输出最小长度循环节循环次数
kmp,考虑 kmp 的 \(fail\) 数组,\(fail_i\) 表示的是如果 \(i\) 失配,那么要从哪一位开始继续匹配
更本质的,他是 \([1,i]\) 中,最长的前缀等于后缀的长度
例如 \(fail_i=m\),那么字串 \([1,m]=[i-m+1,i]\),如果这两个区间的并集等于整个区间,那么 \(i-fail_i\) 则 可能 是这个区间的循环节长度
因为并集等于整个区间,所以 \(\forall 1\le j \le fail_i,s(j)=s(j+i-fail_i)\)
对于这个式子也可以理解为前缀向后移动了 \(i-fail_i\) 位得到了后缀,所以相当于每一个区间内字符都与向后移动 \(i-fail_i\) 位的字符相同
因此,只要 \(i \bmod i-fail_i=0\),就说明了 \(i-fail_i\) 是 \([1,i]\) 最小循环节长度(这也是为什么刚才说可能)
POJ2406
char s[6000005];
int fail[6000005];
int main(){
scanf("%s",s+1);
while(s[1]!='.'){
int n=std::strlen(s+1);
fail[1]=0;
for(reg int j,i=2;i<=n;i++){
j=fail[i-1];
while(j&&s[j+1]!=s[i]) j=fail[j];
if(j) fail[i]=j+1;
else fail[i]=(s[1]==s[i]);
}
printf("%d\n",(n%(n-fail[n]))?1:(n/(n-fail[n])));
scanf("%s",s+1);
}
return 0;
}
POJ1961
char s[6000005];
int fail[6000005];
int n;
int main(){
scanf("%d",&n);
int T=0;
while(n){
scanf("%s",s+1);
fail[1]=0;
for(reg int j,i=2;i<=n;i++){
j=fail[i-1];
while(j&&s[j+1]!=s[i]) j=fail[j];
if(j) fail[i]=j+1;
else fail[i]=(s[1]==s[i]);
}
printf("Test case #%d\n",++T);
for(reg int i=2;i<=n;i++)
if(!(i%(i-fail[i]))&&i/(i-fail[i])>1) printf("%d %d\n",i,i/(i-fail[i]));
EN;
scanf("%d",&n);
}
return 0;
}