Description
给出一个由小写英文字母组成的字符串S,再给出q个询问,要求回答S某个子串的最短循环节。
如果字符串B是字符串A的循环节,那么A可以由B重复若干次得到。
一个正整数n (n<=500,000),表示S的长度。
一个正整数q (q<=2,000,000),表示询问个数。
Solution
暴力屮过系列
不难想到枚举循环长度,判断是否循环可以用hash取巧
容易发现循环的长度一定是b-a+1的因数可以根号枚举,但这显然不够
又可以发现记cnt[i]表示字符i在[a,b]出现的次数,那么分成的段数一定不大于他们全部的gcd
以上两个小优化加起来就能过了,只跑了一半不到的时间
Code
#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
typedef long long LL;
const int MOD1=1000000007;
const int MOD2=1000000009;
const int N=500005;
char str[N];
int s[N][27];
LL h1[N],h2[N];
LL p1[N],p2[N];
int read() {
int x=0,v=1; char ch=getchar();
for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
for (;ch>='0'&&ch<='9';x=x*10+ch-'0',ch=getchar());
return x*v;
}
LL get_hash1(int x,int y) {
return (h1[y]%MOD1-h1[x-1]*p1[y-x+1]%MOD1+MOD1)%MOD1;
}
LL get_hash2(int x,int y) {
return (h2[y]%MOD2-h2[x-1]*p2[y-x+1]%MOD2+MOD2)%MOD2;
}
bool check(int l,int r,int len) {
LL ret1=get_hash1(l,l+len-1);
LL ret2=get_hash2(l,l+len-1);
// for (int pos=l+len;pos<=r;pos+=len) {
LL tmp1=get_hash1(r,r+len-1);
if (tmp1!=ret1) return false;
LL tmp2=get_hash2(r,r+len-1);
if (tmp2!=ret2) return false;
// }
return true;
}
int gcd(int x,int y) {
return !y?x:gcd(y,x%y);
}
int main(void) {
int n,T; scanf("%d%s%d",&n,str+1,&T);
rep(i,1,n) rep(j,0,25) s[i][j]=s[i-1][j]+(str[i]-'a'==j);
rep(i,1,n) h1[i]=(h1[i-1]*27%MOD1+str[i]-'a')%MOD1;
rep(i,1,n) h2[i]=(h2[i-1]*27%MOD2+str[i]-'a')%MOD2;
p1[0]=1; rep(i,1,n) p1[i]=p1[i-1]*27%MOD1;
p2[0]=1; rep(i,1,n) p2[i]=p2[i-1]*27%MOD2;
while (T--) {
int l=read(),r=read(); int lim=r-l+1;
rep(i,0,25) if (s[r][i]-s[l-1][i]) lim=gcd(lim,s[r][i]-s[l-1][i]);
int ans=r-l+1;
for (int lam=1;lam*lam<=lim;lam++) {
if (lim%lam) continue;
int len1=(r-l+1)/lam;
int len2=(r-l+1)/(lim/lam);
if (check(l,l+len1,r-l+1-len1)) ans=std:: min(len1,ans);
if (check(l,l+len2,r-l+1-len2)) ans=std:: min(len2,ans);
}
printf("%d\n", ans);
}
return 0;
}