bzoj2795 [Poi2012]A Horrible Poem hash+乱搞

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;
}