【BZOJ】2795: [Poi2012]A Horrible Poem

题意

一个长度为\(n(n \le 500000)\)的字符串\(s\)，给\(q(q \le 2000000)\)个询问，每个询问给一个区间\([l, r]\)，求这个区间内最短的循环节。

分析

分析以下可以知道：

假设循环节长度为\(len\)，则\(s[l, r-len]=s[l+len, r]\)。

\(len|(r-l+1)\)

如果\(len\)是循环节，则\(len * p\)也是循环节\((len * p|(r-l+1))\)

题解

首先判两个串是否相等用hash即可。

根据\(1, 2\)我们很容易得到\(O(qn^{0.5})\)的做法。

可是由于性质\(3\)的存在，我们先分解\((r-l+1)\)的质因数，然后依次考虑删去质因数即可。

复杂度\(O(qlogn)\)

#include <bits/stdc++.h>
using namespace std;
inline int getint() {
int x=0, c=getchar();
for(; c<48||c>57; c=getchar());
for(; c>47&&c<58; x=x*10+c-48, c=getchar());
return x;
}
const unsigned int N=500005, mo=1e9+7;
int p
, cnt, m
, bg
;
unsigned int a
, po
;
void init(int n) {
for(int i=2; i<=n; ++i) {
if(!bg[i]) {
bg[i]=i;
p[cnt++]=i;
}
for(int j=0, t; j<cnt; ++j) {
t=p[j]*i;
if(t>n) {
break;
}
bg[t]=p[j];
if(i%p[j]==0) {
break;
}
}
}
}
inline bool check(int a, int b, int c, int d) {
return m[b]-m[a-1]*po[b-a+1]==m[d]-m[c-1]*po[d-c+1];
}
int main() {
int n=getint();
init(n);
po[0]=1;
for(int i=1; i<=n; ++i) {
m[i]=m[i-1]*mo+getchar();
po[i]=po[i-1]*mo;
}
int q=getint();
for(; q--; ) {
int l=getint(), r=getint(), len=r-l+1;
for(int x=len; x>1; ) {
int y=bg[x];
for(; len%y==0 && check(l, r-len/y, l+len/y, r); len/=y);
for(; x%y==0; x/=y);
}
printf("%d\n", len);
}
return 0;
}