2795: [Poi2012]A Horrible Poem

最新推荐文章于 2025-06-27 16:48:46 发布

原创最新推荐文章于 2025-06-27 16:48:46 发布 · 319 阅读

0 ·

CC 4.0 BY-SA版权

2751-3000 专栏收录该内容

32 篇文章

订阅专栏

本文介绍了一种解决字符串中最短循环节问题的有效算法。针对给定字符串及多个查询区间，该算法能快速找到每个查询区间的最短循环节长度。通过质数筛法预处理与哈希技巧实现高效判断。

题目链接

题目大意：给出一个长度为N的字符串，有Q次询问；
每次询问给出一个区间，求区间最短循环节长度

题解：膜一发GXZ大爷

我的收获：循环节……e

#include<bits/stdc++.h>
using namespace std;

typedef unsigned long long ull;

const int N=500009;
const ull B=1000173169;

int n,q;
int p[N],lp[N];
char S[N];
bool vis[N];
ull Base[N],H[N];

bool judge(int a,int b,int L){return H[a+L-1]-H[a-1]*Base[L]==H[b]-H[b-L]*Base[L];}

void euler()
{
    for(int i=2;i<=n;i++){
        if(!vis[i]) lp[i]=i,p[++p[0]]=i;
        for(int j=1;j<=p[0]&&i*p[j]<=n;j++)
        {
            vis[i*p[j]]=1;lp[i*p[j]]=p[j];
            if(i%p[j]==0) break;
        }
    }
    Base[0]=1;for(int i=1;i<=n;i++) Base[i]=Base[i-1]*B,H[i]=H[i-1]*B+S[i];
}

void work()
{
    for(scanf("%d",&q);q;q--){
        int l,r;scanf("%d%d",&l,&r);
        int now=r-l+1,v=now;
        while(v!=1){
            if(judge(l,r,r-l+1-now/lp[v])) now/=lp[v];
            v/=lp[v];
        }
        printf("%d\n",now);
    }
}

void init()
{
    scanf("%d%s",&n,S+1);
    euler();
}

int main()
{
    init();
    work();
    return 0;
}