bzoj 2795 [Poi2012]A Horrible Poem hash

最新推荐文章于 2019-03-13 16:12:16 发布

原创最新推荐文章于 2019-03-13 16:12:16 发布 · 292 阅读

0 ·

CC 4.0 BY-SA版权

字符串hash 专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一种利用质因数优化的字符串匹配算法，通过筛选素数并判断每个质因子是否可行，来减少O(qn√)的时间复杂度。算法适用于寻找字符串中重复子串的问题，并给出了完整的实现代码。

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

如果直接枚举长度的约数，hash判断s[l~r-len]和s[l+len~r]是否相等是 $O(q\sqrt{n})$ 的。
由于最后答案的段数对于每个质因子是独立的，所以先把素数筛出来，然后判断每个质因子是否可行，如果可行就将答案的段数乘上这个质因子，就可以把复杂度除log

#include <bits/stdc++.h>
using namespace std;
#define N 510000
#define seed 133
#define ull unsigned long long
int n,cnt,q,l,r;
char s[N];
ull has[N],p[N];
int prime[N],ip[N];
void init()
{
    for(int i=2;i<=n;i++)
    {
        if(!ip[i])prime[++cnt]=i;
        for(int j=1;j<=cnt&&i*prime[j]<=n;j++)
        {
            ip[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}
ull get(int l,int r){return has[r]-has[l-1]*p[r-l+1];}
int check(int x)
{
    int len=(r-l+1)/x;
    return get(l,r-len)==get(l+len,r);
}
int main()
{
    //freopen("tt.in","r",stdin);
    scanf("%d%s",&n,s+1);
    init();p[0]=1;
    for(int i=1;i<=n;i++)p[i]=p[i-1]*seed;
    for(int i=1;i<=n;i++)
        has[i]=has[i-1]*seed+s[i];
    scanf("%d",&q);
    for(;q--;)
    {
        scanf("%d%d",&l,&r);
        int t=r-l+1,ans=1;
        for(int i=1;prime[i]*prime[i]<=t&&i<=cnt;i++)
        {
            while(t%prime[i]==0)
            {
                if(check(ans*prime[i]))ans*=prime[i];
                t/=prime[i];
            }
        }
        if(t!=1&&check(ans*t))ans*=t;
        printf("%d\n",(r-l+1)/ans);
    }
    return 0;
}