spoj 687 repeats

转载声明:http://blog.csdn.net/kongming_acm/article/details/6321363
论文转载:
http://wenku.baidu.com/link?url=xBx1136joEtqVhvjmGnFgsJ9c5LLr9aw2AsupDH8u2UpucAZB2J9WAFc48o6ZU-xphWysH8HAJZLwvImlfg36cJAQJfP3VejDonakPI3_Ey】

题意:简单明了,就是找一个子串重复最多次而形成的.

个人感想:
刚看到这道题的时候,我也是一脸懵逼啊..我还想用自己的能力去想想该怎么做,我是毫无想法,看了论文,我是一头雾水了,还是自己的能力还没达到这样的程度啊.我是看这转载 然后硬刚究竟什么回事.这道题也是出乎我的感觉..我想通过我的硬钢想法来描绘一下这代码.

首先我们知道,循环的次数1是必然的.我们现在要找2次以上的.

代码的意思是先枚举循环节
红色代表长度为1的循环节.
橙色代表长度为2的循环节.
绿色代表长度为3的循环节.
蓝色代表长度为4的循环节.(4的循环节肯定是不存在的,超出了最长长度了)
假设长度为L.
假设当前的左点为i,右点为j

我每条横线上都有一个小点,代表枚举的时候,我要知道这两个位置上的最长公共前缀lcp.
这有什么用了? 通过这个我们就可以之后两个串能过延伸最长长度k,通过 k/L.我们就知道这两个串的 循环节个数 可是只往右边枚举不行啊. 我们得尝试的往左一步求循环节个数 因为我们知道 k%L是不一定能整除的,但是当k>L的时候,其实（j-i）,这块就已经是循环节了.
可以试这画画究竟是不是.

但是多出来的部分k%L！=0 多出来的又是怎么一回事呢?我们又得怎么处理呢?

我们看看这个图,(i,j)的公共前缀=5.
但是我们的循环节长度为3!这超过了2.我先跟我一步步来
（）代表着匹配的长度

（1）
（2）
（3）
这时候我们知道循环节已经=3,那么这时候的循环节是BAC,因为刚好k满足3.

这是我们就可以往左匹配一下,如果（i-（l-k%l ) ,j-（l-k%l ))>=3,那么我们的循环节必然就可以再+1

注意！还没完,因为(i,j)的公共前缀是5,接着继续!
（4）
是不是发现我的j点的颜色加深了,没错这时候的循环节就可以发生改变了为ACB.

这我们就可以比较（i-2，j-2）的公共前缀是不是>=2，我们就可以多+1循环节啦！
（5）
同样道理就是,这时就可以直接比较（i-1，j-1）是不是>=1,我们又可以+1循环节

这就是为什么结果往左的时候要用到（i-（l-k%l ) ,j-（l-k%l ))，这是一个微妙的贪心策略.

这时看到代码的时候,我们会发现一个问题,就是为什么枚举长度的时候, 只需要枚举**（i,L）(L,2*L),
(2*L，3*L）..这样的情况,例如**L=2,它绝对不需要枚举到(1,3),
因为这样枚举的时候,如果循环节满足长度为2,那么不管L=2,的区间怎么枚举,他都会触碰到

不伦检测（1，3）,还是（2，4）那么他整个循环节我都可以检测出来,往,通过往左修正一样能是得到同样结果..这是我感觉出来的,可能这是循环节的特性吧..我也不太懂.

代码:

/* Author:GavinjouElephant * Title: * Number: * main meanning： * * * *///#define OUT#include <iostream>using namespace std;#include <cstdio>#include <cmath>#include <cstring>#include <algorithm>#include <sstream>#include <cctype>#include <vector>#include <set>#include <cstdlib>#include <map>#include <queue>//#include<initializer_list>//#include <windows.h>//#include <fstream>//#include <conio.h>#define MaxN 0x7fffffff#define MinN -0x7fffffff#define Clear(x) memset(x,0,sizeof(x))const int INF=0x3f3f3f3f;int T;const int maxn=500000;int t1[maxn];int t2[maxn];int c[maxn];int N;char in[2];char s[maxn];int t[maxn];bool cmp(int * r,int a,int b,int l){    return r[a]==r[b] &&r[a+l]==r[b+l];}void da(char str[],int sa[],int Rank[],int height[],int n,int m){    n++;    int i,j,p;    int* x=t1;    int* y=t2;    for(i=0;i<m;i++) c[i]=0;    for(i=0;i<n;i++) c[x[i]=str[i]]++;    for(i=1;i<m;i++) c[i]+=c[i-1];    for(i=n-1;i>=0;i--)sa[--c[x[i]]]=i;    for(j=1;j<=n;j<<=1)    {        p=0;        for(i=n-j;i<n;i++)y[p++]=i;        for(i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;        for(i=0;i<m;i++)c[i]=0;        for(i=0;i<n;i++)c[x[y[i]]]++;        for(i=1;i<m;i++) c[i]+=c[i-1];        for(i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];        swap(x,y);        p=1;        x[sa[0]]=0;        for(i=1;i<n;i++)        x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;        if(p>=n)break;        m=p;    }    int k=0;    n--;    for(i=0;i<=n;i++)Rank[sa[i]]=i;    for(i=0;i<n;i++)    {        if(k)k--;        j=sa[Rank[i]-1];        while(str[i+k]==str[j+k])k++;        height[Rank[i]]=k;    }}int sa[maxn];int rk[maxn];int height[maxn];int dpmin[maxn][25];void cdpmin(int n){    int i,j;    for(i=1;i<=n;i++) dpmin[i][0]=height[i];    for(j=1;j<=log(double(n+1))/log(2.0);j++)    {        for(i=1;i+(1<<j)-1<=n;i++)        {            dpmin[i][j]=min(dpmin[i][j-1],dpmin[i+(1<<(j-1))][j-1]);        }    }}int get_min(int x,int y){    int k=(int)(log(double(y-x+1))/log(2.0));    return min(dpmin[x][k],dpmin[y-(1<<k)+1][k]);}int get_min_lcp(int x,int y){    x=rk[x],y=rk[y];    if(x>y) swap(x,y);    x++;    int k=(int)(log(double(y-x+1))/log(2.0));    return min(dpmin[x][k],dpmin[y-(1<<k)+1][k]);}int main(){#ifdef OUT    freopen("coco.txt","r",stdin);    freopen("lala.txt","w",stdout);#endif    scanf("%d",&T);    while(T--)    {        scanf("%d",&N);        for(int i=0;i<N;i++)        {            scanf("%s",in);            s[i]=in[0];        }        s[N]='\0';        da(s,sa,rk,height,N,256);        cdpmin(N);        int ans=0;        for(int i=1;i<N;i++)//长度为i的子串        {            for(int j=0;j+i<N;j+=i)            {                int k=get_min_lcp(j,j+i);                int now=k/i;                int tj=j-(i-k%i);                if(tj>=0)                    if(get_min_lcp(tj,tj+i)>=(i-k%i)) now++;                if(now+1>ans) ans=now+1;            }        }        printf("%d\n",ans);    }    return 0;}