【POI2012】【BZOJ2795】A Horrible Poem

Description

给出一个由小写英文字母组成的字符串S，再给出q个询问，要求回答S某个子串的最短循环节。
如果字符串B是字符串A的循环节，那么A可以由B重复若干次得到。
Input

第一行一个正整数n (n<=500,000)，表示S的长度。
第二行n个小写英文字母，表示字符串S。
第三行一个正整数q (q<=2,000,000)，表示询问个数。
下面q行每行两个正整数a,b (1<=a<=b<=n)，表示询问字符串S[a..b]的最短循环节长度。

Output

依次输出q行正整数，第i行的正整数对应第i个询问的答案。

Sample Input
8
aaabcabc
3
1 3
3 8
4 8

Sample Output
1
3
5
HINT

Source

鸣谢 jiangzoi&oimaster

如果只有一个询问,显然直接KMP就可以
多个询问就hash吧
枚举循环长度,循环长度为询问长度的约数

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#define MAXN 500010#define GET (ch>='0'&&ch<='9')#define ULL unsigned long long#define base 313using namespace std;char ch[MAXN];int n,q;int l,r,maxn;int cnt[MAXN][30];ULL Pow[MAXN],hash[MAXN];int gcd(int a,int b){    return !b?a:gcd(b,a%b);}void in(int &x){    char ch=getchar();x=0;    while (!GET)    ch=getchar();    while (GET) x=x*10+ch-'0',ch=getchar();}void check(int l,int r,int num){    int len=(r-l+1)-(r-l+1)/num;    ULL A=hash[l+len-1]-Pow[len]*hash[l-1],B=hash[r]-hash[r-len]*Pow[len];    if (A==B)   maxn=max(maxn,num);}void solve(int l,int r){    int divisor=r-l+1;maxn=0;    for (int i=0;i<26;i++)  divisor=gcd(divisor,cnt[r][i]-cnt[l-1][i]);    for (int i=1;i*i<=divisor;i++)          if (divisor%i==0)   check(l,r,i),check(l,r,divisor/i);}int main(){    scanf("%d%s",&n,ch+1);Pow[0]=1;    for (int i=1;i<=n;i++)  hash[i]=hash[i-1]*base+ch[i]-'a'+1,Pow[i]=Pow[i-1]*base;    for (int i=1;i<=n;i++)        for (int j=0;j<26;j++)  cnt[i][j]=cnt[i-1][j]+(ch[i]-'a'==j);    in(q);    for (int i=1;i<=q;i++)  in(l),in(r),solve(l,r),printf("%d\n",(r-l+1)/maxn);}