BZOJ2795/POI2012 A horrible poem

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

Solution
假如最短循环节为$len$,设$t=n/len$,现在的问题就是找到最大的t.
t肯定是n的因子,那就可以直接枚举n的所有素因子k,若当前的串可以分成k份,那么把n/k,把串变为原来的1/k,继续操作,直到枚举所有的素因子.
如何判断一个串$[l,r]$是否可以分成相同的k份?哈希!!
直接比较$[l,r-k]$和$[l+k,r]$这两段是否完全相同即可.
若相等:
保证$[l,k-1]=[l+k,l+2k-1]=[l+3k,l+4k-1]…=[r-k,r]$,满足了循环的性质.
一个数最多有$logn$个因子,最后复杂度为$O(qlogn)$.

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#define ll long long
#define ull unsigned long long 
using namespace std;
inline void rd(int &res){
    res=0;char c;
    while(c=getchar(),c<48);
    do res=(res<<1)+(res<<3)+(c^48);
    while(c=getchar(),c>=48);
}
inline void print(int x){
    if(!x)return ;
    print(x/10);
    putchar((x%10)^48);
}
inline void sc(int x){
    if(x<0){x=-x;putchar('-');}
    print(x);
    if(!x)putchar('0');
    putchar('\n');
}
const int M=500005;
const int P=1e9+7,Lay=233;
int n,mark[M],son[M][10];
char A[M];
ull v[M],p[M];
void Init(){
    for(int i=2;i*i<M;i++){
        if(!mark[i]){
            for(int j=i;j<M;j+=i){
                son[j][++son[j][0]]=i;
                mark[j]=1;
            }
        }
    }
}
int get(int l,int r){
    return v[r]-v[l-1]*p[r-l+1];
}
bool chk(int st,int en,int len){//以st为起点,len为循环节 
    if(st==en&&len==1)return true;
    return get(st,en-len)==get(st+len,en);
}
int main(){
    int i,j,k,a,b,Q;
    scanf("%d",&n);
    scanf("%s",A+1);
    p[0]=1;
    for(i=1;i<=n;i++){
        v[i]=v[i-1]*Lay+A[i];
        p[i]=p[i-1]*Lay;
    }
    Init();
    scanf("%d",&Q);
    for(i=1;i<=Q;i++){
        rd(a);rd(b);
        int x=b-a+1,y=b-a+1,z=b-a+1;
        for(j=1;j<=son[x][0];j++){
            while(y%son[x][j]==0){
                z/=son[x][j];
                if(chk(a,a+y-1,y/son[x][j]))y/=son[x][j];
                else break;
            }
            while(z%son[x][j]==0)z/=son[x][j];
        }
        if(z!=1)if(chk(a,a+y-1,y/z))y/=z;
        sc(y);
    }
    return 0;
}