原题是SPOJ - REPEATS
这道题需要求给出的串中重复次数最多的连续重复子串的重复次数。我们可以枚举循环节长度l,然后查询rk[i],rk[i+l],rk[i+2*l]…的lcp,这个可以通过height数组得到,但是查询过程是o(n)的。但我们知道查询height数组时是一个连续取min的过程,所以这个过程可以通过rmq加速。
int dp[maxn][20], mm[maxn];
void initRMQ(int n, int b[])
{
mm[0] = -1;
for(int i = 1; i<=n; i++)
dp[i][0] = b[i], mm[i] = ((i&(i-1))==0)?mm[i-1]+1:mm[i-1];
for(int j = 1; j<=mm[n]; j++)
for(int i = 1; i+(1<<j)-1<=n; i++)
dp[i][j] = min(dp[i][j-1], dp[i+(1<<(j-1))][j-1]);
}
int RMQ(int x, int y)
{
if(x>y) swap(x, y);
x++;
int k = mm[y-x+1];
return min(dp[x][k], dp[y-(1<<k)+1][k]);
}
这里的mm[i]=log2(i)
需要注意的是我们在查询时是查询rk[i],rk[i+l]这样的区间,但前者的排名未必比后者小,所以在查询时得check一下。这样每次查询就变成了o(1)操作。