洛谷 P3538 OKR-A可怕的诗

目录：

---

题目：

传送门

---

分析：

题目要求给出一个字符串和$q$次询问，并求出每次询问的区间内的最短循环节长度
首先有几个性质
1、循环节一定是长度的约数
2、如果$n$是一个循环节，那么$k*n$也必定是一个循环节（关键所在）
3、$n$是$[l,r]$这一段的循环节的充要条件是$[l,r-n]$和$[l+n,r]$相同（利用这个性质我们在判断是否为循环节是可以做到$O(1)$）
所以我们可以在求出这个区间的长度之后，判断它的每个约数是否是循环节（应用性质3），并且因为性质1，它的约数是循环节，原串一定也是。
所以只要不断除以质因数（相当于从大到小枚举约数），缩小L的长度，最后就是最小的长度。
而一个重要的点在于，用上面方法来枚举约数是为了避免$TLE$
在求它的质因数的时候，可以通过线性筛的过程求得，将时间法度由根号$n$ 降为$logn$。
注意这里枚举到全长。

int k=0;
void ss()
{
    for(int i=2;i<=500000;i++)
    {
        if(f[i]==0) p[++k]=i,nex[i]=i;
        for(int j=1;j<=k&&(long long)i*p[j]<=500000;j++)
          {
            f[i*p[j]]=1;
            nex[i*p[j]]=p[j];
            if(i%p[j]==0) break;
          }
    }
    return;
}

这样通过$nxt$数组不断回跳，就可以找出所有质因数（多次乘方的也可以）
判定循环节的时候可以使用$hash$。

---

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>  
#include<cstdlib>
#include<algorithm>
#include<set>
#include<map>
#include<list>
#include<ctime>
#include<iomanip>
#include<string>
#include<bitset>
#define LL long long
#define ull unsigned long long
using namespace std;
const ull bace=13131;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
bool f[500001];int p[500001],nex[500001];
ull h[500001],se[500001];
int k=0;
void ss()
{
    for(int i=2;i<=500000;i++)
    {
        if(f[i]==0) p[++k]=i,nex[i]=i;
        for(int j=1;j<=k&&(long long)i*p[j]<=500000;j++)
          {
            f[i*p[j]]=1;
            nex[i*p[j]]=p[j];
            if(i%p[j]==0) break;
          }
    }
    return;
}
char s[500001];
int a,b;
int yue[500001];
bool check(int l1,int r1,int l2,int r2)
{
    return h[r1]-h[l1-1]*se[r1-l1+1]==h[r2]-h[l2-1]*se[r2-l2+1];
}
int main()
{
    int n=read();
    scanf("%s",s+1);
    for(int i=1;i<=n;i++)
      h[i]=h[i-1]*bace+s[i]-'a'+1;
    se[0]=1;
    for(int i=1;i<=n;i++)
      se[i]=se[i-1]*bace;
    int m=read();
    ss();
    for(int i=1;i<=m;i++)
    {
        int t=0;
        a=read();b=read();
        int len=b-a+1;
        while(len!=1)
        {
            yue[++t]=nex[len];
            len=len/nex[len];
        }
        len=b-a+1;
        for(int j=1;j<=t;j++)
        {
            int w=len/yue[j];
            if(check(a,b-w,a+w,b)) len=w;
        }
        printf("%d\n",len);
    }
    return 0;
}