bzoj2790 [Poi2012]A Horrible Poem（hash+枚举）

6.20upd: 枚举质因数才是科学的复杂度。大概$qlogn$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
int n,m,prime[N<<2],tot=0,minPrime[N];
bool notprime[N];
char s[N];
struct Hash{
	int bin[N],hs[N],base,mod;
	inline int add(int x,int y){x+=y;return x>=mod?x-mod:x;}
	inline int sub(int x,int y){x-=y;return x<0?x+mod:x;}
	inline int mul(int x,int y){return (ll)x*y%mod;}
	inline void init(int _base,int _mod){
		base=_base;mod=_mod;hs[0]=0;bin[0]=1;
		for(int i=1;i<=n;++i) hs[i]=add(mul(hs[i-1],base),s[i]),bin[i]=mul(bin[i-1],base);
	}inline int cal(int x,int y){return sub(hs[y],mul(hs[x-1],bin[y-x+1]));}
	inline bool issame(int l1,int l2,int len){return cal(l1,l1+len-1)==cal(l2,l2+len-1);}
}hs1,hs2;
inline void getPrime(){
	notprime[1]=1;minPrime[1]=1;
	for(int i=2;i<=500000;++i){
		if(!notprime[i]) prime[++tot]=i,minPrime[i]=i;
		for(int j=1;prime[j]*i<=500000;++j){
			notprime[prime[j]*i]=1;minPrime[prime[j]*i]=prime[j];
			if(i%prime[j]==0) break;
		}
	}
}
inline bool jud(int l1,int l2,int len){return hs1.issame(l1,l2,len)&&hs2.issame(l1,l2,len);}
int main(){
//	freopen("a.in","r",stdin);
	scanf("%d%s%d",&n,s+1,&m);
	hs1.init(11113,1e9+7);hs2.init(11119,1e9+9);getPrime();
	while(m--){
		int l,r,ans,rst;scanf("%d%d",&l,&r);ans=rst=r-l+1;
		while(rst!=1){
			int p=minPrime[rst];
			while(ans%p==0&&jud(l,l+ans/p,r-l+1-ans/p)) ans/=p;
			while(rst%p==0) rst/=p;
		}printf("%d\n",ans);
	}return 0;
}

多次询问子串最小循环节。
子串x…y循环节为l的充要条件是x…y-l与x+l…y相同。这个可以由hash来O(1)判断。
循环节长度肯定是y-x+1的因数，很好理解。于是我们可以枚举长度来求最小循环节。
然而这样会T。我们知道循环次数还一定是每个字母出现次数的因数，因此我们对这些取个gcd，枚举gcd的因数即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <map>
using namespace std;
#define ull unsigned long long
#define N 500010
#define k1 11113
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*10+ch-'0',ch=getchar();
	return x*f;
}
int n,cnt[27][N],ans;
ull hs[N],bin[N];char s[N];
inline ull cal(int l,int r){
	return hs[r]-hs[l-1]*bin[r-l+1];
}
inline int gcd(int x,int y){return y?gcd(y,x%y):x;}
inline void solve(int l,int r){
	int len=r-l+1,x=len;
	for(int i=1;i<=26;++i) x=gcd(x,cnt[i][r]-cnt[i][l-1]);
	for(int i=1;i*i<=x;++i){
		if(x%i) continue;int y=len/i;
		if(cal(l,r-y)==cal(l+y,r)) ans=min(ans,y);y=len/(x/i);
		if(cal(l,r-y)==cal(l+y,r)) ans=min(ans,y); 
	}
}
int main(){
//	freopen("a.in","r",stdin);
	n=read();scanf("%s",s+1);bin[0]=1;
	for(int i=1;i<=n;++i)
		for(int j=1;j<=26;++j) cnt[j][i]=cnt[j][i-1]+(s[i]-'a'+1==j);
	for(int i=1;i<=n;++i) hs[i]=hs[i-1]*k1+s[i],bin[i]=bin[i-1]*k1;
	int owo=read();
	while(owo--){
		int x=read(),y=read();ans=y-x+1;solve(x,y);
		printf("%d\n",ans);
	}return 0;
}