本文搬运自本人初中博客园博客,若图片加载不出来,可到原文查看:https://www.cnblogs.com/zhangtingxi/p/15810272.html

题目

原题来自:POI 2012

给出一个由小写英文字母组成的字符串 S,再给出 q 个询问,要求回答 S 某个子串的最短循环节。

如果字符串 B 是字符串 A 的循环节,那么 A 可以由 B 重复若干次得到。

思路

首先,我们如果有三点:

  1. 一个字符串的循环节必然是字符串长度的约数
  2. 循环节的倍数如果长度还是字符串长度的约数,那么他也是循环节
  3. 如果一个长度 ii 是字符串循环节长度,那么 [l,ri][l, r-i][l+i,r][l+i, r] 必然相等

有了这三点,我们只需要打一个素数筛+哈希即可。

总结

这题是一道很综合的题,要求掌握两个知识点,是对循环节的综合考验:

  1. 如下图中,对于字符串 SS,若红色部分长度为 SS 长度的因数,且黄色部分等于橙色部分,则红色部分为 SS 的循环节。
    image
  2. 循环节的倍数如果长度还是字符串长度的约数,那么他也是循环节,于是可以从大到小试除判断

Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
// Problem: P3538 [POI2012]OKR-A Horrible Poem
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P3538
// Memory Limit: 125 MB
// Time Limit: 1000 ms
//
// Powered by CP Editor (https://cpeditor.org)

#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int read(){int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;
ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<1)+
(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define N 500010
//#define M
#define mo 998244353
// #define mo 1000000009
struct node
{
int x, y, n;
}d[N*4];
int n, m, i, j, k;
int a[N], f[N], l, r;
int len, h[N], b[N], g;
char s[N];

inline int MOD(int x)
{
return (x%mo+mo)%mo;
}

inline int qiu(int l, int r)
{
return MOD(a[r]-a[l-1]*f[r-l+1]);
}

void cun(int x, int y)
{
d[++k].x=x; d[k].y=y;
d[k].n=h[x]; h[x]=k;
}

int check(int x)
{
return qiu(l, r-x)==qiu(l+x, r);
}

signed main()
{
// freopen("tiaoshi.in","r",stdin);
// freopen("tiaoshi.out","w",stdout);
n=read(); scanf("%s", s+1);
for(i=f[0]=1; i<=n; ++i)
{
a[i]=MOD(a[i-1]*29+s[i]-'a'+1);
f[i]=MOD(f[i-1]*29);
}
for(i=2; i<=n; ++i)
if(!b[i])
for(j=i; j<=n; j+=i) cun(j, i), b[j]=1;
m=read();
while(m--)
{
l=read(); r=read();
len=r-l+1;
// printf("> %lld\n", len);
for(g=h[len]; g; g=d[g].n)
{
i=d[g].y;
while(len%i==0 && check(len/i)) len/=i;
}
printf("%lld\n", len);
}
return 0;
}