利用rmq加速查询后缀数组中两个不同排名的后缀的lcp

本文介绍了一种解决SPOJ-REPEATS问题的有效算法,通过枚举循环节长度并利用RMQ预处理height数组,实现快速查询连续重复子串的最大重复次数。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

原题是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)操作。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值