题目大意：给出一个字符串s，则存在子串a,a重复k次后得到s。求k的{zd0}值。

分析：这个题目数据规模很大，用后缀数组做的话，倍增算法超时..- -dc3算法效率也不敢恭维...比较好的做法是KMP。不过这不是我论述的重点，毕竟我是在学习后缀数组ing。

用后缀数组做的话，先枚举子串的长度，判断能否被n整除，再判断该串能否通过复制变成原串。判断方法是求suffix(1)和suffix(i+1)的最长公共前缀k，若k=n-i，则说明可行。

为什么呢？

首先，suffix(i+1)的长度是n-i。那么，lcp(1,i+1)=n-i的时候，说明整个串可以由s[1..i]和s[i+1..n]组成，并且s[i+1..n]是由若干个s[1..i]组成的。其实，lcp(1,i+1)的值，表示了从字符串的开端，s[1..i]所延伸（重复出现）的最长长度。所以当这个值等于剩余字符串的长度(n-i)时，表明s[1..i]可以延伸至整个串。如果不理解，可以举个具体的例子细细体会。这是后缀的性质所决定的。

因为只需要求lcp(1,i)，所以可以用一个O(n)的扫描进行预处理。时间消耗主要在后缀数组的构建上，1000000的数据规模用倍增算法只能超时（NlogN），dc3也很勉强..不过能过就好。

算是一次ac吧，毕竟TLE是我阻止不了的...

notes：

1、写程序之前，应该对算法流程有充分的认识。比如说那个地方应该用i，那个地方应该用rank[i]，这些都应该想清楚再敲代码，不然调试的时间会比较漫长。（我就是- -）

2、一开始忘了lcp(i,i)=n-i+1....

codes（dc3算法贴了模板...暂时还不打算学习- -）:

#include<stdio.h>
#include<cstring>
using namespace std;
const int maxn=3000010;
int ws[maxn],wa[maxn],wb[maxn],wv[maxn],sa[maxn],rank[maxn],height[maxn],a[maxn],f[maxn];
char s[maxn];

//dc3
#define F(x) ((x)/3+((x)%3==1?0:tb))
#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)
int c0(int *r,int a,int b)
{return r[a]==r[b]&&r[a+1]==r[b+1]&&r[a+2]==r[b+2];}
int c12(int k,int *r,int a,int b)
{if(k==2) return r[a]<r[b]||r[a]==r[b]&&c12(1,r,a+1,b+1);
else return r[a]<r[b]||r[a]==r[b]&&wv[a+1]<wv[b+1];}
void sort(int *r,int *a,int *b,int n,int m)
{
    int i;
    for(i=0;i<n;i++) wv[i]=r[a[i]];
    for(i=0;i<m;i++) ws[i]=0;
    for(i=0;i<n;i++) ws[wv[i]]++;
    for(i=1;i<m;i++) ws[i]+=ws[i-1];
    for(i=n-1;i>=0;i--) b[--ws[wv[i]]]=a[i];
    return;
}
void dc3(int *r,int *sa,int n,int m)
{
     int i,j,*rn=r+n,*san=sa+n,ta=0,tb=(n+1)/3,tbc=0,p;
     r[n]=r[n+1]=0;
     for(i=0;i<n;i++) if(i%3!=0) wa[tbc++]=i;
     sort(r+2,wa,wb,tbc,m);
     sort(r+1,wb,wa,tbc,m);
     sort(r,wa,wb,tbc,m);
     for(p=1,rn[F(wb[0])]=0,i=1;i<tbc;i++)
     rn[F(wb[i])]=c0(r,wb[i-1],wb[i])?p-1:p++;
     if(p<tbc) dc3(rn,san,tbc,p);
     else for(i=0;i<tbc;i++) san[rn[i]]=i;
     for(i=0;i<tbc;i++) if(san[i]<tb) wb[ta++]=san[i]*3;
     if(n%3==1) wb[ta++]=n-1;
     sort(r,wb,wa,ta,m);
     for(i=0;i<tbc;i++) wv[wb[i]=G(san[i])]=i;
     for(i=0,j=0,p=0;i<ta && j<tbc;p++)
     sa[p]=c12(wb[j]%3,r,wa[i],wb[j])?wa[i++]:wb[j++];
     for(;i<ta;p++) sa[p]=wa[i++];
     for(;j<tbc;p++) sa[p]=wb[j++];
     return;
}
void cal(int *r,int *sa,int n){
     int i,j,k=0;
     for (i=1;i<=n;i++) rank[sa[i]]=i;
     for (i=0;i<n;height[rank[i++]]=k)
         for (k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
     return;
}
void rmq(int *height,int n){
     int i,j=rank[0],min=n-sa[j];
     for (i=j;i>=1;i--){
         f[i]=min;
         min=min<height[i]?min:height[i];
     }
     for (i=j+1,min=1000000;i<=n;i++){
         min=min<height[i]?min:height[i];
         f[i]=min;
     }
     return;
}
int work(int n){
    int i;
    for (i=0;i<n;i++)
        if (n%(i+1)==0)
           if (f[rank[i+1]]==n-i-1) return n/(i+1);
    return 0;
}
int main(){
    while (scanf("%s",s)!=EOF){
          if (strcmp(s,".")==0) return 0;
          int n=strlen(s),i;
          for (i=0;i<n;i++) a[i]=static_cast<int>(s[i]);
          a[n]=0;
          dc3(a,sa,n+1,123);
          cal(a,sa,n);
          rmq(height,n);
          printf("%d\n",work(n));
          }
    return 0;
}