POJ 2406 Power Strings(后缀数组做法超时)
http://poj.org/problem?id=2406
题意:
一个串确定是由一个或多个循环节构成的,现在要你求出这个串由最多几个循环节构成的.(即用串长除以最短循环节的长度).
分析:
本题可以用KMP算法求解,简单方便代码段只有30行,O(n)复杂度:
http://blog.csdn.net/u013480600/article/details/22944853
下面用后缀数组来解.首先对原始串建立后缀数组求出sa和height.
然后我们从小到大一一枚举最小循环节的长度k.
首先k肯定要整除字符串长度len的.如果len%k=0.那么就判断后缀0和后缀0+k的LCP是不是等于len-k.如果是那么k就是最短循环节长度,否则不是.
求后缀0和后缀0+k的LCP可以用RMQ求即可.
RMQ维护height值区间为[1,len].因为height[0]表示的是尾0字符串,没有实际意义可以不管.
最后由于我们每次询问都是问rank[0],到rank[j]之间的height最小值,所以RMQ算法可以用O(n)预处理出dmin[j]表示height数组中第j个数到第rank[0]之间的最小height值.
超时代码:我用的da倍增法算的后缀数组(用的O(n)的RMQ),超时了,我没有试是否能用d3c构造的后缀数组,不过总的来说算循环节还是KMP算法来到快简单(一共个才30行代码).网上找了两个用da倍增法做的这题都是超时代码…
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1000000+100;
struct SuffixArray
{
char s[maxn];
int sa[maxn],rank[maxn],height[maxn];
int t1[maxn],t2[maxn],c[maxn],n;
int dmin[maxn];
void build_sa(int m)
{
int i,*x=t1,*y=t2;
for(i=0;i<m;i++) c[i]=0;
for(i=0;i<n;i++) c[x[i]=s[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(int k=1;k<=n;k<<=1)
{
int p=0;
for(i=n-k;i<n;i++) y[p++]=i;
for(i=0;i<n;i++)if(sa[i]>=k) y[p++]=sa[i]-k;
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]]= y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]? p-1:p++;
if(p>=n) break;
m=p;
}
}
void build_height()
{
int i,j,k=0;
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(s[i+k]==s[j+k])k++;
height[rank[i]]=k;
}
}
void initMin()
{
int min_v=1e8;
for(int i=rank[0];i>1;i--)//求1到rank[0] 之间的最小值height
{
min_v=min(min_v,height[i]);
dmin[i-1]=min_v;
}
min_v=1e8;
for(int i=rank[0]+1;i<n;i++)//求rank[0]到n-1 之间的最小值height
{
min_v=min(min_v,height[i]);
dmin[i]=min_v;
}
}
int LCP(int j)//求后缀i和j的LCP最长公共前缀
{
int R=rank[j];
return dmin[R];
}
}sa;
int main()
{
while(scanf("%s",sa.s)==1)
{
int len = strlen(sa.s);
if(len==1 && strcmp(sa.s,".")==0) break;
sa.n=len+1;
sa.s[len]=0;
sa.build_sa(128);
sa.build_height();
sa.initMin();
bool ok=false;
for(int i=1;i<=len/2;i++)if(len%i==0)
{
int lcp=sa.LCP(i);//其实应该是LCP(0,i)
if(lcp==len-i)
{
printf("%d\n",len/i);
ok=true;
break;
}
}
if(ok==false) printf("1\n");
}
return 0;
}