poj2406 Power Strings kmp

      一个字符串可以由他的一个子串a重复n次得到，现在给一个长度不大于10^6的字符串，求最长n。

      跟白书上的period （训练之南P213，例13）基本一样甚至还简单一点。核心就是利用KMP中的失配函数求最短循环节。可以得到对于字符串的第i位，若有i % (i-f[i])==0，则f[i]到i之间的部分及为这个串的最短循环节，具体画个图，或者参考白书上的图，很容易就能看出来。

      另外这题后缀数组应该也能写，不过10^6的数据必须要写DC3了，倍增的SA目测要T...

      上一个KMP的代码

#include <iostream>
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <algorithm>
#include <stack>
#include <queue>
#include <map>
#include <string>
#include <cstring>
#include <string>
using namespace std;
typedef long long ll;
const int maxn=1000000+100;
char str[maxn],s[maxn],s1[maxn];
int f[maxn];
void getFail(char* P,int* f)
{
    int m=strlen(P);
    f[0]=0;
    f[1]=0;
    for (int i=1; i<m; i++)
    {
        int j=f[i];
        while(j && P[i]!=P[j]) j=f[j];
        f[i+1]=P[i]==P[j]?j+1:0;
    }
}
int find(char* T,char* P,int* f)
{
    int n=strlen(T);
    int m=strlen(P);
    getFail(P,f);
    int j=0;
    for (int i=0; i<n; i++)
    {
        while(j && P[j]!=T[i]) j=f[j];
        if (P[j]==T[i]) j++;
        if (j==m) return i-m+1;
    }
}


int main()
{
//    freopen("in.txt","r",stdin);
    int tt;
    while(~scanf("%s",s))
    {
        if (s[0]=='.') break;
        getFail(s,f);
        int len=strlen(s);
        if (len % (len-f[len])==0)
        {
            printf("%d\n",len/(len-f[len]));
        }
        else printf("1\n");
    }
    return 0;
}