POJ 2406 kmp求解循环节

照着论文题在刷，结果到这的时候看不懂用后缀数组怎么搞，那就用kmp咯o(∩_∩)o

题目大意：

已知一个字符串S有其中某一个字串x循环R次得到，求问R的最大值为多少

结论：

int ans=1; if(len%(len-next[len])==0) ans=len/(len-next[len])

证明：

设：p为循环节，next同kmp中next定义，s为p所构成的串，k=len/p(k>1)

1：len-next[i]|len是具有循环节的充要条件：
由next[len]的定义可知：s[1]…….s[next[i]]==s[n]……s[n-next[i]]
next[len]=前k-1个循环节==后k-1个循环节=len-p;
p=len-next[len]
由 p|len =>len-next[len]|len
故证得： len-next[i]|len 是具有循环节的必要条件

2：用B=A^n 表示B=AAA….AAA(A有n个)
由next[len]得:
①： s[1]……s[next[len]]== s[len-next[len]]……s[len]
令q=len-next[len] 由①得
②： s[q]……s[q+next[len]]=s[q-len+next[len]+1]……s[q-len+2*next[len]]
由②可得：
③：s[len-2*q+1]…s[len]=s[len-q+1]……s[len]^2
由①②③推广可以得到：
④：s[len-n*q+1]…….s[len]=s[len-q+1]……s[len]^n
若满足：len=nq 可得s[len-q+1]……s[len]为一个循环节 长度为q
故q=len-next[len]=p
不难推出len>p=q 充分性得证

综上所述，len-next[len]|len 为p存在的充要条件
PS:上述所有p都的是最小长度的循环节且（p>len）

Code：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e6+10;
char a[N];
int p[N];
int main()
{
    while(scanf("%s",a+1)) {
        if(a[1]=='.') break;
        int n=strlen(a+1);
        int j=0,ans=1;
        p[1]=0;
        for(int i=2;i<=n;i++) {
            while(j>0&&a[j+1]!=a[i])j=p[j];
            if(a[j+1]==a[i])j++;
            p[i]=j;
        }
//      for(int i=1;i<=n;i++) printf("p[%d]=%d\n",i,p[i]);
        if(n%(n-p[n])==0) ans=n/(n-p[n]);
        cout<<ans<<endl;
    }
    return 0;
}