题解:

首先判断长度len是否是l-r的循环节可以做到O(1)

只需判断 l-l+len-1 和 r-len+1-r 是否相同即可,这一步可以用hash

然后循环节一定是总长度的约数。

这样可以做到n√n,但是会T。

一个鬼畜的优化是 枚举在每块内每个字母的个数,这个个数必须是总个数的约数,然后我们取一下gcd,再枚举约数。

实测可过。

另一种做法可以戳这里:http://blog.csdn.net/zeyu_king/article/details/41989279

代码:

 1 #include<cstdio>
 2 #include<cstdlib>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #include<iostream>
 7 #include<vector>
 8 #include<map>
 9 #include<set>
10 #include<queue>
11 #include<string>
12 #define inf 1000000000
13 #define maxn 600000
14 #define maxm 500+100
15 #define eps 1e-10
16 #define ull unsigned long long
17 #define pa pair<int,int>
18 #define for0(i,n) for(int i=0;i<=(n);i++)
19 #define for1(i,n) for(int i=1;i<=(n);i++)
20 #define for2(i,x,y) for(int i=(x);i<=(y);i++)
21 #define for3(i,x,y) for(int i=(x);i>=(y);i--)
22 #define mod 1000000007
23 using namespace std;
24 inline int read()
25 {
26     int x=0,f=1;char ch=getchar();
27     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
28     while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
29     return x*f;
30 }
31 int n,m,a[maxn],c[maxn],s[maxn][30];
32 ull b[maxn],power[maxn];
33 const ull base=13131;
34 inline ull get(int x,int y){return b[x]-b[x+y]*power[y];}
35 inline int gcd(int x,int y){return y?gcd(y,x%y):x;}
36 inline int solve()
37 {
38     int x=read(),y=read(),z=y-x+1,len=y-x+1,t=0;
39     for1(i,26)z=gcd(z,s[y][i]-s[x-1][i]);
40     if(z==1)return len;
41     for(int i=1;i*i<=z;i++)if(z%i==0)
42     {
43        int tmp=len/z*i;
44        if(get(x,len-tmp)==get(y-(len-tmp)+1,len-tmp))return tmp;
45        c[++t]=len/i;
46     }
47     for3(i,t,1)if(get(x,len-c[i])==get(y-(len-c[i])+1,len-c[i]))return c[i];
48     return len;
49 }
50 int main()
51 {
52     freopen("input.txt","r",stdin);
53     freopen("output.txt","w",stdout);
54     n=read();n=0;
55     char ch=getchar();
56     while(ch<'a'||ch>'z')ch=getchar();
57     for(;ch>='a'&&ch<='z';ch=getchar())a[++n]=ch-'a'+1;
58     for3(i,n,1)b[i]=b[i+1]*base+a[i];
59     for1(i,n)for1(j,26)s[i][j]=s[i-1][j]+(a[i]==j);
60     power[0]=1;
61     for1(i,n)power[i]=power[i-1]*base;
62     m=read();
63     while(m--)printf("%d\n",solve());
64     return 0;
65 }
View Code

相关文章: